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
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)=0
virtual void	Initialise ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4Fragment *	EmittedFragment (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 opt)
void	UseSICB (G4bool use)

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

Detailed Description

Definition at line 45 of file G4EvaporationChannel.hh.

Constructor & Destructor Documentation

G4EvaporationChannel()

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

explicit

Definition at line 53 of file G4EvaporationChannel.cc.

                                                       :
  G4VEvaporationChannel(),
  theA(anA),
  theZ(aZ),
  theProbability(aprob),
  theCoulombBarrier(new G4CoulombBarrier(anA, aZ))
{ 
  resA = resZ = 0;
  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();
}

~G4EvaporationChannel()

G4EvaporationChannel::~G4EvaporationChannel ( )

override

Definition at line 70 of file G4EvaporationChannel.cc.

{
  delete theCoulombBarrier;
}

Member Function Documentation

EmittedFragment()

G4Fragment * G4EvaporationChannel::EmittedFragment ( G4Fragment *  theNucleus )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 153 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);
  lv0 -= lv;
  theNucleus->SetZandA_asInt(resZ, resA);
  theNucleus->SetMomentum(lv0);
 
  //G4cout << "Residual: Z= " << resZ << " A= " << resA << " Eex= " 
  //     << theNucleus->GetExcitationEnergy() << G4endl;
  return evFragment; 
} 

GetEmissionProbability()

G4double G4EvaporationChannel::GetEmissionProbability ( G4Fragment *  fragment )

overridevirtual

Implements G4VEvaporationChannel.

Definition at line 81 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;
}

Initialise()

void G4EvaporationChannel::Initialise ( )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 75 of file G4EvaporationChannel.cc.

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

---

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

• G4EvaporationChannel.hh
• G4EvaporationChannel.cc