G4DiscreteGammaTransition Class Reference

#include <G4DiscreteGammaTransition.hh>

Inheritance diagram for G4DiscreteGammaTransition:

Collaboration diagram for G4DiscreteGammaTransition:

Public Member Functions
	G4DiscreteGammaTransition (const G4NuclearLevel *level, G4int Z, G4int ver)
virtual	~G4DiscreteGammaTransition ()
virtual void	SetEnergyFrom (G4double energy)
virtual G4double	GetGammaEnergy ()
virtual G4double	GetGammaCreationTime ()
virtual void	SelectGamma ()
void	SetICM (G4bool ic)
G4bool	GetICM () const
G4double	GetBondEnergy ()
G4int	GetOrbitNumber ()
G4bool	IsAGamma ()
void	Update (const G4NuclearLevel *level, G4int Z)
Public Member Functions inherited from G4VGammaTransition
	G4VGammaTransition ()
virtual	~G4VGammaTransition ()

Private Attributes
G4int	nucleusZ
G4int	orbitE
G4int	verbose
G4double	bondE
G4double	gammaEnergy
G4double	excitation
G4double	gammaCreationTime
G4bool	aGamma
G4bool	icm
const G4NuclearLevel *	aLevel

Detailed Description

Definition at line 78 of file G4DiscreteGammaTransition.hh.

Constructor & Destructor Documentation

G4DiscreteGammaTransition()

G4DiscreteGammaTransition::G4DiscreteGammaTransition	(	const G4NuclearLevel *	level,
		G4int	Z,
		G4int	ver
	)

Definition at line 82 of file G4DiscreteGammaTransition.cc.

  : orbitE(-1), bondE(0.),  gammaEnergy(0.),  excitation(0.), 
    gammaCreationTime(0.), aGamma(true), icm(false)
{
  verbose = verb;
  Update(level, Z);
}

Here is the call graph for this function:

~G4DiscreteGammaTransition()

G4DiscreteGammaTransition::~G4DiscreteGammaTransition ( )

virtual

Definition at line 91 of file G4DiscreteGammaTransition.cc.

{}

Member Function Documentation

GetBondEnergy()

G4double G4DiscreteGammaTransition::GetBondEnergy ( )

inline

Definition at line 94 of file G4DiscreteGammaTransition.hh.

{ return bondE;};

Here is the caller graph for this function:

GetGammaCreationTime()

G4double G4DiscreteGammaTransition::GetGammaCreationTime ( )

virtual

Implements G4VGammaTransition.

Definition at line 241 of file G4DiscreteGammaTransition.cc.

{
  return gammaCreationTime;
}

GetGammaEnergy()

G4double G4DiscreteGammaTransition::GetGammaEnergy ( )

virtual

Implements G4VGammaTransition.

Definition at line 236 of file G4DiscreteGammaTransition.cc.

{
  return gammaEnergy;
}

GetICM()

G4bool G4DiscreteGammaTransition::GetICM ( ) const

inline

Definition at line 93 of file G4DiscreteGammaTransition.hh.

{ return icm;};

GetOrbitNumber()

G4int G4DiscreteGammaTransition::GetOrbitNumber ( )

inline

Definition at line 95 of file G4DiscreteGammaTransition.hh.

{ return orbitE;};

Here is the caller graph for this function:

IsAGamma()

G4bool G4DiscreteGammaTransition::IsAGamma ( )

inline

Definition at line 96 of file G4DiscreteGammaTransition.hh.

{ return aGamma;};

Here is the caller graph for this function:

SelectGamma()

void G4DiscreteGammaTransition::SelectGamma ( )

virtual

Implements G4VGammaTransition.

Definition at line 94 of file G4DiscreteGammaTransition.cc.

{
  // default gamma 
  aGamma = true;    
  gammaEnergy = 0.;
  
  G4int nGammas = aLevel->NumberOfGammas();
  if (nGammas > 0) {
    G4int iGamma = 0;
    if(1 < nGammas) {
      G4double random = G4UniformRand();
      
      //G4cout << "G4DiscreteGammaTransition::SelectGamma  N= " 
      //         << nGammas << " rand= " << random << G4endl;
      for(iGamma=0; iGamma<nGammas; ++iGamma) {
    //G4cout << iGamma << "  prob= " 
    //   << (aLevel->GammaCumulativeProbabilities())[iGamma] << G4endl;
    if(random <= (aLevel->GammaCumulativeProbabilities())[iGamma])
      { break; }
      }
    }
    /*
    G4cout << "Elevel(MeV)= " << aLevel->Energy()
       << " Etran(MeV)= " << (aLevel->GammaEnergies())[iGamma]
       << " Eexc(MeV)= " << excitation << G4endl;
    */
    // VI 2014: initial excitation energy may be not exactly energy of the level
    //          final excitation energy is always energy of some level or zero
    //          transition to the ground state should be always equal to
    //          the excitation energy 
    gammaEnergy = (aLevel->GammaEnergies())[iGamma] 
      + excitation - aLevel->Energy();

    // this check is needed to remove cases when nucleaus is left in 
    // slightly excited state which will require very low energy
    // gamma emission 
    if(excitation <= gammaEnergy + tolerance) { gammaEnergy = excitation; }
    //JMQ: 
    //1)If chosen gamma energy is close enough to excitation energy, 
    //  the later is used instead for gamma dacey to gs (it guarantees 
    //  energy conservation)
    //2)For energy conservation, level energy differences instead of  
    //  tabulated gamma energies must be used (origin of final fake photons)
      
    // VI: remove fake photons - applied only for the last transition
    //     do not applied on each transition
    //if(std::fabs(excitation - gammaEnergy) < tolerance) { 
    //  gammaEnergy = excitation;
    //}

    //  JMQ: Warning: the following check is needed to avoid loops:
    //  Due essentially to missing nuclear levels in data files, it is
    //  possible that gammaEnergy is so low as the nucleus doesn't change
    //  its level after the transition.
    //  When such case is found, force the full deexcitation of the nucleus.
    //
    //    NOTE: you should force the transition to the next lower level,
    //          but this change needs a more complex revision of actual 
    //          design.
    //          I leave this for a later revision.

    // VI: the check is needed to remove very low-energy gamma
    if (gammaEnergy <= tolerance) { gammaEnergy = excitation; }
    /*    
    G4cout << "G4DiscreteGammaTransition::SelectGamma: " << gammaEnergy 
       << " Eexc= " << excitation
       << " icm: " << icm << G4endl;
    */
    // now decide whether Internal Coversion electron should be emitted instead
    if (icm) {
      G4double random = G4UniformRand();
      if ( random <= (aLevel->TotalConvertionProbabilities())[iGamma]
       *(aLevel->GammaWeights())[iGamma]
       /((aLevel->TotalConvertionProbabilities())[iGamma]
         *(aLevel->GammaWeights())[iGamma]
         +(aLevel->GammaWeights())[iGamma])) 
    {
      G4int iShell = 9;
      random = G4UniformRand() ;
      if ( random <= (aLevel->KConvertionProbabilities())[iGamma]) 
        { iShell = 0;}
      else if ( random <= (aLevel->L1ConvertionProbabilities())[iGamma]) 
        { iShell = 1;}
      else if ( random <= (aLevel->L2ConvertionProbabilities())[iGamma]) 
        { iShell = 2;}
      else if ( random <= (aLevel->L3ConvertionProbabilities())[iGamma]) 
        { iShell = 3;}  
      else if ( random <= (aLevel->M1ConvertionProbabilities())[iGamma]) 
        { iShell = 4;}
      else if ( random <= (aLevel->M2ConvertionProbabilities())[iGamma]) 
        { iShell = 5;}
      else if ( random <= (aLevel->M3ConvertionProbabilities())[iGamma]) 
        { iShell = 6;}
      else if ( random <= (aLevel->M4ConvertionProbabilities())[iGamma]) 
        { iShell = 7;}
      else if ( random <= (aLevel->M5ConvertionProbabilities())[iGamma]) 
        { iShell = 8;}
      // the following is needed to match the ishell to that used in  
      // G4AtomicShells
      if ( iShell == 9) {
        if ( (nucleusZ < 28) && (nucleusZ > 20)) {
          iShell--;
        } else if ( nucleusZ == 20 || nucleusZ == 19 ) {
          iShell = iShell -2;
        }
      }
      //L.Desorgher 02/11/2011
      //Atomic shell information is available in Geant4 only up top Z=100
      //To extend the photo evaporation code to Z>100  the call 
      // to G4AtomicShells::GetBindingEnergy should be forbidden for Z>100
      bondE = 0.;
      if (nucleusZ <=100) {
        bondE = G4AtomicShells::GetBindingEnergy(nucleusZ, iShell);
      }
      if (verbose > 1) {
        G4cout << "G4DiscreteGammaTransition: nucleusZ = " <<nucleusZ 
           << " , iShell = " << iShell  
           << " , Shell binding energy = " << bondE/keV
           << " keV " << G4endl;
      }

      // last check on energy
      if(gammaEnergy >  bondE + tolerance) {
        orbitE = iShell;      
        aGamma = false ;   // emitted is not a gamma now 
        gammaEnergy -= bondE; 
      }
      //G4cout << "gammaEnergy = " << gammaEnergy << G4endl;
    }
    }
      
    G4double tau = aLevel->HalfLife() / G4Pow::GetInstance()->logZ(2);

    //09.05.2010 VI rewrite samling of decay time 
    //              assuming ordinary exponential low
    gammaCreationTime = 0.;      
    if(tau > 0.0) {  gammaCreationTime = -tau*G4Log(G4UniformRand()); }
  }
  //G4cout << "G4DiscreteGammaTransition end nGamma= " << nGammas
  //     << "  Egamma= " << gammaEnergy << G4endl;
}

Here is the call graph for this function:

SetEnergyFrom()

void G4DiscreteGammaTransition::SetEnergyFrom ( G4double  energy )

virtual

Implements G4VGammaTransition.

Definition at line 246 of file G4DiscreteGammaTransition.cc.

{
  excitation = energy;
}

Here is the caller graph for this function:

SetICM()

void G4DiscreteGammaTransition::SetICM ( G4bool  ic )

inline

Definition at line 92 of file G4DiscreteGammaTransition.hh.

{ icm = ic; };

Here is the caller graph for this function:

Update()

void G4DiscreteGammaTransition::Update ( const G4NuclearLevel *  level, G4int  Z  )

inline

Definition at line 98 of file G4DiscreteGammaTransition.hh.

  {
    aLevel = level; 
    nucleusZ = Z; 
  }

Here is the caller graph for this function:

Member Data Documentation

aGamma

G4bool G4DiscreteGammaTransition::aGamma

private

Definition at line 115 of file G4DiscreteGammaTransition.hh.

aLevel

const G4NuclearLevel* G4DiscreteGammaTransition::aLevel

private

Definition at line 118 of file G4DiscreteGammaTransition.hh.

bondE

G4double G4DiscreteGammaTransition::bondE

private

Definition at line 110 of file G4DiscreteGammaTransition.hh.

excitation

G4double G4DiscreteGammaTransition::excitation

private

Definition at line 112 of file G4DiscreteGammaTransition.hh.

gammaCreationTime

G4double G4DiscreteGammaTransition::gammaCreationTime

private

Definition at line 113 of file G4DiscreteGammaTransition.hh.

gammaEnergy

G4double G4DiscreteGammaTransition::gammaEnergy

private

Definition at line 111 of file G4DiscreteGammaTransition.hh.

icm

G4bool G4DiscreteGammaTransition::icm

private

Definition at line 116 of file G4DiscreteGammaTransition.hh.

nucleusZ

G4int G4DiscreteGammaTransition::nucleusZ

private

Definition at line 106 of file G4DiscreteGammaTransition.hh.

orbitE

G4int G4DiscreteGammaTransition::orbitE

private

Definition at line 107 of file G4DiscreteGammaTransition.hh.

verbose

G4int G4DiscreteGammaTransition::verbose

private

Definition at line 108 of file G4DiscreteGammaTransition.hh.

---

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

• G4DiscreteGammaTransition.hh
• G4DiscreteGammaTransition.cc