#include <G4PhotonEvaporation.hh>

Inheritance diagram for G4PhotonEvaporation:

Public Member Functions
G4FragmentVector *	BreakItUp (const G4Fragment &theNucleus)

virtual G4bool	BreakUpChain (G4FragmentVector theResult, G4Fragment theNucleus) final

G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

virtual void	Dump () const

virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus) final

	G4PhotonEvaporation (const G4PhotonEvaporation &right)=delete

	G4PhotonEvaporation (G4GammaTransition *ptr=nullptr)

virtual G4double	GetEmissionProbability (G4Fragment *theNucleus) final

virtual G4double	GetFinalLevelEnergy (G4int Z, G4int A, G4double energy) final

virtual G4double	GetLifeTime (G4Fragment *theNucleus)

virtual G4double	GetUpperLevelEnergy (G4int Z, G4int A) final

G4int	GetVacantShellNumber () const

virtual void	Initialise () final

const G4PhotonEvaporation &	operator= (const G4PhotonEvaporation &right)=delete

virtual void	RDMForced (G4bool)

void	SetGammaTransition (G4GammaTransition *)

virtual void	SetICM (G4bool)

void	SetOPTxs (G4int opt)

void	SetVerboseLevel (G4int verbose)

void	UseSICB (G4bool use)

virtual	~G4PhotonEvaporation ()

Protected Attributes
G4int	OPTxs

G4bool	useSICB

Private Member Functions
G4Fragment *	GenerateGamma (G4Fragment *nucleus)

void	InitialiseGRData ()

void	InitialiseLevelManager (G4int Z, G4int A)

G4bool	operator!= (const G4VEvaporationChannel &right) const

G4bool	operator== (const G4VEvaporationChannel &right) const

Private Attributes
G4int	fCode

G4bool	fCorrelatedGamma

G4double	fCummProbability [MAXDEPOINT]

G4double	fExcEnergy

G4bool	fICM

size_t	fIndex

G4bool	fIsomerFlag

G4double	fLevelEnergyMax

const G4LevelManager *	fLevelManager

G4double	fMaxLifeTime

G4NuclearLevelData *	fNuclearLevelData

G4int	fPoints

G4NuclearPolarization *	fPolarization

G4double	fProbability

G4bool	fRDM

G4bool	fSampleTime

G4int	fSecID

G4double	fStep

G4GammaTransition *	fTransition

G4int	fVerbose

G4bool	isInitialised

G4int	theA

G4int	theZ

G4double	Tolerance

G4int	vShellNumber

Static Private Attributes
static G4float	GREnergy [MAXGRDATA] = {0.0f}

static G4float	GRWidth [MAXGRDATA] = {0.0f}

Detailed Description

Definition at line 63 of file G4PhotonEvaporation.hh.

Constructor & Destructor Documentation

◆ G4PhotonEvaporation() [1/2]

G4PhotonEvaporation::G4PhotonEvaporation ( G4GammaTransition * ptr = nullptr )

explicit

Definition at line 64 of file G4PhotonEvaporation.cc.

  : fLevelManager(nullptr), fTransition(p), fPolarization(nullptr),
    fVerbose(1), fPoints(0), vShellNumber(-1), fIndex(0), fSecID(-1),
    fMaxLifeTime(DBL_MAX), 
    fICM(true), fRDM(false), fSampleTime(true), 
    fCorrelatedGamma(false), fIsomerFlag(false), isInitialised(false)
{
  //G4cout << "### New G4PhotonEvaporation() " << this << G4endl;   
  fNuclearLevelData = G4NuclearLevelData::GetInstance(); 
  Tolerance = 20*CLHEP::eV;
 
  if(!fTransition) { fTransition = new G4GammaTransition(); }
 
  theA = theZ = fCode = 0;
  fSecID = G4PhysicsModelCatalog::GetModelID("model_G4PhotonEvaporation");
  fLevelEnergyMax = fStep = fExcEnergy = fProbability = 0.0;
 
  for(G4int i=0; i<MAXDEPOINT; ++i) { fCummProbability[i] = 0.0; }
  if(0.0f == GREnergy[1]) { InitialiseGRData(); }
}

References CLHEP::eV, fCode, fCummProbability, fExcEnergy, fLevelEnergyMax, fNuclearLevelData, fProbability, fSecID, fStep, fTransition, G4NuclearLevelData::GetInstance(), G4PhysicsModelCatalog::GetModelID(), GREnergy, InitialiseGRData(), MAXDEPOINT, theA, theZ, and Tolerance.

◆ ~G4PhotonEvaporation()

G4PhotonEvaporation::~G4PhotonEvaporation ( )

virtual

Definition at line 85 of file G4PhotonEvaporation.cc.

{ 
  delete fTransition;
}

References fTransition.

◆ G4PhotonEvaporation() [2/2]

G4PhotonEvaporation::G4PhotonEvaporation ( const G4PhotonEvaporation & right )

delete

Member Function Documentation

◆ BreakItUp()

G4FragmentVector * G4PhotonEvaporation::BreakItUp ( const G4Fragment & theNucleus )

Definition at line 182 of file G4PhotonEvaporation.cc.

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

References BreakUpChain(), fSecID, fVerbose, G4cout, G4endl, and G4Fragment::SetCreatorModelID().

Referenced by G4LENDCapture::ApplyYourself(), and G4ParticleHPCaptureFS::ApplyYourself().

◆ BreakUpChain()

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

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 195 of file G4PhotonEvaporation.cc.

{
  if(!isInitialised) { Initialise(); }
  if(fVerbose > 1) {
    G4cout << "G4PhotonEvaporation::BreakUpChain RDM= " << fRDM << " "
           << *nucleus << G4endl;
  }
  G4Fragment* gamma = nullptr;
  fSampleTime = !fRDM;
 
  // start decay chain from unpolarized state
  if(fCorrelatedGamma) {
    fPolarization = new G4NuclearPolarization(nucleus->GetZ_asInt(),
                                              nucleus->GetA_asInt(),
                                              nucleus->GetExcitationEnergy());
    nucleus->SetNuclearPolarization(fPolarization);
  }
 
  do {
    gamma = GenerateGamma(nucleus);
    if(gamma) {
      gamma->SetCreatorModelID(fSecID);
      products->push_back(gamma);
      if(fVerbose > 2) {
        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);
 
  // clear nuclear polarization end of chain
  if(fPolarization) {
    delete fPolarization;
    fPolarization = nullptr;
    nucleus->SetNuclearPolarization(fPolarization);
  }  
  return false;
}

References fCorrelatedGamma, fPolarization, fRDM, fSampleTime, fSecID, fVerbose, G4cout, G4endl, GenerateGamma(), G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetZ_asInt(), Initialise(), isInitialised, G4Fragment::SetCreatorModelID(), and G4Fragment::SetNuclearPolarization().

Referenced by BreakItUp().

◆ BreakUpFragment()

G4FragmentVector * G4VEvaporationChannel::BreakUpFragment ( G4Fragment * theNucleus )

inlineinherited

Definition at line 106 of file G4VEvaporationChannel.hh.

{
  G4FragmentVector* results = new G4FragmentVector();
  BreakUpChain(results, theNucleus);
  return results;
}

References G4VEvaporationChannel::BreakUpChain().

Referenced by G4NeutronRadCapture::ApplyYourself().

◆ Dump()

void G4VEvaporationChannel::Dump ( ) const

virtualinherited

Reimplemented in G4GEMChannel, and G4GEMChannelVI.

Definition at line 71 of file G4VEvaporationChannel.cc.

72{}

◆ EmittedFragment()

G4Fragment * G4PhotonEvaporation::EmittedFragment ( G4Fragment * theNucleus )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 131 of file G4PhotonEvaporation.cc.

{
  if(!isInitialised) { Initialise(); }
  fSampleTime = (fRDM) ? false : true;
 
  // potentially external code may set initial polarization
  // but only for radioactive decay nuclear polarization is considered
  G4NuclearPolarizationStore* fNucPStore = nullptr;
  if(fCorrelatedGamma && fRDM) {
    fNucPStore = G4NuclearPolarizationStore::GetInstance();
    if(nucleus->GetNuclearPolarization()) { 
      fNucPStore->RemoveMe(nucleus->GetNuclearPolarization());
      delete nucleus->GetNuclearPolarization(); 
    } 
    fPolarization = fNucPStore->FindOrBuild(nucleus->GetZ_asInt(),
                                            nucleus->GetA_asInt(),
                                            nucleus->GetExcitationEnergy());
    nucleus->SetNuclearPolarization(fPolarization);
  }
  if(fVerbose > 2) { 
    G4cout << "G4PhotonEvaporation::EmittedFragment: " 
           << *nucleus << G4endl;
    if(fPolarization) { G4cout << "NucPolar: " << fPolarization << G4endl; }
    G4cout << " CorrGamma: " << fCorrelatedGamma << " RDM: " << fRDM
           << " fPolarization: " << fPolarization << G4endl;
  }
  G4Fragment* gamma = GenerateGamma(nucleus);
  
  if(gamma != nullptr) { gamma->SetCreatorModelID(fSecID); }
  
  // remove G4NuclearPolarizaton when reach ground state
  if(fNucPStore && fPolarization && 0 == fIndex) {
    if(fVerbose > 3) { 
      G4cout << "G4PhotonEvaporation::EmittedFragment: remove " 
             << fPolarization << G4endl;
    }
    fNucPStore->RemoveMe(fPolarization);
    fPolarization = nullptr;
    nucleus->SetNuclearPolarization(fPolarization);
  }
 
  if(fVerbose > 2) {
    G4cout << "G4PhotonEvaporation::EmittedFragment: RDM= " 
           << fRDM << " done:" << G4endl; 
    if(gamma) { G4cout << *gamma << G4endl; }
    G4cout << "   Residual: " << *nucleus << G4endl;
  }
  return gamma; 
}

References fCorrelatedGamma, fIndex, G4NuclearPolarizationStore::FindOrBuild(), fPolarization, fRDM, fSampleTime, fSecID, fVerbose, G4cout, G4endl, GenerateGamma(), G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4NuclearPolarizationStore::GetInstance(), G4Fragment::GetNuclearPolarization(), G4Fragment::GetZ_asInt(), Initialise(), isInitialised, G4NuclearPolarizationStore::RemoveMe(), G4Fragment::SetCreatorModelID(), and G4Fragment::SetNuclearPolarization().

Referenced by G4ITDecay::DecayIt().

◆ GenerateGamma()

G4Fragment * G4PhotonEvaporation::GenerateGamma ( G4Fragment * nucleus )

private

Definition at line 343 of file G4PhotonEvaporation.cc.

{
  if(!isInitialised) { Initialise(); }
  G4Fragment* result = nullptr;
  G4double eexc = nucleus->GetExcitationEnergy();
  if(eexc <= Tolerance) { return result; }
 
  InitialiseLevelManager(nucleus->GetZ_asInt(), nucleus->GetA_asInt());
 
  G4double time = nucleus->GetCreationTime();
 
  G4double efinal = 0.0;
  G4double ratio  = 0.0;
  vShellNumber    = -1;
  G4int  JP1      = 0;
  G4int  JP2      = 0;
  G4int  multiP   = 0;
  G4bool isGamma  = true;
  G4bool isDiscrete = false;
 
  const G4NucLevel* level = nullptr;
  size_t ntrans = 0;
 
  if(fVerbose > 2) {
    G4cout << "GenerateGamma: " << " Eex= " << eexc
           << " Eexmax= " << fLevelEnergyMax << G4endl;
  }
  // initial discrete state
  if(fLevelManager && eexc <= fLevelEnergyMax + Tolerance) {
    fIndex = fLevelManager->NearestLevelIndex(eexc, fIndex);
    // initial state below 1st level
    if(0 == fIndex && eexc >= Tolerance 
       && fLevelManager->NumberOfTransitions() > 0) { fIndex = 1; }
    isDiscrete = true; 
    if(fVerbose > 2) {
      G4cout << "              index= " << fIndex 
             << " lTime= " << fLevelManager->LifeTime(fIndex) << G4endl;
    }
    if(0 < fIndex) {
      // for discrete transition  
      level = fLevelManager->GetLevel(fIndex);
      if(level) { 
        ntrans = level->NumberOfTransitions();
        JP1 = fLevelManager->SpinTwo(fIndex); 
        if(fVerbose > 2) {
          G4cout << "              ntrans= " << ntrans << " JP= " << JP1
                 << " RDM: " << fRDM << G4endl;
        }
        if(0 == ntrans && fLevelManager->FloatingLevel(fIndex) > 0) {
          --fIndex;
          level = fLevelManager->GetLevel(fIndex);
          ntrans = level->NumberOfTransitions();
          JP1 = fLevelManager->SpinTwo(fIndex); 
        }
      }
    }
  }
  // if a level has no defined transitions
  if(0 == ntrans) { isDiscrete = false; }
  if(fVerbose > 2) {
    G4int prec = G4cout.precision(4);
    G4cout << "GenerateGamma: Z= " << nucleus->GetZ_asInt()
           << " A= " << nucleus->GetA_asInt() 
           << " Exc= " << eexc << " Emax= " 
           << fLevelEnergyMax << " idx= " << fIndex
           << " fCode= " << fCode << " fPoints= " << fPoints
           << " Ntr= " << ntrans << " discrete: " << isDiscrete
           << " fProb= " << fProbability << G4endl;
    G4cout.precision(prec);
  }
 
  // continues part
  if(!isDiscrete) {
    // we compare current excitation versus value used for probability 
    // computation and also Z and A used for probability computation 
    if(fCode != 1000*theZ + theA || eexc != fExcEnergy) { 
      GetEmissionProbability(nucleus); 
    }
    if(fProbability == 0.0) { 
      fPoints = 1; 
      efinal = 0.0; 
    } else {
      G4double y = fCummProbability[fPoints-1]*G4UniformRand();
      for(G4int i=1; i<fPoints; ++i) {
        if(fVerbose > 3) {
          G4cout << "y= " << y << " cummProb= " << fCummProbability[i] 
                 << " fPoints= " << fPoints << " fStep= " << fStep << G4endl;
        }
        if(y <= fCummProbability[i]) {
          efinal = fStep*((i - 1) + (y - fCummProbability[i-1])
                          /(fCummProbability[i] - fCummProbability[i-1]));
          break;
        }
      }
    }
    // final discrete level
    if(fVerbose > 2) {
      G4cout << "Continues proposes Efinal= " << efinal << G4endl;
    }
    if(fLevelManager) {
      if(efinal < fLevelEnergyMax) {
        fIndex = fLevelManager->NearestLevelIndex(efinal, fIndex);
        efinal = fLevelManager->LevelEnergy(fIndex);
        // protection - take level below
        if(efinal >= eexc && 0 < fIndex) {
          --fIndex;
          efinal = fLevelManager->LevelEnergy(fIndex);
        } 
        nucleus->SetFloatingLevelNumber(fLevelManager->FloatingLevel(fIndex));
 
        // not allowed to have final energy above max energy
        // if G4LevelManager exist
      } else {
        efinal = fLevelEnergyMax;
        fIndex = fLevelManager->NearestLevelIndex(efinal, fIndex);
      }
    }
    if(fVerbose > 2) {
      G4cout << "Continues emission efinal(MeV)= " << efinal << G4endl; 
    }
    //discrete part ground state
  } else if(0 == fIndex) {
    return result;
 
    //discrete part
  } else {
 
    if(fVerbose > 2) {
      G4cout << "Discrete emission from level Index= " << fIndex 
             << " Elevel= " << fLevelManager->LevelEnergy(fIndex)
             << " Ltime= " << fLevelManager->LifeTime(fIndex)
             << " LtimeMax= " << fMaxLifeTime
             << "  RDM= " << fRDM << "  ICM= " << fICM << G4endl;
    }
 
    // stable fragment has life time -1 or above the limit
    // if is called from the radioactive decay the life time is not checked
    G4double ltime = fLevelManager->LifeTime(fIndex);
    if(ltime < 0.0 || (!fRDM && ltime > fMaxLifeTime)) { return result; }
 
    size_t idx = 0;
    if(1 < ntrans) {
      idx = level->SampleGammaTransition(G4UniformRand());
    }
    if(fVerbose > 2) {
      G4cout << "Ntrans= " << ntrans << " idx= " << idx
         << " ICM= " << fICM << "  JP1= " << JP1 << G4endl;
    }
    G4double prob = (G4double)level->GammaProbability(idx);
    // prob = 0 means that there is only internal conversion
    if(fICM && prob < 1.0) {
      G4double rndm = G4UniformRand();
      if(rndm > prob) {
    isGamma = false;
    rndm = (rndm - prob)/(1.0 - prob);
    vShellNumber = level->SampleShell(idx, rndm);
      }
    }
    // it is discrete transition with possible gamma correlation
    ratio  = level->MultipolarityRatio(idx);
    multiP = level->TransitionType(idx);
    fIndex = level->FinalExcitationIndex(idx);
    JP2 = fLevelManager->SpinTwo(fIndex); 
 
    // final energy and time
    efinal = fLevelManager->LevelEnergy(fIndex);
    // time is sampled if decay not prompt and this class called not 
    // from radioactive decay and isomer production is enabled 
    if(fSampleTime && fIsomerFlag && ltime > 0.0) { 
      time -= ltime*G4Log(G4UniformRand()); 
    }
    nucleus->SetFloatingLevelNumber(fLevelManager->FloatingLevel(fIndex));
  }
  // protection for floating levels
  if(std::abs(efinal - eexc) <= Tolerance) { return result; }
 
  result = fTransition->SampleTransition(nucleus, efinal, ratio, JP1,
                                         JP2, multiP, vShellNumber, 
                                         isDiscrete, isGamma);
  if(result) { result->SetCreationTime(time); }
 
  // updated residual nucleus
  nucleus->SetCreationTime(time);
  nucleus->SetSpin(0.5*JP2);
  if(fPolarization) { fPolarization->SetExcitationEnergy(efinal); }
 
  // ignore the floating levels with zero energy and create ground state
  if(efinal == 0.0 && fIndex > 0) {
    fIndex = 0;
    nucleus->SetFloatingLevelNumber(fLevelManager->FloatingLevel(fIndex));
  }
      
  if(fVerbose > 2) { 
    G4cout << "Final level E= " << efinal << " time= " << time 
           << " idxFinal= " << fIndex << " isDiscrete: " << isDiscrete
           << " isGamma: " << isGamma << " multiP= " << multiP 
           << " shell= " << vShellNumber 
           << " JP1= " << JP1 << " JP2= " << JP2 << G4endl;
  }
  return result;
}

References fCode, fCummProbability, fExcEnergy, fICM, G4NucLevel::FinalExcitationIndex(), fIndex, fIsomerFlag, fLevelEnergyMax, fLevelManager, G4LevelManager::FloatingLevel(), fMaxLifeTime, fPoints, fPolarization, fProbability, fRDM, fSampleTime, fStep, fTransition, fVerbose, G4cout, G4endl, G4Log(), G4UniformRand, G4NucLevel::GammaProbability(), G4Fragment::GetA_asInt(), G4Fragment::GetCreationTime(), GetEmissionProbability(), G4Fragment::GetExcitationEnergy(), G4LevelManager::GetLevel(), G4Fragment::GetZ_asInt(), Initialise(), InitialiseLevelManager(), isInitialised, G4LevelManager::LevelEnergy(), G4LevelManager::LifeTime(), G4NucLevel::MultipolarityRatio(), G4LevelManager::NearestLevelIndex(), G4LevelManager::NumberOfTransitions(), G4NucLevel::NumberOfTransitions(), CLHEP::prec, G4NucLevel::SampleGammaTransition(), G4NucLevel::SampleShell(), G4GammaTransition::SampleTransition(), G4Fragment::SetCreationTime(), G4NuclearPolarization::SetExcitationEnergy(), G4Fragment::SetFloatingLevelNumber(), G4Fragment::SetSpin(), G4LevelManager::SpinTwo(), theA, theZ, Tolerance, G4NucLevel::TransitionType(), and vShellNumber.

Referenced by BreakUpChain(), and EmittedFragment().

◆ GetEmissionProbability()

G4double G4PhotonEvaporation::GetEmissionProbability ( G4Fragment * theNucleus )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 240 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 > 2) {
    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;
  const G4double MaxDeltaEnergy = CLHEP::MeV;
  fPoints = std::min((G4int)(fStep/MaxDeltaEnergy) + 2, MAXDEPOINT);
  fStep /= ((G4double)(fPoints - 1));
  if(fVerbose > 2) {
    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 levelDensity = fNuclearLevelData->GetLevelDensity(Z,A,fExcEnergy);
  G4double xsqr = std::sqrt(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(levelDensity*std::abs(fExcEnergy - egam)) - xsqr))
      *gammaR2*gammaE2/(egdp2*egdp2 + gammaR2);
    fProbability += (p1 + p0);
    fCummProbability[i] = fProbability;
    if(fVerbose > 3) {
      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;
}

References A, G4Fragment::ComputeGroundStateMass(), emax, fCode, fCummProbability, fExcEnergy, fNuclearLevelData, fPoints, fProbability, fStep, fVerbose, G4cout, G4endl, G4Exp(), G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4NuclearLevelData::GetLevelDensity(), G4Fragment::GetZ_asInt(), GREnergy, GRWidth, CLHEP::hbarc, Initialise(), isInitialised, G4INCL::Math::max(), MAXDEPOINT, MAXGRDATA, CLHEP::MeV, CLHEP::millibarn, G4INCL::Math::min(), CLHEP::neutron_mass_c2, CLHEP::pi2, Tolerance, and Z.

Referenced by GenerateGamma().

◆ GetFinalLevelEnergy()

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

finalvirtual

Definition at line 325 of file G4PhotonEvaporation.cc.

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

References A, G4INCL::KinematicsUtils::energy(), fIndex, fLevelEnergyMax, fLevelManager, InitialiseLevelManager(), G4LevelManager::NearestLevelEnergy(), Tolerance, and Z.

◆ GetLifeTime()

G4double G4VEvaporationChannel::GetLifeTime ( G4Fragment * theNucleus )

virtualinherited

Definition at line 50 of file G4VEvaporationChannel.cc.

{
  return 0.0;
}

◆ GetUpperLevelEnergy()

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

finalvirtual

Definition at line 336 of file G4PhotonEvaporation.cc.

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

References A, fLevelEnergyMax, InitialiseLevelManager(), and Z.

◆ GetVacantShellNumber()

G4int G4PhotonEvaporation::GetVacantShellNumber ( ) const

inline

Definition at line 173 of file G4PhotonEvaporation.hh.

{ 
  return vShellNumber;
}

References vShellNumber.

Referenced by G4ITDecay::DecayIt().

◆ Initialise()

void G4PhotonEvaporation::Initialise ( )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 90 of file G4PhotonEvaporation.cc.

{
  if(isInitialised) { return; }
  isInitialised = true;
 
  G4DeexPrecoParameters* param = fNuclearLevelData->GetParameters();
  Tolerance = param->GetMinExcitation();
  fMaxLifeTime = param->GetMaxLifeTime();
  fCorrelatedGamma = param->CorrelatedGamma();
  fICM = param->GetInternalConversionFlag();
  fIsomerFlag = param->IsomerProduction();
  if(fRDM) { fIsomerFlag = true; }
  fVerbose = param->GetVerbose();
 
  fTransition->SetPolarizationFlag(fCorrelatedGamma);
  fTransition->SetTwoJMAX(param->GetTwoJMAX());
  fTransition->SetVerbose(fVerbose);
  if(fVerbose > 1) {
    G4cout << "### G4PhotonEvaporation is initialized " << this << G4endl;   
  }
}

References G4DeexPrecoParameters::CorrelatedGamma(), fCorrelatedGamma, fICM, fIsomerFlag, fMaxLifeTime, fNuclearLevelData, fRDM, fTransition, fVerbose, G4cout, G4endl, G4DeexPrecoParameters::GetInternalConversionFlag(), G4DeexPrecoParameters::GetMaxLifeTime(), G4DeexPrecoParameters::GetMinExcitation(), G4NuclearLevelData::GetParameters(), G4DeexPrecoParameters::GetTwoJMAX(), G4DeexPrecoParameters::GetVerbose(), isInitialised, G4DeexPrecoParameters::IsomerProduction(), G4GammaTransition::SetPolarizationFlag(), G4GammaTransition::SetTwoJMAX(), G4GammaTransition::SetVerbose(), and Tolerance.

Referenced by BreakUpChain(), EmittedFragment(), GenerateGamma(), and GetEmissionProbability().

◆ InitialiseGRData()

void G4PhotonEvaporation::InitialiseGRData ( )

private

Definition at line 112 of file G4PhotonEvaporation.cc.

{
#ifdef G4MULTITHREADED
  G4MUTEXLOCK(&G4PhotonEvaporation::PhotonEvaporationMutex);
#endif
  if(0.0f == GREnergy[1]) { 
    G4Pow* g4calc = G4Pow::GetInstance();
    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];
    }
  } 
#ifdef G4MULTITHREADED
  G4MUTEXUNLOCK(&G4PhotonEvaporation::PhotonEvaporationMutex);
#endif
}

References A, G4MUTEXLOCK, G4MUTEXUNLOCK, G4Pow::GetInstance(), GREnergy, GRWidth, MAXGRDATA, CLHEP::MeV, and G4Pow::powZ().

Referenced by G4PhotonEvaporation().

◆ InitialiseLevelManager()

void G4PhotonEvaporation::InitialiseLevelManager	(	G4int	Z,
		G4int	A
	)

inlineprivate

Definition at line 162 of file G4PhotonEvaporation.hh.

{
  if(Z != theZ || A != theA) {
    theZ = Z;
    theA = A;
    fIndex = 0;
    fLevelManager = fNuclearLevelData->GetLevelManager(theZ, theA);
    fLevelEnergyMax = fLevelManager ? fLevelManager->MaxLevelEnergy() : 0.0;
  }
}

References A, fIndex, fLevelEnergyMax, fLevelManager, fNuclearLevelData, G4NuclearLevelData::GetLevelManager(), G4LevelManager::MaxLevelEnergy(), theA, theZ, and Z.

Referenced by GenerateGamma(), GetFinalLevelEnergy(), and GetUpperLevelEnergy().

◆ operator!=()

G4bool G4VEvaporationChannel::operator!= ( const G4VEvaporationChannel & right ) const

privateinherited

◆ operator=()

const G4PhotonEvaporation & G4PhotonEvaporation::operator= ( const G4PhotonEvaporation & right )

delete

◆ operator==()

G4bool G4VEvaporationChannel::operator== ( const G4VEvaporationChannel & right ) const

privateinherited

◆ RDMForced()

void G4PhotonEvaporation::RDMForced ( G4bool val )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 558 of file G4PhotonEvaporation.cc.

{
  fRDM = val;
}

References fRDM.

Referenced by G4RadioactiveDecay::G4RadioactiveDecay().

◆ SetGammaTransition()

void G4PhotonEvaporation::SetGammaTransition ( G4GammaTransition * p )

Definition at line 545 of file G4PhotonEvaporation.cc.

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

References fTransition.

◆ SetICM()

void G4PhotonEvaporation::SetICM ( G4bool val )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 553 of file G4PhotonEvaporation.cc.

{
  fICM = val;
}

References fICM.

Referenced by G4LENDCapture::ApplyYourself(), G4ParticleHPCaptureFS::ApplyYourself(), and G4RadioactiveDecay::G4RadioactiveDecay().

◆ SetOPTxs()

void G4VEvaporationChannel::SetOPTxs ( G4int opt )

inlineinherited

Definition at line 113 of file G4VEvaporationChannel.hh.

114{}

◆ SetVerboseLevel()

void G4PhotonEvaporation::SetVerboseLevel ( G4int verbose )

inline

Definition at line 156 of file G4PhotonEvaporation.hh.

{
  fVerbose = verbose;
}

References fVerbose.

◆ UseSICB()

void G4VEvaporationChannel::UseSICB ( G4bool use )

inlineinherited

Definition at line 116 of file G4VEvaporationChannel.hh.

117{}

Field Documentation

◆ fCode

G4int G4PhotonEvaporation::fCode

private

Definition at line 125 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), and GetEmissionProbability().

◆ fCorrelatedGamma

G4bool G4PhotonEvaporation::fCorrelatedGamma

private

Definition at line 147 of file G4PhotonEvaporation.hh.

Referenced by BreakUpChain(), EmittedFragment(), and Initialise().

◆ fCummProbability

G4double G4PhotonEvaporation::fCummProbability[MAXDEPOINT]

private

Definition at line 134 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), and GetEmissionProbability().

◆ fExcEnergy

G4double G4PhotonEvaporation::fExcEnergy

private

Definition at line 137 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), and GetEmissionProbability().

◆ fICM

G4bool G4PhotonEvaporation::fICM

private

Definition at line 144 of file G4PhotonEvaporation.hh.

Referenced by GenerateGamma(), Initialise(), and SetICM().

◆ fIndex

size_t G4PhotonEvaporation::fIndex

private

Definition at line 127 of file G4PhotonEvaporation.hh.

Referenced by EmittedFragment(), GenerateGamma(), GetFinalLevelEnergy(), and InitialiseLevelManager().

◆ fIsomerFlag

G4bool G4PhotonEvaporation::fIsomerFlag

private

Definition at line 148 of file G4PhotonEvaporation.hh.

Referenced by GenerateGamma(), and Initialise().

◆ fLevelEnergyMax

G4double G4PhotonEvaporation::fLevelEnergyMax

private

Definition at line 136 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), GetFinalLevelEnergy(), GetUpperLevelEnergy(), and InitialiseLevelManager().

◆ fLevelManager

const G4LevelManager* G4PhotonEvaporation::fLevelManager

private

Definition at line 114 of file G4PhotonEvaporation.hh.

Referenced by GenerateGamma(), GetFinalLevelEnergy(), and InitialiseLevelManager().

◆ fMaxLifeTime

G4double G4PhotonEvaporation::fMaxLifeTime

private

Definition at line 140 of file G4PhotonEvaporation.hh.

Referenced by GenerateGamma(), and Initialise().

◆ fNuclearLevelData

G4NuclearLevelData* G4PhotonEvaporation::fNuclearLevelData

private

Definition at line 113 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GetEmissionProbability(), Initialise(), and InitialiseLevelManager().

◆ fPoints

G4int G4PhotonEvaporation::fPoints

private

Definition at line 124 of file G4PhotonEvaporation.hh.

Referenced by GenerateGamma(), and GetEmissionProbability().

◆ fPolarization

G4NuclearPolarization* G4PhotonEvaporation::fPolarization

private

Definition at line 119 of file G4PhotonEvaporation.hh.

Referenced by BreakUpChain(), EmittedFragment(), and GenerateGamma().

◆ fProbability

G4double G4PhotonEvaporation::fProbability

private

Definition at line 138 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), and GetEmissionProbability().

◆ fRDM

G4bool G4PhotonEvaporation::fRDM

private

Definition at line 145 of file G4PhotonEvaporation.hh.

Referenced by BreakUpChain(), EmittedFragment(), GenerateGamma(), Initialise(), and RDMForced().

◆ fSampleTime

G4bool G4PhotonEvaporation::fSampleTime

private

Definition at line 146 of file G4PhotonEvaporation.hh.

Referenced by BreakUpChain(), EmittedFragment(), and GenerateGamma().

◆ fSecID

G4int G4PhotonEvaporation::fSecID

private

Definition at line 129 of file G4PhotonEvaporation.hh.

Referenced by BreakItUp(), BreakUpChain(), EmittedFragment(), and G4PhotonEvaporation().

◆ fStep

G4double G4PhotonEvaporation::fStep

private

Definition at line 139 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), and GetEmissionProbability().

◆ fTransition

G4GammaTransition* G4PhotonEvaporation::fTransition

private

Definition at line 115 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), Initialise(), SetGammaTransition(), and ~G4PhotonEvaporation().

◆ fVerbose

G4int G4PhotonEvaporation::fVerbose

private

Definition at line 121 of file G4PhotonEvaporation.hh.

Referenced by BreakItUp(), BreakUpChain(), EmittedFragment(), GenerateGamma(), GetEmissionProbability(), Initialise(), and SetVerboseLevel().

◆ GREnergy

G4float G4PhotonEvaporation::GREnergy = {0.0f}

staticprivate

Definition at line 131 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GetEmissionProbability(), and InitialiseGRData().

◆ GRWidth

G4float G4PhotonEvaporation::GRWidth = {0.0f}

staticprivate

Definition at line 132 of file G4PhotonEvaporation.hh.

Referenced by GetEmissionProbability(), and InitialiseGRData().

◆ isInitialised

G4bool G4PhotonEvaporation::isInitialised

private

Definition at line 149 of file G4PhotonEvaporation.hh.

Referenced by BreakUpChain(), EmittedFragment(), GenerateGamma(), GetEmissionProbability(), and Initialise().

◆ OPTxs

G4int G4VEvaporationChannel::OPTxs

protectedinherited

Definition at line 93 of file G4VEvaporationChannel.hh.

Referenced by G4EvaporationChannel::GetEmissionProbability().

◆ theA

G4int G4PhotonEvaporation::theA

private

Definition at line 123 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), and InitialiseLevelManager().

◆ theZ

G4int G4PhotonEvaporation::theZ

private

Definition at line 122 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), and InitialiseLevelManager().

◆ Tolerance

G4double G4PhotonEvaporation::Tolerance

private

Definition at line 142 of file G4PhotonEvaporation.hh.

Referenced by G4PhotonEvaporation(), GenerateGamma(), GetEmissionProbability(), GetFinalLevelEnergy(), and Initialise().

◆ useSICB

G4bool G4VEvaporationChannel::useSICB

protectedinherited

Definition at line 94 of file G4VEvaporationChannel.hh.

◆ vShellNumber

G4int G4PhotonEvaporation::vShellNumber

private

Definition at line 126 of file G4PhotonEvaporation.hh.

Referenced by GenerateGamma(), and GetVacantShellNumber().

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

source/processes/hadronic/models/de_excitation/photon_evaporation/include/G4PhotonEvaporation.hh
source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4PhotonEvaporation.cc

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4PhotonEvaporation() [1/2]

◆ ~G4PhotonEvaporation()

◆ G4PhotonEvaporation() [2/2]

Member Function Documentation

◆ BreakItUp()

◆ BreakUpChain()

◆ BreakUpFragment()

◆ Dump()

◆ EmittedFragment()

◆ GenerateGamma()

◆ GetEmissionProbability()

◆ GetFinalLevelEnergy()

◆ GetLifeTime()

◆ GetUpperLevelEnergy()

◆ GetVacantShellNumber()

◆ Initialise()

◆ InitialiseGRData()

◆ InitialiseLevelManager()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ RDMForced()

◆ SetGammaTransition()

◆ SetICM()

◆ SetOPTxs()

◆ SetVerboseLevel()

◆ UseSICB()

Field Documentation

◆ fCode

◆ fCorrelatedGamma

◆ fCummProbability

◆ fExcEnergy

◆ fICM

◆ fIndex

◆ fIsomerFlag

◆ fLevelEnergyMax

◆ fLevelManager

◆ fMaxLifeTime

◆ fNuclearLevelData

◆ fPoints

◆ fPolarization

◆ fProbability

◆ fRDM

◆ fSampleTime

◆ fSecID

◆ fStep

◆ fTransition

◆ fVerbose

◆ GREnergy

◆ GRWidth

◆ isInitialised

◆ OPTxs

◆ theA

◆ theZ

◆ Tolerance

◆ useSICB

◆ vShellNumber