Garfield::TrackPAI Class Reference

#include <TrackPAI.hh>

Inheritance diagram for Garfield::TrackPAI:

Public Member Functions
	TrackPAI ()
	~TrackPAI ()
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 &ncls, double &ecls, double &extra)
double	GetClusterDensity ()
double	GetStoppingPower ()
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 ()

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)
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

Detailed Description

Definition at line 13 of file TrackPAI.hh.

Constructor & Destructor Documentation

TrackPAI()

Garfield::TrackPAI::TrackPAI

(

)

Definition at line 13 of file TrackPAI.cc.

    : ready(false),
      x(0.),
      y(0.),
      z(0.),
      t(0.),
      dx(0.),
      dy(0),
      dz(1.),
      e(0.),
      speed(0.),
      emax(0.),
      imfp(0.),
      dedx(0.),
      nSteps(1000),
      mediumName(""),
      mediumDensity(0.),
      electronDensity(0.) {
 
  className = "TrackPAI";
 
  electrons.clear();
  holes.clear();
}

~TrackPAI()

Garfield::TrackPAI::~TrackPAI ( )

inline

Definition at line 19 of file TrackPAI.hh.

{}

Member Function Documentation

GetCluster()

bool Garfield::TrackPAI::GetCluster ( double &  xcls, double &  ycls, double &  zcls, double &  tcls, int &  ncls, double &  ecls, double &  extra  )

virtual

Implements Garfield::Track.

Definition at line 115 of file TrackPAI.cc.

                                         {
 
  ncls = 0;
  edep = extra = 0.;
 
  // Clear the stack.
  electrons.clear();
  holes.clear();
 
  if (!ready) {
    std::cerr << className << "::GetCluster:\n";
    std::cerr << "    Track not initialized.\n";
    std::cerr << "    Call NewTrack first.\n";
    return false;
  }
 
  if (isChanged) {
    if (SetupCrossSectionTable()) {
      isChanged = false;
    } else {
      std::cerr << className << "::GetCluster:\n";
      std::cerr << "    Calculation of ionisation cross-section failed.\n";
      return false;
    }
  }
 
  // Draw a step length and propagate the particle.
  const double d = -imfp * log(RndmUniformPos());
  x += d * dx;
  y += d * dy;
  z += d * dz;
  t += d / speed;
 
  // Check the medium at this location.
  Medium* medium = 0;
  if (!sensor->GetMedium(x, y, z, medium)) {
    ready = false;
    return false;
  }
  if (medium->GetName() != mediumName ||
      medium->GetNumberDensity() != mediumDensity || !medium->IsIonisable()) {
    ready = false;
    return false;
  }
 
  // Check if the particle is still inside the drift area.
  if (!sensor->IsInArea(x, y, z)) {
    ready = false;
    return false;
  }
 
  xcls = x;
  ycls = y;
  zcls = z;
  tcls = t;
 
  // Sample the energy deposition.
  double f = 0.;
  const double u = RndmUniform();
  if (u < cdf.back()) {
    if (u <= cdf[0]) {
      edep = energies[0];
    } else if (u >= 1.) {
      edep = energies.back();
    } else {
      // Find the energy loss by interpolation
      // from the cumulative distribution table
      int iLow = 0, iUp = cdf.size(), iMid;
      while (iUp - iLow > 1) {
        iMid = (iUp + iLow) >> 1;
        if (u >= cdf[iMid]) {
          iLow = iMid;
        } else {
          iUp = iMid;
        }
      }
      if (edep < 100.) {
        edep = energies[iLow] + (u - cdf[iLow]) *
                                    (energies[iUp] - energies[iLow]) /
                                    (cdf[iUp] - cdf[iLow]);
        f = rutherford[iLow] + (edep - energies[iLow]) *
                                   (rutherford[iUp] - rutherford[iLow]) /
                                   (energies[iUp] - energies[iLow]);
      } else {
        edep = log(energies[iLow]) +
               (log(u) - log(cdf[iLow])) *
                   (log(energies[iUp]) - log(energies[iLow])) /
                   (log(cdf[iUp]) - log(cdf[iLow]));
        edep = exp(edep);
        f = rutherford[iLow] + (log(edep) - log(energies[iLow])) *
                                   (rutherford[iUp] - rutherford[iLow]) /
                                   (log(energies[iUp]) - log(energies[iLow]));
      }
    }
  } else {
    // Use the free-electron differential cross-section.
    f = 1.;
    edep = SampleAsymptoticCs(u);
  }
  // Update the particle energy.
  e -= edep;
 
  // Number of electron/hole (or electron/ion pairs) produced.
  ncls = 1;
 
  if (debug) {
    std::cout << className << "::GetCluster:\n";
    std::cout << "   Fraction of Rutherford scattering: " << f << "\n";
  }
  return true;
}

GetClusterDensity()

double Garfield::TrackPAI::GetClusterDensity ( )

virtual

Reimplemented from Garfield::Track.

Definition at line 309 of file TrackPAI.cc.

                                   {
 
  if (!ready) {
    std::cerr << className << "::GetClusterDensity:\n";
    std::cerr << "    Track has not been initialized.\n";
    return 0.;
  }
 
  if (isChanged) {
    if (SetupCrossSectionTable()) {
      isChanged = false;
    } else {
      std::cerr << className << "::GetClusterDensity:\n";
      std::cerr << "    Ionisation cross-section could not be calculated.\n";
      return 0.;
    }
  }
 
  return 1. / imfp;
}

GetStoppingPower()

double Garfield::TrackPAI::GetStoppingPower ( )

virtual

Reimplemented from Garfield::Track.

Definition at line 330 of file TrackPAI.cc.

                                  {
 
  if (!ready) {
    std::cerr << className << "::GetStoppingPower:\n";
    std::cerr << "    Track has not been initialised.\n";
    return 0.;
  }
 
  if (isChanged) {
    if (SetupCrossSectionTable()) {
      isChanged = false;
    } else {
      std::cerr << className << "::GetStoppingPower:\n";
      std::cerr << "    Ionisation cross-section could not be calculated.\n";
      return 0.;
    }
  }
 
  return dedx;
}

NewTrack()

bool Garfield::TrackPAI::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 38 of file TrackPAI.cc.

                                          {
 
  ready = false;
 
  // Make sure the sensor has been set.
  if (sensor == 0) {
    std::cerr << className << "::NewTrack:\n";
    std::cerr << "    Sensor is not defined.\n";
    return false;
  }
 
  // Get the medium at this location and check if it is "ionisable".
  Medium* medium = 0;
  if (!sensor->GetMedium(x0, y0, z0, medium)) {
    std::cerr << className << "::NewTrack:\n";
    std::cerr << "    No medium at initial position.\n";
    return false;
  }
  if (!medium->IsIonisable()) {
    std::cerr << className << "::NewTrack:\n";
    std::cerr << "    Medium at initial position is not ionisable.\n";
    return false;
  }
 
  if (medium->GetName() != mediumName ||
      medium->GetNumberDensity() != mediumDensity) {
    isChanged = true;
    if (!SetupMedium(medium)) {
      std::cerr << className << "::NewTrack:\n";
      std::cerr << "    Properties of medium " << medium->GetName()
                << " are not available.\n";
      return false;
    }
    mediumName = medium->GetName();
    mediumDensity = medium->GetNumberDensity();
  }
 
  ready = true;
 
  if (isChanged) {
    if (!SetupCrossSectionTable()) {
      std::cerr << className << "::NewTrack:\n";
      std::cerr << "    Calculation of ionisation cross-section failed.\n";
      ready = false;
      return false;
    }
    isChanged = false;
  }
 
  x = x0;
  y = y0;
  z = z0;
  t = t0;
  const double d = sqrt(dx0 * dx0 + dy0 * dy0 + dz0 * dz0);
  if (d < Small) {
    if (debug) {
      std::cout << className << "::NewTrack:\n";
      std::cout << "    Direction vector has zero norm.\n";
      std::cout << "    Initial direction is randomized.\n";
    }
    const double ctheta = 1. - 2. * RndmUniform();
    const double stheta = sqrt(1. - ctheta * ctheta);
    const double phi = TwoPi * RndmUniform();
    dx = cos(phi) * stheta;
    dy = sin(phi) * stheta;
    dz = ctheta;
  } else {
    // Normalize the direction vector.
    dx = dx0 / d;
    dy = dy0 / d;
    dz = dz0 / d;
  }
  return true;
}

---

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

• TrackPAI.hh
• TrackPAI.cc