#include <TrackSimple.hh>

Inheritance diagram for Garfield::TrackSimple:

Public Member Functions
	TrackSimple ()

	~TrackSimple ()

void	SetEqualSpacing ()

void	SetExponentialSpacing ()

void	SetClusterDensity (const double d)

double	GetClusterDensity ()

void	SetStoppingPower (const double dedx)

double	GetStoppingPower ()

bool	NewTrack (const double x0, const double y0, const double z0, const double t0, const double dx0, const double dy0, const double dz0)

bool	GetCluster (double &xcls, double &ycls, double &zcls, double &tcls, int &n, double &e, double &extra)

Public Member Functions inherited from Garfield::Track
	Track ()

virtual	~Track ()

virtual void	SetParticle (std::string part)

void	SetEnergy (const double e)

void	SetBetaGamma (const double bg)

void	SetBeta (const double beta)

void	SetGamma (const double gamma)

void	SetMomentum (const double p)

void	SetKineticEnergy (const double ekin)

double	GetEnergy () const

double	GetBetaGamma () const

double	GetBeta () const

double	GetGamma () const

double	GetMomentum () const

double	GetKineticEnergy () const

void	SetSensor (Sensor *s)

virtual bool	NewTrack (const double x0, const double y0, const double z0, const double t0, const double dx0, const double dy0, const double dz0)=0

virtual bool	GetCluster (double &xcls, double &ycls, double &zcls, double &tcls, int &n, double &e, double &extra)=0

virtual double	GetClusterDensity ()

virtual double	GetStoppingPower ()

void	EnablePlotting (ViewDrift *viewer)

void	DisablePlotting ()

void	EnableDebugging ()

void	DisableDebugging ()

Protected Attributes
bool	isReady

double	x

double	y

double	z

double	t

double	dx

double	dy

double	dz

double	mfp

double	eloss

bool	useEqualSpacing

Protected Attributes inherited from Garfield::Track
std::string	className

double	q

int	spin

double	mass

double	energy

double	beta2

bool	isElectron

std::string	particleName

Sensor *	sensor

bool	isChanged

bool	usePlotting

ViewDrift *	viewer

bool	debug

int	plotId

Additional Inherited Members
Protected Member Functions inherited from Garfield::Track
void	PlotNewTrack (const double x0, const double y0, const double z0)

void	PlotCluster (const double x0, const double y0, const double z0)

Detailed Description

Definition at line 9 of file TrackSimple.hh.

Constructor & Destructor Documentation

◆ TrackSimple()

Garfield::TrackSimple::TrackSimple ( )

Definition at line 11 of file TrackSimple.cc.

    : isReady(false),
      x(0.),
      y(0.),
      z(0.),
      t(0.),
      dx(0.),
      dy(0.),
      dz(1.),
      mfp(0.04),
      eloss(2530.),
      useEqualSpacing(false) {
 
  className = "TrackSimple";
}

◆ ~TrackSimple()

Garfield::TrackSimple::~TrackSimple ( )

inline

Definition at line 15 of file TrackSimple.hh.

15{}

Member Function Documentation

◆ GetCluster()

bool Garfield::TrackSimple::GetCluster	(	double &	xcls,
		double &	ycls,
		double &	zcls,
		double &	tcls,
		int &	n,
		double &	e,
		double &	extra
	)

virtual

Implements Garfield::Track.

Definition at line 101 of file TrackSimple.cc.

                                                                             {
 
  extra = 0.;
  if (!isReady) return false;
 
  if (useEqualSpacing) {
    x += dx * mfp;
    y += dy * mfp;
    z += dz * mfp;
  } else {
    const double d = -mfp * log(RndmUniformPos());
    x += dx * d;
    y += dy * d;
    z += dz * d;
  }
 
  xcls = x;
  ycls = y;
  zcls = z;
  tcls = t;
 
  n = 1;
  e = eloss * mfp;
 
  Medium* medium;
  if (!sensor->GetMedium(x, y, z, medium)) {
    isReady = false;
    if (debug) {
      std::cout << className << "::GetCluster:\n";
      std::cout << "    Particle left the medium.\n";
    }
    return false;
  }
 
  return true;
}

◆ GetClusterDensity()

double Garfield::TrackSimple::GetClusterDensity ( )

virtual

Reimplemented from Garfield::Track.

Definition at line 39 of file TrackSimple.cc.

39{ return 1. / mfp; }

◆ GetStoppingPower()

double Garfield::TrackSimple::GetStoppingPower ( )

virtual

Reimplemented from Garfield::Track.

Definition at line 53 of file TrackSimple.cc.

53{ return eloss; }

◆ NewTrack()

bool Garfield::TrackSimple::NewTrack	(	const double	x0,
		const double	y0,
		const double	z0,
		const double	t0,
		const double	dx0,
		const double	dy0,
		const double	dz0
	)

virtual

Implements Garfield::Track.

Definition at line 55 of file TrackSimple.cc.

                                             {
 
  // Check if a sensor has been defined
  if (sensor == 0) {
    std::cerr << className << "::NewTrack:\n";
    std::cerr << "    Sensor is not defined.\n";
    isReady = false;
    return false;
  }
 
  // Make sure we are inside a medium
  Medium* medium;
  if (!sensor->GetMedium(x0, y0, z0, medium)) {
    std::cerr << className << "::NewTrack:\n";
    std::cerr << "    No medium at initial position.\n";
    isReady = false;
    return false;
  }
 
  isReady = true;
 
  x = x0;
  y = y0;
  z = z0;
  t = t0;
 
  // Normalise the direction
  const double d = sqrt(dx0 * dx0 + dy0 * dy0 + dz0 * dz0);
  if (d < Small) {
    // Choose random direction
    double phi = TwoPi * RndmUniform();
    double ctheta = 1. - 2. * RndmUniform();
    double stheta = sqrt(1. - ctheta * ctheta);
    dx = cos(phi) * stheta;
    dy = sin(phi) * stheta;
    dz = ctheta;
  } else {
    dx = dx0 / d;
    dy = dy0 / d;
    dz = dz0 / d;
  }
  return true;
}

◆ SetClusterDensity()

void Garfield::TrackSimple::SetClusterDensity ( const double d )

Definition at line 27 of file TrackSimple.cc.

                                                  {
 
  if (d < Small) {
    std::cerr << className << "::SetClusterDensity:\n";
    std::cerr << "    Cluster density (number of clusters per cm)"
              << " must be positive.\n" << std::endl;
    return;
  }
 
  mfp = 1. / d;
}

◆ SetEqualSpacing()

void Garfield::TrackSimple::SetEqualSpacing ( )

inline

Definition at line 17 of file TrackSimple.hh.

17{ useEqualSpacing = true; }

◆ SetExponentialSpacing()

void Garfield::TrackSimple::SetExponentialSpacing ( )

inline

Definition at line 18 of file TrackSimple.hh.

18{ useEqualSpacing = false; }

◆ SetStoppingPower()

void Garfield::TrackSimple::SetStoppingPower ( const double dedx )

Definition at line 41 of file TrackSimple.cc.

                                                    {
 
  if (dedx < Small) {
    std::cerr << className << "::SetStoppingPower:\n";
    std::cerr << "    Stopping power (average energy loss [eV] per cm)"
              << " must be positive.\n";
    return;
  }
 
  eloss = dedx;
}