65 #ifdef G4MULTITHREADED
75 : verbose(1),
name(modname), isActive(false), flagAuger(false),
76 flagAugerCascade(false), flagPIXE(false), ignoreCuts(false),
77 isActiveLocked(false), isAugerLocked(false),
78 isAugerCascadeLocked(false), isPIXELocked(false)
86 #ifdef G4MULTITHREADED
92 #ifdef G4MULTITHREADED
133 if(!
isActive && 0 == nRegions) {
return; }
143 G4cout <<
"### === Deexcitation model " <<
name
144 <<
" is activated for " << nRegions;
151 for(
size_t j=0; j<nRegions; ++j) {
153 if(reg && 0 < numOfCouples) {
160 for(
G4int i=0; i<numOfCouples; ++i) {
173 for(
G4int k=0; k<nelm; ++k) {
175 if(Z > 5 && Z < 93) {
193 G4cout <<
"### === PIXE model for hadrons: "
196 G4cout <<
"### === PIXE model for e+-: "
211 if(rname ==
"DefaultRegionForParallelWorld") {
return; }
219 if(ss ==
"world" || ss ==
"World" || ss ==
"WORLD") {
220 ss =
"DefaultRegionForTheWorld";
223 for(
size_t i=0; i<
n; ++i) {
241 if(ss ==
"DefaultRegionForTheWorld") {
245 if(ss == (*regions)[i]->
GetName()) {
continue; }
263 if(gCut < as->BindingEnergy()) {
286 if(eLossMax <= 0.0 || truelength <= 0.0) {
return; }
314 const G4double* theAtomNumDensityVector =
319 for(
G4int i=0; i<nelm; ++i) {
320 G4int Z = (*theElementVector)[i]->GetZasInt();
323 std::min(9,(*theElementVector)[i]->GetNbOfAtomicShells());
324 G4double rho = truelength*theAtomNumDensityVector[i];
326 for(
G4int ii=0; ii<nshells; ++ii) {
331 if(gCut > bindingEnergy) {
break; }
333 if(eLossMax > bindingEnergy) {
345 if( stot > 1.0 || eLossMax < bindingEnergy) {
break; }
353 for(
G4int j=0; j<nsec; ++j) {
376 }
while (stot < 1.0);
T max(const T t1, const T t2)
brief Return the largest of the two arguments
std::vector< G4bool > AugerRegions
G4EmParameters * theParameters
virtual const G4AtomicShell * GetAtomicShell(G4int Z, G4AtomicShellEnumerator shell)=0
std::vector< G4bool > deRegions
G4bool isAugerCascadeLocked
G4StepPoint * GetPreStepPoint() const
G4ProductionCutsTable * theCoupleTable
void AlongStepDeexcitation(std::vector< G4Track * > &tracks, const G4Step &step, G4double &eLoss, G4int coupleIndex)
const G4String & GetName() const
std::vector< G4bool > activeZ
void SetCreatorModelIndex(G4int idx)
void SetDeexcitationActiveRegion(const G4String &rname, G4bool valDeexcitation, G4bool valAuger, G4bool valPIXE)
void InitialiseAtomicDeexcitation()
const G4ElementVector * GetElementVector() const
std::vector< G4String > activeRegions
const G4double * GetVecNbOfAtomsPerVolume() const
G4ParticleDefinition * GetDefinition() const
const std::vector< G4double > * GetEnergyCutsVector(size_t pcIdx) const
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
std::vector< G4bool > PIXERegions
#define G4MUTEX_INITIALIZER
G4bool CheckAugerActiveRegion(G4int coupleIndex)
G4ParticleDefinition * GetDefinition() const
G4Material * GetMaterial() const
void GenerateParticles(std::vector< G4DynamicParticle * > *secVect, const G4AtomicShell *, G4int Z, G4int coupleIndex)
G4double GetStepLength() const
const G4ThreeVector & GetPosition() const
std::vector< G4DynamicParticle * > vdyn
std::vector< G4bool > activePIXEMedia
std::vector< G4bool > activeAugerMedia
G4Track * GetTrack() const
G4double GetGlobalTime() const
const G4ParticleDefinition * gamma
G4StepPoint * GetPostStepPoint() const
G4double GetKineticEnergy() const
size_t GetTableSize() const
#define G4MUTEXUNLOCK(mutex)
G4ProductionCuts * GetProductionCuts() const
G4bool AugerCascade() const
const G4String & PIXECrossSectionModel()
static G4ProductionCutsTable * GetProductionCutsTable()
G4VAtomDeexcitation(const G4String &modname="Deexcitation")
static G4RegionStore * GetInstance()
std::vector< G4Element * > G4ElementVector
G4bool DeexcitationIgnoreCut() const
const G4String & PIXEElectronCrossSectionModel()
#define G4MUTEXLOCK(mutex)
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
void DefineRegParamForDeex(G4VAtomDeexcitation *) const
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
G4ProductionCuts * GetProductionCuts() const
static size_t GetNumberOfElements()
static const G4double reg
std::vector< G4bool > activeDeexcitationMedia
G4double bindingEnergy(G4int A, G4int Z)
virtual G4double GetShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=nullptr)=0
virtual void InitialiseForNewRun()=0
virtual ~G4VAtomDeexcitation()
static G4EmParameters * Instance()
size_t GetNumberOfElements() const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static G4ElementTable * GetElementTable()
G4double BindingEnergy() const
static G4int Register(const G4String &)