#include <G4EvaporationChannel.hh>

Inheritance diagram for G4EvaporationChannel:

Public Member Functions
virtual G4bool	BreakUpChain (G4FragmentVector theResult, G4Fragment theNucleus)

G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

virtual void	Dump () const

G4Fragment *	EmittedFragment (G4Fragment *theNucleus) override

	G4EvaporationChannel (G4int A, G4int Z, G4EvaporationProbability *)

G4double	GetEmissionProbability (G4Fragment *fragment) override

virtual G4double	GetLifeTime (G4Fragment *theNucleus)

void	Initialise () override

virtual void	RDMForced (G4bool)

virtual void	SetICM (G4bool)

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

	~G4EvaporationChannel () override

Protected Attributes
G4int	OPTxs

G4bool	useSICB

Private Member Functions
	G4EvaporationChannel (const G4EvaporationChannel &right)

G4bool	operator!= (const G4EvaporationChannel &right) const

const G4EvaporationChannel &	operator= (const G4EvaporationChannel &right)

G4bool	operator== (const G4EvaporationChannel &right) const

Private Attributes
G4double	evapMass

G4double	evapMass2

G4double	mass

G4int	resA

G4double	resMass

G4int	resZ

G4int	secID

G4int	theA

G4CoulombBarrier *	theCoulombBarrier

G4NuclearLevelData *	theLevelData

G4EvaporationProbability *	theProbability

G4int	theZ

Detailed Description

Definition at line 45 of file G4EvaporationChannel.hh.

Constructor & Destructor Documentation

◆ G4EvaporationChannel() [1/2]

G4EvaporationChannel::G4EvaporationChannel	(	G4int	A,
		G4int	Z,
		G4EvaporationProbability *	aprob
	)

explicit

Definition at line 54 of file G4EvaporationChannel.cc.

                                                       :
  G4VEvaporationChannel(),
  theA(anA),
  theZ(aZ),
  secID(-1),
  theProbability(aprob),
  theCoulombBarrier(new G4CoulombBarrier(anA, aZ))
{ 
  resA = resZ = 0;
  secID = G4PhysicsModelCatalog::GetModelID("model_G4EvaporationChannel");
  mass = resMass = 0.0; 
  evapMass = G4NucleiProperties::GetNuclearMass(theA, theZ);
  //G4cout << "G4EvaporationChannel: Z= " << theZ << " A= " << theA 
  //      << " M(GeV)= " << evapMass/GeV << G4endl;
  evapMass2 = evapMass*evapMass;
  theLevelData = G4NuclearLevelData::GetInstance();
}

References evapMass, evapMass2, G4NuclearLevelData::GetInstance(), G4PhysicsModelCatalog::GetModelID(), G4NucleiProperties::GetNuclearMass(), mass, resA, resMass, resZ, secID, theA, theLevelData, and theZ.

◆ ~G4EvaporationChannel()

G4EvaporationChannel::~G4EvaporationChannel ( )

override

Definition at line 73 of file G4EvaporationChannel.cc.

{
  delete theCoulombBarrier;
}

References theCoulombBarrier.

◆ G4EvaporationChannel() [2/2]

G4EvaporationChannel::G4EvaporationChannel ( const G4EvaporationChannel & right )

private

Member Function Documentation

◆ BreakUpChain()

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

virtualinherited

Reimplemented in G4UnstableFragmentBreakUp, and G4PhotonEvaporation.

Definition at line 66 of file G4VEvaporationChannel.cc.

{
  return false;
}

Referenced by G4VEvaporationChannel::BreakUpFragment().

◆ 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 * G4EvaporationChannel::EmittedFragment ( G4Fragment * theNucleus )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 156 of file G4EvaporationChannel.cc.

{
  G4double ekin;
  // assumed, that TotalProbability(...) was already called
  // if value iz zero no possiblity to sample final state
  if(resA <= 4 || theProbability->GetProbability() == 0.0) {
    ekin = 0.5*(mass*mass - resMass*resMass + evapMass2)/mass - evapMass;
  } else {
    ekin = theProbability->SampleEnergy();
  }
  ekin = std::max(ekin, 0.0);
  G4LorentzVector lv0 = theNucleus->GetMomentum();
  G4LorentzVector lv(std::sqrt(ekin*(ekin + 2.0*evapMass))*G4RandomDirection(), 
                     ekin + evapMass);
  lv.boost(lv0.boostVector());
 
  G4Fragment* evFragment = new G4Fragment(theA, theZ, lv);
  if(evFragment != nullptr) { evFragment->SetCreatorModelID(secID); }
  lv0 -= lv;
  theNucleus->SetZandA_asInt(resZ, resA);
  theNucleus->SetMomentum(lv0);
  theNucleus->SetCreatorModelID(secID);
  
  //G4cout << "Residual: Z= " << resZ << " A= " << resA << " Eex= " 
  //     << theNucleus->GetExcitationEnergy() << G4endl;
  return evFragment; 
} 

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), evapMass, evapMass2, G4RandomDirection(), G4Fragment::GetMomentum(), mass, G4INCL::Math::max(), resA, resMass, resZ, G4VEmissionProbability::SampleEnergy(), secID, G4Fragment::SetCreatorModelID(), G4Fragment::SetMomentum(), G4Fragment::SetZandA_asInt(), theA, theProbability, and theZ.

◆ GetEmissionProbability()

G4double G4EvaporationChannel::GetEmissionProbability ( G4Fragment * fragment )

overridevirtual

Implements G4VEvaporationChannel.

Definition at line 84 of file G4EvaporationChannel.cc.

{
  theProbability->ResetProbability();
  G4int fragA = fragment->GetA_asInt();
  G4int fragZ = fragment->GetZ_asInt();
  resA = fragA - theA;
  resZ = fragZ - theZ;
 
  // Only channels which are physically allowed are taken into account 
  if(resA < theA || resA < resZ || resZ < 0 || (resA == theA && resZ < theZ)
     || ((resA > 1) && (resA == resZ || resZ == 0)))
    { return 0.0; }
 
  G4double exEnergy = fragment->GetExcitationEnergy();
  G4double delta0 = theLevelData->GetPairingCorrection(fragZ,fragA);
  /*
  G4cout << "G4EvaporationChannel::Initialize Z= "<<theZ<<" A= "<<theA 
     << " FragZ= " << fragZ << " FragA= " << fragA 
     << " exEnergy= " << exEnergy << " d0= " << delta0 << G4endl;
  */
  if(exEnergy < delta0) { return 0.0; }
 
  G4double fragMass = fragment->GetGroundStateMass();
  mass = fragMass + exEnergy;
 
  resMass = G4NucleiProperties::GetNuclearMass(resA, resZ);
  G4double bCoulomb = 0.0;
  G4double elim = 0.0;
  if(theZ > 0) {
    bCoulomb = theCoulombBarrier->GetCoulombBarrier(resA,resZ,exEnergy);
 
    // for OPTxs >0 penetration under the barrier is taken into account
    const G4double dCB = 3.5*CLHEP::MeV;
    elim = (0 != OPTxs) ? 
      std::max(bCoulomb*0.5, bCoulomb - dCB*theZ) : bCoulomb;
  }
  /*
  G4cout << "exEnergy= " << exEnergy << " Ec= " << bCoulomb
         << " d0= " << delta0 
     << " Free= " << mass - resMass - evapMass 
     << G4endl;
  */
  if(mass <= resMass + evapMass + elim) { return 0.0; }
 
  G4double twoMass = mass + mass;
  G4double ekinmax = 
    ((mass-resMass)*(mass+resMass) + evapMass2)/twoMass - evapMass;
  G4double ekinmin = 0.0;
  if(elim > 0.0) {
    G4double resM = std::max(mass - evapMass - elim, resMass);
    ekinmin = 
      std::max(((mass-resM)*(mass+resM) + evapMass2)/twoMass - evapMass,0.0);
  }
  /*
  G4cout << "Emin= " <<ekinmin<<" Emax= "<<ekinmax
     << " mass= " << mass << " resM= " << resMass 
     << " evapM= " << evapMass << G4endl;
  */
  if(ekinmax <= ekinmin) { return 0.0; }
 
  theProbability->SetDecayKinematics(resZ, resA, resMass, mass);
  G4double prob = theProbability->TotalProbability(*fragment, ekinmin, 
                                                   ekinmax, bCoulomb,
                                                   exEnergy - delta0);
  /*  
  G4cout<<"G4EvaporationChannel: prob= "<< prob << " Z= " << theZ 
        << " A= " << theA << " E1= " << ekinmin << " E2= " << ekinmax 
        << G4endl;
  */
  return prob;
}

References ekinmin, evapMass, evapMass2, G4Fragment::GetA_asInt(), G4CoulombBarrier::GetCoulombBarrier(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4NucleiProperties::GetNuclearMass(), G4NuclearLevelData::GetPairingCorrection(), G4Fragment::GetZ_asInt(), mass, G4INCL::Math::max(), CLHEP::MeV, G4VEvaporationChannel::OPTxs, resA, G4VEmissionProbability::ResetProbability(), resMass, resZ, G4VEmissionProbability::SetDecayKinematics(), theA, theCoulombBarrier, theLevelData, theProbability, theZ, and G4EvaporationProbability::TotalProbability().

◆ GetLifeTime()

G4double G4VEvaporationChannel::GetLifeTime ( G4Fragment * theNucleus )

virtualinherited

Definition at line 50 of file G4VEvaporationChannel.cc.

{
  return 0.0;
}

◆ Initialise()

void G4EvaporationChannel::Initialise ( )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 78 of file G4EvaporationChannel.cc.

{
  theProbability->Initialise();
  G4VEvaporationChannel::Initialise();  
}

References G4VEmissionProbability::Initialise(), G4VEvaporationChannel::Initialise(), and theProbability.

◆ operator!=()

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

private

◆ operator=()

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

private

◆ operator==()

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

private

◆ RDMForced()

void G4VEvaporationChannel::RDMForced ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 58 of file G4VEvaporationChannel.cc.

59{}

◆ SetICM()

void G4VEvaporationChannel::SetICM ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 55 of file G4VEvaporationChannel.cc.

56{}

Referenced by G4NeutronRadCapture::InitialiseModel().

◆ SetOPTxs()

void G4VEvaporationChannel::SetOPTxs ( G4int opt )

inlineinherited

Definition at line 113 of file G4VEvaporationChannel.hh.

114{}

◆ UseSICB()

void G4VEvaporationChannel::UseSICB ( G4bool use )

inlineinherited

Definition at line 116 of file G4VEvaporationChannel.hh.

117{}

Field Documentation

◆ evapMass

G4double G4EvaporationChannel::evapMass

private

Definition at line 79 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ evapMass2

G4double G4EvaporationChannel::evapMass2

private

Definition at line 80 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ mass

G4double G4EvaporationChannel::mass

private

Definition at line 76 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ OPTxs

G4int G4VEvaporationChannel::OPTxs

protectedinherited

Definition at line 93 of file G4VEvaporationChannel.hh.

Referenced by GetEmissionProbability().

◆ resA

G4int G4EvaporationChannel::resA

private

Definition at line 71 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ resMass

G4double G4EvaporationChannel::resMass

private

Definition at line 77 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ resZ

G4int G4EvaporationChannel::resZ

private

Definition at line 72 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ secID

G4int G4EvaporationChannel::secID

private

Definition at line 74 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), and G4EvaporationChannel().

◆ theA

G4int G4EvaporationChannel::theA

private

Definition at line 69 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ theCoulombBarrier

G4CoulombBarrier* G4EvaporationChannel::theCoulombBarrier

private

Definition at line 86 of file G4EvaporationChannel.hh.

Referenced by GetEmissionProbability(), and ~G4EvaporationChannel().

◆ theLevelData

G4NuclearLevelData* G4EvaporationChannel::theLevelData

private

Definition at line 88 of file G4EvaporationChannel.hh.

Referenced by G4EvaporationChannel(), and GetEmissionProbability().

◆ theProbability

G4EvaporationProbability* G4EvaporationChannel::theProbability

private

Definition at line 83 of file G4EvaporationChannel.hh.

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

◆ theZ

G4int G4EvaporationChannel::theZ

private

Definition at line 70 of file G4EvaporationChannel.hh.

Referenced by EmittedFragment(), G4EvaporationChannel(), and GetEmissionProbability().

◆ useSICB

G4bool G4VEvaporationChannel::useSICB

protectedinherited

Definition at line 94 of file G4VEvaporationChannel.hh.

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

source/processes/hadronic/models/de_excitation/evaporation/include/G4EvaporationChannel.hh
source/processes/hadronic/models/de_excitation/evaporation/src/G4EvaporationChannel.cc

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4EvaporationChannel() [1/2]

◆ ~G4EvaporationChannel()

◆ G4EvaporationChannel() [2/2]

Member Function Documentation

◆ BreakUpChain()

◆ BreakUpFragment()

◆ Dump()

◆ EmittedFragment()

◆ GetEmissionProbability()

◆ GetLifeTime()

◆ Initialise()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ RDMForced()

◆ SetICM()

◆ SetOPTxs()

◆ UseSICB()

Field Documentation

◆ evapMass

◆ evapMass2

◆ mass

◆ OPTxs

◆ resA

◆ resMass

◆ resZ

◆ secID

◆ theA

◆ theCoulombBarrier

◆ theLevelData

◆ theProbability

◆ theZ

◆ useSICB