#include <G4UAtomicDeexcitation.hh>

Inheritance diagram for G4UAtomicDeexcitation:

Collaboration diagram for G4UAtomicDeexcitation:

Public Member Functions
	G4UAtomicDeexcitation ()

virtual	~G4UAtomicDeexcitation ()

virtual void	InitialiseForNewRun ()

virtual void	InitialiseForExtraAtom (G4int Z)

void	SetCutForSecondaryPhotons (G4double cut)

void	SetCutForAugerElectrons (G4double cut)

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

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

virtual G4double	GetShellIonisationCrossSectionPerAtom (const G4ParticleDefinition , G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material mat=0)

virtual G4double	ComputeShellIonisationCrossSectionPerAtom (const G4ParticleDefinition , G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material mat=0)

Public Member Functions inherited from G4VAtomDeexcitation
	G4VAtomDeexcitation (const G4String &modname="Deexcitation")

virtual	~G4VAtomDeexcitation ()

void	InitialiseAtomicDeexcitation ()

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)

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

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

Detailed Description

Definition at line 61 of file G4UAtomicDeexcitation.hh.

Constructor & Destructor Documentation

G4UAtomicDeexcitation::G4UAtomicDeexcitation ( )

Definition at line 76 of file G4UAtomicDeexcitation.cc.

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

Here is the call graph for this function:

G4UAtomicDeexcitation::~G4UAtomicDeexcitation ( )

virtual

Definition at line 90 of file G4UAtomicDeexcitation.cc.

 {
   delete anaPIXEshellCS;
   delete PIXEshellCS;
   delete ePIXEshellCS;
 }

Member Function Documentation

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 404 of file G4UAtomicDeexcitation.cc.

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

Here is the call graph for this function:

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 176 of file G4UAtomicDeexcitation.cc.

 {
 
   // Defined initial conditions
   G4int givenShellId = atomicShell->ShellId();
   //G4cout << "generating particles for vacancy in shellId: " 
   // << givenShellId << G4endl; // debug
   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);
         //if (aParticle != 0) { 
         // G4cout << "****FLUO!_1**** " 
                 // << aParticle->GetParticleDefinition()->GetParticleType() 
                 // << " " << aParticle->GetKineticEnergy()/keV << G4endl ;}
           }
         else if ( provShellId == -1)
           {
         // G4cout << "Try to generate Auger 1" << G4endl; 
         aParticle = GenerateAuger(Z, givenShellId);
         // if (aParticle != 0) { G4cout << "****AUGER!****" << G4endl;} 
           }
         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);
         //if (aParticle != 0) { G4cout << "****FLUO!_2****" << aParticle->GetParticleDefinition()->GetParticleType() << " " << aParticle->GetKineticEnergy()/keV << G4endl;} //debug
           }
         else if ( provShellId == -1)
           {
         //      G4cout << "Try to generate Auger 2" << G4endl; //debug
         aParticle = GenerateAuger(Z, newShellId);
         //      if (aParticle != 0) { G4cout << "****AUGER!****" << G4endl;} //debug
           }
         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");
     }
   
   //G4cout << "---------FATTO!---------" << G4endl; //debug 
 
 } // 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);
 //      if (aParticle != 0) {
 //       G4cout << "****FLUO!_1**** " 
 //                 << aParticle->GetParticleDefinition()->GetParticleType() 
 //                 << " " << aParticle->GetKineticEnergy()/keV << G4endl ;}
     }
     else if(provShellId == -1){
         aParticle = GenerateAuger(Z, givenShellId);
 //        if (aParticle != 0) { G4cout << "****AUGER!****" << 
 //            aParticle->GetParticleDefinition()->GetParticleType()  
 //                << " " << aParticle->GetKineticEnergy()/keV << G4endl ; }
 //        else G4cout<<G4endl; 
     }
     //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
 
 //ENDSI
 }

Here is the call graph for this function:

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 171 of file G4UAtomicDeexcitation.cc.

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

Here is the call graph for this function:

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 337 of file G4UAtomicDeexcitation.cc.

 {
   // we must put a control on the shell that are passed: 
   // some shells should not pass (line "0" or "2")
   //G4cout << pdef->GetParticleName() << " Z= " << Z << " Shell= " << shellEnum
   //     << "  E= " << kineticEnergy << G4endl;
 
   // 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; }
 
   //G4cout << pdef->GetParticleName() << " Z= " << Z << "  " << PIXEshellCS 
   //     << "  " << ePIXEshellCS << G4endl;
   // 
   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;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4UAtomicDeexcitation::InitialiseForExtraAtom ( G4int Z )

virtual

Implements G4VAtomDeexcitation.

Definition at line 167 of file G4UAtomicDeexcitation.cc.

168 {}

void G4UAtomicDeexcitation::InitialiseForNewRun ( )

virtual

Implements G4VAtomDeexcitation.

Definition at line 97 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();
     
   //G4cout << namePIXExsModel << "  " << namePIXExsElectronModel << G4endl;
 
   // 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 == "Empirical")
       {
     PIXEshellCS = new G4empCrossSection(namePIXExsModel);
       }
   }
   //G4cout << "PIXE is initialised" << G4endl;
 
   // 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();
     }
     } 
   //G4cout << "ePIXE is initialised" << G4endl;
 }

Here is the call graph for this function:

void G4UAtomicDeexcitation::SetCutForAugerElectrons ( G4double cut )

Definition at line 398 of file G4UAtomicDeexcitation.cc.

 {
   minElectronEnergy = cut;
 }

void G4UAtomicDeexcitation::SetCutForSecondaryPhotons ( G4double cut )

Definition at line 393 of file G4UAtomicDeexcitation.cc.

 {
   minGammaEnergy = cut;
 }

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

source/geant4.10.03.p02/source/processes/electromagnetic/lowenergy/include/G4UAtomicDeexcitation.hh
source/geant4.10.03.p02/source/processes/electromagnetic/lowenergy/src/G4UAtomicDeexcitation.cc

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation