#include <G4ContinuumGammaTransition.hh>

Inheritance diagram for G4ContinuumGammaTransition:

Collaboration diagram for G4ContinuumGammaTransition:

Public Member Functions
	G4ContinuumGammaTransition (const G4NuclearLevelManager *manager, G4int Z, G4int A, G4double exc, G4int ver)

virtual	~G4ContinuumGammaTransition ()

virtual void	SetEnergyFrom (G4double energy)

virtual G4double	GetGammaEnergy ()

virtual G4double	GetGammaCreationTime ()

virtual void	SelectGamma ()

void	Update (const G4NuclearLevelManager *manager, G4int Z, G4int A, G4double exc)

Public Member Functions inherited from G4VGammaTransition
	G4VGammaTransition ()

virtual	~G4VGammaTransition ()

Private Member Functions
G4double	E1Pdf (G4double energy)

G4double	GammaTime ()

Private Attributes
G4int	nucleusA

G4int	nucleusZ

G4int	nBins

G4int	verbose

G4double	eMin

G4double	eMax

G4double	minLevelE

G4double	excitation

G4double	eGamma

G4double	gammaCreationTime

G4double	energyGDR2

G4double	widthGDR

G4double	widthGDR2

std::vector< G4double >	sampleArray

const G4NuclearLevelManager *	levelManager

G4Pow *	g4pow

G4ConstantLevelDensityParameter	ldPar

Detailed Description

Definition at line 62 of file G4ContinuumGammaTransition.hh.

Constructor & Destructor Documentation

◆ G4ContinuumGammaTransition()

G4ContinuumGammaTransition::G4ContinuumGammaTransition	(	const G4NuclearLevelManager *	manager,
		G4int	Z,
		G4int	A,
		G4double	exc,
		G4int	ver
	)

Definition at line 66 of file G4ContinuumGammaTransition.cc.

 {
   nBins = 100;
   verbose = verb;
   eMin = keV;
   eMax = minLevelE = eGamma = gammaCreationTime = 0.0;
   sampleArray.resize(nBins+1,0.0);
   g4pow = G4Pow::GetInstance();
   Update(manager, Z, A, exc);
 }

Here is the call graph for this function:

◆ ~G4ContinuumGammaTransition()

G4ContinuumGammaTransition::~G4ContinuumGammaTransition ( )

virtual

Definition at line 79 of file G4ContinuumGammaTransition.cc.

80 {}

Member Function Documentation

◆ E1Pdf()

G4double G4ContinuumGammaTransition::E1Pdf ( G4double energy )

private

Definition at line 182 of file G4ContinuumGammaTransition.cc.

 {
   G4double theProb = 0.0;
   G4double U = excitation - e;
 
   if(U < 0.0) { return theProb; }
 
   G4double aLevelDensityParam = 
     ldPar.LevelDensityParameter(nucleusA,nucleusZ,excitation);
 
   //G4double levelDensBef = G4Exp(2.0*std::sqrt(aLevelDensityParam*excitation));
   //G4double levelDensAft = G4Exp(2.0*std::sqrt(aLevelDensityParam*(excitation - e)));
   G4double coeff = G4Exp(2.0*(std::sqrt(aLevelDensityParam*U) 
                   - std::sqrt(aLevelDensityParam*excitation)));
 
   //if(verbose > 20)
   //  G4cout << nucleusA << " LevelDensityParameter = " <<  aLevelDensityParam
   //       << " Bef Aft " << levelDensBef << " " << levelDensAft << G4endl;
   
   // Now form the probability density
   // Define constants for the photoabsorption cross-section (the reverse
   // process of our de-excitation)
 
   G4double sigma0 = 2.5 * nucleusA;  
  
   G4double e2 = e*e;
   G4double numerator = sigma0 * e2 * widthGDR2;
   G4double denominator = (e2 - energyGDR2)* (e2 - energyGDR2) + widthGDR2*e2;
   //  if (denominator < 1.0e-9) denominator = 1.0e-9;
 
   G4double sigmaAbs = numerator/denominator; 
 
   if(verbose > 2) {
     G4cout << "E_GDR(MeV)= " << std::sqrt(energyGDR2) << " W_GDR(MeV)= " << widthGDR 
        << " sigAbs= " << sigmaAbs  
        << " E(MeV)= " << e << " coeff= " << coeff
        << G4endl;
   }
 
   theProb = sigmaAbs * e2 * coeff;
 
   return theProb;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GammaTime()

G4double G4ContinuumGammaTransition::GammaTime ( )

private

Definition at line 226 of file G4ContinuumGammaTransition.cc.

 {
   G4double tau = hbar_Planck/widthGDR;
   G4double creationTime = -tau*G4Log(G4UniformRand());
 
   return creationTime;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetGammaCreationTime()

G4double G4ContinuumGammaTransition::GetGammaCreationTime ( )

virtual

Implements G4VGammaTransition.

Definition at line 172 of file G4ContinuumGammaTransition.cc.

 {
   return gammaCreationTime;
 }

◆ GetGammaEnergy()

G4double G4ContinuumGammaTransition::GetGammaEnergy ( )

virtual

Implements G4VGammaTransition.

Definition at line 167 of file G4ContinuumGammaTransition.cc.

 {
   return eGamma;
 }

◆ SelectGamma()

void G4ContinuumGammaTransition::SelectGamma ( )

virtual

Implements G4VGammaTransition.

Definition at line 109 of file G4ContinuumGammaTransition.cc.

 {
   eGamma = 0.;
   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 > 1) {
       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 > 1) {
     G4cout << "*---*---* G4ContinuumTransition: eGamma = " << eGamma
        << "   finalExcitation = " << finalExcitation << G4endl;
   }
 
   if(finalExcitation <= minLevelE) {
     eGamma = excitation;
 
   } else {
     finalExcitation = levelManager->NearestLevel(finalExcitation)->Energy();
     eGamma = excitation - finalExcitation;
   }  
 
   gammaCreationTime = GammaTime();
 
   if(verbose > 1) {
     G4cout << "*---*---* G4ContinuumTransition: gammaCreationTime = "
        << gammaCreationTime/second << G4endl;
   }
 }

Here is the call graph for this function:

◆ SetEnergyFrom()

void G4ContinuumGammaTransition::SetEnergyFrom ( G4double energy )

virtual

Implements G4VGammaTransition.

Definition at line 177 of file G4ContinuumGammaTransition.cc.

 {
   excitation = energy;
 }

◆ Update()

void G4ContinuumGammaTransition::Update	(	const G4NuclearLevelManager *	manager,
		G4int	Z,
		G4int	A,
		G4double	exc
	)

Definition at line 82 of file G4ContinuumGammaTransition.cc.

 {
   levelManager = manager;
   nucleusZ = Z;
   nucleusA = A;
   excitation = exc;
   eGamma = 0.;
   gammaCreationTime = 0.;
 
   minLevelE = DBL_MAX;
   if(levelManager) { 
     minLevelE = std::max(0.5*levelManager->MinLevelEnergy(), tolerance);
   }
   // Energy range for photon generation; 
   // upper limit is defined 5*Gamma(GDR) from GDR peak
   // Giant Dipole Resonance energy
   G4double energyGDR = (40.3 / g4pow->powZ(nucleusA,0.2) ) * MeV;
   energyGDR2 = energyGDR*energyGDR;
   // Giant Dipole Resonance width
   widthGDR = 0.30 * energyGDR;
   widthGDR2 = widthGDR*widthGDR;
   // Extend 
   eMax = energyGDR + factor * widthGDR;
   if (eMax > excitation) { eMax = excitation; }
 }