G4ContinuumGammaTransition Class Reference

#include <G4ContinuumGammaTransition.hh>

Inheritance diagram for G4ContinuumGammaTransition:

Public Member Functions
	G4ContinuumGammaTransition (const G4NuclearLevelManager *levelManager, G4int Z, G4int A, G4double excitation, G4int verbose)
virtual	~G4ContinuumGammaTransition ()
virtual void	SetEnergyFrom (G4double energy)
virtual G4double	GetGammaEnergy ()
virtual G4double	GetGammaCreationTime ()
virtual void	SelectGamma ()

---

Detailed Description

Definition at line 59 of file G4ContinuumGammaTransition.hh.

---

Constructor & Destructor Documentation

G4ContinuumGammaTransition::G4ContinuumGammaTransition	(	const G4NuclearLevelManager *	levelManager,
		G4int	Z,
		G4int	A,
		G4double	excitation,
		G4int	verbose
	)

Definition at line 66 of file G4ContinuumGammaTransition.cc.

References G4VGammaTransition::_verbose, G4NuclearLevel::Energy(), G4Pow::GetInstance(), G4NuclearLevelManager::GetLevel(), G4NuclearLevelManager::MaxLevelEnergy(), G4NuclearLevelManager::MinLevelEnergy(), G4NuclearLevelManager::NumberOfLevels(), and G4Pow::powZ().

                                          :
  _nucleusA(A), _nucleusZ(Z), _excitation(excitation), _levelManager(levelManager) 
{
  G4double eTolerance = 0.;
  G4int lastButOne = _levelManager->NumberOfLevels() - 2;
  if (lastButOne >= 0)
    {
      eTolerance = (_levelManager->MaxLevelEnergy() -
                    _levelManager->GetLevel(lastButOne)->Energy());
      if (eTolerance < 0.) eTolerance = 0.;
    }
  

  _verbose = verbose;
  _eGamma = 0.;
  _gammaCreationTime = 0.;

  _maxLevelE = _levelManager->MaxLevelEnergy() + eTolerance;
  _minLevelE = _levelManager->MinLevelEnergy();

  // Energy range for photon generation; upper limit is defined 5*Gamma(GDR) from GDR peak
  _eMin = 0.001 * MeV;
  // Giant Dipole Resonance energy
  G4double energyGDR = (40.3 / G4Pow::GetInstance()->powZ(_nucleusA,0.2) ) * MeV;
  // Giant Dipole Resonance width
  G4double widthGDR = 0.30 * energyGDR;
  // Extend 
  G4double factor = 5;
  _eMax = energyGDR + factor * widthGDR;
  if (_eMax > excitation) _eMax = _excitation;

}

G4ContinuumGammaTransition::~G4ContinuumGammaTransition

(

)

[virtual]

Definition at line 107 of file G4ContinuumGammaTransition.cc.

{}

---

Member Function Documentation

G4double G4ContinuumGammaTransition::GetGammaCreationTime

(

)

[virtual]

Implements G4VGammaTransition.

Definition at line 167 of file G4ContinuumGammaTransition.cc.

{
  return _gammaCreationTime;
}

G4double G4ContinuumGammaTransition::GetGammaEnergy

(

)

[virtual]

Implements G4VGammaTransition.

Definition at line 162 of file G4ContinuumGammaTransition.cc.

{
  return _eGamma;
}

void G4ContinuumGammaTransition::SelectGamma

(

)

[virtual]

Implements G4VGammaTransition.

Definition at line 110 of file G4ContinuumGammaTransition.cc.

References G4VGammaTransition::_verbose, G4NuclearLevel::Energy(), G4cout, G4endl, and G4NuclearLevelManager::NearestLevel().

{

  _eGamma = 0.;

  G4int nBins = 200;
  G4double sampleArray[200];
  G4int i;
  for (i=0; i<nBins; i++)
    {
      G4double e = _eMin + ( (_eMax - _eMin) / nBins) * i;
      sampleArray[i] = E1Pdf(e);

      if(_verbose > 10)
        G4cout << "*---* G4ContinuumTransition: e = " << e 
               << " pdf = " << sampleArray[i] << G4endl;
    }
  G4RandGeneralTmp randGeneral(sampleArray, nBins);
  G4double random = randGeneral.shoot();
  
  _eGamma = _eMin + (_eMax - _eMin) * random;
  
  G4double finalExcitation = _excitation - _eGamma;
  
  if(_verbose > 10) {
    G4cout << "*---*---* G4ContinuumTransition: eGamma = " << _eGamma
           << "   finalExcitation = " << finalExcitation 
           << " random = " << random << G4endl;
  }
  //  if (finalExcitation < 0)
  if(finalExcitation < _minLevelE/2.)
    {
      _eGamma = _excitation;
      finalExcitation = 0.;
    }
  
  if (finalExcitation < _maxLevelE && finalExcitation > 0.) 
    {
      G4double levelE = _levelManager->NearestLevel(finalExcitation)->Energy();
      G4double diff = finalExcitation - levelE;
      _eGamma = _eGamma + diff;
    }  

  _gammaCreationTime = GammaTime();

  if(_verbose > 10) {
    G4cout << "*---*---* G4ContinuumTransition: _gammaCreationTime = "
           << _gammaCreationTime/second << G4endl;
  }
  return;  
}

void G4ContinuumGammaTransition::SetEnergyFrom

(

G4double

energy

)

[virtual]

Implements G4VGammaTransition.

Definition at line 173 of file G4ContinuumGammaTransition.cc.

{
  if (energy > 0.) _excitation = energy;
}

---

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

• G4ContinuumGammaTransition.hh
• G4ContinuumGammaTransition.cc

---