71 : verbose(1),
name(modname), isActive(false),
72 flagAuger(false),flagAugerCascade(false),
73 flagPIXE(false), ignoreCuts(false)
125 G4cout <<
"### === Deexcitation model " <<
name
126 <<
" is activated for " << nRegions;
133 for(
size_t j=0; j<nRegions; ++j) {
135 if(reg && 0 < numOfCouples) {
140 for(
G4int i=0; i<numOfCouples; ++i) {
153 for(
G4int k=0; k<nelm; ++k) {
155 if(Z > 5 && Z < 93) {
165 G4cout <<
"### === PIXE model for hadrons: "
168 G4cout <<
"### === PIXE model for e+-: "
183 if(rname ==
"DefaultRegionForParallelWorld") {
return; }
189 if(ss ==
"world" || ss ==
"World" || ss ==
"WORLD") {
190 ss =
"DefaultRegionForTheWorld";
194 for(
size_t i=0; i<
n; ++i) {
213 if(ss ==
"DefaultRegionForTheWorld") {
234 if(eLossMax <= 0.0 || truelength <= 0.0) {
return; }
262 const G4double* theAtomNumDensityVector =
267 for(
G4int i=0; i<nelm; ++i) {
271 std::min(9,(*theElementVector)[i]->GetNbOfAtomicShells());
272 G4double rho = truelength*theAtomNumDensityVector[i];
274 for(
G4int ii=0; ii<nshells; ++ii) {
279 if(gCut > bindingEnergy) {
break; }
281 if(eLossMax > bindingEnergy) {
293 if( stot > 1.0 || eLossMax < bindingEnergy) {
break; }
301 for(
G4int j=0; j<nsec; ++j) {
324 }
while (stot < 1.0);
G4ParticleDefinition * GetDefinition() const
G4ProductionCuts * GetProductionCuts() const
G4Region * GetRegion(const G4String &name, G4bool verbose=true) const
const std::vector< G4double > * GetEnergyCutsVector(size_t pcIdx) const
void InitialiseAtomicDeexcitation()
void AlongStepDeexcitation(std::vector< G4Track * > &tracks, const G4Step &step, G4double &eLoss, G4int coupleIndex)
std::vector< G4String > activeRegions
std::vector< G4Element * > G4ElementVector
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
G4double GetStepLength() const
G4bool AugerCascade() const
std::vector< G4DynamicParticle * > vdyn
G4Material * GetMaterial() const
virtual G4double GetShellIonisationCrossSectionPerAtom(const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=0)=0
void SetDeexcitationActiveRegion(const G4String &rname, G4bool valDeexcitation, G4bool valAuger, G4bool valPIXE)
G4ParticleDefinition * GetDefinition() const
G4double BindingEnergy() const
const G4ElementVector * GetElementVector() const
G4bool DeexcitationIgnoreCut() const
std::vector< G4bool > activeAugerMedia
std::vector< G4bool > AugerRegions
void SetCreatorModelIndex(G4int idx)
static const G4double reg
static G4RegionStore * GetInstance()
G4StepPoint * GetPreStepPoint() const
virtual const G4AtomicShell * GetAtomicShell(G4int Z, G4AtomicShellEnumerator shell)=0
const G4double * GetVecNbOfAtomsPerVolume() const
G4GLOB_DLL std::ostream G4cout
static size_t GetNumberOfElements()
size_t GetTableSize() const
const G4ThreeVector & GetPosition() const
const G4String & PIXECrossSectionModel()
const G4ParticleDefinition * gamma
G4EmParameters * theParameters
G4ProductionCutsTable * theCoupleTable
static G4ProductionCutsTable * GetProductionCutsTable()
G4bool CheckAugerActiveRegion(G4int coupleIndex)
std::vector< G4bool > activeDeexcitationMedia
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
static G4int Register(const G4String &)
const G4String & GetName() const
G4StepPoint * GetPostStepPoint() const
static G4EmParameters * Instance()
virtual ~G4VAtomDeexcitation()
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
virtual void InitialiseForNewRun()=0
std::vector< G4bool > deRegions
G4double GetGlobalTime() const
size_t GetNumberOfElements() const
std::vector< G4bool > activeZ
G4double GetKineticEnergy() const
void GenerateParticles(std::vector< G4DynamicParticle * > *secVect, const G4AtomicShell *, G4int Z, G4int coupleIndex)
G4Track * GetTrack() const
static G4ElementTable * GetElementTable()
G4double bindingEnergy(G4int A, G4int Z)
G4ProductionCuts * GetProductionCuts() const
G4VAtomDeexcitation(const G4String &modname="Deexcitation")
const G4String & PIXEElectronCrossSectionModel()
std::vector< G4bool > PIXERegions
std::vector< G4bool > activePIXEMedia