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

#include <G4eBremParametrizedModel.hh>

+ Inheritance diagram for G4eBremParametrizedModel:

Public Member Functions

 G4eBremParametrizedModel (const G4ParticleDefinition *p=nullptr, const G4String &nam="eBremParam")
 
 ~G4eBremParametrizedModel () override
 
void Initialise (const G4ParticleDefinition *, const G4DataVector &) override
 
void InitialiseLocal (const G4ParticleDefinition *, G4VEmModel *masterModel) override
 
G4double MinEnergyCut (const G4ParticleDefinition *, const G4MaterialCutsCouple *) override
 
G4double ComputeDEDXPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy) override
 
G4double ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double tkin, G4double Z, G4double, G4double cutEnergy, G4double maxEnergy=DBL_MAX) override
 
void SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double cutEnergy, G4double maxEnergy) override
 
void SetupForMaterial (const G4ParticleDefinition *, const G4Material *, G4double) override
 
G4eBremParametrizedModel & operator= (const G4eBremParametrizedModel &right)=delete
 
 G4eBremParametrizedModel (const G4eBremParametrizedModel &)=delete
 
- Public Member Functions inherited from G4VEmModel
 G4VEmModel (const G4String &nam)
 
virtual ~G4VEmModel ()
 
virtual void InitialiseForMaterial (const G4ParticleDefinition *, const G4Material *)
 
virtual void InitialiseForElement (const G4ParticleDefinition *, G4int Z)
 
virtual G4double CrossSectionPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
virtual G4double GetPartialCrossSection (const G4Material *, G4int level, const G4ParticleDefinition *, G4double kineticEnergy)
 
virtual G4double ComputeCrossSectionPerShell (const G4ParticleDefinition *, G4int Z, G4int shellIdx, G4double kinEnergy, 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 *, const G4double &length, G4double &eloss)
 
virtual G4double Value (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy)
 
virtual G4double MinPrimaryEnergy (const G4Material *, const G4ParticleDefinition *, G4double cut=0.0)
 
virtual void DefineForRegion (const G4Region *)
 
virtual void FillNumberOfSecondaries (G4int &numberOfTriplets, G4int &numberOfRecoil)
 
virtual void ModelDescription (std::ostream &outFile) const
 
void InitialiseElementSelectors (const G4ParticleDefinition *, const G4DataVector &)
 
std::vector< G4EmElementSelector * > * GetElementSelectors ()
 
void SetElementSelectors (std::vector< G4EmElementSelector * > *)
 
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)
 
const G4Element * SelectRandomAtom (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
const G4Element * SelectTargetAtom (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double logKineticEnergy, 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)
 
const G4Element * GetCurrentElement (const G4Material *mat=nullptr) const
 
G4int SelectRandomAtomNumber (const G4Material *) const
 
const G4Isotope * GetCurrentIsotope (const G4Element *elm=nullptr) const
 
G4int SelectIsotopeNumber (const G4Element *) const
 
void SetParticleChange (G4VParticleChange *, G4VEmFluctuationModel *f=nullptr)
 
void SetCrossSectionTable (G4PhysicsTable *, G4bool isLocal)
 
G4ElementData * GetElementData ()
 
G4PhysicsTable * GetCrossSectionTable ()
 
G4VEmFluctuationModel * GetModelOfFluctuations ()
 
G4VEmAngularDistribution * GetAngularDistribution ()
 
G4VEmModel * GetTripletModel ()
 
void SetTripletModel (G4VEmModel *)
 
void SetAngularDistribution (G4VEmAngularDistribution *)
 
G4double HighEnergyLimit () const
 
G4double LowEnergyLimit () const
 
G4double HighEnergyActivationLimit () const
 
G4double LowEnergyActivationLimit () const
 
G4double PolarAngleLimit () const
 
G4double SecondaryThreshold () const
 
G4bool DeexcitationFlag () const
 
G4bool ForceBuildTableFlag () const
 
G4bool UseAngularGeneratorFlag () const
 
void SetAngularGeneratorFlag (G4bool)
 
void SetHighEnergyLimit (G4double)
 
void SetLowEnergyLimit (G4double)
 
void SetActivationHighEnergyLimit (G4double)
 
void SetActivationLowEnergyLimit (G4double)
 
G4bool IsActive (G4double kinEnergy) const
 
void SetPolarAngleLimit (G4double)
 
void SetSecondaryThreshold (G4double)
 
void SetDeexcitationFlag (G4bool val)
 
void SetForceBuildTable (G4bool val)
 
void SetFluctuationFlag (G4bool val)
 
void SetMasterThread (G4bool val)
 
G4bool IsMaster () const
 
void SetUseBaseMaterials (G4bool val)
 
G4bool UseBaseMaterials () const
 
G4double MaxSecondaryKinEnergy (const G4DynamicParticle *dynParticle)
 
const G4String & GetName () const
 
void SetCurrentCouple (const G4MaterialCutsCouple *)
 
G4bool IsLocked () const
 
void SetLocked (G4bool)
 
void SetLPMFlag (G4bool)
 
G4VEmModel & operator= (const G4VEmModel &right)=delete
 
 G4VEmModel (const G4VEmModel &)=delete
 

Protected Attributes

G4NistManager * nist
 
const G4ParticleDefinition * particle
 
G4ParticleDefinition * theGamma
 
G4ParticleChangeForLoss * fParticleChange
 
G4double minThreshold
 
G4double particleMass
 
G4double kinEnergy
 
G4double totalEnergy
 
G4double currentZ
 
G4double z13
 
G4double z23
 
G4double lnZ
 
G4double densityFactor
 
G4double densityCorr
 
G4double Fel
 
G4double Finel
 
G4double facFel
 
G4double facFinel
 
G4double fMax
 
G4double fCoulomb
 
G4bool isElectron
 
- Protected Attributes inherited from G4VEmModel
G4ElementData * fElementData = nullptr
 
G4VParticleChange * pParticleChange = nullptr
 
G4PhysicsTable * xSectionTable = nullptr
 
const G4Material * pBaseMaterial = nullptr
 
const std::vector< G4double > * theDensityFactor = nullptr
 
const std::vector< G4int > * theDensityIdx = nullptr
 
G4double inveplus
 
G4double pFactor = 1.0
 
std::size_t currentCoupleIndex = 0
 
std::size_t basedCoupleIndex = 0
 
G4bool lossFlucFlag = true
 

Static Protected Attributes

static const G4double xgi [8]
 
static const G4double wgi [8]
 

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 G4eBremParametrizedModel.hh.

Constructor & Destructor Documentation

◆ G4eBremParametrizedModel() [1/2]

G4eBremParametrizedModel::G4eBremParametrizedModel ( const G4ParticleDefinition * p = nullptr,
const G4String & nam = "eBremParam" )
explicit

Definition at line 98 of file G4eBremParametrizedModel.cc.

100 : G4VEmModel(nam),
101 particle(nullptr),
102 fMigdalConstant(classic_electr_radius*electron_Compton_length*electron_Compton_length*4.0*pi),
103 bremFactor(fine_structure_const*classic_electr_radius*classic_electr_radius*16./3.),
104 isInitialised(false),
105 isElectron(true)
106{
107 theGamma = G4Gamma::Gamma();
108
109 minThreshold = 0.1*keV;
110 lowKinEnergy = 10.*MeV;
111 SetLowEnergyLimit(lowKinEnergy);
112
113 nist = G4NistManager::Instance();
114
115 SetAngularDistribution(new G4ModifiedTsai());
116
117 particleMass = kinEnergy = totalEnergy = currentZ = z13 = z23 = lnZ = Fel = Finel
118 = densityFactor = densityCorr = fMax = fCoulomb = 0.;
119
120 InitialiseConstants();
121 if(nullptr != p) { SetParticle(p); }
122}
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition G4Gamma.cc:81
G4NistManager::Instance
static G4NistManager * Instance()
Definition G4NistManager.cc:59
G4VEmModel::SetLowEnergyLimit
void SetLowEnergyLimit(G4double)
Definition G4VEmModel.hh:745
G4VEmModel::G4VEmModel
G4VEmModel(const G4String &nam)
Definition G4VEmModel.cc:67
G4VEmModel::SetAngularDistribution
void SetAngularDistribution(G4VEmAngularDistribution *)
Definition G4VEmModel.hh:606
G4eBremParametrizedModel::particle
const G4ParticleDefinition * particle
Definition G4eBremParametrizedModel.hh:126

◆ ~G4eBremParametrizedModel()

G4eBremParametrizedModel::~G4eBremParametrizedModel ( )
overridedefault

◆ G4eBremParametrizedModel() [2/2]

G4eBremParametrizedModel::G4eBremParametrizedModel ( const G4eBremParametrizedModel & )
delete

Member Function Documentation

◆ ComputeCrossSectionPerAtom()

G4double G4eBremParametrizedModel::ComputeCrossSectionPerAtom ( const G4ParticleDefinition * p,
G4double tkin,
G4double Z,
G4double ,
G4double cutEnergy,
G4double maxEnergy = DBL_MAX )
overridevirtual

Reimplemented from G4VEmModel.

Definition at line 263 of file G4eBremParametrizedModel.cc.

269{
270 if(!particle) { SetParticle(p); }
271 if(kineticEnergy < lowKinEnergy) { return 0.0; }
272 G4double cut = std::min(cutEnergy, kineticEnergy);
273 G4double tmax = std::min(maxEnergy, kineticEnergy);
274
275 if(cut >= tmax) { return 0.0; }
276
277 SetCurrentElement(Z);
278
279 G4double cross = ComputeXSectionPerAtom(cut);
280
281 // allow partial integration
282 if(tmax < kinEnergy) { cross -= ComputeXSectionPerAtom(tmax); }
283
284 cross *= Z*Z*bremFactor;
285
286 return cross;
287}
G4double
double G4double
Definition G4Types.hh:83

◆ ComputeDEDXPerVolume()

G4double G4eBremParametrizedModel::ComputeDEDXPerVolume ( const G4Material * material,
const G4ParticleDefinition * p,
G4double kineticEnergy,
G4double cutEnergy )
overridevirtual

Reimplemented from G4VEmModel.

Definition at line 197 of file G4eBremParametrizedModel.cc.

202{
203 if(!particle) { SetParticle(p); }
204 if(kineticEnergy < lowKinEnergy) { return 0.0; }
205 G4double cut = std::min(cutEnergy, kineticEnergy);
206 if(cut == 0.0) { return 0.0; }
207
208 SetupForMaterial(particle, material,kineticEnergy);
209
210 const G4ElementVector* theElementVector = material->GetElementVector();
211 const G4double* theAtomicNumDensityVector = material->GetAtomicNumDensityVector();
212
213 G4double dedx = 0.0;
214
215 // loop for elements in the material
216 for (size_t i=0; i<material->GetNumberOfElements(); i++) {
217
218 G4VEmModel::SetCurrentElement((*theElementVector)[i]);
219 SetCurrentElement((*theElementVector)[i]->GetZ());
220
221 dedx += theAtomicNumDensityVector[i]*currentZ*currentZ*ComputeBremLoss(cut);
222 }
223 dedx *= bremFactor;
224
225 return dedx;
226}
G4ElementVector
std::vector< const G4Element * > G4ElementVector
Definition G4ElementVector.hh:33
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition G4Material.hh:183
G4Material::GetAtomicNumDensityVector
const G4double * GetAtomicNumDensityVector() const
Definition G4Material.hh:202
G4Material::GetNumberOfElements
std::size_t GetNumberOfElements() const
Definition G4Material.hh:180
G4VEmModel::SetCurrentElement
void SetCurrentElement(const G4Element *)
Definition G4VEmModel.hh:492
G4eBremParametrizedModel::SetupForMaterial
void SetupForMaterial(const G4ParticleDefinition *, const G4Material *, G4double) override
Definition G4eBremParametrizedModel.cc:156

◆ Initialise()

void G4eBremParametrizedModel::Initialise ( const G4ParticleDefinition * p,
const G4DataVector & cuts )
overridevirtual

Implements G4VEmModel.

Definition at line 171 of file G4eBremParametrizedModel.cc.

173{
174 if(p) { SetParticle(p); }
175
176 lowKinEnergy = LowEnergyLimit();
177
178 currentZ = 0.;
179
180 if(IsMaster()) { InitialiseElementSelectors(p, cuts); }
181
182 if(isInitialised) { return; }
183 fParticleChange = GetParticleChangeForLoss();
184 isInitialised = true;
185}
G4VEmModel::LowEnergyLimit
G4double LowEnergyLimit() const
Definition G4VEmModel.hh:640
G4VEmModel::IsMaster
G4bool IsMaster() const
Definition G4VEmModel.hh:717
G4VEmModel::InitialiseElementSelectors
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
Definition G4VEmModel.cc:135
G4VEmModel::GetParticleChangeForLoss
G4ParticleChangeForLoss * GetParticleChangeForLoss()
Definition G4VEmModel.cc:105
G4eBremParametrizedModel::fParticleChange
G4ParticleChangeForLoss * fParticleChange
Definition G4eBremParametrizedModel.hh:128

◆ InitialiseLocal()

void G4eBremParametrizedModel::InitialiseLocal ( const G4ParticleDefinition * ,
G4VEmModel * masterModel )
overridevirtual

Reimplemented from G4VEmModel.

Definition at line 189 of file G4eBremParametrizedModel.cc.

191{
192 SetElementSelectors(masterModel->GetElementSelectors());
193}
G4VEmModel::SetElementSelectors
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
Definition G4VEmModel.hh:816
G4VEmModel::GetElementSelectors
std::vector< G4EmElementSelector * > * GetElementSelectors()
Definition G4VEmModel.hh:808

◆ MinEnergyCut()

G4double G4eBremParametrizedModel::MinEnergyCut ( const G4ParticleDefinition * ,
const G4MaterialCutsCouple *  )
overridevirtual

Reimplemented from G4VEmModel.

Definition at line 148 of file G4eBremParametrizedModel.cc.

150{
151 return minThreshold;
152}

◆ operator=()

G4eBremParametrizedModel & G4eBremParametrizedModel::operator= ( const G4eBremParametrizedModel & right)
delete

◆ SampleSecondaries()

void G4eBremParametrizedModel::SampleSecondaries ( std::vector< G4DynamicParticle * > * vdp,
const G4MaterialCutsCouple * couple,
const G4DynamicParticle * dp,
G4double cutEnergy,
G4double maxEnergy )
overridevirtual

Implements G4VEmModel.

Definition at line 470 of file G4eBremParametrizedModel.cc.

476{
477 G4double kineticEnergy = dp->GetKineticEnergy();
478 if(kineticEnergy < lowKinEnergy) { return; }
479 G4double cut = std::min(cutEnergy, kineticEnergy);
480 G4double emax = std::min(maxEnergy, kineticEnergy);
481 if(cut >= emax) { return; }
482
483 SetupForMaterial(particle, couple->GetMaterial(),kineticEnergy);
484
485 const G4Element* elm = SelectTargetAtom(couple,particle,kineticEnergy,
486 dp->GetLogKineticEnergy(),cut,emax);
487 SetCurrentElement(elm->GetZ());
488
489 kinEnergy = kineticEnergy;
490 totalEnergy = kineticEnergy + particleMass;
491 densityCorr = densityFactor*totalEnergy*totalEnergy;
492
493 G4double xmin = G4Log(cut*cut + densityCorr);
494 G4double xmax = G4Log(emax*emax + densityCorr);
495 G4double gammaEnergy, f, x;
496
497 CLHEP::HepRandomEngine* rndmEngine = G4Random::getTheEngine();
498
499 do {
500 x = G4Exp(xmin + rndmEngine->flat()*(xmax - xmin)) - densityCorr;
501 if(x < 0.0) x = 0.0;
502 gammaEnergy = sqrt(x);
503 f = ComputeDXSectionPerAtom(gammaEnergy);
504
505 if ( f > fMax ) {
506 G4cout << "### G4eBremParametrizedModel Warning: Majoranta exceeded! "
507 << f << " > " << fMax
508 << " Egamma(MeV)= " << gammaEnergy
509 << " E(mEV)= " << kineticEnergy
510 << G4endl;
511 }
512
513 // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
514 } while (f < fMax*rndmEngine->flat());
515
516 //
517 // angles of the emitted gamma. ( Z - axis along the parent particle)
518 // use general interface
519 //
520 G4ThreeVector gammaDirection =
521 GetAngularDistribution()->SampleDirection(dp, totalEnergy-gammaEnergy,
522 G4lrint(currentZ),
523 couple->GetMaterial());
524
525 // create G4DynamicParticle object for the Gamma
526 auto gamma = new G4DynamicParticle(theGamma,gammaDirection, gammaEnergy);
527 vdp->push_back(gamma);
528
529 G4double totMomentum = sqrt(kineticEnergy*(totalEnergy + electron_mass_c2));
530 G4ThreeVector direction = (totMomentum*dp->GetMomentumDirection()
531 - gammaEnergy*gammaDirection).unit();
532
533 // energy of primary
534 G4double finalE = kineticEnergy - gammaEnergy;
535
536 // stop tracking and create new secondary instead of primary
537 if(gammaEnergy > SecondaryThreshold()) {
538 fParticleChange->ProposeTrackStatus(fStopAndKill);
539 fParticleChange->SetProposedKineticEnergy(0.0);
540 auto el =
541 new G4DynamicParticle(const_cast<G4ParticleDefinition*>(particle),
542 direction, finalE);
543 vdp->push_back(el);
544
545 // continue tracking
546 } else {
547 fParticleChange->SetProposedMomentumDirection(direction);
548 fParticleChange->SetProposedKineticEnergy(finalE);
549 }
550}
G4Exp
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition G4Exp.hh:180
G4Log
G4double G4Log(G4double x)
Definition G4Log.hh:227
fStopAndKill
@ fStopAndKill
Definition G4TrackStatus.hh:46
G4endl
#define G4endl
Definition G4ios.hh:67
G4cout
G4GLOB_DLL std::ostream G4cout
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::HepRandomEngine::flat
virtual double flat()=0
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetLogKineticEnergy
G4double GetLogKineticEnergy() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Element::GetZ
G4double GetZ() const
Definition G4Element.hh:119
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition G4MaterialCutsCouple.hh:166
G4ParticleChangeForLoss::SetProposedKineticEnergy
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition G4ParticleChangeForLoss.hh:108
G4ParticleChangeForLoss::SetProposedMomentumDirection
void SetProposedMomentumDirection(const G4ThreeVector &dir)
Definition G4ParticleChangeForLoss.hh:139
G4VEmAngularDistribution::SampleDirection
virtual G4ThreeVector & SampleDirection(const G4DynamicParticle *dp, G4double finalTotalEnergy, G4int Z, const G4Material *)=0
G4VEmModel::GetAngularDistribution
G4VEmAngularDistribution * GetAngularDistribution()
Definition G4VEmModel.hh:599
G4VEmModel::SelectTargetAtom
const G4Element * SelectTargetAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double logKineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition G4VEmModel.hh:576
G4VEmModel::SecondaryThreshold
G4double SecondaryThreshold() const
Definition G4VEmModel.hh:668
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4AblaRandom::flat
double flat()
Definition G4AblaRandom.cc:48
G4lrint
int G4lrint(double ad)
Definition templates.hh:134

◆ SetupForMaterial()

void G4eBremParametrizedModel::SetupForMaterial ( const G4ParticleDefinition * ,
const G4Material * mat,
G4double kineticEnergy )
overridevirtual

Reimplemented from G4VEmModel.

Definition at line 156 of file G4eBremParametrizedModel.cc.

159{
160 densityFactor = mat->GetElectronDensity()*fMigdalConstant;
161
162 // calculate threshold for density effect
163 kinEnergy = kineticEnergy;
164 totalEnergy = kineticEnergy + particleMass;
165 densityCorr = densityFactor*totalEnergy*totalEnergy;
166}
G4Material::GetElectronDensity
G4double GetElectronDensity() const
Definition G4Material.hh:203

Referenced by ComputeDEDXPerVolume(), and SampleSecondaries().

Member Data Documentation

◆ currentZ

G4double G4eBremParametrizedModel::currentZ
protected

Definition at line 138 of file G4eBremParametrizedModel.hh.

Referenced by ComputeDEDXPerVolume(), G4eBremParametrizedModel(), Initialise(), and SampleSecondaries().

◆ densityCorr

G4double G4eBremParametrizedModel::densityCorr
protected

Definition at line 141 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel(), SampleSecondaries(), and SetupForMaterial().

◆ densityFactor

G4double G4eBremParametrizedModel::densityFactor
protected

Definition at line 140 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel(), SampleSecondaries(), and SetupForMaterial().

◆ facFel

G4double G4eBremParametrizedModel::facFel
protected

Definition at line 143 of file G4eBremParametrizedModel.hh.

◆ facFinel

G4double G4eBremParametrizedModel::facFinel
protected

Definition at line 143 of file G4eBremParametrizedModel.hh.

◆ fCoulomb

G4double G4eBremParametrizedModel::fCoulomb
protected

Definition at line 144 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel().

◆ Fel

G4double G4eBremParametrizedModel::Fel
protected

Definition at line 142 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel().

◆ Finel

G4double G4eBremParametrizedModel::Finel
protected

Definition at line 142 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel().

◆ fMax

G4double G4eBremParametrizedModel::fMax
protected

Definition at line 144 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel(), and SampleSecondaries().

◆ fParticleChange

G4ParticleChangeForLoss* G4eBremParametrizedModel::fParticleChange
protected

Definition at line 128 of file G4eBremParametrizedModel.hh.

Referenced by Initialise(), and SampleSecondaries().

◆ isElectron

G4bool G4eBremParametrizedModel::isElectron
protected

Definition at line 156 of file G4eBremParametrizedModel.hh.

◆ kinEnergy

G4double G4eBremParametrizedModel::kinEnergy
protected

Definition at line 136 of file G4eBremParametrizedModel.hh.

Referenced by ComputeCrossSectionPerAtom(), G4eBremParametrizedModel(), SampleSecondaries(), and SetupForMaterial().

◆ lnZ

G4double G4eBremParametrizedModel::lnZ
protected

Definition at line 139 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel().

◆ minThreshold

G4double G4eBremParametrizedModel::minThreshold
protected

Definition at line 132 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel(), and MinEnergyCut().

◆ nist

G4NistManager* G4eBremParametrizedModel::nist
protected

Definition at line 125 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel().

◆ particle

const G4ParticleDefinition* G4eBremParametrizedModel::particle
protected

Definition at line 126 of file G4eBremParametrizedModel.hh.

Referenced by ComputeCrossSectionPerAtom(), ComputeDEDXPerVolume(), and SampleSecondaries().

◆ particleMass

G4double G4eBremParametrizedModel::particleMass
protected

Definition at line 135 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel(), SampleSecondaries(), and SetupForMaterial().

◆ theGamma

G4ParticleDefinition* G4eBremParametrizedModel::theGamma
protected

Definition at line 127 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel(), and SampleSecondaries().

◆ totalEnergy

G4double G4eBremParametrizedModel::totalEnergy
protected

Definition at line 137 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel(), SampleSecondaries(), and SetupForMaterial().

◆ wgi

const G4double G4eBremParametrizedModel::wgi
staticprotected
Initial value:
={ 0.0506, 0.1112, 0.1569, 0.1813,
0.1813, 0.1569, 0.1112, 0.0506 }

Definition at line 68 of file G4eBremParametrizedModel.hh.

100 :
101
102 void InitialiseConstants();
103
104 G4double ComputeBremLoss(G4double cutEnergy);
105
106 G4double ComputeXSectionPerAtom(G4double cutEnergy);
107
108 G4double ComputeDXSectionPerAtom(G4double gammaEnergy);
109
110 G4double ComputeParametrizedDXSectionPerAtom(G4double kineticEnergy,
111 G4double gammaEnergy,
112 G4double Z);
113
114 void SetParticle(const G4ParticleDefinition* p);
115
116 G4double ScreenFunction1(G4double ScreenVariable);
117
118 G4double ScreenFunction2(G4double ScreenVariable);
119
120 inline void SetCurrentElement(const G4double);
121
122protected:
123
124 G4NistManager* nist;
125 const G4ParticleDefinition* particle;
126 G4ParticleDefinition* theGamma;
127 G4ParticleChangeForLoss* fParticleChange;
128
129 static const G4double xgi[8], wgi[8];
130
131 G4double minThreshold;
132
133 // cash
134 G4double particleMass;
135 G4double kinEnergy;
136 G4double totalEnergy;
137 G4double currentZ;
138 G4double z13, z23, lnZ;
139 G4double densityFactor;
140 G4double densityCorr;
141 G4double Fel, Finel;
142 G4double facFel, facFinel;
143 G4double fMax,fCoulomb;
144
145private:
146
147 G4double lowKinEnergy;
148 G4double fMigdalConstant;
149 G4double bremFactor;
150
151 G4bool isInitialised;
152
153protected:
154
155 G4bool isElectron;
156
157};
158
159//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
160
161inline void G4eBremParametrizedModel::SetCurrentElement(const G4double Z)
162{
163 if(Z != currentZ) {
164 currentZ = Z;
165
166 G4int iz = G4lrint(Z);
167 z13 = nist->GetZ13(iz);
168 z23 = z13*z13;
169 lnZ = nist->GetLOGZ(iz);
170
171 Fel = facFel - lnZ/3. ;
172 Finel = facFinel - 2.*lnZ/3. ;
173
174 fCoulomb = GetCurrentElement()->GetfCoulomb();
175 fMax = Fel-fCoulomb + Finel/currentZ + (1.+1./currentZ)/12.;
176 }
177}
178
179//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
180
181
182#endif
G4bool
bool G4bool
Definition G4Types.hh:86
G4int
int G4int
Definition G4Types.hh:85
G4Element::GetfCoulomb
G4double GetfCoulomb() const
Definition G4Element.hh:165
G4NistManager::GetZ13
G4double GetZ13(G4double Z) const
Definition G4NistManager.hh:490
G4NistManager::GetLOGZ
G4double GetLOGZ(G4int Z) const
Definition G4NistManager.hh:505
G4VEmModel::GetCurrentElement
const G4Element * GetCurrentElement(const G4Material *mat=nullptr) const
Definition G4VEmModel.cc:239

◆ xgi

const G4double G4eBremParametrizedModel::xgi
staticprotected
Initial value:
={ 0.0199, 0.1017, 0.2372, 0.4083,
0.5917, 0.7628, 0.8983, 0.9801 }

Definition at line 66 of file G4eBremParametrizedModel.hh.

◆ z13

G4double G4eBremParametrizedModel::z13
protected

Definition at line 139 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel().

◆ z23

G4double G4eBremParametrizedModel::z23
protected

Definition at line 139 of file G4eBremParametrizedModel.hh.

Referenced by G4eBremParametrizedModel().

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