G4PhotonEvaporation Class Reference

#include <G4PhotonEvaporation.hh>

Inheritance diagram for G4PhotonEvaporation:

Collaboration diagram for G4PhotonEvaporation:

Public Member Functions
	G4PhotonEvaporation (G4GammaTransition *ptr=nullptr)
virtual	~G4PhotonEvaporation ()
virtual void	Initialise () final
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus) final
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus) final
G4FragmentVector *	BreakItUp (const G4Fragment &theNucleus)
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus) final
virtual G4double	GetFinalLevelEnergy (G4int Z, G4int A, G4double energy) final
virtual G4double	GetUpperLevelEnergy (G4int Z, G4int A) final
void	SetGammaTransition (G4GammaTransition *)
void	SetMaxHalfLife (G4double)
virtual void	SetICM (G4bool)
virtual void	RDMForced (G4bool)
void	SetVerboseLevel (G4int verbose)
G4int	GetVacantShellNumber () const
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual void	Dump () const
G4double	GetMaxLevelEnergy (G4int Z, G4int A)
G4double	GetNearestLevelEnergy (G4int Z, G4int A, G4double energy)
void	SetPhotonEvaporation (G4VEvaporationChannel *p)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Additional Inherited Members
Protected Attributes inherited from G4VEvaporationChannel
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 63 of file G4PhotonEvaporation.hh.

Constructor & Destructor Documentation

G4PhotonEvaporation::G4PhotonEvaporation ( G4GammaTransition *  ptr = nullptr )

explicit

Definition at line 60 of file G4PhotonEvaporation.cc.

  : fLevelManager(nullptr), fTransition(p), fVerbose(0), fPoints(0), 
    vShellNumber(-1), fIndex(0), fMaxLifeTime(DBL_MAX), 
    fICM(false), fRDM(false), fSampleTime(true), isInitialised(false)
{
  //G4cout << "### New G4PhotonEvaporation() " << this << G4endl;   
  fNuclearLevelData = G4NuclearLevelData::GetInstance(); 
  LevelDensity = 0.125/CLHEP::MeV;
  Tolerance = 20*CLHEP::eV;

  if(!fTransition) { fTransition = new G4GammaTransition(); }

  theA = theZ = fCode = 0;
  fLevelEnergyMax = fStep = fExcEnergy = fProbability = 0.0;

  for(G4int i=0; i<MAXDEPOINT; ++i) { fCummProbability[i] = 0.0; }
  if(0.0f == GREnergy[1]) {
    G4Pow* g4calc = G4Pow::GetInstance();
    static const G4float GRWfactor = 0.3f;
    for (G4int A=1; A<MAXGRDATA; ++A) {
      GREnergy[A] = (G4float)(40.3*CLHEP::MeV/g4calc->powZ(A,0.2));
      GRWidth[A] = GRWfactor*GREnergy[A];
    }
  } 
}

Here is the call graph for this function:

G4PhotonEvaporation::~G4PhotonEvaporation ( )

virtual

Definition at line 86 of file G4PhotonEvaporation.cc.

{ 
  delete fTransition;
}

Member Function Documentation

G4FragmentVector * G4PhotonEvaporation::BreakItUp

(

const G4Fragment &

theNucleus

)

Definition at line 124 of file G4PhotonEvaporation.cc.

{
  if(fVerbose > 0) {
    G4cout << "G4PhotonEvaporation::BreakItUp" << G4endl;
  }
  G4Fragment* aNucleus = new G4Fragment(nucleus);
  G4FragmentVector* products = new G4FragmentVector();
  BreakUpChain(products, aNucleus);
  products->push_back(aNucleus);
  return products;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4PhotonEvaporation::BreakUpChain ( G4FragmentVector *  theResult, G4Fragment *  theNucleus  )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 136 of file G4PhotonEvaporation.cc.

{
  if(!isInitialised) { Initialise(); }
  if(fVerbose > 0) {
    G4cout << "G4PhotonEvaporation::BreakUpChain RDM= " << fRDM << " "
       << *nucleus << G4endl;
  }
  G4Fragment* gamma = nullptr;
  if(fRDM) { fSampleTime = false; }
  else     { fSampleTime = true; }

  do {
    gamma = GenerateGamma(nucleus);
    if(gamma) { 
      products->push_back(gamma);
      if(fVerbose > 0) {
    G4cout << "G4PhotonEvaporation::BreakUpChain: "   
           << *gamma << G4endl;
    G4cout << "   Residual: " << *nucleus << G4endl;
      }
      // for next decays in the chain always sample time
      fSampleTime = true;
    } 
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  } while(gamma);
  return false;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4Fragment * G4PhotonEvaporation::EmittedFragment ( G4Fragment *  theNucleus )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 109 of file G4PhotonEvaporation.cc.

{
  if(fRDM) { fSampleTime = false; }
  else     { fSampleTime = true; }

  G4Fragment* gamma = GenerateGamma(nucleus);
  if(fVerbose > 0) {
    G4cout << "G4PhotonEvaporation::EmittedFragment: RDM= " << fRDM << G4endl; 
    if(gamma) { G4cout << *gamma << G4endl; }
    G4cout << "   Residual: " << *nucleus << G4endl;
  }
  return gamma; 
}

Here is the caller graph for this function:

G4double G4PhotonEvaporation::GetEmissionProbability ( G4Fragment *  theNucleus )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 166 of file G4PhotonEvaporation.cc.

{
  if(!isInitialised) { Initialise(); }
  fProbability = 0.0;
  fExcEnergy = nucleus->GetExcitationEnergy();
  G4int Z = nucleus->GetZ_asInt();
  G4int A = nucleus->GetA_asInt();
  fCode   = 1000*Z + A; 
  if(fVerbose > 1) {
    G4cout << "G4PhotonEvaporation::GetEmissionProbability: Z=" 
       << Z << " A=" << A << " Eexc(MeV)= " << fExcEnergy << G4endl; 
  }
  // ignore gamma de-excitation for exotic fragments 
  // and for very low excitations
  if(0 >= Z || 1 >= A || Z == A || Tolerance >= fExcEnergy)
    { return fProbability; }

  // ignore gamma de-excitation for highly excited levels
  if(A >= MAXGRDATA) { A =  MAXGRDATA-1; }
  //G4cout<<" GREnergy= "<< GREnergy[A]<<" GRWidth= "<<GRWidth[A]<<G4endl; 

  static const G4float GREfactor = 5.0f;
  if(fExcEnergy >= (G4double)(GREfactor*GRWidth[A] + GREnergy[A])) { 
    return fProbability; 
  }
  // probability computed assuming continium transitions
  // VI: continium transition are limited only to final states
  //     below Fermi energy (this approach needs further evaluation)
  G4double emax = std::max(0.0, nucleus->ComputeGroundStateMass(Z, A-1) 
    + CLHEP::neutron_mass_c2 - nucleus->GetGroundStateMass());

  // max energy level for continues transition
  emax = std::min(emax, fExcEnergy);
  const G4double eexcfac = 0.99;
  if(0.0 == emax || fExcEnergy*eexcfac <= emax) { emax = fExcEnergy*eexcfac; }

  fStep = emax;
  static const G4double MaxDeltaEnergy = CLHEP::MeV;
  fPoints = std::min((G4int)(fStep/MaxDeltaEnergy) + 2, MAXDEPOINT);
  fStep /= ((G4double)(fPoints - 1));
  if(fVerbose > 1) {
    G4cout << "Emax= " << emax << " Npoints= " << fPoints 
       << "  Eex= " << fExcEnergy << G4endl;
  }
  // integrate probabilities
  G4double eres = (G4double)GREnergy[A];
  G4double wres = (G4double)GRWidth[A];
  G4double eres2= eres*eres;
  G4double wres2= wres*wres;
  G4double xsqr = std::sqrt(A*LevelDensity*fExcEnergy);

  G4double egam    = fExcEnergy;
  G4double gammaE2 = egam*egam;
  G4double gammaR2 = gammaE2*wres2;
  G4double egdp2   = gammaE2 - eres2;

  G4double p0 = G4Exp(-2.0*xsqr)*gammaR2*gammaE2/(egdp2*egdp2 + gammaR2);
  G4double p1(0.0);

  for(G4int i=1; i<fPoints; ++i) {
    egam -= fStep;
    gammaE2 = egam*egam;
    gammaR2 = gammaE2*wres2;
    egdp2   = gammaE2 - eres2;
    p1 = G4Exp(2.0*(std::sqrt(A*LevelDensity*std::abs(fExcEnergy - egam)) - xsqr))
      *gammaR2*gammaE2/(egdp2*egdp2 + gammaR2);
    fProbability += (p1 + p0);
    fCummProbability[i] = fProbability;
    //G4cout << "Egamma= " << egam << "  Eex= " << fExcEnergy
    //<< "  p0= " << p0 << " p1= " << p1 << " sum= " << fCummProbability[i] <<G4endl;
    p0 = p1;
  }

  static const G4double NormC = 1.25*CLHEP::millibarn
    /(CLHEP::pi2*CLHEP::hbarc*CLHEP::hbarc);
  fProbability *= fStep*NormC*A;
  if(fVerbose > 1) { G4cout << "prob= " << fProbability << G4endl; }
  return fProbability;
}

Here is the call graph for this function:

G4double G4PhotonEvaporation::GetFinalLevelEnergy ( G4int  Z, G4int  A, G4double  energy  )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 247 of file G4PhotonEvaporation.cc.

{
  G4double E = energy;
  InitialiseLevelManager(Z, A);
  if(fLevelManager) { 
    E = (G4double)fLevelManager->NearestLevelEnergy(energy, fIndex); 
    if(E > fLevelEnergyMax + Tolerance) { E = energy; }
  }
  return E;
}

Here is the call graph for this function:

G4double G4PhotonEvaporation::GetUpperLevelEnergy ( G4int  Z, G4int  A  )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 258 of file G4PhotonEvaporation.cc.

{
  InitialiseLevelManager(Z, A);
  return fLevelEnergyMax;
}

G4int G4PhotonEvaporation::GetVacantShellNumber ( ) const

inline

Definition at line 161 of file G4PhotonEvaporation.hh.

{ 
  return vShellNumber;
}

Here is the caller graph for this function:

void G4PhotonEvaporation::Initialise ( )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 91 of file G4PhotonEvaporation.cc.

{
  if(isInitialised) { return; }
  isInitialised = true;

  if(fVerbose > 0) {
    G4cout << "### G4PhotonEvaporation is initialized " << this << G4endl;   
  }
  G4DeexPrecoParameters* param = fNuclearLevelData->GetParameters();
  LevelDensity = param->GetLevelDensity();
  Tolerance = param->GetMinExcitation();
  fMaxLifeTime = param->GetMaxLifeTime();

  fTransition->SetPolarizationFlag(param->CorrelatedGamma());
  fTransition->SetVerbose(fVerbose);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void G4PhotonEvaporation::RDMForced ( G4bool  val )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 448 of file G4PhotonEvaporation.cc.

{
  fRDM = val;
}

Here is the caller graph for this function:

void G4PhotonEvaporation::SetGammaTransition

(

G4GammaTransition *

p

)

Definition at line 435 of file G4PhotonEvaporation.cc.

{
  if(p != fTransition) {
    delete fTransition;
    fTransition = p;
  }
}

void G4PhotonEvaporation::SetICM ( G4bool  val )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 443 of file G4PhotonEvaporation.cc.

{
  fICM = val;
}

Here is the caller graph for this function:

void G4PhotonEvaporation::SetMaxHalfLife

(

G4double

)

Definition at line 432 of file G4PhotonEvaporation.cc.

{}

void G4PhotonEvaporation::SetVerboseLevel ( G4int  verbose )

inline

Definition at line 144 of file G4PhotonEvaporation.hh.

{
  fVerbose = verbose;
}

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/photon_evaporation/include/G4PhotonEvaporation.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4PhotonEvaporation.cc