G4NeutronEvaporationChannel Class Reference

#include <G4NeutronEvaporationChannel.hh>

Inheritance diagram for G4NeutronEvaporationChannel:

Public Member Functions
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus)
G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual void	Dump () const
G4Fragment *	EmittedFragment (G4Fragment *theNucleus) override
	G4NeutronEvaporationChannel ()
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)
	~G4NeutronEvaporationChannel () override

Protected Attributes
G4int	OPTxs
G4bool	useSICB

Private Member Functions
	G4NeutronEvaporationChannel (const G4NeutronEvaporationChannel &right)
G4bool	operator!= (const G4NeutronEvaporationChannel &right) const
const G4NeutronEvaporationChannel &	operator= (const G4NeutronEvaporationChannel &right)
G4bool	operator== (const G4NeutronEvaporationChannel &right) const

Private Attributes
G4double	evapMass
G4double	evapMass2
G4double	mass
G4NeutronEvaporationProbability	pr
G4int	resA
G4double	resMass
G4int	resZ
G4int	secID
G4int	theA
G4CoulombBarrier *	theCoulombBarrier
G4NuclearLevelData *	theLevelData
G4EvaporationProbability *	theProbability
G4int	theZ

Detailed Description

Definition at line 38 of file G4NeutronEvaporationChannel.hh.

Constructor & Destructor Documentation

G4NeutronEvaporationChannel() [1/2]

G4NeutronEvaporationChannel::G4NeutronEvaporationChannel ( )

explicit

Definition at line 34 of file G4NeutronEvaporationChannel.cc.

  : G4EvaporationChannel(1,0,&pr) 
{}

~G4NeutronEvaporationChannel()

G4NeutronEvaporationChannel::~G4NeutronEvaporationChannel ( )

override

Definition at line 38 of file G4NeutronEvaporationChannel.cc.

{}

G4NeutronEvaporationChannel() [2/2]

G4NeutronEvaporationChannel::G4NeutronEvaporationChannel ( const G4NeutronEvaporationChannel &  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.

{}

EmittedFragment()

G4Fragment * G4EvaporationChannel::EmittedFragment ( G4Fragment *  theNucleus )

overridevirtualinherited

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(), G4EvaporationChannel::evapMass, G4EvaporationChannel::evapMass2, G4RandomDirection(), G4Fragment::GetMomentum(), G4EvaporationChannel::mass, G4INCL::Math::max(), G4EvaporationChannel::resA, G4EvaporationChannel::resMass, G4EvaporationChannel::resZ, G4VEmissionProbability::SampleEnergy(), G4EvaporationChannel::secID, G4Fragment::SetCreatorModelID(), G4Fragment::SetMomentum(), G4Fragment::SetZandA_asInt(), G4EvaporationChannel::theA, G4EvaporationChannel::theProbability, and G4EvaporationChannel::theZ.

GetEmissionProbability()

G4double G4EvaporationChannel::GetEmissionProbability ( G4Fragment *  fragment )

overridevirtualinherited

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, G4EvaporationChannel::evapMass, G4EvaporationChannel::evapMass2, G4Fragment::GetA_asInt(), G4CoulombBarrier::GetCoulombBarrier(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4NucleiProperties::GetNuclearMass(), G4NuclearLevelData::GetPairingCorrection(), G4Fragment::GetZ_asInt(), G4EvaporationChannel::mass, G4INCL::Math::max(), CLHEP::MeV, G4VEvaporationChannel::OPTxs, G4EvaporationChannel::resA, G4VEmissionProbability::ResetProbability(), G4EvaporationChannel::resMass, G4EvaporationChannel::resZ, G4VEmissionProbability::SetDecayKinematics(), G4EvaporationChannel::theA, G4EvaporationChannel::theCoulombBarrier, G4EvaporationChannel::theLevelData, G4EvaporationChannel::theProbability, G4EvaporationChannel::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 ( )

overridevirtualinherited

Reimplemented from G4VEvaporationChannel.

Definition at line 78 of file G4EvaporationChannel.cc.

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

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

operator!=()

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

private

operator=()

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

private

operator==()

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

private

RDMForced()

void G4VEvaporationChannel::RDMForced ( G4bool  )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 58 of file G4VEvaporationChannel.cc.

{}

SetICM()

void G4VEvaporationChannel::SetICM ( G4bool  )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 55 of file G4VEvaporationChannel.cc.

{}

Referenced by G4NeutronRadCapture::InitialiseModel().

SetOPTxs()

void G4VEvaporationChannel::SetOPTxs ( G4int  opt )

inlineinherited

Definition at line 113 of file G4VEvaporationChannel.hh.

{}

UseSICB()

void G4VEvaporationChannel::UseSICB ( G4bool  use )

inlineinherited

Definition at line 116 of file G4VEvaporationChannel.hh.

{}

Field Documentation

evapMass

G4double G4EvaporationChannel::evapMass

privateinherited

Definition at line 79 of file G4EvaporationChannel.hh.

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

evapMass2

G4double G4EvaporationChannel::evapMass2

privateinherited

Definition at line 80 of file G4EvaporationChannel.hh.

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

mass

G4double G4EvaporationChannel::mass

privateinherited

Definition at line 76 of file G4EvaporationChannel.hh.

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

OPTxs

G4int G4VEvaporationChannel::OPTxs

protectedinherited

Definition at line 93 of file G4VEvaporationChannel.hh.

Referenced by G4EvaporationChannel::GetEmissionProbability().

pr

G4NeutronEvaporationProbability G4NeutronEvaporationChannel::pr

private

Definition at line 55 of file G4NeutronEvaporationChannel.hh.

resA

G4int G4EvaporationChannel::resA

privateinherited

Definition at line 71 of file G4EvaporationChannel.hh.

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

resMass

G4double G4EvaporationChannel::resMass

privateinherited

Definition at line 77 of file G4EvaporationChannel.hh.

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

resZ

G4int G4EvaporationChannel::resZ

privateinherited

Definition at line 72 of file G4EvaporationChannel.hh.

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

secID

G4int G4EvaporationChannel::secID

privateinherited

Definition at line 74 of file G4EvaporationChannel.hh.

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

theA

G4int G4EvaporationChannel::theA

privateinherited

Definition at line 69 of file G4EvaporationChannel.hh.

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

theCoulombBarrier

G4CoulombBarrier* G4EvaporationChannel::theCoulombBarrier

privateinherited

Definition at line 86 of file G4EvaporationChannel.hh.

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

theLevelData

G4NuclearLevelData* G4EvaporationChannel::theLevelData

privateinherited

Definition at line 88 of file G4EvaporationChannel.hh.

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

theProbability

G4EvaporationProbability* G4EvaporationChannel::theProbability

privateinherited

Definition at line 83 of file G4EvaporationChannel.hh.

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

theZ

G4int G4EvaporationChannel::theZ

privateinherited

Definition at line 70 of file G4EvaporationChannel.hh.

Referenced by G4EvaporationChannel::EmittedFragment(), G4EvaporationChannel::G4EvaporationChannel(), and G4EvaporationChannel::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/G4NeutronEvaporationChannel.hh
• source/processes/hadronic/models/de_excitation/evaporation/src/G4NeutronEvaporationChannel.cc