Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4RToEConvForElectron Class Reference

#include <G4RToEConvForElectron.hh>

Inheritance diagram for G4RToEConvForElectron:
Collaboration diagram for G4RToEConvForElectron:

Public Member Functions

 G4RToEConvForElectron ()
 
virtual ~G4RToEConvForElectron ()
 
- Public Member Functions inherited from G4VRangeToEnergyConverter
 G4VRangeToEnergyConverter ()
 
 G4VRangeToEnergyConverter (const G4VRangeToEnergyConverter &right)
 
G4VRangeToEnergyConverteroperator= (const G4VRangeToEnergyConverter &right)
 
virtual ~G4VRangeToEnergyConverter ()
 
G4int operator== (const G4VRangeToEnergyConverter &right) const
 
G4int operator!= (const G4VRangeToEnergyConverter &right) const
 
virtual G4double Convert (G4double rangeCut, const G4Material *material)
 
const G4ParticleDefinitionGetParticleType () const
 
const G4PhysicsTableGetLossTable () const
 
virtual void Reset ()
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 

Protected Member Functions

virtual G4double ComputeLoss (G4double AtomicNumber, G4double KineticEnergy)
 
- Protected Member Functions inherited from G4VRangeToEnergyConverter
virtual void BuildLossTable ()
 
virtual void BuildRangeVector (const G4Material *aMaterial, G4RangeVector *rangeVector)
 
G4double ConvertCutToKineticEnergy (G4RangeVector *theRangeVector, G4double theCutInLength, size_t materialIndex) const
 

Protected Attributes

G4double Mass
 
G4double Z
 
G4double taul
 
G4double ionpot
 
G4double ionpotlog
 
G4double bremfactor
 
- Protected Attributes inherited from G4VRangeToEnergyConverter
G4double fMaxEnergyCut
 
const G4ParticleDefinitiontheParticle
 
G4LossTabletheLossTable
 
G4int NumberOfElements
 
const G4int TotBin
 
std::vector< G4RangeVector * > fRangeVectorStore
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VRangeToEnergyConverter
static void SetEnergyRange (G4double lowedge, G4double highedge)
 
static G4double GetLowEdgeEnergy ()
 
static G4double GetHighEdgeEnergy ()
 
static G4double GetMaxEnergyCut ()
 
static void SetMaxEnergyCut (G4double value)
 
- Protected Types inherited from G4VRangeToEnergyConverter
typedef G4PhysicsTable G4LossTable
 
typedef G4PhysicsLogVector G4LossVector
 
typedef G4PhysicsLogVector G4RangeVector
 
- Static Protected Attributes inherited from G4VRangeToEnergyConverter
static G4double LowestEnergy = 0.99e-3*MeV
 
static G4double HighestEnergy = 100.0e6*MeV
 
static G4double MaxEnergyCut = 10.0*GeV
 

Detailed Description

Definition at line 51 of file G4RToEConvForElectron.hh.

Constructor & Destructor Documentation

G4RToEConvForElectron::G4RToEConvForElectron ( )

Definition at line 45 of file G4RToEConvForElectron.cc.

47  Mass(0.0),
48  Z(-1.),
49  taul(0.0),
50  ionpot(0.0),
51  ionpotlog(-1.0e-10),
52  bremfactor(0.1)
53 {
55  if (theParticle ==0) {
56 #ifdef G4VERBOSE
57  if (GetVerboseLevel()>0) {
58  G4cout << " G4RToEConvForElectron::G4RToEConvForElectron() ";
59  G4cout << " Electron is not defined !!" << G4endl;
60  }
61 #endif
62  } else {
64  }
65 }
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
G4GLOB_DLL std::ostream G4cout
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
#define G4endl
Definition: G4ios.hh:61
const G4ParticleDefinition * theParticle

Here is the call graph for this function:

G4RToEConvForElectron::~G4RToEConvForElectron ( )
virtual

Definition at line 67 of file G4RToEConvForElectron.cc.

68 {
69 }

Member Function Documentation

G4double G4RToEConvForElectron::ComputeLoss ( G4double  AtomicNumber,
G4double  KineticEnergy 
)
protectedvirtual

Implements G4VRangeToEnergyConverter.

Definition at line 75 of file G4RToEConvForElectron.cc.

77 {
78  const G4double cbr1=0.02, cbr2=-5.7e-5, cbr3=1., cbr4=0.072;
79  const G4double Tlow=10.*keV, Thigh=1.*GeV;
80 
81  // calculate dE/dx for electrons
82  if( std::fabs(AtomicNumber-Z)>0.1 ) {
83  Z = AtomicNumber;
84  taul = Tlow/Mass;
85  ionpot = 1.6e-5*MeV*std::exp(0.9*std::log(Z))/Mass;
86  ionpotlog = std::log(ionpot);
87  }
88 
89 
90  G4double tau = KineticEnergy/Mass;
91  G4double dEdx;
92 
93  if(tau<taul) {
94  G4double t1 = taul+1.;
95  G4double t2 = taul+2.;
96  G4double tsq = taul*taul;
97  G4double beta2 = taul*t2/(t1*t1);
98  G4double f = 1.-beta2+std::log(tsq/2.)
99  +(0.5+0.25*tsq+(1.+2.*taul)*std::log(0.5))/(t1*t1);
100  dEdx = (std::log(2.*taul+4.)-2.*ionpotlog+f)/beta2;
101  dEdx = twopi_mc2_rcl2*Z*dEdx;
102  G4double clow = dEdx*std::sqrt(taul);
103  dEdx = clow/std::sqrt(KineticEnergy/Mass);
104 
105  } else {
106  G4double t1 = tau+1.;
107  G4double t2 = tau+2.;
108  G4double tsq = tau*tau;
109  G4double beta2 = tau*t2/(t1*t1);
110  G4double f = 1.-beta2+std::log(tsq/2.)
111  +(0.5+0.25*tsq+(1.+2.*tau)*std::log(0.5))/(t1*t1);
112  dEdx = (std::log(2.*tau+4.)-2.*ionpotlog+f)/beta2;
113  dEdx = twopi_mc2_rcl2*Z*dEdx;
114 
115  // loss from bremsstrahlung follows
116  G4double cbrem = (cbr1+cbr2*Z)
117  *(cbr3+cbr4*std::log(KineticEnergy/Thigh));
118  cbrem = Z*(Z+1.)*cbrem*tau/beta2;
119 
120  cbrem *= bremfactor ;
121 
122  dEdx += twopi_mc2_rcl2*cbrem;
123  }
124 
125  return dEdx;
126 }
tuple t1
Definition: plottest35.py:33
static constexpr double GeV
Definition: G4SIunits.hh:217
static constexpr double MeV
Definition: G4SIunits.hh:214
double G4double
Definition: G4Types.hh:76
static constexpr double keV
Definition: G4SIunits.hh:216

Member Data Documentation

G4double G4RToEConvForElectron::bremfactor
protected

Definition at line 72 of file G4RToEConvForElectron.hh.

G4double G4RToEConvForElectron::ionpot
protected

Definition at line 70 of file G4RToEConvForElectron.hh.

G4double G4RToEConvForElectron::ionpotlog
protected

Definition at line 71 of file G4RToEConvForElectron.hh.

G4double G4RToEConvForElectron::Mass
protected

Definition at line 67 of file G4RToEConvForElectron.hh.

G4double G4RToEConvForElectron::taul
protected

Definition at line 69 of file G4RToEConvForElectron.hh.

G4double G4RToEConvForElectron::Z
protected

Definition at line 68 of file G4RToEConvForElectron.hh.


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