G4UAtomicDeexcitation Class Reference

#include <G4UAtomicDeexcitation.hh>

Inheritance 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", const G4String &pixename="")
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	SetPIXE (G4bool)
G4bool	IsPIXEActive () const
const G4String &	GetName () const
void	SetPIXECrossSectionModel (const G4String &)
const G4String &	PIXECrossSectionModel () const
void	SetPIXEElectronCrossSectionModel (const G4String &)
const G4String &	PIXEElectronCrossSectionModel () 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=0)=0
virtual G4double	ComputeShellIonisationCrossSectionPerAtom (const G4ParticleDefinition *, G4int Z, G4AtomicShellEnumerator shell, G4double kinE, const G4Material *mat=0)=0
void	AlongStepDeexcitation (std::vector< G4Track * > &tracks, const G4Step &step, G4double &eLoss, G4int coupleIndex)

Detailed Description

Definition at line 62 of file G4UAtomicDeexcitation.hh.

Constructor & Destructor Documentation

G4UAtomicDeexcitation()

G4UAtomicDeexcitation::G4UAtomicDeexcitation

(

)

Definition at line 76 of file G4UAtomicDeexcitation.cc.

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

~G4UAtomicDeexcitation()

G4UAtomicDeexcitation::~G4UAtomicDeexcitation ( )

virtual

Definition at line 91 of file G4UAtomicDeexcitation.cc.

{
  delete PIXEshellCS;
  delete anaPIXEshellCS;
  delete ePIXEshellCS;
}

Member Function Documentation

ComputeShellIonisationCrossSectionPerAtom()

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 389 of file G4UAtomicDeexcitation.cc.

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

GenerateParticles()

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 234 of file G4UAtomicDeexcitation.cc.

{
 
  // Defined initial conditions
  G4int givenShellId = atomicShell->ShellId();
  //  G4cout << "generating particles for vacancy in shellId: " << givenShellId << G4endl;
  minGammaEnergy = gammaCut;
  minElectronEnergy = eCut;
 
  // V.I. short-cut
  if(!IsAugerActive()) {  minElectronEnergy = DBL_MAX; }
 
  // generation secondaries
  G4DynamicParticle* aParticle=0;
  G4int provShellId = 0;
  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 ;} //debug  
          }
        else if ( provShellId == -1)
          {
        aParticle = GenerateAuger(Z, givenShellId);
        //if (aParticle != 0) { G4cout << "****AUGER!****" << G4endl;} //debug
          }
        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)
          {
        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 
 
}

GetAtomicShell()

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 229 of file G4UAtomicDeexcitation.cc.

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

GetShellIonisationCrossSectionPerAtom()

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

virtual

Implements G4VAtomDeexcitation.

Definition at line 325 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;
}

Referenced by ComputeShellIonisationCrossSectionPerAtom().

InitialiseForExtraAtom()

void G4UAtomicDeexcitation::InitialiseForExtraAtom ( G4int  Z )

virtual

Implements G4VAtomDeexcitation.

Definition at line 226 of file G4UAtomicDeexcitation.cc.

{}

InitialiseForNewRun()

void G4UAtomicDeexcitation::InitialiseForNewRun ( )

virtual

Implements G4VAtomDeexcitation.

Definition at line 98 of file G4UAtomicDeexcitation.cc.

{
  if(!IsFluoActive()) { return; }
  transitionManager = G4AtomicTransitionManager::Instance();
  if(IsPIXEActive()) {
    G4cout << G4endl;
    G4cout << "### === G4UAtomicDeexcitation::InitialiseForNewRun()" << G4endl;
    anaPIXEshellCS = new G4teoCrossSection("Analytical");
 
  }
  else  {return;}
  // initializing PIXE x-section name
  // 
  if (PIXECrossSectionModel() == "" ||
      PIXECrossSectionModel() == "Empirical" ||
      PIXECrossSectionModel() == "empirical") 
    {
      SetPIXECrossSectionModel("Empirical");
    }
  else if (PIXECrossSectionModel() == "ECPSSR_Analytical" ||
       PIXECrossSectionModel() == "Analytical" || 
       PIXECrossSectionModel() == "analytical") 
    {
      SetPIXECrossSectionModel("Analytical");
    }
  else if (PIXECrossSectionModel() == "ECPSSR_FormFactor" ||
       PIXECrossSectionModel() == "ECPSSR_Tabulated" ||
       PIXECrossSectionModel() == "Analytical_Tabulated") 
    {
      SetPIXECrossSectionModel("ECPSSR_FormFactor");
    }
  else 
    {
      G4cout << "### G4UAtomicDeexcitation::InitialiseForNewRun WARNING "
         << G4endl;
      G4cout << "    PIXE cross section name " << PIXECrossSectionModel()
         << " is unknown, Analytical cross section will be used" << G4endl; 
      SetPIXECrossSectionModel("Analytical");
    }
    
  // Check if old model should be deleted 
  if(PIXEshellCS) 
    {
      if(PIXECrossSectionModel() != PIXEshellCS->GetName()) 
    {
      delete PIXEshellCS;
          PIXEshellCS = 0;
    }
    }
 
  // Instantiate empirical model
  if(!PIXEshellCS) {
    if (PIXECrossSectionModel() == "Empirical")
      {
    PIXEshellCS = new G4empCrossSection("Empirical");
      }
 
    if (PIXECrossSectionModel() == "ECPSSR_FormFactor")
      {
    PIXEshellCS = new G4teoCrossSection("ECPSSR_FormFactor");
      }
  }
 
  // Electron cross section
  // initializing PIXE x-section name
  // 
  if (PIXEElectronCrossSectionModel() == "" ||
      PIXEElectronCrossSectionModel() == "Livermore")
    {
      SetPIXEElectronCrossSectionModel("Livermore");
    }
  else if (PIXEElectronCrossSectionModel() == "ProtonAnalytical" ||
       PIXEElectronCrossSectionModel() == "Analytical" || 
       PIXEElectronCrossSectionModel() == "analytical") 
    {
      SetPIXEElectronCrossSectionModel("ProtonAnalytical");
    }
  else if (PIXEElectronCrossSectionModel() == "ProtonEmpirical" ||
       PIXEElectronCrossSectionModel() == "Empirical" || 
       PIXEElectronCrossSectionModel() == "empirical") 
    {
      SetPIXEElectronCrossSectionModel("ProtonEmpirical");
    }
  else if (PIXEElectronCrossSectionModel() == "Penelope")
    SetPIXEElectronCrossSectionModel("Penelope");
  else 
    {
      G4cout << "### G4UAtomicDeexcitation::InitialiseForNewRun WARNING "
         << G4endl;
      G4cout << "    PIXE e- cross section name " << PIXEElectronCrossSectionModel()
         << " is unknown, PIXE is disabled" << G4endl; 
      SetPIXEElectronCrossSectionModel("Livermore");
    }
    
  // Check if old model should be deleted 
  if(ePIXEshellCS) 
    {
      if(PIXEElectronCrossSectionModel() != ePIXEshellCS->GetName()) 
    {
      delete ePIXEshellCS;
          ePIXEshellCS = 0;
    }
    }
 
  // Instantiate empirical model
  if(!ePIXEshellCS) 
    {
      if(PIXEElectronCrossSectionModel() == "Empirical")
    {
      ePIXEshellCS = new G4empCrossSection("Empirical");
    }
 
      else if(PIXEElectronCrossSectionModel() == "Analytical") 
    {
      ePIXEshellCS = new G4teoCrossSection("Analytical");
    }
 
      else if(PIXEElectronCrossSectionModel() == "Livermore")
    {
      ePIXEshellCS = new G4LivermoreIonisationCrossSection();
    }
      else if (PIXEElectronCrossSectionModel() == "Penelope")
    {
      ePIXEshellCS = new G4PenelopeIonisationCrossSection();
    }
    } 
}

SetCutForAugerElectrons()

void G4UAtomicDeexcitation::SetCutForAugerElectrons

(

G4double

cut

)

Definition at line 383 of file G4UAtomicDeexcitation.cc.

{
  minElectronEnergy = cut;
}

SetCutForSecondaryPhotons()

void G4UAtomicDeexcitation::SetCutForSecondaryPhotons

(

G4double

cut

)

Definition at line 378 of file G4UAtomicDeexcitation.cc.

{
  minGammaEnergy = cut;
}

---

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

• G4UAtomicDeexcitation.hh
• G4UAtomicDeexcitation.cc