Generate tracks based on a cluster density given by the user. More...

#include <TrackSimple.hh>

Inheritance diagram for Garfield::TrackSimple:

Public Member Functions
	TrackSimple ()
	Constructor.

virtual	~TrackSimple ()
	Destructor.

void	SetEqualSpacing ()
	Constant distance between clusters.

void	SetExponentialSpacing ()
	Exponentially distributed distance between clusters.

void	SetClusterDensity (const double d)
	Set the cluster density (inverse mean free path).

virtual double	GetClusterDensity ()

void	SetStoppingPower (const double dedx)
	Set the stopping power (dE/dx).

virtual double	GetStoppingPower ()
	Get the stopping power (mean energy loss [eV] per cm).

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

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

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

virtual	~Track ()
	Destructor.

virtual void	SetParticle (const std::string &part)
	Set the type of particle.

void	SetEnergy (const double e)
	Set the particle energy.

void	SetBetaGamma (const double bg)
	Set the relative momentum of the particle.

void	SetBeta (const double beta)
	Set the speed ( $\beta = v/c$ ) of the particle.

void	SetGamma (const double gamma)
	Set the Lorentz factor of the particle.

void	SetMomentum (const double p)
	Set the particle momentum.

void	SetKineticEnergy (const double ekin)
	Set the kinetic energy of the particle.

double	GetEnergy () const

double	GetBetaGamma () const

double	GetBeta () const

double	GetGamma () const

double	GetMomentum () const

double	GetKineticEnergy () const

double	GetCharge () const
	Get the charge of the projectile.

double	GetMass () const
	Get the mass [eV / c2] of the projectile.

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 ()
	Get the stopping power (mean energy loss [eV] per cm).

void	EnablePlotting (ViewDrift *viewer)

void	DisablePlotting ()

void	EnableDebugging ()

void	DisableDebugging ()

Protected Attributes
bool	m_isReady

double	m_x

double	m_y

double	m_z

double	m_t

double	m_dx

double	m_dy

double	m_dz

double	m_mfp

double	m_eloss

bool	m_useEqualSpacing

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

double	m_q

int	m_spin

double	m_mass

double	m_energy

double	m_beta2

bool	m_isElectron

std::string	m_particleName

Sensor *	m_sensor

bool	m_isChanged

bool	m_usePlotting

ViewDrift *	m_viewer

bool	m_debug

int	m_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

Generate tracks based on a cluster density given by the user.

Definition at line 10 of file TrackSimple.hh.

Constructor & Destructor Documentation

◆ TrackSimple()

Garfield::TrackSimple::TrackSimple ( )

Constructor.

Definition at line 11 of file TrackSimple.cc.

    : m_isReady(false),
      m_x(0.),
      m_y(0.),
      m_z(0.),
      m_t(0.),
      m_dx(0.),
      m_dy(0.),
      m_dz(1.),
      m_mfp(0.04),
      m_eloss(2530.),
      m_useEqualSpacing(false) {
 
  m_className = "TrackSimple";
}

◆ ~TrackSimple()

virtual Garfield::TrackSimple::~TrackSimple ( )

inlinevirtual

Destructor.

Definition at line 16 of file TrackSimple.hh.

16{}

Member Function Documentation

◆ GetCluster()

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

virtual

Get the next "cluster" (ionising collision of the charged particle).

Parameters

xcls,ycls,zcls	coordinates of the collision
tcls	time of the collision
n	number of electrons produced
e	deposited energy
extra	additional information (not always implemented)

Implements Garfield::Track.

Definition at line 96 of file TrackSimple.cc.

                                                                             {
 
  extra = 0.;
  if (!m_isReady) return false;
 
  if (m_useEqualSpacing) {
    m_x += m_dx * m_mfp;
    m_y += m_dy * m_mfp;
    m_z += m_dz * m_mfp;
  } else {
    const double d = -m_mfp * log(RndmUniformPos());
    m_x += m_dx * d;
    m_y += m_dy * d;
    m_z += m_dz * d;
  }
 
  xcls = m_x;
  ycls = m_y;
  zcls = m_z;
  tcls = m_t;
 
  n = 1;
  e = m_eloss * m_mfp;
 
  Medium* medium = NULL;
  if (!m_sensor->GetMedium(m_x, m_y, m_z, medium)) {
    m_isReady = false;
    if (m_debug) {
      std::cout << m_className << "::GetCluster: Particle left the medium.\n";
    }
    return false;
  }
 
  return true;
}

◆ GetClusterDensity()

double Garfield::TrackSimple::GetClusterDensity ( )

virtual

Get the cluster density (number of ionizing collisions per cm or inverse mean free path for ionization).

Reimplemented from Garfield::Track.

Definition at line 39 of file TrackSimple.cc.

39{ return 1. / m_mfp; }

◆ GetStoppingPower()

double Garfield::TrackSimple::GetStoppingPower ( )

virtual

Get the stopping power (mean energy loss [eV] per cm).

Reimplemented from Garfield::Track.

Definition at line 53 of file TrackSimple.cc.

53{ return m_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

Calculate a new track starting from (x0, y0, z0) at time t0 in direction (dx0, dy0, dz0).

Implements Garfield::Track.

Definition at line 55 of file TrackSimple.cc.

                                             {
 
  // Check if a sensor has been defined
  if (!m_sensor) {
    std::cerr << m_className << "::NewTrack:\n"
              << "    Sensor is not defined.\n";
    m_isReady = false;
    return false;
  }
 
  // Make sure we are inside a medium
  Medium* medium = NULL;
  if (!m_sensor->GetMedium(x0, y0, z0, medium)) {
    std::cerr << m_className << "::NewTrack:\n";
    std::cerr << "    No medium at initial position.\n";
    m_isReady = false;
    return false;
  }
 
  m_isReady = true;
 
  m_x = x0;
  m_y = y0;
  m_z = z0;
  m_t = t0;
 
  // Normalise the direction.
  const double d = sqrt(dx0 * dx0 + dy0 * dy0 + dz0 * dz0);
  if (d < Small) {
    // Choose random direction.
    RndmDirection(m_dx, m_dy, m_dz);
  } else {
    m_dx = dx0 / d;
    m_dy = dy0 / d;
    m_dz = dz0 / d;
  }
  return true;
}

◆ SetClusterDensity()

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

Set the cluster density (inverse mean free path).

Definition at line 27 of file TrackSimple.cc.

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

◆ SetEqualSpacing()

void Garfield::TrackSimple::SetEqualSpacing ( )

inline

Constant distance between clusters.

Definition at line 19 of file TrackSimple.hh.

19{ m_useEqualSpacing = true; }

◆ SetExponentialSpacing()

void Garfield::TrackSimple::SetExponentialSpacing ( )

inline

Exponentially distributed distance between clusters.

Definition at line 21 of file TrackSimple.hh.

21{ m_useEqualSpacing = false; }

◆ SetStoppingPower()

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

Set the stopping power (dE/dx).

Definition at line 41 of file TrackSimple.cc.

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

Member Data Documentation

◆ m_dx

double Garfield::TrackSimple::m_dx

protected

Definition at line 41 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

◆ m_dy

double Garfield::TrackSimple::m_dy

protected

Definition at line 41 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

◆ m_dz

double Garfield::TrackSimple::m_dz

protected

Definition at line 41 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

◆ m_eloss

double Garfield::TrackSimple::m_eloss

protected

Definition at line 45 of file TrackSimple.hh.

Referenced by GetCluster(), GetStoppingPower(), and SetStoppingPower().

◆ m_isReady

bool Garfield::TrackSimple::m_isReady

protected

Definition at line 37 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

◆ m_mfp

double Garfield::TrackSimple::m_mfp

protected

Definition at line 43 of file TrackSimple.hh.

Referenced by GetCluster(), GetClusterDensity(), and SetClusterDensity().

◆ m_t

double Garfield::TrackSimple::m_t

protected

Definition at line 40 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

◆ m_useEqualSpacing

bool Garfield::TrackSimple::m_useEqualSpacing

protected

Definition at line 47 of file TrackSimple.hh.

Referenced by GetCluster(), SetEqualSpacing(), and SetExponentialSpacing().

◆ m_x

double Garfield::TrackSimple::m_x

protected

Definition at line 40 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

◆ m_y

double Garfield::TrackSimple::m_y

protected

Definition at line 40 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

◆ m_z

double Garfield::TrackSimple::m_z

protected

Definition at line 40 of file TrackSimple.hh.

Referenced by GetCluster(), and NewTrack().

The documentation for this class was generated from the following files:

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TrackSimple()

◆ ~TrackSimple()

Member Function Documentation

◆ GetCluster()

◆ GetClusterDensity()

◆ GetStoppingPower()

◆ NewTrack()

◆ SetClusterDensity()

◆ SetEqualSpacing()

◆ SetExponentialSpacing()

◆ SetStoppingPower()

Member Data Documentation

◆ m_dx

◆ m_dy

◆ m_dz

◆ m_eloss

◆ m_isReady

◆ m_mfp

◆ m_t

◆ m_useEqualSpacing

◆ m_x

◆ m_y

◆ m_z