Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | Protected Attributes | List of all members
G4HeatedKleinNishinaCompton Class Reference

#include <G4HeatedKleinNishinaCompton.hh>

+ Inheritance diagram for G4HeatedKleinNishinaCompton:

Public Member Functions

 G4HeatedKleinNishinaCompton (const G4ParticleDefinition *p=0, const G4String &nam="Heated-Klein-Nishina")
 
virtual ~G4HeatedKleinNishinaCompton ()
 
virtual void Initialise (const G4ParticleDefinition *, const G4DataVector &)
 
virtual G4double ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double cut, G4double emax)
 
virtual void SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
 
void SetElectronTemperature (G4double t)
 
G4double GetElectronTemperature ()
 
- Public Member Functions inherited from G4VEmModel
 G4VEmModel (const G4String &nam)
 
virtual ~G4VEmModel ()
 
virtual void Initialise (const G4ParticleDefinition *, const G4DataVector &)=0
 
virtual void SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double tmax=DBL_MAX)=0
 
virtual G4double ComputeDEDXPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
 
virtual G4double CrossSectionPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
virtual G4double ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
virtual G4double ChargeSquareRatio (const G4Track &)
 
virtual G4double GetChargeSquareRatio (const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
 
virtual G4double GetParticleCharge (const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
 
virtual void StartTracking (G4Track *)
 
virtual void CorrectionsAlongStep (const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
 
virtual G4double Value (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy)
 
virtual G4double MinPrimaryEnergy (const G4Material *, const G4ParticleDefinition *)
 
virtual void SetupForMaterial (const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
 
virtual void DefineForRegion (const G4Region *)
 
void InitialiseElementSelectors (const G4ParticleDefinition *, const G4DataVector &)
 
G4double ComputeDEDX (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
 
G4double CrossSection (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
G4double ComputeMeanFreePath (const G4ParticleDefinition *, G4double kineticEnergy, const G4Material *, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
G4double ComputeCrossSectionPerAtom (const G4ParticleDefinition *, const G4Element *, G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
G4int SelectIsotopeNumber (const G4Element *)
 
const G4Element * SelectRandomAtom (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
const G4Element * SelectRandomAtom (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
void SetParticleChange (G4VParticleChange *, G4VEmFluctuationModel *f=0)
 
void SetCrossSectionTable (G4PhysicsTable *)
 
G4PhysicsTable * GetCrossSectionTable ()
 
G4VEmFluctuationModel * GetModelOfFluctuations ()
 
G4VEmAngularDistribution * GetAngularDistribution ()
 
void SetAngularDistribution (G4VEmAngularDistribution *)
 
G4double HighEnergyLimit () const
 
G4double LowEnergyLimit () const
 
G4double HighEnergyActivationLimit () const
 
G4double LowEnergyActivationLimit () const
 
G4double PolarAngleLimit () const
 
G4double SecondaryThreshold () const
 
G4bool LPMFlag () const
 
G4bool DeexcitationFlag () const
 
G4bool ForceBuildTableFlag () const
 
void SetHighEnergyLimit (G4double)
 
void SetLowEnergyLimit (G4double)
 
void SetActivationHighEnergyLimit (G4double)
 
void SetActivationLowEnergyLimit (G4double)
 
G4bool IsActive (G4double kinEnergy)
 
void SetPolarAngleLimit (G4double)
 
void SetSecondaryThreshold (G4double)
 
void SetLPMFlag (G4bool val)
 
void SetDeexcitationFlag (G4bool val)
 
void ForceBuildTable (G4bool val)
 
G4double MaxSecondaryKinEnergy (const G4DynamicParticle *dynParticle)
 
const G4String & GetName () const
 
void SetCurrentCouple (const G4MaterialCutsCouple *)
 
const G4Element * GetCurrentElement () const
 

Protected Attributes

G4ParticleDefinition * theGamma
 
G4ParticleDefinition * theElectron
 
G4ParticleChangeForGamma * fParticleChange
 
G4double lowestGammaEnergy
 
G4double fTemperature
 
- Protected Attributes inherited from G4VEmModel
G4VParticleChange * pParticleChange
 
G4PhysicsTable * xSectionTable
 
const std::vector< G4double > * theDensityFactor
 
const std::vector< G4int > * theDensityIdx
 

Additional Inherited Members

- Protected Member Functions inherited from G4VEmModel
G4ParticleChangeForLoss * GetParticleChangeForLoss ()
 
G4ParticleChangeForGamma * GetParticleChangeForGamma ()
 
virtual G4double MaxSecondaryEnergy (const G4ParticleDefinition *, G4double kineticEnergy)
 
const G4MaterialCutsCouple * CurrentCouple () const
 
void SetCurrentElement (const G4Element *)
 

Detailed Description

Definition at line 59 of file G4HeatedKleinNishinaCompton.hh.

Constructor & Destructor Documentation

◆ G4HeatedKleinNishinaCompton()

G4HeatedKleinNishinaCompton::G4HeatedKleinNishinaCompton ( const G4ParticleDefinition *  p = 0,
const G4String &  nam = "Heated-Klein-Nishina" 
)

Definition at line 67 of file G4HeatedKleinNishinaCompton.cc.

69 : G4VEmModel(nam)
70{
71 theGamma = G4Gamma::Gamma();
72 theElectron = G4Electron::Electron();
73 lowestGammaEnergy = 1.0*eV;
74 fTemperature = 1.0*keV;
75 fParticleChange = 0;
76}
G4Electron::Electron
static G4Electron * Electron()
Definition: G4Electron.cc:94
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:86
G4HeatedKleinNishinaCompton::fParticleChange
G4ParticleChangeForGamma * fParticleChange
Definition: G4HeatedKleinNishinaCompton.hh:92

◆ ~G4HeatedKleinNishinaCompton()

G4HeatedKleinNishinaCompton::~G4HeatedKleinNishinaCompton ( )
virtual

Definition at line 80 of file G4HeatedKleinNishinaCompton.cc.

81{}

Member Function Documentation

◆ ComputeCrossSectionPerAtom()

G4double G4HeatedKleinNishinaCompton::ComputeCrossSectionPerAtom ( const G4ParticleDefinition *  ,
G4double  kinEnergy,
G4double  Z,
G4double  A,
G4double  cut,
G4double  emax 
)
virtual

Reimplemented from G4VEmModel.

Definition at line 95 of file G4HeatedKleinNishinaCompton.cc.

100{
101 G4double xSection = 0.0 ;
102 if ( Z < 0.9999 ) return xSection;
103 if ( GammaEnergy < 0.01*keV ) return xSection;
104 // if ( GammaEnergy > (100.*GeV/Z) ) return xSection;
105
106 static const G4double a = 20.0 , b = 230.0 , c = 440.0;
107
108 static const G4double
109 d1= 2.7965e-1*barn, d2=-1.8300e-1*barn, d3= 6.7527 *barn, d4=-1.9798e+1*barn,
110 e1= 1.9756e-5*barn, e2=-1.0205e-2*barn, e3=-7.3913e-2*barn, e4= 2.7079e-2*barn,
111 f1=-3.9178e-7*barn, f2= 6.8241e-5*barn, f3= 6.0480e-5*barn, f4= 3.0274e-4*barn;
112
113 G4double p1Z = Z*(d1 + e1*Z + f1*Z*Z), p2Z = Z*(d2 + e2*Z + f2*Z*Z),
114 p3Z = Z*(d3 + e3*Z + f3*Z*Z), p4Z = Z*(d4 + e4*Z + f4*Z*Z);
115
116 G4double T0 = 15.0*keV;
117 if (Z < 1.5) T0 = 40.0*keV;
118
119 G4double X = max(GammaEnergy, T0) / electron_mass_c2;
120 xSection = p1Z*std::log(1.+2.*X)/X
121 + (p2Z + p3Z*X + p4Z*X*X)/(1. + a*X + b*X*X + c*X*X*X);
122
123 // modification for low energy. (special case for Hydrogen)
124 if (GammaEnergy < T0) {
125 G4double dT0 = 1.*keV;
126 X = (T0+dT0) / electron_mass_c2 ;
127 G4double sigma = p1Z*log(1.+2*X)/X
128 + (p2Z + p3Z*X + p4Z*X*X)/(1. + a*X + b*X*X + c*X*X*X);
129 G4double c1 = -T0*(sigma-xSection)/(xSection*dT0);
130 G4double c2 = 0.150;
131 if (Z > 1.5) c2 = 0.375-0.0556*log(Z);
132 G4double y = log(GammaEnergy/T0);
133 xSection *= exp(-y*(c1+c2*y));
134 }
135 // G4cout << "e= " << GammaEnergy << " Z= " << Z << " cross= " << xSection << G4endl;
136 return xSection;
137}
G4double
double G4double
Definition: G4Types.hh:64

◆ GetElectronTemperature()

G4double G4HeatedKleinNishinaCompton::GetElectronTemperature ( )
inline

Definition at line 86 of file G4HeatedKleinNishinaCompton.hh.

86{return fTemperature;};

◆ Initialise()

void G4HeatedKleinNishinaCompton::Initialise ( const G4ParticleDefinition *  ,
const G4DataVector &   
)
virtual

Implements G4VEmModel.

Definition at line 85 of file G4HeatedKleinNishinaCompton.cc.

87{
88 if(!fParticleChange) fParticleChange = GetParticleChangeForGamma();
89}
G4VEmModel::GetParticleChangeForGamma
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition: G4VEmModel.cc:109

◆ SampleSecondaries()

void G4HeatedKleinNishinaCompton::SampleSecondaries ( std::vector< G4DynamicParticle * > *  fvect,
const G4MaterialCutsCouple *  ,
const G4DynamicParticle *  aDynamicGamma,
G4double  tmin,
G4double  maxEnergy 
)
virtual

Implements G4VEmModel.

Definition at line 143 of file G4HeatedKleinNishinaCompton.cc.

148{
149 // The scattered gamma energy is sampled according to Klein - Nishina formula.
150 // The random number techniques of Butcher & Messel are used
151 // (Nuc Phys 20(1960),15).
152 // Note : Effects due to binding of atomic electrons are negliged.
153
154 // We start to prepare a heated electron from Maxwell distribution.
155 // Then we try to boost to the electron rest frame and make scattering.
156 // The final step is to recover new gamma 4momentum in the lab frame
157
158 G4double eMomentumC2 = CLHEP::RandGamma::shoot(1.5,1.);
159 eMomentumC2 *= 2*electron_mass_c2*fTemperature; // electron (pc)^2
160 G4ThreeVector eMomDir = G4RandomDirection();
161 eMomDir *= std::sqrt(eMomentumC2);
162 G4double eEnergy = std::sqrt(eMomentumC2+electron_mass_c2*electron_mass_c2);
163 G4LorentzVector electron4v = G4LorentzVector(eMomDir,eEnergy);
164 G4ThreeVector bst = electron4v.boostVector();
165
166 G4LorentzVector gamma4v = aDynamicGamma->Get4Momentum();
167 gamma4v.boost(-bst);
168
169 G4ThreeVector gammaMomV = gamma4v.vect();
170 G4double gamEnergy0 = gammaMomV.mag();
171
172
173 // G4double gamEnergy0 = aDynamicGamma->GetKineticEnergy();
174 G4double E0_m = gamEnergy0 / electron_mass_c2 ;
175
176 // G4ThreeVector gamDirection0 = /aDynamicGamma->GetMomentumDirection();
177
178 G4ThreeVector gamDirection0 = gammaMomV/gamEnergy0;
179
180 // sample the energy rate of the scattered gamma in the electron rest frame
181 //
182
183 G4double epsilon, epsilonsq, onecost, sint2, greject ;
184
185 G4double eps0 = 1./(1. + 2.*E0_m);
186 G4double epsilon0sq = eps0*eps0;
187 G4double alpha1 = - log(eps0);
188 G4double alpha2 = 0.5*(1.- epsilon0sq);
189
190 do
191 {
192 if ( alpha1/(alpha1+alpha2) > G4UniformRand() )
193 {
194 epsilon = exp(-alpha1*G4UniformRand()); // eps0**r
195 epsilonsq = epsilon*epsilon;
196
197 }
198 else
199 {
200 epsilonsq = epsilon0sq + (1.- epsilon0sq)*G4UniformRand();
201 epsilon = sqrt(epsilonsq);
202 };
203
204 onecost = (1.- epsilon)/(epsilon*E0_m);
205 sint2 = onecost*(2.-onecost);
206 greject = 1. - epsilon*sint2/(1.+ epsilonsq);
207
208 } while (greject < G4UniformRand());
209
210 //
211 // scattered gamma angles. ( Z - axis along the parent gamma)
212 //
213
214 G4double cosTeta = 1. - onecost;
215 G4double sinTeta = sqrt (sint2);
216 G4double Phi = twopi * G4UniformRand();
217 G4double dirx = sinTeta*cos(Phi), diry = sinTeta*sin(Phi), dirz = cosTeta;
218
219 //
220 // update G4VParticleChange for the scattered gamma
221 //
222
223 G4ThreeVector gamDirection1 ( dirx,diry,dirz );
224 gamDirection1.rotateUz(gamDirection0);
225 G4double gamEnergy1 = epsilon*gamEnergy0;
226 gamDirection1 *= gamEnergy1;
227
228 G4LorentzVector gamma4vfinal = G4LorentzVector(gamDirection1,gamEnergy1);
229
230
231 // kinematic of the scattered electron
232 //
233
234 G4double eKinEnergy = gamEnergy0 - gamEnergy1;
235 G4ThreeVector eDirection = gamEnergy0*gamDirection0 - gamEnergy1*gamDirection1;
236 eDirection = eDirection.unit();
237 G4double eFinalMom = std::sqrt(eKinEnergy*(eKinEnergy+2*electron_mass_c2));
238 eDirection *= eFinalMom;
239 G4LorentzVector e4vfinal = G4LorentzVector(eDirection,gamEnergy1+electron_mass_c2);
240
241 gamma4vfinal.boost(bst);
242 e4vfinal.boost(bst);
243
244 gamDirection1 = gamma4vfinal.vect();
245 gamEnergy1 = gamDirection1.mag();
246 gamDirection1 /= gamEnergy1;
247
248
249
250
251 fParticleChange->SetProposedKineticEnergy(gamEnergy1);
252
253 if( gamEnergy1 > lowestGammaEnergy )
254 {
255 gamDirection1 /= gamEnergy1;
256 fParticleChange->ProposeMomentumDirection(gamDirection1);
257 }
258 else
259 {
260 fParticleChange->ProposeTrackStatus(fStopAndKill);
261 gamEnergy1 += fParticleChange->GetLocalEnergyDeposit();
262 fParticleChange->ProposeLocalEnergyDeposit(gamEnergy1);
263 }
264
265 eKinEnergy = e4vfinal.t()-electron_mass_c2;
266
267 if( eKinEnergy > DBL_MIN )
268 {
269 // create G4DynamicParticle object for the electron.
270 eDirection = e4vfinal.vect();
271 G4double eFinMomMag = eDirection.mag();
272 eDirection /= eFinMomMag;
273 G4DynamicParticle* dp = new G4DynamicParticle(theElectron,eDirection,eKinEnergy);
274 fvect->push_back(dp);
275 }
276}
G4LorentzVector
CLHEP::HepLorentzVector G4LorentzVector
Definition: G4LorentzVector.hh:35
G4RandomDirection
G4ThreeVector G4RandomDirection()
Definition: G4RandomDirection.hh:55
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:56
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:53
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::Hep3Vector::mag
double mag() const
CLHEP::HepLorentzVector::boostVector
Hep3Vector boostVector() const
Definition: LorentzVector.cc:177
CLHEP::HepLorentzVector::boost
HepLorentzVector & boost(double, double, double)
Definition: LorentzVector.cc:59
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
CLHEP::RandGamma::shoot
static double shoot()
G4DynamicParticle::Get4Momentum
G4LorentzVector Get4Momentum() const
G4ParticleChangeForGamma::SetProposedKineticEnergy
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition: G4ParticleChangeForGamma.hh:129
G4ParticleChangeForGamma::ProposeMomentumDirection
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
Definition: G4ParticleChangeForGamma.hh:147
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4VParticleChange::GetLocalEnergyDeposit
G4double GetLocalEnergyDeposit() const
G4VParticleChange::ProposeLocalEnergyDeposit
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
DBL_MIN
#define DBL_MIN
Definition: templates.hh:75

◆ SetElectronTemperature()

void G4HeatedKleinNishinaCompton::SetElectronTemperature ( G4double  t)
inline

Definition at line 85 of file G4HeatedKleinNishinaCompton.hh.

85{fTemperature=t;};

Member Data Documentation

◆ fParticleChange

G4ParticleChangeForGamma* G4HeatedKleinNishinaCompton::fParticleChange
protected

Definition at line 92 of file G4HeatedKleinNishinaCompton.hh.

Referenced by G4HeatedKleinNishinaCompton(), Initialise(), and SampleSecondaries().

◆ fTemperature

G4double G4HeatedKleinNishinaCompton::fTemperature
protected

Definition at line 94 of file G4HeatedKleinNishinaCompton.hh.

Referenced by G4HeatedKleinNishinaCompton(), GetElectronTemperature(), SampleSecondaries(), and SetElectronTemperature().

◆ lowestGammaEnergy

G4double G4HeatedKleinNishinaCompton::lowestGammaEnergy
protected

Definition at line 93 of file G4HeatedKleinNishinaCompton.hh.

Referenced by G4HeatedKleinNishinaCompton(), and SampleSecondaries().

◆ theElectron

G4ParticleDefinition* G4HeatedKleinNishinaCompton::theElectron
protected

Definition at line 91 of file G4HeatedKleinNishinaCompton.hh.

Referenced by G4HeatedKleinNishinaCompton(), and SampleSecondaries().

◆ theGamma

G4ParticleDefinition* G4HeatedKleinNishinaCompton::theGamma
protected

Definition at line 90 of file G4HeatedKleinNishinaCompton.hh.

Referenced by G4HeatedKleinNishinaCompton().

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