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 ()
Public Member Functions inherited from G4VGammaTransition
	G4VGammaTransition ()
virtual	~G4VGammaTransition ()
virtual void	SelectGamma ()=0
virtual G4double	GetGammaEnergy ()=0
virtual G4double	GetGammaCreationTime ()=0
virtual void	SetEnergyFrom (G4double energy)=0

Additional Inherited Members
Protected Attributes inherited from G4VGammaTransition
G4int	_verbose

Detailed Description

Definition at line 59 of file G4ContinuumGammaTransition.hh.

Constructor & Destructor Documentation

G4ContinuumGammaTransition()

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

Definition at line 66 of file G4ContinuumGammaTransition.cc.

                              :
  _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::~G4ContinuumGammaTransition ( )

virtual

Definition at line 107 of file G4ContinuumGammaTransition.cc.

{}

Member Function Documentation

GetGammaCreationTime()

G4double G4ContinuumGammaTransition::GetGammaCreationTime ( )

virtual

Implements G4VGammaTransition.

Definition at line 167 of file G4ContinuumGammaTransition.cc.

{
  return _gammaCreationTime;
}

GetGammaEnergy()

G4double G4ContinuumGammaTransition::GetGammaEnergy ( )

virtual

Implements G4VGammaTransition.

Definition at line 162 of file G4ContinuumGammaTransition.cc.

{
  return _eGamma;
}

SelectGamma()

void G4ContinuumGammaTransition::SelectGamma ( )

virtual

Implements G4VGammaTransition.

Definition at line 110 of file G4ContinuumGammaTransition.cc.

{
 
  _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;  
}

SetEnergyFrom()

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