G4EvaporationChannel Class Reference

#include <G4EvaporationChannel.hh>

Inheritance diagram for G4EvaporationChannel:

Public Member Functions
	G4EvaporationChannel (G4int A, G4int Z, G4EvaporationProbability *)
	~G4EvaporationChannel () override
void	Initialise () override
G4double	GetEmissionProbability (G4Fragment *fragment) override
G4Fragment *	EmittedFragment (G4Fragment *theNucleus) override
G4double	ComputeInverseXSection (G4Fragment *, G4double kinEnergy) override
G4double	ComputeProbability (G4Fragment *, G4double kinEnergy) override
G4int	GetZ () const
G4int	GetA () const
G4EvaporationProbability *	GetEvaporationProbability ()
	G4EvaporationChannel (const G4EvaporationChannel &right)=delete
const G4EvaporationChannel &	operator= (const G4EvaporationChannel &right)=delete
G4bool	operator== (const G4EvaporationChannel &right) const =delete
G4bool	operator!= (const G4EvaporationChannel &right) const =delete
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()=default
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus)
G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual void	Dump () const
virtual void	SetICM (G4bool)
virtual void	RDMForced (G4bool)
void	SetOPTxs (G4int)
void	UseSICB (G4bool)
	G4VEvaporationChannel (const G4VEvaporationChannel &right)=delete
const G4VEvaporationChannel &	operator= (const G4VEvaporationChannel &right)=delete
G4bool	operator== (const G4VEvaporationChannel &right) const =delete
G4bool	operator!= (const G4VEvaporationChannel &right) const =delete

Additional Inherited Members
Protected Attributes inherited from G4VEvaporationChannel
G4int	OPTxs {3}
G4bool	useSICB {true}

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 53 of file G4EvaporationChannel.cc.

                                                       :
  G4VEvaporationChannel(),
  theProbability(aprob),
  theCoulombBarrier(new G4CoulombBarrier(anA, aZ)),
  theA(anA), theZ(aZ)
{ 
  secID = G4PhysicsModelCatalog::GetModelID("model_G4EvaporationChannel");
  evapMass = G4NucleiProperties::GetNuclearMass(theA, theZ);
  evapMass2 = evapMass*evapMass;
  theLevelData = G4NuclearLevelData::GetInstance();
}

Referenced by G4AlphaEvaporationChannel::G4AlphaEvaporationChannel(), G4DeuteronEvaporationChannel::G4DeuteronEvaporationChannel(), G4EvaporationChannel(), G4He3EvaporationChannel::G4He3EvaporationChannel(), G4NeutronEvaporationChannel::G4NeutronEvaporationChannel(), G4ProtonEvaporationChannel::G4ProtonEvaporationChannel(), G4TritonEvaporationChannel::G4TritonEvaporationChannel(), operator!=(), operator=(), and operator==().

~G4EvaporationChannel()

G4EvaporationChannel::~G4EvaporationChannel ( )

override

Definition at line 66 of file G4EvaporationChannel.cc.

{
  delete theCoulombBarrier;
}

G4EvaporationChannel() [2/2]

G4EvaporationChannel::G4EvaporationChannel ( const G4EvaporationChannel & right )

delete

Member Function Documentation

ComputeInverseXSection()

G4double G4EvaporationChannel::ComputeInverseXSection ( G4Fragment * frag, G4double kinEnergy )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 154 of file G4EvaporationChannel.cc.

{
  G4double p = ComputeProbability(frag, kinEnergy);
  return (p > 0.0) ? theProbability->RecentXS() : 0.0;
}

ComputeProbability()

G4double G4EvaporationChannel::ComputeProbability ( G4Fragment * frag, G4double kinEnergy )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 161 of file G4EvaporationChannel.cc.

{
  G4double prob = GetEmissionProbability(frag);
  if (prob <= 0.0) { return 0.0; }
 
  bCoulomb = (theZ > 0) ? theCoulombBarrier->GetCoulombBarrier(resA, resZ, 0.0) : 0.0;
  G4double p = theProbability->ComputeProbability(kinEnergy, bCoulomb);
  return p;
}

Referenced by ComputeInverseXSection().

EmittedFragment()

G4Fragment * G4EvaporationChannel::EmittedFragment ( G4Fragment * theNucleus )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 132 of file G4EvaporationChannel.cc.

{
  G4double ekin = ekinmax;
  // assumed, that TotalProbability(...) was already called
  // if value iz zero no possiblity to sample final state
  if(resA > 4 && theProbability->GetProbability() > 0.0) {
    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);
  evFragment->SetCreatorModelID(secID);
  lv0 -= lv;
  theNucleus->SetZAandMomentum(lv0, resZ, resA);
  theNucleus->SetCreatorModelID(secID);
  return evFragment; 
} 

GetA()

G4int G4EvaporationChannel::GetA ( ) const

inline

Definition at line 66 of file G4EvaporationChannel.hh.

{ return theA; };

GetEmissionProbability()

G4double G4EvaporationChannel::GetEmissionProbability ( G4Fragment * fragment )

overridevirtual

Implements G4VEvaporationChannel.

Definition at line 77 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 fragMass = fragment->GetGroundStateMass();
  mass = fragMass + exEnergy;
  resMass = G4NucleiProperties::GetNuclearMass(resA, resZ);
  if (mass <= evapMass + resMass) { return 0.0; } 
 
  ekinmax = 0.5*((mass-resMass)*(mass+resMass) + evapMass2)/mass - evapMass;
 
  // for OPTxs=1 elim=0 for all fragments - x-section include the CoulombBarrier
  G4double elim = 0.0;
  if(theZ > 0) {
    bCoulomb = theCoulombBarrier->GetCoulombBarrier(resA, resZ, 0.0);
 
    // for OPTxs >0 penetration under the barrier is taken into account
    elim = (0 < OPTxs) ? bCoulomb*0.5 : bCoulomb;
  }
  /*
  G4cout << "G4EvaporationChannel::Initialize Z=" << theZ <<" A=" << theA 
     << " FragZ=" << fragZ << " FragA=" << fragA << G4endl;
  G4cout << "      Eex=" << exEnergy << " CB=" << bCoulomb
         << " Elim=" << elim << " Efree=" << mass - resMass - evapMass 
     << G4endl;
  */
  // Coulomb barrier compound at rest
  G4double resM = mass - evapMass - elim;
  if (resM < resMass) { return 0.0; }
  G4double ekinmin = 0.5*((mass-resM)*(mass+resM) + evapMass2)/mass - evapMass;
 
  /*  
  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);
  return prob;
}

Referenced by ComputeProbability().

GetEvaporationProbability()

G4EvaporationProbability * G4EvaporationChannel::GetEvaporationProbability ( )

inline

Definition at line 68 of file G4EvaporationChannel.hh.

  { return theProbability; }

GetZ()

G4int G4EvaporationChannel::GetZ ( ) const

inline

Definition at line 64 of file G4EvaporationChannel.hh.

{ return theZ; };

Initialise()

void G4EvaporationChannel::Initialise ( )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 71 of file G4EvaporationChannel.cc.

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

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

---

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

• G4EvaporationChannel.hh
• G4EvaporationChannel.cc