Garfield::TrackBichsel Class Reference

#include <TrackBichsel.hh>

Inheritance diagram for Garfield::TrackBichsel:

Public Member Functions
	TrackBichsel ()
	~TrackBichsel ()
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)
double	GetClusterDensity ()
double	GetStoppingPower ()
void	SetDataFile (const std::string filename)
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 TrackBichsel.hh.

Constructor & Destructor Documentation

TrackBichsel()

Garfield::TrackBichsel::TrackBichsel

(

)

Definition at line 15 of file TrackBichsel.cc.

    : bg(3.16228),
      speed(SpeedOfLight * bg / sqrt(1. + bg * bg)),
      x(0.),
      y(0.),
      z(0.),
      t(0.),
      dx(0.),
      dy(0.),
      dz(1.),
      imfp(4.05090e4),
      datafile("SiM0invw.inv"),
      iCdf(2),
      nCdfEntries(-1),
      isInitialised(false),
      isInMedium(false) {
 
  className = "TrackBichsel";
}

~TrackBichsel()

Garfield::TrackBichsel::~TrackBichsel ( )

inline

Definition at line 19 of file TrackBichsel.hh.

{}

Member Function Documentation

GetCluster()

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

virtual

Implements Garfield::Track.

Definition at line 116 of file TrackBichsel.cc.

                                                                              {
 
  if (!isInitialised || !isInMedium) return false;
 
  double d = -log(RndmUniformPos()) / imfp;
  x += dx * d;
  y += dy * d;
  z += dz * d;
  t += d / speed;
 
  xcls = x;
  ycls = y;
  zcls = z;
  tcls = t;
  n = 0;
  e = 0.;
  extra = 0.;
 
  Medium* medium;
  if (!sensor->GetMedium(x, y, z, medium)) {
    isInMedium = false;
    if (debug) {
      std::cout << className << "::GetCluster:\n";
      std::cout << "    Particle left the medium.\n";
    }
    return false;
  }
 
  if (medium->GetName() != "Si" || !medium->IsIonisable()) {
    isInMedium = false;
    if (debug) {
      std::cout << className << "::GetCluster:\n";
      std::cout << "    Particle left the medium.\n";
    }
    return false;
  }
 
  const double u = nCdfEntries * RndmUniform();
  const int j = int(u);
  if (j == 0) {
    e = 0. + u * cdf[0][iCdf];
  } else if (j >= nCdfEntries) {
    e = cdf[nCdfEntries - 1][iCdf];
  } else {
    e = cdf[j - 1][iCdf] + (u - j) * (cdf[j][iCdf] - cdf[j - 1][iCdf]);
  }
 
  return true;
}

GetClusterDensity()

double Garfield::TrackBichsel::GetClusterDensity ( )

virtual

Reimplemented from Garfield::Track.

Definition at line 167 of file TrackBichsel.cc.

                                       {
 
  const int nEntries = 38;
 
  const double tabBg[nEntries] = {
      0.316,   0.398,   0.501,   0.631,    0.794,    1.000,   1.259,   1.585,
      1.995,   2.512,   3.162,   3.981,    5.012,    6.310,   7.943,   10.000,
      12.589,  15.849,  19.953,  25.119,   31.623,   39.811,  50.119,  63.096,
      79.433,  100.000, 125.893, 158.489,  199.526,  251.189, 316.228, 398.107,
      501.187, 630.958, 794.329, 1000.000, 1258.926, 1584.894};
 
  const double tabImfp[nEntries] = {
      30.32496, 21.14965, 15.06555, 11.05635, 8.43259, 6.72876, 5.63184,
      4.93252,  4.49174,  4.21786,  4.05090,  3.95186, 3.89531, 3.86471,
      3.84930,  3.84226,  3.83952,  3.83887,  3.83912, 3.83970, 3.84035,
      3.84095,  3.84147,  3.84189,  3.84223,  3.84249, 3.84269, 3.84283,
      3.84293,  3.84300,  3.84304,  3.84308,  3.84310, 3.84311, 3.84312,
      3.84313,  3.84313,  3.84314};
 
  if (isChanged) bg = GetBetaGamma();
 
  if (bg < tabBg[0]) {
    if (debug) {
      std::cerr << className << "::GetClusterDensity:\n";
      std::cerr << "    Bg is below the tabulated range.\n";
    }
    return tabImfp[0] * 1.e4;
  } else if (bg > tabBg[nEntries - 1]) {
    return tabImfp[nEntries - 1] * 1.e4;
  }
 
  // Locate the requested energy in the table
  int iLow = 0;
  int iUp = nEntries - 1;
  int iM;
  while (iUp - iLow > 1) {
    iM = (iUp + iLow) >> 1;
    if (bg >= tabBg[iM]) {
      iLow = iM;
    } else {
      iUp = iM;
    }
  }
 
  if (fabs(bg - tabBg[iLow]) < 1.e-6 * (tabBg[iUp] - tabBg[iLow])) {
    return tabImfp[iLow] * 1.e4;
  }
  if (fabs(bg - tabBg[iUp]) < 1.e-6 * (tabBg[iUp] - tabBg[iLow])) {
    return tabImfp[iUp] * 1.e4;
  }
 
  // Log-log interpolation
  const double logX0 = log(tabBg[iLow]);
  const double logX1 = log(tabBg[iUp]);
  const double logY0 = log(tabImfp[iLow]);
  const double logY1 = log(tabImfp[iUp]);
  double d = logY0 + (log(bg) - logX0) * (logY1 - logY0) / (logX1 - logX0);
  return 1.e4 * exp(d);
}

Referenced by NewTrack().

GetStoppingPower()

double Garfield::TrackBichsel::GetStoppingPower ( )

virtual

Reimplemented from Garfield::Track.

Definition at line 227 of file TrackBichsel.cc.

                                      {
 
  const int nEntries = 51;
 
  const double tabBg[nEntries] = {
      0.316,     0.398,    0.501,    0.631,     0.794,     1.000,     1.259,
      1.585,     1.995,    2.512,    3.162,     3.981,     5.012,     6.310,
      7.943,     10.000,   12.589,   15.849,    19.953,    25.119,    31.623,
      39.811,    50.119,   63.096,   79.433,    100.000,   125.893,   158.489,
      199.526,   251.189,  316.228,  398.107,   501.187,   630.958,   794.329,
      1000.000,  1258.926, 1584.894, 1995.263,  2511.888,  3162.280,  3981.074,
      5011.875,  6309.578, 7943.287, 10000.010, 12589.260, 15848.940, 19952.640,
      25118.880, 31622.800};
 
  const double tabdEdx[nEntries] = {
      2443.71800, 1731.65600, 1250.93400, 928.69920, 716.37140, 578.28850,
      490.83670,  437.33820,  406.58490,  390.95170, 385.29000, 386.12000,
      391.07730,  398.53930,  407.39420,  416.90860, 426.63010, 436.30240,
      445.78980,  455.02530,  463.97370,  472.61410, 480.92980, 488.90240,
      496.51900,  503.77130,  510.65970,  517.19570, 523.39830, 529.29120,
      534.90670,  540.27590,  545.42880,  550.39890, 555.20800, 559.88820,
      564.45780,  568.93850,  573.34700,  577.69140, 581.99010, 586.25090,
      590.47720,  594.68660,  598.86880,  603.03510, 607.18890, 611.33250,
      615.46810,  619.59740,  623.72150};
 
  if (isChanged) bg = GetBetaGamma();
 
  if (bg < tabBg[0]) {
    if (debug) {
      std::cerr << className << "::GetStoppingPower:\n";
      std::cerr << "    Bg is below the tabulated range.\n";
    }
    return tabdEdx[0] * 1.e4;
  } else if (bg > tabBg[nEntries - 1]) {
    return tabdEdx[nEntries - 1] * 1.e4;
  }
 
  // Locate the requested energy in the table
  int iLow = 0;
  int iUp = nEntries - 1;
  int iM;
  while (iUp - iLow > 1) {
    iM = (iUp + iLow) >> 1;
    if (bg >= tabBg[iM]) {
      iLow = iM;
    } else {
      iUp = iM;
    }
  }
 
  if (debug) {
    std::cout << className << "::GetStoppingPower:\n";
    std::cout << "    Bg = " << bg << "\n";
    std::cout << "    Interpolating between " << tabBg[iLow] << " and "
              << tabBg[iUp] << "\n";
  }
 
  if (fabs(bg - tabBg[iLow]) < 1.e-6 * (tabBg[iUp] - tabBg[iLow])) {
    return tabdEdx[iLow] * 1.e4;
  }
  if (fabs(bg - tabBg[iUp]) < 1.e-6 * (tabBg[iUp] - tabBg[iLow])) {
    return tabdEdx[iUp] * 1.e4;
  }
 
  // Log-log interpolation
  const double logX0 = log(tabBg[iLow]);
  const double logX1 = log(tabBg[iUp]);
  const double logY0 = log(tabdEdx[iLow]);
  const double logY1 = log(tabdEdx[iUp]);
  const double dedx =
      logY0 + (log(bg) - logX0) * (logY1 - logY0) / (logX1 - logX0);
  return 1.e4 * exp(dedx);
}

NewTrack()

bool Garfield::TrackBichsel::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 35 of file TrackBichsel.cc.

                                              {
 
  // Make sure a sensor has been defined.
  if (sensor == 0) {
    std::cerr << className << "::NewTrack:\n";
    std::cerr << "    Sensor is not defined.\n";
    isInMedium = false;
    return false;
  }
 
  // If not yet done, load the cross-section table from file.
  if (!isInitialised) {
    if (!LoadCrossSectionTable(datafile)) {
      std::cerr << className << "::NewTrack:\n";
      std::cerr << "    Cross-section table could not be loaded.\n";
      return false;
    }
    isInitialised = true;
  }
 
  // 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";
    isInMedium = false;
    return false;
  }
 
  // Check if the medium is silicon.
  if (medium->GetName() != "Si") {
    std::cerr << className << "::NewTrack:" << std::endl;
    std::cerr << "    Medium at initial position is not silicon.\n";
    isInMedium = false;
    return false;
  }
 
  // Check if primary ionisation has been enabled.
  if (!medium->IsIonisable()) {
    std::cerr << className << "::NewTrack:\n";
    std::cerr << "    Medium at initial position is not ionisable.\n";
    isInMedium = false;
    return false;
  }
 
  isInMedium = true;
  x = x0;
  y = y0;
  z = z0;
  t = t0;
 
  // Normalise the direction vector.
  const double d = sqrt(dx0 * dx0 + dy0 * dy0 + dz0 * dz0);
  if (d < Small) {
    // In case of a null vector, choose a random direction.
    const double phi = TwoPi * RndmUniform();
    const double ctheta = 1. - 2. * RndmUniform();
    const 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;
  }
 
  // If the particle properties have changed, update the cross-section table.
  if (isChanged) {
    bg = GetBetaGamma();
    imfp = GetClusterDensity();
    speed = SpeedOfLight * GetBeta();
    SelectCrossSectionTable();
    isChanged = false;
  }
 
  return true;
}

SetDataFile()

void Garfield::TrackBichsel::SetDataFile ( const std::string  filename )

inline

Definition at line 30 of file TrackBichsel.hh.

{ datafile = filename; }

---

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

• TrackBichsel.hh
• TrackBichsel.cc