#include <G4UAtomicDeexcitation.hh>

Inheritance diagram for G4UAtomicDeexcitation:

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

G4bool	CheckAugerActiveRegion (G4int coupleIndex)

G4bool	CheckDeexcitationActiveRegion (G4int coupleIndex)

G4double	ComputeShellIonisationCrossSectionPerAtom (const G4ParticleDefinition , G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material mat=nullptr) override
	access or compute PIXE cross section More...

	G4UAtomicDeexcitation ()

	G4UAtomicDeexcitation (G4UAtomicDeexcitation &)=delete

void	GenerateParticles (std::vector< G4DynamicParticle * > secVect, const G4AtomicShell , G4int Z, G4double gammaCut, G4double eCut) override
	generation of deexcitation for given atom, shell vacancy and cuts More...

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

const G4AtomicShell *	GetAtomicShell (G4int Z, G4AtomicShellEnumerator shell) override

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

const G4String &	GetName () const

G4double	GetShellIonisationCrossSectionPerAtom (const G4ParticleDefinition , G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material mat=nullptr) override
	access or compute PIXE cross section More...

G4int	GetVerboseLevel () const

void	InitialiseAtomicDeexcitation ()

void	InitialiseForExtraAtom (G4int Z) override

void	InitialiseForNewRun () override
	initialisation methods More...

G4bool	IsAugerActive () const

G4bool	IsAugerCascadeActive () const

G4bool	IsFluoActive () const

G4bool	IsPIXEActive () const

G4UAtomicDeexcitation &	operator= (const G4UAtomicDeexcitation &right)=delete

void	SetAuger (G4bool)

void	SetAugerCascade (G4bool)

void	SetCutForAugerElectrons (G4double cut)
	Set threshold energy for Auger electron production. More...

void	SetCutForSecondaryPhotons (G4double cut)
	Set threshold energy for fluorescence. More...

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

void	SetFluo (G4bool)

void	SetPIXE (G4bool)

void	SetVerboseLevel (G4int)

virtual	~G4UAtomicDeexcitation ()

Private Member Functions
G4DynamicParticle *	GenerateAuger (G4int Z, G4int shellId)
	Generates a particle from a non-radiative transition and returns it. More...

G4DynamicParticle *	GenerateAuger (G4int Z, G4int shellId, G4int &newAugerShellId)

G4DynamicParticle *	GenerateFluorescence (G4int Z, G4int shellId, G4int provShellId)
	Generates a particle from a radiative transition and returns it. More...

G4int	SelectTypeOfTransition (G4int Z, G4int shellId)

Private Attributes
std::vector< G4bool >	activeAugerMedia

std::vector< G4bool >	activeDeexcitationMedia

std::vector< G4bool >	activePIXEMedia

std::vector< G4String >	activeRegions

std::vector< G4bool >	activeZ

G4VhShellCrossSection *	anaPIXEshellCS

std::vector< G4bool >	AugerRegions

std::vector< G4bool >	deRegions

G4EmCorrections *	emcorr

G4VhShellCrossSection *	ePIXEshellCS

G4bool	flagAuger = false

G4bool	flagPIXE = false

const G4ParticleDefinition *	gamma

G4bool	ignoreCuts = false

G4bool	isActive = false

G4bool	isActiveLocked = false

G4bool	isAugerLocked = false

G4bool	isPIXELocked = false

G4double	minElectronEnergy

G4double	minGammaEnergy

G4String	name

G4int	nCouples = 0

G4int	newShellId

std::vector< G4bool >	PIXERegions

G4VhShellCrossSection *	PIXEshellCS
	Data member for the calculation of the proton and alpha ionisation XS. More...

const G4ProductionCutsTable *	theCoupleTable = nullptr

const G4ParticleDefinition *	theElectron

const G4ParticleDefinition *	thePositron

G4AtomicTransitionManager *	transitionManager

std::vector< int >	vacancyArray

std::vector< G4DynamicParticle * >	vdyn

G4int	verbose = 1

Detailed Description

Definition at line 60 of file G4UAtomicDeexcitation.hh.

Constructor & Destructor Documentation

◆ G4UAtomicDeexcitation() [1/2]

G4UAtomicDeexcitation::G4UAtomicDeexcitation ( )

explicit

Definition at line 77 of file G4UAtomicDeexcitation.cc.

                                            :
  G4VAtomDeexcitation("UAtomDeexcitation"),
  minGammaEnergy(DBL_MAX), 
  minElectronEnergy(DBL_MAX),
  newShellId(-1)
{
  anaPIXEshellCS = nullptr;
  PIXEshellCS    = nullptr;
  ePIXEshellCS   = nullptr;
  emcorr = G4LossTableManager::Instance()->EmCorrections();
  theElectron = G4Electron::Electron();
  thePositron = G4Positron::Positron();
  transitionManager = G4AtomicTransitionManager::Instance();
}

References anaPIXEshellCS, G4Electron::Electron(), emcorr, G4LossTableManager::EmCorrections(), ePIXEshellCS, G4AtomicTransitionManager::Instance(), G4LossTableManager::Instance(), PIXEshellCS, G4Positron::Positron(), theElectron, thePositron, and transitionManager.

◆ ~G4UAtomicDeexcitation()

G4UAtomicDeexcitation::~G4UAtomicDeexcitation ( )

virtual

Definition at line 94 of file G4UAtomicDeexcitation.cc.

{
  delete anaPIXEshellCS;
  delete PIXEshellCS;
  delete ePIXEshellCS;
}

References anaPIXEshellCS, ePIXEshellCS, and PIXEshellCS.

◆ G4UAtomicDeexcitation() [2/2]

G4UAtomicDeexcitation::G4UAtomicDeexcitation ( G4UAtomicDeexcitation & )

delete

Member Function Documentation

◆ AlongStepDeexcitation()

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

inherited

Definition at line 261 of file G4VAtomDeexcitation.cc.

{
  G4double truelength = step.GetStepLength();
  if(!flagPIXE && !activePIXEMedia[coupleIndex]) { return; }
  if(eLossMax <= 0.0 || truelength <= 0.0)       { return; }
 
  // step parameters
  const G4StepPoint* preStep = step.GetPreStepPoint();
  const G4ThreeVector prePos = preStep->GetPosition();
  const G4ThreeVector delta = step.GetPostStepPoint()->GetPosition() - prePos;
  const G4double preTime = preStep->GetGlobalTime();
  const G4double dt = step.GetPostStepPoint()->GetGlobalTime() - preTime;
 
  // particle parameters
  const G4Track* track = step.GetTrack();
  const G4ParticleDefinition* part = track->GetDefinition();
  G4double ekin = preStep->GetKineticEnergy(); 
 
  // media parameters
  G4double gCut = (*theCoupleTable->GetEnergyCutsVector(0))[coupleIndex];
  if(ignoreCuts) { gCut = 0.0; }
  G4double eCut = DBL_MAX;
  if(CheckAugerActiveRegion(coupleIndex)) { 
    eCut = (*theCoupleTable->GetEnergyCutsVector(1))[coupleIndex];
    if(ignoreCuts) { eCut = 0.0; }    
  }
 
  //G4cout<<"!Sample PIXE gCut(MeV)= "<<gCut<<"  eCut(MeV)= "<<eCut
  //        <<" Ekin(MeV)= " << ekin/MeV << G4endl;
 
  const G4Material* material = preStep->GetMaterial();
  const G4ElementVector* theElementVector = material->GetElementVector();
  const G4double* theAtomNumDensityVector = 
    material->GetVecNbOfAtomsPerVolume();
  const G4int nelm = material->GetNumberOfElements();
 
  // loop over deexcitations
  for(G4int i=0; i<nelm; ++i) {
    G4int Z = (*theElementVector)[i]->GetZasInt();
    if(activeZ[Z] && Z < 93) {  
      G4int nshells = 
        std::min(9,(*theElementVector)[i]->GetNbOfAtomicShells());
      G4double rho = truelength*theAtomNumDensityVector[i];
      //G4cout<<"   Z "<< Z <<" is active  x(mm)= " << truelength/mm << G4endl;
      for(G4int ii=0; ii<nshells; ++ii) {
        G4AtomicShellEnumerator as = G4AtomicShellEnumerator(ii);
        const G4AtomicShell* shell = GetAtomicShell(Z, as);
        const G4double bindingEnergy = shell->BindingEnergy();
 
        if(gCut > bindingEnergy) { break; }
 
        if(eLossMax > bindingEnergy) { 
          G4double sig = rho*
            GetShellIonisationCrossSectionPerAtom(part, Z, as, ekin, material);
 
          // mfp is mean free path in units of step size
          if(sig > 0.0) {
            G4double mfp = 1.0/sig;
            G4double stot = 0.0;
            //G4cout << " Shell " << ii << " mfp(mm)= " << mfp/mm << G4endl;
            // sample ionisation points
            do {
              stot -= mfp*G4Log(G4UniformRand());
              if( stot > 1.0 || eLossMax < bindingEnergy) { break; }
              // sample deexcitation
              vdyn.clear();
              GenerateParticles(&vdyn, shell, Z, gCut, eCut); 
              G4int nsec = vdyn.size();
              if(nsec > 0) {
                G4ThreeVector r = prePos  + stot*delta;
                G4double time   = preTime + stot*dt;
                for(G4int j=0; j<nsec; ++j) {
                  G4DynamicParticle* dp = vdyn[j];
                  G4double e = dp->GetKineticEnergy();
 
                  // save new secondary if there is enough energy
                  if(eLossMax >= e) {
                    eLossMax -= e;                    
                    G4Track* t = new G4Track(dp, time, r);
 
                    // defined secondary type
                    if(dp->GetDefinition() == gamma) { 
                      t->SetCreatorModelID(_GammaPIXE);
                    } else {
                      t->SetCreatorModelID(_ePIXE);
                    }
                    tracks.push_back(t);
                  } else {
                    delete dp;
                  }
                }
              }
              // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
            } while (stot < 1.0);
          }
        }
      }
    } 
  }
  return;
}

Referenced by G4VEnergyLossProcess::AlongStepDoIt().

◆ CheckAugerActiveRegion()

G4bool G4VAtomDeexcitation::CheckAugerActiveRegion ( G4int coupleIndex )

inlineinherited

Definition at line 266 of file G4VAtomDeexcitation.hh.

{
  return (idx < nCouples) ? activeAugerMedia[idx] : false;
}

References G4VAtomDeexcitation::activeAugerMedia, and G4VAtomDeexcitation::nCouples.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), and G4VAtomDeexcitation::GenerateParticles().

◆ CheckDeexcitationActiveRegion()

G4bool G4VAtomDeexcitation::CheckDeexcitationActiveRegion ( G4int coupleIndex )

inlineinherited

Definition at line 260 of file G4VAtomDeexcitation.hh.

{
  return (idx < nCouples) ? activeDeexcitationMedia[idx] : false;
}

References G4VAtomDeexcitation::activeDeexcitationMedia, and G4VAtomDeexcitation::nCouples.

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

◆ ComputeShellIonisationCrossSectionPerAtom()

G4double G4UAtomicDeexcitation::ComputeShellIonisationCrossSectionPerAtom	(	const G4ParticleDefinition *	p,
		G4int	Z,
		G4AtomicShellEnumerator	shell,
		G4double	kinE,
		const G4Material *	mat = `nullptr`
	)

overridevirtual

access or compute PIXE cross section

Implements G4VAtomDeexcitation.

Definition at line 386 of file G4UAtomicDeexcitation.cc.

{
  return GetShellIonisationCrossSectionPerAtom(p,Z,shell,kinE,mat);
}

References GetShellIonisationCrossSectionPerAtom(), and Z.

◆ GenerateAuger() [1/2]

G4DynamicParticle * G4UAtomicDeexcitation::GenerateAuger	(	G4int	Z,
		G4int	shellId
	)

private

Generates a particle from a non-radiative transition and returns it.

Definition at line 537 of file G4UAtomicDeexcitation.cc.

{
  if(!IsAugerActive()) { 
    //    G4cout << "auger inactive!" << G4endl; //debug
    return nullptr; 
  }
  
  if (shellId <=0 ) {
    //G4Exception("G4UAtomicDeexcitation::GenerateAuger()","de0002",
    //      JustWarning, "Energy deposited locally");
    return nullptr;
  }
 
  G4int maxNumOfShells = transitionManager->NumberOfReachableAugerShells(Z);  
  
  const G4AugerTransition* refAugerTransition = 
    transitionManager->ReachableAugerShell(Z,maxNumOfShells-1);
 
  // This loop gives to shellNum the value of the index of shellId
  // in the vector storing the list of the vacancies in the variuos shells 
  // that can originate a NON-radiative transition
  G4int shellNum = 0;
    
  if ( shellId <= refAugerTransition->FinalShellId() ) 
    // "FinalShellId" is final from the point of view of the electron 
    // who makes the transition, 
    // being the Id of the shell in which there is a vacancy
    {
      G4int pippo = transitionManager->ReachableAugerShell(Z,shellNum)->FinalShellId();
      if (shellId  != pippo ) {
    do { 
      shellNum++;
      if(shellNum == maxNumOfShells)
        {
          // G4cout << "No Auger transition found" << G4endl; //debug
          return 0;
        }
    }
    while (shellId != (transitionManager->ReachableAugerShell(Z,shellNum)->FinalShellId()) );
      }
 
      // Now we have that shellnum is the shellIndex of the shell named ShellId
      //      G4cout << " the index of the shell is: "<<shellNum<<G4endl;
      // But we have now to select two shells: one for the transition, 
      // and another for the auger emission.
      G4int transitionLoopShellIndex = 0;      
      G4double partSum = 0;
      const G4AugerTransition* anAugerTransition = 
    transitionManager->ReachableAugerShell(Z,shellNum);
 
      G4int transitionSize = 
    (anAugerTransition->TransitionOriginatingShellIds())->size();
      while (transitionLoopShellIndex < transitionSize) {
 
        std::vector<G4int>::const_iterator pos = 
      anAugerTransition->TransitionOriginatingShellIds()->begin();
 
        G4int transitionLoopShellId = *(pos+transitionLoopShellIndex);
        G4int numberOfPossibleAuger = 
      (anAugerTransition->AugerTransitionProbabilities(transitionLoopShellId))->size();
        G4int augerIndex = 0;
      
    if (augerIndex < numberOfPossibleAuger) {
      do 
        {
          G4double thisProb = anAugerTransition->AugerTransitionProbability(augerIndex, 
                                        transitionLoopShellId);
          partSum += thisProb;
          augerIndex++;
          
        } while (augerIndex < numberOfPossibleAuger);
    }
        transitionLoopShellIndex++;
      }
     
      G4double totalVacancyAugerProbability = partSum;
 
      //And now we start to select the right auger transition and emission
      G4int transitionRandomShellIndex = 0;
      G4int transitionRandomShellId = 1;
      G4int augerIndex = 0;
      partSum = 0; 
      G4double partialProb = G4UniformRand();
      
      G4int numberOfPossibleAuger = 0;      
      G4bool foundFlag = false;
 
      while (transitionRandomShellIndex < transitionSize) {
 
        std::vector<G4int>::const_iterator pos = 
      anAugerTransition->TransitionOriginatingShellIds()->begin();
 
        transitionRandomShellId = *(pos+transitionRandomShellIndex);
        
    augerIndex = 0;
    numberOfPossibleAuger = (anAugerTransition-> 
                 AugerTransitionProbabilities(transitionRandomShellId))->size();
 
        while (augerIndex < numberOfPossibleAuger) {
      G4double thisProb =anAugerTransition->AugerTransitionProbability(augerIndex, 
                                       transitionRandomShellId);
 
          partSum += thisProb;
          
          if (partSum >= (partialProb*totalVacancyAugerProbability) ) { // was /
        foundFlag = true;
        break;
      }
          augerIndex++;
        }
        if (partSum >= (partialProb*totalVacancyAugerProbability) ) {break;} // was /
        transitionRandomShellIndex++;
      }
 
      // Now we have the index of the shell from wich comes the auger electron (augerIndex), 
      // and the id of the shell, from which the transition e- come (transitionRandomShellid)
      // If no Transition has been found, 0 is returned.  
      if (!foundFlag) {
    return nullptr;
      } 
      
      // Isotropic angular distribution for the outcoming e-
      G4double newcosTh = 1.-2.*G4UniformRand();
      G4double  newsinTh = std::sqrt(1.-newcosTh*newcosTh);
      G4double newPhi = twopi*G4UniformRand();
      
      G4double xDir =  newsinTh*std::sin(newPhi);
      G4double yDir = newsinTh*std::cos(newPhi);
      G4double zDir = newcosTh;
      
      G4ThreeVector newElectronDirection(xDir,yDir,zDir);
      
      // energy of the auger electron emitted            
      G4double transitionEnergy = 
    anAugerTransition->AugerTransitionEnergy(augerIndex, transitionRandomShellId);
      
      if (transitionEnergy < minElectronEnergy) {
    return nullptr;
      }
 
      // This is the shell where the new vacancy is: it is the same
      // shell where the electron came from
      newShellId = transitionRandomShellId;
      
      //Auger cascade by Burkhant Suerfu on March 24 2015 (Bugzilla 1727)
      if (IsAugerCascadeActive())
    {
      vacancyArray.push_back(newShellId);
      vacancyArray.push_back(anAugerTransition->AugerOriginatingShellId(augerIndex,transitionRandomShellId));
    }
     
      return new G4DynamicParticle(G4Electron::Electron(), 
                   newElectronDirection,
                   transitionEnergy);
    }
  else 
    {
      return nullptr;
    }
}

References G4AugerTransition::AugerOriginatingShellId(), G4AugerTransition::AugerTransitionEnergy(), G4AugerTransition::AugerTransitionProbabilities(), G4AugerTransition::AugerTransitionProbability(), G4Electron::Electron(), G4AugerTransition::FinalShellId(), G4UniformRand, G4VAtomDeexcitation::IsAugerActive(), G4VAtomDeexcitation::IsAugerCascadeActive(), minElectronEnergy, newShellId, G4AtomicTransitionManager::NumberOfReachableAugerShells(), pos, G4AtomicTransitionManager::ReachableAugerShell(), transitionManager, G4AugerTransition::TransitionOriginatingShellIds(), twopi, vacancyArray, and Z.

Referenced by GenerateParticles().

◆ GenerateAuger() [2/2]

G4DynamicParticle * G4UAtomicDeexcitation::GenerateAuger	(	G4int	Z,
		G4int	shellId,
		G4int &	newAugerShellId
	)

private

Auger cascade by Burkhant Suerfu on March 24 2015 (Bugzilla 1727) Generates auger electron cascade.

◆ GenerateFluorescence()

G4DynamicParticle * G4UAtomicDeexcitation::GenerateFluorescence	(	G4int	Z,
		G4int	shellId,
		G4int	provShellId
	)

private

Generates a particle from a radiative transition and returns it.

Definition at line 464 of file G4UAtomicDeexcitation.cc.

{ 
  if (shellId <=0 )
    {
      //G4Exception("G4UAtomicDeexcitation::GenerateFluorescence()","de0002",JustWarning, "Energy deposited locally");
      return nullptr;
    }
 
  //isotropic angular distribution for the outcoming photon
  G4double newcosTh = 1.-2.*G4UniformRand();
  G4double newsinTh = std::sqrt((1.-newcosTh)*(1. + newcosTh));
  G4double newPhi = twopi*G4UniformRand();
  
  G4double xDir = newsinTh*std::sin(newPhi);
  G4double yDir = newsinTh*std::cos(newPhi);
  G4double zDir = newcosTh;
  
  G4ThreeVector newGammaDirection(xDir,yDir,zDir);
  
  G4int shellNum = 0;
  G4int maxNumOfShells = transitionManager->NumberOfReachableShells(Z);
  
  // find the index of the shell named shellId
  while (shellId != transitionManager->
     ReachableShell(Z,shellNum)->FinalShellId())
    {
      if(shellNum == maxNumOfShells-1)
    {
      break;
    }
      shellNum++;
    }
  // number of shell from wich an electron can reach shellId
  size_t transitionSize = transitionManager->
    ReachableShell(Z,shellNum)->OriginatingShellIds().size();
  
  size_t index = 0;
  
  // find the index of the shell named provShellId in the vector
  // storing the shells from which shellId can be reached 
  while (provShellId != transitionManager->
     ReachableShell(Z,shellNum)->OriginatingShellId(index))
    {
      if(index ==  transitionSize-1)
    {
      break;
    }
      index++;
    }
  // energy of the gamma leaving provShellId for shellId
  G4double transitionEnergy = transitionManager->
    ReachableShell(Z,shellNum)->TransitionEnergy(index);
  
  if (transitionEnergy < minGammaEnergy) return nullptr;
 
  // This is the shell where the new vacancy is: it is the same
  // shell where the electron came from
  newShellId = transitionManager->
    ReachableShell(Z,shellNum)->OriginatingShellId(index);
    
  G4DynamicParticle* newPart = new G4DynamicParticle(G4Gamma::Gamma(), 
                             newGammaDirection,
                             transitionEnergy);
 
  //Auger cascade by Burkhant Suerfu on March 24 2015 (Bugzilla 1727)
  if (IsAugerCascadeActive()) vacancyArray.push_back(newShellId);
 
  return newPart;
}

References G4UniformRand, G4Gamma::Gamma(), G4VAtomDeexcitation::IsAugerCascadeActive(), minGammaEnergy, newShellId, G4AtomicTransitionManager::NumberOfReachableShells(), transitionManager, twopi, vacancyArray, and Z.

Referenced by GenerateParticles().

◆ GenerateParticles() [1/2]

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

overridevirtual

generation of deexcitation for given atom, shell vacancy and cuts

Implements G4VAtomDeexcitation.

Definition at line 186 of file G4UAtomicDeexcitation.cc.

{
  // Defined initial conditions
  G4int givenShellId = atomicShell->ShellId();
  minGammaEnergy = gammaCut;
  minElectronEnergy = eCut;
 
  // generation secondaries
  G4DynamicParticle* aParticle=0;
  G4int provShellId = 0;
  
  //ORIGINAL METHOD BY ALFONSO MANTERO
  if (!IsAugerCascadeActive())
    {
      //----------------------------  
      G4int counter = 0;
      
      // let's check that 5<Z<100
      
      if (Z>5 && Z<100) {
 
    // The aim of this loop is to generate more than one fluorecence photon 
    // from the same ionizing event 
    do
      {
        if (counter == 0) 
          // First call to GenerateParticles(...):
          // givenShellId is given by the process
          {
        provShellId = SelectTypeOfTransition(Z, givenShellId);
        
        if  ( provShellId >0) 
          {
            aParticle = GenerateFluorescence(Z,givenShellId,provShellId);           
          }
        else if ( provShellId == -1)
          {
            aParticle = GenerateAuger(Z, givenShellId);
          }
        else
          {
            //G4Exception("G4UAtomicDeexcitation::GenerateParticles()",
            //      "de0002",JustWarning, "Energy deposited locally");
          }
          }
        else 
          // Following calls to GenerateParticles(...):
          // newShellId is given by GenerateFluorescence(...)
          {
        provShellId = SelectTypeOfTransition(Z,newShellId);
        if  (provShellId >0)
          {
            aParticle = GenerateFluorescence(Z,newShellId,provShellId);
          }
        else if ( provShellId == -1)
          {
            aParticle = GenerateAuger(Z, newShellId);       
          }
        else
          {
            //G4Exception("G4UAtomicDeexcitation::GenerateParticles()","de0002",JustWarning, "Energy deposited locally");
          }
          }
        counter++;
        if (aParticle != 0) 
          {
        vectorOfParticles->push_back(aParticle);
        //G4cout << "Deexcitation Occurred!" << G4endl; //debug
          }
        else {provShellId = -2;}
      }  
    while (provShellId > -2); 
      }
      else
    {
      //G4Exception("G4UAtomicDeexcitation::GenerateParticles()","de0001",JustWarning, "Energy deposited locally");
    }
 
    } // Auger cascade is not active
 
  //END OF ORIGINAL METHOD BY ALFONSO MANTERO
  //----------------------
 
  // NEW METHOD
  // Auger cascade by Burkhant Suerfu on March 24 2015 (Bugzilla 1727)
  if (IsAugerCascadeActive())
    {
      //----------------------
      vacancyArray.push_back(givenShellId);
 
      // let's check that 5<Z<100
      if (Z<6 || Z>99){
    //G4Exception("G4UAtomicDeexcitation::GenerateParticles()","de0001",JustWarning, "Energy deposited locally");
    return;
      }
 
      // as long as there is vacancy to be filled by either fluo or auger, stay in the loop.
      while(!vacancyArray.empty()){
    //  prepare to process the last element, and then delete it from the vector.
    givenShellId = vacancyArray[0];
    provShellId = SelectTypeOfTransition(Z,givenShellId);
 
    //G4cout<<"\n------ Atom Transition with Z: "<<Z<<"\tbetween current:"
    //      <<givenShellId<<" & target:"<<provShellId<<G4endl;
    if(provShellId>0){
      aParticle = GenerateFluorescence(Z,givenShellId,provShellId);
    }
    else if(provShellId == -1){
      aParticle = GenerateAuger(Z, givenShellId);
    }
    //else
    //  G4Exception("G4UAtomicDeexcitation::GenerateParticles()","de0002",JustWarning, "Energy deposited locally");
 
    //  if a particle is created, put it in the vector of new particles
    if(aParticle!=0)
      vectorOfParticles->push_back(aParticle);
    else{;}
    //  one vacancy has been processed. Erase it.
    vacancyArray.erase(vacancyArray.begin());
      }
      //----------------------
      //End of Auger cascade by Burkhant Suerfu on March 24 2015 (Bugzilla 1727)
 
    } // Auger cascade is active
}

References GenerateAuger(), GenerateFluorescence(), G4VAtomDeexcitation::IsAugerCascadeActive(), minElectronEnergy, minGammaEnergy, newShellId, SelectTypeOfTransition(), G4AtomicShell::ShellId(), vacancyArray, and Z.

◆ GenerateParticles() [2/2]

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

inherited

Definition at line 235 of file G4VAtomDeexcitation.cc.

{
  G4double gCut = DBL_MAX;
  if(ignoreCuts) {
    gCut = 0.0;
  } else if (nullptr != theCoupleTable) {
    gCut = (*(theCoupleTable->GetEnergyCutsVector(0)))[idx];
  }
  if(gCut < as->BindingEnergy()) {
    G4double eCut = DBL_MAX;
    if(CheckAugerActiveRegion(idx)) {
      if(ignoreCuts) {
        eCut = 0.0;
      } else if (nullptr != theCoupleTable) {
        eCut = (*(theCoupleTable->GetEnergyCutsVector(1)))[idx];
      }
    }
    GenerateParticles(v, as, Z, gCut, eCut);
  }
}

References G4VAtomDeexcitation::CheckAugerActiveRegion(), DBL_MAX, G4VAtomDeexcitation::GenerateParticles(), G4ProductionCutsTable::GetEnergyCutsVector(), G4VAtomDeexcitation::ignoreCuts, G4VAtomDeexcitation::theCoupleTable, and Z.

◆ GetAtomicShell()

const G4AtomicShell * G4UAtomicDeexcitation::GetAtomicShell	(	G4int	Z,
		G4AtomicShellEnumerator	shell
	)

overridevirtual

Get atomic shell by shell index, used by discrete processes (for example, photoelectric), when shell vacancy sampled by the model

Implements G4VAtomDeexcitation.

Definition at line 179 of file G4UAtomicDeexcitation.cc.

{
  return transitionManager->Shell(Z, size_t(shell));
}

References G4AtomicTransitionManager::Shell(), transitionManager, and Z.

◆ GetListOfActiveAtoms()

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

inlineinherited

Definition at line 244 of file G4VAtomDeexcitation.hh.

{
  return activeZ;
}

References G4VAtomDeexcitation::activeZ.

◆ GetName()

const G4String & G4VAtomDeexcitation::GetName ( ) const

inlineinherited

Definition at line 238 of file G4VAtomDeexcitation.hh.

{
  return name;
}

References G4VAtomDeexcitation::name.

Referenced by G4VAtomDeexcitation::SetDeexcitationActiveRegion().

◆ GetShellIonisationCrossSectionPerAtom()

G4double G4UAtomicDeexcitation::GetShellIonisationCrossSectionPerAtom	(	const G4ParticleDefinition *	pdef,
		G4int	Z,
		G4AtomicShellEnumerator	shell,
		G4double	kinE,
		const G4Material *	mat = `nullptr`
	)

overridevirtual

access or compute PIXE cross section

Implements G4VAtomDeexcitation.

Definition at line 320 of file G4UAtomicDeexcitation.cc.

{
  // we must put a control on the shell that are passed: 
  // some shells should not pass (line "0" or "2")
 
  // check atomic number
  G4double xsec = 0.0;
  if(Z > 93 || Z < 6 ) { return xsec; } //corrected by alf - Z<6 missing
  G4int idx = G4int(shellEnum);
  if(idx >= G4AtomicShells::GetNumberOfShells(Z)) { return xsec; }
 
  if(pdef == theElectron || pdef == thePositron) {
    xsec = ePIXEshellCS->CrossSection(Z,shellEnum,kineticEnergy,0.0,mat);
    return xsec;
  }
 
  G4double mass = pdef->GetPDGMass();
  G4double escaled = kineticEnergy;
  G4double q2 = 0.0;
 
  // scaling to protons
  if ((pdef->GetParticleName() != "proton" && pdef->GetParticleName() != "alpha" ) )
    {
      mass = proton_mass_c2;
      escaled = kineticEnergy*mass/(pdef->GetPDGMass());
 
      if(mat) {
    q2 = emcorr->EffectiveChargeSquareRatio(pdef,mat,kineticEnergy);
      } else {
    G4double q = pdef->GetPDGCharge()/eplus;
    q2 = q*q;
      }
    }
  
  if(PIXEshellCS) { xsec = PIXEshellCS->CrossSection(Z,shellEnum,escaled,mass,mat); }
  if(xsec < 1e-100) {     
    xsec = anaPIXEshellCS->CrossSection(Z,shellEnum,escaled,mass,mat); 
  }
 
  if (q2)  {xsec *= q2;}
 
  return xsec;
}

References anaPIXEshellCS, G4VhShellCrossSection::CrossSection(), G4EmCorrections::EffectiveChargeSquareRatio(), emcorr, ePIXEshellCS, eplus, G4AtomicShells::GetNumberOfShells(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), PIXEshellCS, source.hepunit::proton_mass_c2, theElectron, thePositron, and Z.

Referenced by ComputeShellIonisationCrossSectionPerAtom().

◆ GetVerboseLevel()

G4int G4VAtomDeexcitation::GetVerboseLevel ( ) const

inlineinherited

Definition at line 254 of file G4VAtomDeexcitation.hh.

{
  return verbose;
}

References G4VAtomDeexcitation::verbose.

◆ InitialiseAtomicDeexcitation()

void G4VAtomDeexcitation::InitialiseAtomicDeexcitation ( )

inherited

Definition at line 87 of file G4VAtomDeexcitation.cc.

{
  G4EmParameters* theParameters = G4EmParameters::Instance();
  theParameters->DefineRegParamForDeex(this);
 
  // Define list of couples
  theCoupleTable = G4ProductionCutsTable::GetProductionCutsTable();
  nCouples = theCoupleTable->GetTableSize();
 
  // needed for unit tests
  size_t nn = std::max(nCouples, 1);
  if(activeDeexcitationMedia.size() != nn) {
    activeDeexcitationMedia.resize(nn, false);
    activeAugerMedia.resize(nn, false);
    activePIXEMedia.resize(nn, false);
  }
  if(activeZ.size() != 93) { activeZ.resize(93, false); }
 
  // initialisation of flags and options
  // normally there is no locksed flags
  if(!isActiveLocked) { isActive  = theParameters->Fluo(); }
  if(!isAugerLocked)  { flagAuger = theParameters->Auger(); }
  if(!isPIXELocked)   { flagPIXE  = theParameters->Pixe(); }
  ignoreCuts = theParameters->DeexcitationIgnoreCut();
 
  // Define list of regions
  size_t nRegions = deRegions.size();
  // check if deexcitation is active for the given run
  if(!isActive && 0 == nRegions) { return; }
 
  // if no active regions add a world
  if(0 == nRegions) {
    SetDeexcitationActiveRegion("World",isActive,flagAuger,flagPIXE);
    nRegions = deRegions.size();
  }
 
  if(0 < verbose) {
    G4cout << G4endl;
    G4cout << "### ===  Deexcitation model " << name 
           << " is activated for " << nRegions;
    if(1 == nRegions) { G4cout << " region:" << G4endl; }
    else              { G4cout << " regions:" << G4endl;}
  }
 
  // Identify active media
  const G4RegionStore* regionStore = G4RegionStore::GetInstance();
  for(size_t j=0; j<nRegions; ++j) {
    const G4Region* reg = regionStore->GetRegion(activeRegions[j], false);
    if(nullptr != reg && 0 < nCouples) {
      const G4ProductionCuts* rpcuts = reg->GetProductionCuts();
      if(0 < verbose) {
        G4cout << "          " << activeRegions[j]
               << "  " << deRegions[j]  << "  " << AugerRegions[j]
               << "  " << PIXERegions[j] << G4endl;  
      }
      for(G4int i=0; i<nCouples; ++i) {
        const G4MaterialCutsCouple* couple =
          theCoupleTable->GetMaterialCutsCouple(i);
        if (couple->GetProductionCuts() == rpcuts) {
          activeDeexcitationMedia[i] = deRegions[j];
          activeAugerMedia[i] = AugerRegions[j];
          activePIXEMedia[i] = PIXERegions[j];
        }
      }
    }
  }
  G4int nelm = G4Element::GetNumberOfElements();
  //G4cout << nelm << G4endl;
  for(G4int k=0; k<nelm; ++k) {
    G4int Z = (*(G4Element::GetElementTable()))[k]->GetZasInt();
    if(Z > 5 && Z < 93) { 
      activeZ[Z] = true;
      //G4cout << "!!! Active de-excitation Z= " << Z << G4endl;  
    }
  }
 
  // Initialise derived class
  InitialiseForNewRun();
 
  if(0 < verbose && flagAuger) {
    G4cout << "### ===  Auger flag: " << flagAuger 
           << G4endl;
  }
  if(0 < verbose) {
    G4cout << "### ===  Ignore cuts flag:   " << ignoreCuts
           << G4endl;
  }
  if(0 < verbose && flagPIXE) {
    G4cout << "### ===  PIXE model for hadrons: " 
           << theParameters->PIXECrossSectionModel()
           << G4endl;  
    G4cout << "### ===  PIXE model for e+-:     " 
           << theParameters->PIXEElectronCrossSectionModel()
           << G4endl;  
  }
}

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

◆ InitialiseForExtraAtom()

void G4UAtomicDeexcitation::InitialiseForExtraAtom ( G4int Z )

overridevirtual

Implements G4VAtomDeexcitation.

Definition at line 173 of file G4UAtomicDeexcitation.cc.

174{}

◆ InitialiseForNewRun()

void G4UAtomicDeexcitation::InitialiseForNewRun ( )

overridevirtual

initialisation methods

Implements G4VAtomDeexcitation.

Definition at line 103 of file G4UAtomicDeexcitation.cc.

{
  if(!IsFluoActive()) { return; }
  transitionManager->Initialise();
  if(!IsPIXEActive()) { return; }
 
  if(!anaPIXEshellCS) {
    anaPIXEshellCS = new G4teoCrossSection("ECPSSR_Analytical");
  }
  G4cout << G4endl;
  G4cout << "### === G4UAtomicDeexcitation::InitialiseForNewRun()" << G4endl;
 
  G4EmParameters* param = G4EmParameters::Instance();
  G4String namePIXExsModel = param->PIXECrossSectionModel();
  G4String namePIXExsElectronModel = param->PIXEElectronCrossSectionModel();
 
  // Check if old cross section for p/ion should be deleted 
  if(PIXEshellCS && namePIXExsModel != PIXEshellCS->GetName()) 
    {
      delete PIXEshellCS;
      PIXEshellCS = nullptr;
    }
 
  // Instantiate new proton/ion cross section
  if(!PIXEshellCS) {
    if (namePIXExsModel == "ECPSSR_FormFactor")
      {
    PIXEshellCS = new G4teoCrossSection(namePIXExsModel);
      }
    else if(namePIXExsModel == "ECPSSR_ANSTO")
      {
    PIXEshellCS = new G4teoCrossSection(namePIXExsModel);
      }    
    else if(namePIXExsModel == "Empirical")
      {
    PIXEshellCS = new G4empCrossSection(namePIXExsModel);
      }
  }
 
  // Check if old cross section for e+- should be deleted 
  if(ePIXEshellCS && namePIXExsElectronModel != ePIXEshellCS->GetName()) 
    {
      delete ePIXEshellCS;
      ePIXEshellCS = nullptr;
    } 
 
  // Instantiate new e+- cross section
  if(!ePIXEshellCS) 
    {
      if(namePIXExsElectronModel == "Empirical")
    {
      ePIXEshellCS = new G4empCrossSection("Empirical");
    }
      else if(namePIXExsElectronModel == "ECPSSR_Analytical") 
    {
      ePIXEshellCS = new G4teoCrossSection("ECPSSR_Analytical");
    }
      else if (namePIXExsElectronModel == "Penelope")
    {
      ePIXEshellCS = new G4PenelopeIonisationCrossSection();
    }
      else 
    {
      ePIXEshellCS = new G4LivermoreIonisationCrossSection();
    }
    } 
}

References anaPIXEshellCS, ePIXEshellCS, G4cout, G4endl, G4VhShellCrossSection::GetName(), G4AtomicTransitionManager::Initialise(), G4EmParameters::Instance(), G4VAtomDeexcitation::IsFluoActive(), G4VAtomDeexcitation::IsPIXEActive(), G4EmParameters::PIXECrossSectionModel(), G4EmParameters::PIXEElectronCrossSectionModel(), PIXEshellCS, and transitionManager.

◆ IsAugerActive()

G4bool G4VAtomDeexcitation::IsAugerActive ( ) const

inlineinherited

Definition at line 213 of file G4VAtomDeexcitation.hh.

{
  return flagAuger;
}

References G4VAtomDeexcitation::flagAuger.

Referenced by GenerateAuger().

◆ IsAugerCascadeActive()

G4bool G4VAtomDeexcitation::IsAugerCascadeActive ( ) const

inlineinherited

Definition at line 223 of file G4VAtomDeexcitation.hh.

{
  return flagAuger;
}

References G4VAtomDeexcitation::flagAuger.

Referenced by GenerateAuger(), GenerateFluorescence(), and GenerateParticles().

◆ IsFluoActive()

G4bool G4VAtomDeexcitation::IsFluoActive ( ) const

inlineinherited

Definition at line 203 of file G4VAtomDeexcitation.hh.

{
  return isActive;
}

References G4VAtomDeexcitation::isActive.

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

◆ IsPIXEActive()

G4bool G4VAtomDeexcitation::IsPIXEActive ( ) const

inlineinherited

Definition at line 233 of file G4VAtomDeexcitation.hh.

{
  return flagPIXE;
}

References G4VAtomDeexcitation::flagPIXE.

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

◆ operator=()

G4UAtomicDeexcitation & G4UAtomicDeexcitation::operator= ( const G4UAtomicDeexcitation & right )

delete

◆ SelectTypeOfTransition()

G4int G4UAtomicDeexcitation::SelectTypeOfTransition	(	G4int	Z,
		G4int	shellId
	)

private

Decides wether a radiative transition is possible and, if it is, returns the identity of the starting shell for the transition

Definition at line 398 of file G4UAtomicDeexcitation.cc.

{
  if (shellId <=0 ) {
    //G4Exception("G4UAtomicDeexcitation::SelecttypeOfTransition()","de0002",
    //      JustWarning, "Energy deposited locally");
    return 0;
  }
  
  G4int provShellId = -1;
  G4int shellNum = 0;
  G4int maxNumOfShells = transitionManager->NumberOfReachableShells(Z);  
  
  const G4FluoTransition* refShell = 
    transitionManager->ReachableShell(Z,maxNumOfShells-1);
 
  // This loop gives shellNum the value of the index of shellId
  // in the vector storing the list of the shells reachable through
  // a radiative transition
  if ( shellId <= refShell->FinalShellId())
    {
      while (shellId != transitionManager->ReachableShell(Z,shellNum)->FinalShellId())
    {
      if(shellNum ==maxNumOfShells-1)
        {
          break;
        }
      shellNum++;
    }
      G4int transProb = 0; //AM change 29/6/07 was 1
   
      G4double partialProb = G4UniformRand();      
      G4double partSum = 0;
      const G4FluoTransition* aShell = transitionManager->ReachableShell(Z,shellNum);
      G4int trSize =  (aShell->TransitionProbabilities()).size();
    
      // Loop over the shells wich can provide an electron for a 
      // radiative transition towards shellId:
      // in every loop the partial sum of the first transProb shells
      // is calculated and compared with a random number [0,1].
      // If the partial sum is greater, the shell whose index is transProb
      // is chosen as the starting shell for a radiative transition
      // and its identity is returned
      // Else, terminateded the loop, -1 is returned
      while(transProb < trSize){
    partSum += aShell->TransitionProbability(transProb);
 
    if(partialProb <= partSum)
      {
        provShellId = aShell->OriginatingShellId(transProb);
        break;
      }
    transProb++;
      }
      // here provShellId is the right one or is -1.
      // if -1, the control is passed to the Auger generation part of the package 
    }
  else 
    {
      provShellId = -1;
    }
  return provShellId;
}

References G4FluoTransition::FinalShellId(), G4UniformRand, G4AtomicTransitionManager::NumberOfReachableShells(), G4FluoTransition::OriginatingShellId(), G4AtomicTransitionManager::ReachableShell(), transitionManager, G4FluoTransition::TransitionProbabilities(), G4FluoTransition::TransitionProbability(), and Z.

Referenced by GenerateParticles().

◆ SetAuger()

void G4VAtomDeexcitation::SetAuger ( G4bool val )

inlineinherited

Definition at line 208 of file G4VAtomDeexcitation.hh.

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

References G4VAtomDeexcitation::flagAuger, and G4VAtomDeexcitation::isAugerLocked.

Referenced by G4VAtomDeexcitation::SetAugerCascade().

◆ SetAugerCascade()

void G4VAtomDeexcitation::SetAugerCascade ( G4bool val )

inlineinherited

Definition at line 218 of file G4VAtomDeexcitation.hh.

{
  SetAuger(val);
}

References G4VAtomDeexcitation::SetAuger().

◆ SetCutForAugerElectrons()

void G4UAtomicDeexcitation::SetCutForAugerElectrons ( G4double cut )

Set threshold energy for Auger electron production.

Definition at line 378 of file G4UAtomicDeexcitation.cc.

{
  minElectronEnergy = cut;
}

References minElectronEnergy.

◆ SetCutForSecondaryPhotons()

void G4UAtomicDeexcitation::SetCutForSecondaryPhotons ( G4double cut )

Set threshold energy for fluorescence.

Definition at line 371 of file G4UAtomicDeexcitation.cc.

{
  minGammaEnergy = cut;
}

References minGammaEnergy.

◆ SetDeexcitationActiveRegion()

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

inherited

Definition at line 187 of file G4VAtomDeexcitation.cc.

{
  // no PIXE in parallel world
  if(rname == "DefaultRegionForParallelWorld") { return; }
 
  G4String ss = rname;
  /*  
  G4cout << "### G4VAtomDeexcitation::SetDeexcitationActiveRegion " << ss 
         << "  " << valDeexcitation << "  " << valAuger
         << "  " << valPIXE << G4endl;
  */
  if(ss == "world" || ss == "World" || ss == "WORLD") {
    ss = "DefaultRegionForTheWorld";
  }
  size_t n = deRegions.size();
  for(size_t i=0; i<n; ++i) {
 
    // Region already exist
    if(ss == activeRegions[i]) {
      deRegions[i] = valDeexcitation;
      AugerRegions[i] = valAuger;
      PIXERegions[i] = valPIXE;
      return;  
    }
  }
  // New region
  activeRegions.push_back(ss);
  deRegions.push_back(valDeexcitation);
  AugerRegions.push_back(valAuger);
  PIXERegions.push_back(valPIXE);
 
  // if de-excitation defined for the world volume 
  // it should be active for all G4Regions
  if(ss == "DefaultRegionForTheWorld") {
    G4RegionStore* regions = G4RegionStore::GetInstance();
    G4int nn = regions->size();
    for(G4int i=0; i<nn; ++i) {
      if(ss == (*regions)[i]->GetName()) { continue; }
      SetDeexcitationActiveRegion((*regions)[i]->GetName(), valDeexcitation,
                                  valAuger, valPIXE);
                                  
    }
  }
}

References G4VAtomDeexcitation::activeRegions, G4VAtomDeexcitation::AugerRegions, G4VAtomDeexcitation::deRegions, G4RegionStore::GetInstance(), G4VAtomDeexcitation::GetName(), CLHEP::detail::n, G4InuclParticleNames::nn, G4VAtomDeexcitation::PIXERegions, and G4VAtomDeexcitation::SetDeexcitationActiveRegion().

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

◆ SetFluo()

void G4VAtomDeexcitation::SetFluo ( G4bool val )

inlineinherited

Definition at line 198 of file G4VAtomDeexcitation.hh.

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

References G4VAtomDeexcitation::isActive, and G4VAtomDeexcitation::isActiveLocked.

Referenced by G4EmDNAPhysics_stationary_option2::ConstructProcess().

◆ SetPIXE()

void G4VAtomDeexcitation::SetPIXE ( G4bool val )

inlineinherited

Definition at line 228 of file G4VAtomDeexcitation.hh.

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

References G4VAtomDeexcitation::flagPIXE, and G4VAtomDeexcitation::isPIXELocked.

◆ SetVerboseLevel()

void G4VAtomDeexcitation::SetVerboseLevel ( G4int val )

inlineinherited

Definition at line 249 of file G4VAtomDeexcitation.hh.

{
  verbose = val;
}

References G4VAtomDeexcitation::verbose.

Referenced by G4LossTableManager::ResetParameters().

Field Documentation

◆ activeAugerMedia

std::vector<G4bool> G4VAtomDeexcitation::activeAugerMedia

privateinherited

Definition at line 183 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::CheckAugerActiveRegion(), and G4VAtomDeexcitation::InitialiseAtomicDeexcitation().

◆ activeDeexcitationMedia

std::vector<G4bool> G4VAtomDeexcitation::activeDeexcitationMedia

privateinherited

Definition at line 182 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::CheckDeexcitationActiveRegion(), and G4VAtomDeexcitation::InitialiseAtomicDeexcitation().

◆ activePIXEMedia

std::vector<G4bool> G4VAtomDeexcitation::activePIXEMedia

privateinherited

Definition at line 184 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), and G4VAtomDeexcitation::InitialiseAtomicDeexcitation().

◆ activeRegions

std::vector<G4String> G4VAtomDeexcitation::activeRegions

privateinherited

Definition at line 189 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetDeexcitationActiveRegion().

◆ activeZ

std::vector<G4bool> G4VAtomDeexcitation::activeZ

privateinherited

Definition at line 181 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), G4VAtomDeexcitation::GetListOfActiveAtoms(), and G4VAtomDeexcitation::InitialiseAtomicDeexcitation().

◆ anaPIXEshellCS

G4VhShellCrossSection* G4UAtomicDeexcitation::anaPIXEshellCS

private

Definition at line 128 of file G4UAtomicDeexcitation.hh.

Referenced by G4UAtomicDeexcitation(), GetShellIonisationCrossSectionPerAtom(), InitialiseForNewRun(), and ~G4UAtomicDeexcitation().

◆ AugerRegions

std::vector<G4bool> G4VAtomDeexcitation::AugerRegions

privateinherited

Definition at line 186 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetDeexcitationActiveRegion().

◆ deRegions

std::vector<G4bool> G4VAtomDeexcitation::deRegions

privateinherited

Definition at line 185 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetDeexcitationActiveRegion().

◆ emcorr

G4EmCorrections* G4UAtomicDeexcitation::emcorr

private

Definition at line 130 of file G4UAtomicDeexcitation.hh.

Referenced by G4UAtomicDeexcitation(), and GetShellIonisationCrossSectionPerAtom().

◆ ePIXEshellCS

G4VhShellCrossSection* G4UAtomicDeexcitation::ePIXEshellCS

private

Definition at line 129 of file G4UAtomicDeexcitation.hh.

Referenced by G4UAtomicDeexcitation(), GetShellIonisationCrossSectionPerAtom(), InitialiseForNewRun(), and ~G4UAtomicDeexcitation().

◆ flagAuger

G4bool G4VAtomDeexcitation::flagAuger = false

privateinherited

Definition at line 173 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), G4VAtomDeexcitation::IsAugerActive(), G4VAtomDeexcitation::IsAugerCascadeActive(), and G4VAtomDeexcitation::SetAuger().

◆ flagPIXE

G4bool G4VAtomDeexcitation::flagPIXE = false

privateinherited

Definition at line 174 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), G4VAtomDeexcitation::IsPIXEActive(), and G4VAtomDeexcitation::SetPIXE().

◆ gamma

const G4ParticleDefinition* G4VAtomDeexcitation::gamma

privateinherited

Definition at line 166 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), and G4VAtomDeexcitation::G4VAtomDeexcitation().

◆ ignoreCuts

G4bool G4VAtomDeexcitation::ignoreCuts = false

privateinherited

Definition at line 175 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), G4VAtomDeexcitation::GenerateParticles(), and G4VAtomDeexcitation::InitialiseAtomicDeexcitation().

◆ isActive

G4bool G4VAtomDeexcitation::isActive = false

privateinherited

Definition at line 172 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), G4VAtomDeexcitation::IsFluoActive(), and G4VAtomDeexcitation::SetFluo().

◆ isActiveLocked

G4bool G4VAtomDeexcitation::isActiveLocked = false

privateinherited

Definition at line 177 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetFluo().

◆ isAugerLocked

G4bool G4VAtomDeexcitation::isAugerLocked = false

privateinherited

Definition at line 178 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetAuger().

◆ isPIXELocked

G4bool G4VAtomDeexcitation::isPIXELocked = false

privateinherited

Definition at line 179 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetPIXE().

◆ minElectronEnergy

G4double G4UAtomicDeexcitation::minElectronEnergy

private

Definition at line 142 of file G4UAtomicDeexcitation.hh.

Referenced by GenerateAuger(), GenerateParticles(), and SetCutForAugerElectrons().

◆ minGammaEnergy

G4double G4UAtomicDeexcitation::minGammaEnergy

private

Data member which stores the shells to be filled by the radiative transition

Definition at line 141 of file G4UAtomicDeexcitation.hh.

Referenced by GenerateFluorescence(), GenerateParticles(), and SetCutForSecondaryPhotons().

◆ name

G4String G4VAtomDeexcitation::name

privateinherited

Definition at line 191 of file G4VAtomDeexcitation.hh.

Referenced by source.g4viscp.G4Scene::create_scene(), G4VAtomDeexcitation::GetName(), G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), mcscore.MCParticle::printout(), and source.g4viscp.G4Scene::update_scene().

◆ nCouples

G4int G4VAtomDeexcitation::nCouples = 0

privateinherited

Definition at line 169 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::CheckAugerActiveRegion(), G4VAtomDeexcitation::CheckDeexcitationActiveRegion(), and G4VAtomDeexcitation::InitialiseAtomicDeexcitation().

◆ newShellId

G4int G4UAtomicDeexcitation::newShellId

private

Definition at line 143 of file G4UAtomicDeexcitation.hh.

Referenced by GenerateAuger(), GenerateFluorescence(), and GenerateParticles().

◆ PIXERegions

std::vector<G4bool> G4VAtomDeexcitation::PIXERegions

privateinherited

Definition at line 187 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetDeexcitationActiveRegion().

◆ PIXEshellCS

G4VhShellCrossSection* G4UAtomicDeexcitation::PIXEshellCS

private

Data member for the calculation of the proton and alpha ionisation XS.

Definition at line 127 of file G4UAtomicDeexcitation.hh.

Referenced by G4UAtomicDeexcitation(), GetShellIonisationCrossSectionPerAtom(), InitialiseForNewRun(), and ~G4UAtomicDeexcitation().

◆ theCoupleTable

const G4ProductionCutsTable* G4VAtomDeexcitation::theCoupleTable = nullptr

privateinherited

Definition at line 167 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), G4VAtomDeexcitation::G4VAtomDeexcitation(), G4VAtomDeexcitation::GenerateParticles(), and G4VAtomDeexcitation::InitialiseAtomicDeexcitation().

◆ theElectron

const G4ParticleDefinition* G4UAtomicDeexcitation::theElectron

private

Definition at line 132 of file G4UAtomicDeexcitation.hh.

Referenced by G4UAtomicDeexcitation(), and GetShellIonisationCrossSectionPerAtom().

◆ thePositron

const G4ParticleDefinition* G4UAtomicDeexcitation::thePositron

private

Definition at line 133 of file G4UAtomicDeexcitation.hh.

Referenced by G4UAtomicDeexcitation(), and GetShellIonisationCrossSectionPerAtom().

◆ transitionManager

G4AtomicTransitionManager* G4UAtomicDeexcitation::transitionManager

private

Definition at line 124 of file G4UAtomicDeexcitation.hh.

Referenced by G4UAtomicDeexcitation(), GenerateAuger(), GenerateFluorescence(), GetAtomicShell(), InitialiseForNewRun(), and SelectTypeOfTransition().

◆ vacancyArray

std::vector<int> G4UAtomicDeexcitation::vacancyArray

private

Definition at line 137 of file G4UAtomicDeexcitation.hh.

Referenced by GenerateAuger(), GenerateFluorescence(), and GenerateParticles().

◆ vdyn

std::vector<G4DynamicParticle*> G4VAtomDeexcitation::vdyn

privateinherited

Definition at line 188 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::AlongStepDeexcitation(), and G4VAtomDeexcitation::G4VAtomDeexcitation().

◆ verbose

G4int G4VAtomDeexcitation::verbose = 1

privateinherited

Definition at line 170 of file G4VAtomDeexcitation.hh.

Referenced by G4VAtomDeexcitation::GetVerboseLevel(), G4VAtomDeexcitation::InitialiseAtomicDeexcitation(), and G4VAtomDeexcitation::SetVerboseLevel().

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

source/processes/electromagnetic/lowenergy/include/G4UAtomicDeexcitation.hh
source/processes/electromagnetic/lowenergy/src/G4UAtomicDeexcitation.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4UAtomicDeexcitation() [1/2]

◆ ~G4UAtomicDeexcitation()

◆ G4UAtomicDeexcitation() [2/2]

Member Function Documentation

◆ AlongStepDeexcitation()

◆ CheckAugerActiveRegion()

◆ CheckDeexcitationActiveRegion()

◆ ComputeShellIonisationCrossSectionPerAtom()

◆ GenerateAuger() [1/2]

◆ GenerateAuger() [2/2]

◆ GenerateFluorescence()

◆ GenerateParticles() [1/2]

◆ GenerateParticles() [2/2]

◆ GetAtomicShell()

◆ GetListOfActiveAtoms()

◆ GetName()

◆ GetShellIonisationCrossSectionPerAtom()

◆ GetVerboseLevel()

◆ InitialiseAtomicDeexcitation()

◆ InitialiseForExtraAtom()

◆ InitialiseForNewRun()

◆ IsAugerActive()

◆ IsAugerCascadeActive()

◆ IsFluoActive()

◆ IsPIXEActive()

◆ operator=()

◆ SelectTypeOfTransition()

◆ SetAuger()

◆ SetAugerCascade()

◆ SetCutForAugerElectrons()

◆ SetCutForSecondaryPhotons()

◆ SetDeexcitationActiveRegion()

◆ SetFluo()

◆ SetPIXE()

◆ SetVerboseLevel()

Field Documentation

◆ activeAugerMedia

◆ activeDeexcitationMedia

◆ activePIXEMedia

◆ activeRegions

◆ activeZ

◆ anaPIXEshellCS

◆ AugerRegions

◆ deRegions

◆ emcorr

◆ ePIXEshellCS

◆ flagAuger

◆ flagPIXE

◆ gamma

◆ ignoreCuts

◆ isActive

◆ isActiveLocked

◆ isAugerLocked

◆ isPIXELocked

◆ minElectronEnergy

◆ minGammaEnergy

◆ name

◆ nCouples

◆ newShellId

◆ PIXERegions

◆ PIXEshellCS

◆ theCoupleTable

◆ theElectron

◆ thePositron

◆ transitionManager

◆ vacancyArray

◆ vdyn

◆ verbose