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 ()

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	(	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(), python.hepunit::MeV, 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 183 of file G4ContinuumGammaTransition.cc.

{
  return _gammaCreationTime;
}

G4double G4ContinuumGammaTransition::GetGammaEnergy ( )

virtual

Implements G4VGammaTransition.

Definition at line 178 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, G4UniformRand, G4NuclearLevelManager::NearestLevel(), and python.hepunit::second.

{

  _eGamma = 0.;

  G4int nBins = 100;
  G4double sampleArray[101];
  sampleArray[0] = 0.0;
  G4int i;
  G4double del = (_eMax - _eMin) / G4double(nBins);
  G4double sum = 0;
  G4double w1 = E1Pdf(_eMin);
  G4double w2;
  //G4cout << _eMin << "  " << _eMax << "  " << del << G4endl;
  for (i=1; i<=nBins; i++) {
    G4double e = _eMin + del * i;
    w2 = E1Pdf(e);
    sum += 0.5*(w1 + w2);
    w1 = w2;
    sampleArray[i] = sum;
    if(_verbose > 10) {
      G4cout << "*---* G4ContinuumTransition: e = " << e 
         << " pdf = " << sampleArray[i] << G4endl;
    }
  }
  sum *= G4UniformRand();
  _eGamma = _eMax;
  for (i=1; i<=nBins; i++) {
    if(sum <= sampleArray[i]) {
      _eGamma = _eMin + del * i;
      G4double w = sampleArray[i] - sampleArray[i-1];
      //G4cout << _eGamma << "  " << w << G4endl;
      if(w != 0.0) {
    _eGamma -= (sampleArray[i] - sum)*del/w;
      }
      break;
    }
  }   
  
  G4double finalExcitation = _excitation - _eGamma;
  
  if(_verbose > 10) {
    G4cout << "*---*---* G4ContinuumTransition: eGamma = " << _eGamma
       << "   finalExcitation = " << finalExcitation << 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 189 of file G4ContinuumGammaTransition.cc.

References energy().

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

---

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

• G4ContinuumGammaTransition.hh
• G4ContinuumGammaTransition.cc