Garfield Namespace Reference

Namespaces
namespace	HeedInterface
namespace	Magboltz
namespace	Numerics

Classes
class	AvalancheMC
class	AvalancheMicroscopic
class	ComponentAnalyticField
class	ComponentAnsys121
class	ComponentAnsys123
class	ComponentBase
class	ComponentConstant
class	ComponentCST
class	ComponentElmer
class	ComponentFieldMap
class	ComponentNeBem2d
class	ComponentTcad2d
class	ComponentTcad3d
class	ComponentUser
class	ComponentUserMapBase
class	ComponentVoxel
class	DriftLineRKF
class	GeometryBase
class	GeometryRoot
class	GeometrySimple
class	HeedChamber
class	Medium
class	MediumCdTe
class	MediumConductor
class	MediumGaAs
class	MediumGas
class	MediumMagboltz
class	MediumPlastic
class	MediumSilicon
class	OpticalData
class	PlottingEngine
class	PlottingEngineRoot
struct	PolygonInfo
class	RandomEngine
class	RandomEngineRoot
class	Sensor
class	Solid
class	SolidBox
class	SolidSphere
class	SolidTube
class	Track
class	TrackBichsel
class	TrackElectron
class	TrackHeed
class	TrackPAI
class	TrackSimple
class	ViewCell
class	ViewDrift
class	ViewFEMesh
class	ViewField
class	ViewGeometry
class	ViewMedium
class	ViewSignal

Functions
double	RndmUniform ()
double	RndmUniformPos ()
double	RndmGaussian ()
double	RndmGaussian (const double mu, const double sigma)
double	RndmLorentzian (const double mu, const double gamma)
double	RndmVoigt (const double mu, const double sigma, const double gamma)
double	RndmPolya (const double theta)

Variables
PlottingEngineRoot	plottingEngine
RandomEngineRoot	randomEngine

Function Documentation

RndmGaussian() [1/2]

double Garfield::RndmGaussian ( )

inline

Definition at line 27 of file Random.hh.

                             {
 
  static bool cached = false;
  static double u = 0.;
  if (cached) {
    cached = false;
    return u;
  }
  // Box-Muller algorithm
  u = 2. * RndmUniform() - 1.;
  double v = 2. * RndmUniform() - 1.;
  double r2 = u * u + v * v;
  while (r2 > 1.) {
    u = 2. * RndmUniform() - 1.;
    v = 2. * RndmUniform() - 1.;
    r2 = u * u + v * v;
  }
  const double p = sqrt(-2. * log(r2) / r2);
  u *= p;
  cached = true;
  return v * p;
}

Referenced by RndmGaussian(), and RndmVoigt().

RndmGaussian() [2/2]

double Garfield::RndmGaussian ( const double  mu, const double  sigma  )

inline

Definition at line 51 of file Random.hh.

                                                                {
 
  return mu + sigma * RndmGaussian();
}

RndmLorentzian()

double Garfield::RndmLorentzian ( const double  mu, const double  gamma  )

inline

Definition at line 58 of file Random.hh.

                                                                  {
 
  return mu + gamma * tan(Pi * (RndmUniform() - 0.5));
}

Referenced by RndmVoigt().

RndmPolya()

double Garfield::RndmPolya ( const double  theta )

inline

Definition at line 77 of file Random.hh.

                                            {
 
  // Algorithm from Review of Particle Physics
  // C. Amsler et al, Phys. Lett. B 667 (2008)
  if (theta <= 0.) return -log(RndmUniformPos());
  const double c = 3 * (theta + 1.) - 0.75;
  double u1, u2, v1, v2, v3;
  double x;
  while (1) {
    u1 = RndmUniformPos();
    v1 = u1 * (1. - u1);
    if (v1 == 0.) continue;
    v2 = (u1 - 0.5) * sqrt(c / v1);
    x = theta + v2;
    if (x <= 0.) continue;
    u2 = RndmUniformPos();
    v3 = 64 * u2 * u2 * pow(v1, 3);
    if (v3 <= 1. - 2 * v2 * v2 / x ||
        log(v3) <= 2 * (theta * log(x / theta) - v2)) {
      return x / (theta + 1.);
    }
  }
}

RndmUniform()

double Garfield::RndmUniform ( )

inline

Definition at line 16 of file Random.hh.

{ return randomEngine.Draw(); }

Referenced by Garfield::MediumSilicon::ComputeSecondaries(), Garfield::TrackBichsel::GetCluster(), Garfield::TrackElectron::GetCluster(), Garfield::TrackPAI::GetCluster(), Garfield::Medium::GetElectronCollision(), Garfield::MediumMagboltz::GetElectronCollision(), Garfield::MediumSilicon::GetElectronCollision(), Garfield::Medium::GetElectronMomentum(), Garfield::MediumSilicon::GetElectronMomentum(), Garfield::MediumMagboltz::GetPhotonCollision(), Garfield::TrackBichsel::NewTrack(), Garfield::TrackElectron::NewTrack(), Garfield::TrackHeed::NewTrack(), Garfield::TrackPAI::NewTrack(), Garfield::TrackSimple::NewTrack(), RndmGaussian(), RndmLorentzian(), RndmUniformPos(), Garfield::TrackHeed::TransportDeltaElectron(), and Garfield::TrackHeed::TransportPhoton().

RndmUniformPos()

double Garfield::RndmUniformPos ( )

inline

Definition at line 19 of file Random.hh.

                               {
 
  double r = RndmUniform();
  while (r <= 0.) r = RndmUniform();
  return r;
}

Referenced by Garfield::MediumSilicon::ComputeSecondaries(), Garfield::TrackBichsel::GetCluster(), Garfield::TrackSimple::GetCluster(), Garfield::TrackElectron::GetCluster(), Garfield::TrackPAI::GetCluster(), Garfield::MediumMagboltz::GetElectronCollision(), and RndmPolya().

RndmVoigt()

double Garfield::RndmVoigt ( const double  mu, const double  sigma, const double  gamma  )

inline

Definition at line 67 of file Random.hh.

                                            {
 
  if (sigma <= 0.) return RndmLorentzian(mu, gamma);
  const double a = gamma / (Sqrt2 * sigma);
  const double x = RndmLorentzian(0., a) + RndmGaussian(0., 1. / Sqrt2);
  return mu + x * Sqrt2 * sigma;
}

Variable Documentation

plottingEngine

PlottingEngineRoot Garfield::plottingEngine

Definition at line 11 of file PlottingEngineRoot.cc.

Referenced by Garfield::ViewDrift::NewElectronDriftLine(), Garfield::ViewDrift::NewHoleDriftLine(), Garfield::ViewDrift::NewPhotonTrack(), Garfield::ViewSignal::PlotSignal(), Garfield::ViewCell::ViewCell(), Garfield::ViewFEMesh::ViewFEMesh(), Garfield::ViewField::ViewField(), Garfield::ViewGeometry::ViewGeometry(), Garfield::ViewMedium::ViewMedium(), and Garfield::ViewSignal::ViewSignal().

randomEngine

RandomEngineRoot Garfield::randomEngine

Definition at line 6 of file RandomEngineRoot.cc.

Referenced by main(), and RndmUniform().