Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4VAtomDeexcitation Class Referenceabstract

#include <G4VAtomDeexcitation.hh>

+ Inheritance diagram for G4VAtomDeexcitation:

Public Member Functions

 G4VAtomDeexcitation (const G4String &modname="Deexcitation")
 
virtual ~G4VAtomDeexcitation ()
 
void InitialiseAtomicDeexcitation ()
 
virtual void InitialiseForNewRun ()=0
 
virtual void InitialiseForExtraAtom (G4int Z)=0
 
void SetDeexcitationActiveRegion (const G4String &rname, G4bool valDeexcitation, G4bool valAuger, G4bool valPIXE)
 
void SetFluo (G4bool)
 
G4bool IsFluoActive () const
 
void SetAuger (G4bool)
 
G4bool IsAugerActive () const
 
void SetAugerCascade (G4bool)
 
G4bool IsAugerCascadeActive () const
 
void SetPIXE (G4bool)
 
G4bool IsPIXEActive () const
 
const G4String & GetName () const
 
const std::vector< G4bool > & GetListOfActiveAtoms () const
 
void SetVerboseLevel (G4int)
 
G4int GetVerboseLevel () const
 
G4bool CheckDeexcitationActiveRegion (G4int coupleIndex)
 
G4bool CheckAugerActiveRegion (G4int coupleIndex)
 
virtual const G4AtomicShell * GetAtomicShell (G4int Z, G4AtomicShellEnumerator shell)=0
 
void GenerateParticles (std::vector< G4DynamicParticle * > *secVect, const G4AtomicShell *, G4int Z, G4int coupleIndex)
 
virtual void GenerateParticles (std::vector< G4DynamicParticle * > *secVect, const G4AtomicShell *, G4int Z, G4double gammaCut, G4double eCut)=0
 
virtual G4double GetShellIonisationCrossSectionPerAtom (const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=nullptr)=0
 
virtual G4double ComputeShellIonisationCrossSectionPerAtom (const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=nullptr)=0
 
void AlongStepDeexcitation (std::vector< G4Track * > &tracks, const G4Step &step, G4double &eLoss, G4int coupleIndex)
 
 G4VAtomDeexcitation (G4VAtomDeexcitation &)=delete
 
G4VAtomDeexcitation & operator= (const G4VAtomDeexcitation &right)=delete
 

Detailed Description

Definition at line 64 of file G4VAtomDeexcitation.hh.

Constructor & Destructor Documentation

◆ G4VAtomDeexcitation() [1/2]

G4VAtomDeexcitation::G4VAtomDeexcitation ( const G4String &  modname = "Deexcitation")
explicit

Definition at line 73 of file G4VAtomDeexcitation.cc.

74 : verbose(1), name(modname), isActive(false), flagAuger(false),
75 flagAugerCascade(false), flagPIXE(false), ignoreCuts(false),
76 isActiveLocked(false), isAugerLocked(false),
77 isAugerCascadeLocked(false), isPIXELocked(false)
78{
79 theParameters = G4EmParameters::Instance();
80 vdyn.reserve(5);
81 theCoupleTable = nullptr;
82 G4String gg = "gammaPIXE";
83 G4String ee = "e-PIXE";
84 if(pixeIDg < 0) {
85#ifdef G4MULTITHREADED
86 G4MUTEXLOCK(&atomDeexcitationMutex);
87 if(pixeIDg < 0) {
88#endif
89 pixeIDg = G4PhysicsModelCatalog::Register(gg);
90 pixeIDe = G4PhysicsModelCatalog::Register(ee);
91#ifdef G4MULTITHREADED
92 }
93 G4MUTEXUNLOCK(&atomDeexcitationMutex);
94#endif
95 }
96 gamma = G4Gamma::Gamma();
97}
G4MUTEXLOCK
#define G4MUTEXLOCK(mutex)
Definition: G4Threading.hh:251
G4MUTEXUNLOCK
#define G4MUTEXUNLOCK(mutex)
Definition: G4Threading.hh:254
G4EmParameters::Instance
static G4EmParameters * Instance()
Definition: G4EmParameters.cc:67
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:85
G4PhysicsModelCatalog::Register
static G4int Register(const G4String &)
Definition: G4PhysicsModelCatalog.cc:51

◆ ~G4VAtomDeexcitation()

G4VAtomDeexcitation::~G4VAtomDeexcitation ( )
virtual

Definition at line 101 of file G4VAtomDeexcitation.cc.

102{}

◆ G4VAtomDeexcitation() [2/2]

G4VAtomDeexcitation::G4VAtomDeexcitation ( G4VAtomDeexcitation &  )
delete

Member Function Documentation

◆ AlongStepDeexcitation()

void G4VAtomDeexcitation::AlongStepDeexcitation ( std::vector< G4Track * > &  tracks,
const G4Step &  step,
G4double &  eLoss,
G4int  coupleIndex 
)

Definition at line 280 of file G4VAtomDeexcitation.cc.

284{
285 G4double truelength = step.GetStepLength();
286 if(!flagPIXE && !activePIXEMedia[coupleIndex]) { return; }
287 if(eLossMax <= 0.0 || truelength <= 0.0) { return; }
288
289 // step parameters
290 const G4StepPoint* preStep = step.GetPreStepPoint();
291 G4ThreeVector prePos = preStep->GetPosition();
292 G4ThreeVector delta = step.GetPostStepPoint()->GetPosition() - prePos;
293 G4double preTime = preStep->GetGlobalTime();
294 G4double dt = step.GetPostStepPoint()->GetGlobalTime() - preTime;
295
296 // particle parameters
297 const G4Track* track = step.GetTrack();
298 const G4ParticleDefinition* part = track->GetDefinition();
299 G4double ekin = preStep->GetKineticEnergy();
300
301 // media parameters
302 G4double gCut = (*theCoupleTable->GetEnergyCutsVector(0))[coupleIndex];
303 if(ignoreCuts) { gCut = 0.0; }
304 G4double eCut = DBL_MAX;
305 if(CheckAugerActiveRegion(coupleIndex)) {
306 eCut = (*theCoupleTable->GetEnergyCutsVector(1))[coupleIndex];
307 if(ignoreCuts) { eCut = 0.0; }
308 }
309
310 //G4cout<<"!Sample PIXE gCut(MeV)= "<<gCut<<" eCut(MeV)= "<<eCut
311 // <<" Ekin(MeV)= " << ekin/MeV << G4endl;
312
313 const G4Material* material = preStep->GetMaterial();
314 const G4ElementVector* theElementVector = material->GetElementVector();
315 const G4double* theAtomNumDensityVector =
316 material->GetVecNbOfAtomsPerVolume();
317 G4int nelm = material->GetNumberOfElements();
318
319 // loop over deexcitations
320 for(G4int i=0; i<nelm; ++i) {
321 G4int Z = (*theElementVector)[i]->GetZasInt();
322 if(activeZ[Z] && Z < 93) {
323 G4int nshells =
324 std::min(9,(*theElementVector)[i]->GetNbOfAtomicShells());
325 G4double rho = truelength*theAtomNumDensityVector[i];
326 //G4cout<<" Z "<< Z <<" is active x(mm)= " << truelength/mm << G4endl;
327 for(G4int ii=0; ii<nshells; ++ii) {
328 G4AtomicShellEnumerator as = G4AtomicShellEnumerator(ii);
329 const G4AtomicShell* shell = GetAtomicShell(Z, as);
330 G4double bindingEnergy = shell->BindingEnergy();
331
332 if(gCut > bindingEnergy) { break; }
333
334 if(eLossMax > bindingEnergy) {
335 G4double sig = rho*
336 GetShellIonisationCrossSectionPerAtom(part, Z, as, ekin, material);
337
338 // mfp is mean free path in units of step size
339 if(sig > 0.0) {
340 G4double mfp = 1.0/sig;
341 G4double stot = 0.0;
342 //G4cout << " Shell " << ii << " mfp(mm)= " << mfp/mm << G4endl;
343 // sample ionisation points
344 do {
345 stot -= mfp*std::log(G4UniformRand());
346 if( stot > 1.0 || eLossMax < bindingEnergy) { break; }
347 // sample deexcitation
348 vdyn.clear();
349 GenerateParticles(&vdyn, shell, Z, gCut, eCut);
350 G4int nsec = vdyn.size();
351 if(nsec > 0) {
352 G4ThreeVector r = prePos + stot*delta;
353 G4double time = preTime + stot*dt;
354 for(G4int j=0; j<nsec; ++j) {
355 G4DynamicParticle* dp = vdyn[j];
356 G4double e = dp->GetKineticEnergy();
357
358 // save new secondary if there is enough energy
359 if(eLossMax >= e) {
360 eLossMax -= e;
361 G4Track* t = new G4Track(dp, time, r);
362
363 // defined secondary type
364 if(dp->GetDefinition() == gamma) {
365 t->SetCreatorModelIndex(pixeIDg);
366 } else {
367 t->SetCreatorModelIndex(pixeIDe);
368 }
369
370 tracks.push_back(t);
371 } else {
372 delete dp;
373 }
374 }
375 }
376 // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
377 } while (stot < 1.0);
378 }
379 }
380 }
381 }
382 }
383 return;
384}
G4AtomicShellEnumerator
G4AtomicShellEnumerator
Definition: G4AtomicShellEnumerator.hh:50
G4ElementVector
std::vector< G4Element * > G4ElementVector
Definition: G4ElementVector.hh:43
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
G4AtomicShell::BindingEnergy
G4double BindingEnergy() const
Definition: G4AtomicShell.hh:79
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:188
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:184
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:204
G4ProductionCutsTable::GetEnergyCutsVector
const std::vector< G4double > * GetEnergyCutsVector(std::size_t pcIdx) const
Definition: G4ProductionCutsTable.hh:233
G4StepPoint::GetGlobalTime
G4double GetGlobalTime() const
G4StepPoint::GetMaterial
G4Material * GetMaterial() const
G4StepPoint::GetPosition
const G4ThreeVector & GetPosition() const
G4StepPoint::GetKineticEnergy
G4double GetKineticEnergy() const
G4Step::GetTrack
G4Track * GetTrack() const
G4Step::GetPreStepPoint
G4StepPoint * GetPreStepPoint() const
G4Step::GetStepLength
G4double GetStepLength() const
G4Step::GetPostStepPoint
G4StepPoint * GetPostStepPoint() const
G4Track::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4Track::SetCreatorModelIndex
void SetCreatorModelIndex(G4int idx)
G4VAtomDeexcitation::GetAtomicShell
virtual const G4AtomicShell * GetAtomicShell(G4int Z, G4AtomicShellEnumerator shell)=0
G4VAtomDeexcitation::CheckAugerActiveRegion
G4bool CheckAugerActiveRegion(G4int coupleIndex)
Definition: G4VAtomDeexcitation.hh:270
G4VAtomDeexcitation::GenerateParticles
void GenerateParticles(std::vector< G4DynamicParticle * > *secVect, const G4AtomicShell *, G4int Z, G4int coupleIndex)
Definition: G4VAtomDeexcitation.cc:254
G4VAtomDeexcitation::GetShellIonisationCrossSectionPerAtom
virtual G4double GetShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=nullptr)=0
G4InuclSpecialFunctions::bindingEnergy
G4double bindingEnergy(G4int A, G4int Z)
Definition: bindingEnergy.cc:36
DBL_MAX
#define DBL_MAX
Definition: templates.hh:62

Referenced by G4VEnergyLossProcess::AlongStepDoIt().

◆ CheckAugerActiveRegion()

G4bool G4VAtomDeexcitation::CheckAugerActiveRegion ( G4int  coupleIndex)
inline

Definition at line 270 of file G4VAtomDeexcitation.hh.

271{
272 return (activeAugerMedia[coupleIndex]);
273}

Referenced by AlongStepDeexcitation(), and GenerateParticles().

◆ CheckDeexcitationActiveRegion()

G4bool G4VAtomDeexcitation::CheckDeexcitationActiveRegion ( G4int  coupleIndex)
inline

Definition at line 264 of file G4VAtomDeexcitation.hh.

265{
266 return (activeDeexcitationMedia[coupleIndex]);
267}

Referenced by G4LivermoreComptonModel::SampleSecondaries(), G4LivermoreComptonModifiedModel::SampleSecondaries(), G4LivermorePhotoElectricModel::SampleSecondaries(), G4LivermorePolarizedComptonModel::SampleSecondaries(), G4LivermorePolarizedPhotoElectricGDModel::SampleSecondaries(), G4LivermorePolarizedPhotoElectricModel::SampleSecondaries(), G4LowEPComptonModel::SampleSecondaries(), G4LowEPPolarizedComptonModel::SampleSecondaries(), G4PenelopeComptonModel::SampleSecondaries(), G4PenelopeIonisationModel::SampleSecondaries(), G4PenelopePhotoElectricModel::SampleSecondaries(), G4KleinNishinaModel::SampleSecondaries(), and G4PEEffectFluoModel::SampleSecondaries().

◆ ComputeShellIonisationCrossSectionPerAtom()

virtual G4double G4VAtomDeexcitation::ComputeShellIonisationCrossSectionPerAtom ( const G4ParticleDefinition *  ,
G4int  Z,
G4AtomicShellEnumerator  shell,
G4double  kinE,
const G4Material *  mat = nullptr 
)
pure virtual

Implemented in G4UAtomicDeexcitation.

Referenced by G4EmCalculator::ComputeShellIonisationCrossSectionPerAtom().

◆ GenerateParticles() [1/2]

virtual void G4VAtomDeexcitation::GenerateParticles ( std::vector< G4DynamicParticle * > *  secVect,
const G4AtomicShell *  ,
G4int  Z,
G4double  gammaCut,
G4double  eCut 
)
pure virtual

Implemented in G4UAtomicDeexcitation.

◆ GenerateParticles() [2/2]

void G4VAtomDeexcitation::GenerateParticles ( std::vector< G4DynamicParticle * > *  secVect,
const G4AtomicShell *  as,
G4int  Z,
G4int  coupleIndex 
)

Definition at line 254 of file G4VAtomDeexcitation.cc.

257{
258 G4double gCut = DBL_MAX;
259 if(ignoreCuts) {
260 gCut = 0.0;
261 } else if (theCoupleTable) {
262 gCut = (*(theCoupleTable->GetEnergyCutsVector(0)))[idx];
263 }
264 if(gCut < as->BindingEnergy()) {
265 G4double eCut = DBL_MAX;
266 if(CheckAugerActiveRegion(idx)) {
267 if(ignoreCuts) {
268 eCut = 0.0;
269 } else if (theCoupleTable) {
270 eCut = (*(theCoupleTable->GetEnergyCutsVector(1)))[idx];
271 }
272 }
273 GenerateParticles(v, as, Z, gCut, eCut);
274 }
275}

Referenced by AlongStepDeexcitation(), G4ECDecay::DecayIt(), G4ITDecay::DecayIt(), GenerateParticles(), G4DNABornIonisationModel1::SampleSecondaries(), G4DNABornIonisationModel2::SampleSecondaries(), G4DNACPA100IonisationModel::SampleSecondaries(), G4DNAEmfietzoglouIonisationModel::SampleSecondaries(), G4DNARuddIonisationExtendedModel::SampleSecondaries(), G4DNARuddIonisationModel::SampleSecondaries(), G4LivermoreComptonModel::SampleSecondaries(), G4LivermoreComptonModifiedModel::SampleSecondaries(), G4LivermorePhotoElectricModel::SampleSecondaries(), G4LivermorePolarizedComptonModel::SampleSecondaries(), G4LivermorePolarizedPhotoElectricGDModel::SampleSecondaries(), G4LivermorePolarizedPhotoElectricModel::SampleSecondaries(), G4LowEPComptonModel::SampleSecondaries(), G4LowEPPolarizedComptonModel::SampleSecondaries(), G4MicroElecInelasticModel::SampleSecondaries(), G4PenelopeComptonModel::SampleSecondaries(), G4PenelopeIonisationModel::SampleSecondaries(), G4PenelopePhotoElectricModel::SampleSecondaries(), G4MicroElecInelasticModel_new::SampleSecondaries(), G4KleinNishinaModel::SampleSecondaries(), and G4PEEffectFluoModel::SampleSecondaries().

◆ GetAtomicShell()

virtual const G4AtomicShell * G4VAtomDeexcitation::GetAtomicShell ( G4int  Z,
G4AtomicShellEnumerator  shell 
)
pure virtual

Implemented in G4UAtomicDeexcitation.

Referenced by AlongStepDeexcitation(), G4ECDecay::DecayIt(), G4ITDecay::DecayIt(), G4DNABornIonisationModel1::SampleSecondaries(), G4DNABornIonisationModel2::SampleSecondaries(), G4DNACPA100IonisationModel::SampleSecondaries(), G4DNAEmfietzoglouIonisationModel::SampleSecondaries(), G4DNARuddIonisationExtendedModel::SampleSecondaries(), G4DNARuddIonisationModel::SampleSecondaries(), G4LivermoreComptonModel::SampleSecondaries(), G4LivermoreComptonModifiedModel::SampleSecondaries(), G4LivermorePhotoElectricModel::SampleSecondaries(), G4LivermorePolarizedComptonModel::SampleSecondaries(), G4LivermorePolarizedPhotoElectricGDModel::SampleSecondaries(), G4LivermorePolarizedPhotoElectricModel::SampleSecondaries(), G4LowEPComptonModel::SampleSecondaries(), G4LowEPPolarizedComptonModel::SampleSecondaries(), G4MicroElecInelasticModel::SampleSecondaries(), G4MicroElecInelasticModel_new::SampleSecondaries(), G4KleinNishinaModel::SampleSecondaries(), and G4PEEffectFluoModel::SampleSecondaries().

◆ GetListOfActiveAtoms()

const std::vector< G4bool > & G4VAtomDeexcitation::GetListOfActiveAtoms ( ) const
inline

Definition at line 248 of file G4VAtomDeexcitation.hh.

249{
250 return activeZ;
251}

◆ GetName()

const G4String & G4VAtomDeexcitation::GetName ( ) const
inline

Definition at line 242 of file G4VAtomDeexcitation.hh.

243{
244 return name;
245}

Referenced by SetDeexcitationActiveRegion().

◆ GetShellIonisationCrossSectionPerAtom()

virtual G4double G4VAtomDeexcitation::GetShellIonisationCrossSectionPerAtom ( const G4ParticleDefinition *  ,
G4int  Z,
G4AtomicShellEnumerator  shell,
G4double  kinE,
const G4Material *  mat = nullptr 
)
pure virtual

Implemented in G4UAtomicDeexcitation.

Referenced by AlongStepDeexcitation(), and G4EmCalculator::GetShellIonisationCrossSectionPerAtom().

◆ GetVerboseLevel()

G4int G4VAtomDeexcitation::GetVerboseLevel ( ) const
inline

Definition at line 258 of file G4VAtomDeexcitation.hh.

259{
260 return verbose;
261}

◆ InitialiseAtomicDeexcitation()

void G4VAtomDeexcitation::InitialiseAtomicDeexcitation ( )

Definition at line 106 of file G4VAtomDeexcitation.cc.

107{
108 theParameters->DefineRegParamForDeex(this);
109
110 // Define list of couples
111 theCoupleTable = G4ProductionCutsTable::GetProductionCutsTable();
112 G4int numOfCouples = theCoupleTable->GetTableSize();
113
114 // needed for unit tests
115 size_t nn = std::max(numOfCouples, 1);
116 if(activeDeexcitationMedia.size() != nn) {
117 activeDeexcitationMedia.resize(nn, false);
118 activeAugerMedia.resize(nn, false);
119 activePIXEMedia.resize(nn, false);
120 }
121 if(activeZ.size() != 93) { activeZ.resize(93, false); }
122
123 // initialisation of flags and options
124 // normally there is no locksed flags
125 if(!isActiveLocked) { isActive = theParameters->Fluo(); }
126 if(!isAugerLocked) { flagAuger = theParameters->Auger(); }
127 if(!isAugerCascadeLocked) { flagAugerCascade = theParameters->AugerCascade(); }
128 if(!isPIXELocked) { flagPIXE = theParameters->Pixe(); }
129 ignoreCuts = theParameters->DeexcitationIgnoreCut();
130
131 // Define list of regions
132 size_t nRegions = deRegions.size();
133 // check if deexcitation is active for the given run
134 if(!isActive && 0 == nRegions) { return; }
135
136 // if no active regions add a world
137 if(0 == nRegions) {
138 SetDeexcitationActiveRegion("World",isActive,flagAuger,flagPIXE);
139 nRegions = deRegions.size();
140 }
141
142 if(0 < verbose) {
143 G4cout << G4endl;
144 G4cout << "### === Deexcitation model " << name
145 << " is activated for " << nRegions;
146 if(1 == nRegions) { G4cout << " region:" << G4endl; }
147 else { G4cout << " regions:" << G4endl;}
148 }
149
150 // Identify active media
151 G4RegionStore* regionStore = G4RegionStore::GetInstance();
152 for(size_t j=0; j<nRegions; ++j) {
153 const G4Region* reg = regionStore->GetRegion(activeRegions[j], false);
154 if(reg && 0 < numOfCouples) {
155 const G4ProductionCuts* rpcuts = reg->GetProductionCuts();
156 if(0 < verbose) {
157 G4cout << " " << activeRegions[j]
158 << " " << deRegions[j] << " " << AugerRegions[j]
159 << " " << PIXERegions[j] << G4endl;
160 }
161 for(G4int i=0; i<numOfCouples; ++i) {
162 const G4MaterialCutsCouple* couple =
163 theCoupleTable->GetMaterialCutsCouple(i);
164 if (couple->GetProductionCuts() == rpcuts) {
165 activeDeexcitationMedia[i] = deRegions[j];
166 activeAugerMedia[i] = AugerRegions[j];
167 activePIXEMedia[i] = PIXERegions[j];
168 }
169 }
170 }
171 }
172 G4int nelm = G4Element::GetNumberOfElements();
173 //G4cout << nelm << G4endl;
174 for(G4int k=0; k<nelm; ++k) {
175 G4int Z = (*(G4Element::GetElementTable()))[k]->GetZasInt();
176 if(Z > 5 && Z < 93) {
177 activeZ[Z] = true;
178 //G4cout << "!!! Active de-excitation Z= " << Z << G4endl;
179 }
180 }
181
182 // Initialise derived class
183 InitialiseForNewRun();
184
185 if(0 < verbose && flagAuger) {
186 G4cout << "### === Auger cascade flag: " << flagAugerCascade
187 << G4endl;
188 }
189 if(0 < verbose) {
190 G4cout << "### === Ignore cuts flag: " << ignoreCuts
191 << G4endl;
192 }
193 if(0 < verbose && flagPIXE) {
194 G4cout << "### === PIXE model for hadrons: "
195 << theParameters->PIXECrossSectionModel()
196 << G4endl;
197 G4cout << "### === PIXE model for e+-: "
198 << theParameters->PIXEElectronCrossSectionModel()
199 << G4endl;
200 }
201}
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4Element::GetElementTable
static G4ElementTable * GetElementTable()
Definition: G4Element.cc:397
G4Element::GetNumberOfElements
static size_t GetNumberOfElements()
Definition: G4Element.cc:404
G4EmParameters::PIXECrossSectionModel
const G4String & PIXECrossSectionModel()
Definition: G4EmParameters.cc:1065
G4EmParameters::PIXEElectronCrossSectionModel
const G4String & PIXEElectronCrossSectionModel()
Definition: G4EmParameters.cc:1076
G4EmParameters::Fluo
G4bool Fluo() const
Definition: G4EmParameters.cc:248
G4EmParameters::Pixe
G4bool Pixe() const
Definition: G4EmParameters.cc:292
G4EmParameters::DefineRegParamForDeex
void DefineRegParamForDeex(G4VAtomDeexcitation *) const
Definition: G4EmParameters.cc:1259
G4EmParameters::AugerCascade
G4bool AugerCascade() const
Definition: G4EmParameters.cc:281
G4EmParameters::DeexcitationIgnoreCut
G4bool DeexcitationIgnoreCut() const
Definition: G4EmParameters.cc:303
G4EmParameters::Auger
G4bool Auger() const
Definition: G4EmParameters.cc:270
G4MaterialCutsCouple::GetProductionCuts
G4ProductionCuts * GetProductionCuts() const
Definition: G4MaterialCutsCouple.hh:141
G4ProductionCutsTable::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
Definition: G4ProductionCutsTable.hh:246
G4ProductionCutsTable::GetTableSize
std::size_t GetTableSize() const
Definition: G4ProductionCutsTable.hh:239
G4ProductionCutsTable::GetProductionCutsTable
static G4ProductionCutsTable * GetProductionCutsTable()
Definition: G4ProductionCutsTable.cc:58
G4RegionStore::GetInstance
static G4RegionStore * GetInstance()
Definition: G4RegionStore.cc:156
G4RegionStore::GetRegion
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
Definition: G4RegionStore.cc:211
G4Region::GetProductionCuts
G4ProductionCuts * GetProductionCuts() const
G4VAtomDeexcitation::SetDeexcitationActiveRegion
void SetDeexcitationActiveRegion(const G4String &rname, G4bool valDeexcitation, G4bool valAuger, G4bool valPIXE)
Definition: G4VAtomDeexcitation.cc:206
G4VAtomDeexcitation::InitialiseForNewRun
virtual void InitialiseForNewRun()=0

Referenced by LBE::ConstructGeneral(), and G4LossTableManager::ResetParameters().

◆ InitialiseForExtraAtom()

virtual void G4VAtomDeexcitation::InitialiseForExtraAtom ( G4int  Z)
pure virtual

Implemented in G4UAtomicDeexcitation.

◆ InitialiseForNewRun()

virtual void G4VAtomDeexcitation::InitialiseForNewRun ( )
pure virtual

Implemented in G4UAtomicDeexcitation.

Referenced by InitialiseAtomicDeexcitation().

◆ IsAugerActive()

G4bool G4VAtomDeexcitation::IsAugerActive ( ) const
inline

Definition at line 217 of file G4VAtomDeexcitation.hh.

218{
219 return flagAuger;
220}

◆ IsAugerCascadeActive()

G4bool G4VAtomDeexcitation::IsAugerCascadeActive ( ) const
inline

Definition at line 227 of file G4VAtomDeexcitation.hh.

228{
229 return flagAugerCascade;
230}

Referenced by G4UAtomicDeexcitation::GenerateParticles().

◆ IsFluoActive()

G4bool G4VAtomDeexcitation::IsFluoActive ( ) const
inline

Definition at line 207 of file G4VAtomDeexcitation.hh.

208{
209 return isActive;
210}

Referenced by G4ECDecay::DecayIt(), G4ITDecay::DecayIt(), G4LivermorePhotoElectricModel::Initialise(), G4LivermorePolarizedPhotoElectricGDModel::Initialise(), G4LivermorePolarizedPhotoElectricModel::Initialise(), and G4UAtomicDeexcitation::InitialiseForNewRun().

◆ IsPIXEActive()

G4bool G4VAtomDeexcitation::IsPIXEActive ( ) const
inline

Definition at line 237 of file G4VAtomDeexcitation.hh.

238{
239 return flagPIXE;
240}

Referenced by G4VEnergyLossProcess::BuildPhysicsTable(), G4PenelopeIonisationModel::Initialise(), and G4UAtomicDeexcitation::InitialiseForNewRun().

◆ operator=()

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

◆ SetAuger()

void G4VAtomDeexcitation::SetAuger ( G4bool  val)
inline

Definition at line 212 of file G4VAtomDeexcitation.hh.

213{
214 if(!isAugerLocked) { flagAuger = val; isAugerLocked = true; }
215}

◆ SetAugerCascade()

void G4VAtomDeexcitation::SetAugerCascade ( G4bool  val)
inline

Definition at line 222 of file G4VAtomDeexcitation.hh.

223{
224 if(!isAugerCascadeLocked) { flagAugerCascade = val; isAugerCascadeLocked = true; }
225}

◆ SetDeexcitationActiveRegion()

void G4VAtomDeexcitation::SetDeexcitationActiveRegion ( const G4String &  rname,
G4bool  valDeexcitation,
G4bool  valAuger,
G4bool  valPIXE 
)

Definition at line 206 of file G4VAtomDeexcitation.cc.

210{
211 // no PIXE in parallel world
212 if(rname == "DefaultRegionForParallelWorld") { return; }
213
214 G4String ss = rname;
215 /*
216 G4cout << "### G4VAtomDeexcitation::SetDeexcitationActiveRegion " << ss
217 << " " << valDeexcitation << " " << valAuger
218 << " " << valPIXE << G4endl;
219 */
220 if(ss == "world" || ss == "World" || ss == "WORLD") {
221 ss = "DefaultRegionForTheWorld";
222 }
223 size_t n = deRegions.size();
224 for(size_t i=0; i<n; ++i) {
225
226 // Region already exist
227 if(ss == activeRegions[i]) {
228 deRegions[i] = valDeexcitation;
229 AugerRegions[i] = valAuger;
230 PIXERegions[i] = valPIXE;
231 return;
232 }
233 }
234 // New region
235 activeRegions.push_back(ss);
236 deRegions.push_back(valDeexcitation);
237 AugerRegions.push_back(valAuger);
238 PIXERegions.push_back(valPIXE);
239
240 // if de-excitation defined for the world volume
241 // it should be active for all G4Regions
242 if(ss == "DefaultRegionForTheWorld") {
243 G4RegionStore* regions = G4RegionStore::GetInstance();
244 G4int nn = regions->size();
245 for(G4int i=0; i<nn; ++i) {
246 if(ss == (*regions)[i]->GetName()) { continue; }
247 SetDeexcitationActiveRegion((*regions)[i]->GetName(), valDeexcitation,
248 valAuger, valPIXE);
249
250 }
251 }
252}
G4VAtomDeexcitation::GetName
const G4String & GetName() const
Definition: G4VAtomDeexcitation.hh:242

Referenced by G4EmLowEParameters::DefineRegParamForDeex(), InitialiseAtomicDeexcitation(), and SetDeexcitationActiveRegion().

◆ SetFluo()

void G4VAtomDeexcitation::SetFluo ( G4bool  val)
inline

Definition at line 202 of file G4VAtomDeexcitation.hh.

203{
204 if(!isActiveLocked) { isActive = val; isActiveLocked = true; }
205}

Referenced by G4EmDNAPhysics_stationary_option2::ConstructProcess().

◆ SetPIXE()

void G4VAtomDeexcitation::SetPIXE ( G4bool  val)
inline

Definition at line 232 of file G4VAtomDeexcitation.hh.

233{
234 if(!isPIXELocked) { flagPIXE = val; isPIXELocked = true; }
235}

◆ SetVerboseLevel()

void G4VAtomDeexcitation::SetVerboseLevel ( G4int  val)
inline

Definition at line 253 of file G4VAtomDeexcitation.hh.

254{
255 verbose = val;
256}

Referenced by G4LossTableManager::ResetParameters().

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