G4EvaporationChannel Class Reference

#include <G4EvaporationChannel.hh>

Inheritance diagram for G4EvaporationChannel:

Collaboration diagram for G4EvaporationChannel:

Public Member Functions
	G4EvaporationChannel (G4int A, G4int Z, const G4String &aName, G4EvaporationProbability *, G4VCoulombBarrier *)
virtual	~G4EvaporationChannel ()
void	Initialise ()
virtual G4double	GetEmissionProbability (G4Fragment *fragment)
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()
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)
virtual G4double	GetFinalLevelEnergy (G4int Z, G4int A, G4double energy)
virtual G4double	GetUpperLevelEnergy (G4int Z, G4int A)
G4double	GetMaxLevelEnergy (G4int Z, G4int A)
G4double	GetNearestLevelEnergy (G4int Z, G4int A, G4double energy)
void	SetPhotonEvaporation (G4VEvaporationChannel *p)
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 47 of file G4EvaporationChannel.hh.

Constructor & Destructor Documentation

G4EvaporationChannel::G4EvaporationChannel ( G4int  A, G4int  Z, const G4String &  aName, G4EvaporationProbability *  aprob, G4VCoulombBarrier *  barrier  )

explicit

Definition at line 53 of file G4EvaporationChannel.cc.

                                                                      :
    G4VEvaporationChannel(aName),
    theA(anA),
    theZ(aZ),
    theProbability(aprob),
    theCoulombBarrier(barrier)
{ 
  ResA = ResZ = 0;
  Mass = CoulombBarrier = MinKinEnergy = MaxKinEnergy = EmissionProbability = 0.0; 
  EvapMass = G4NucleiProperties::GetNuclearMass(theA, theZ);
  pairingCorrection = G4PairingCorrection::GetInstance();
}

Here is the call graph for this function:

G4EvaporationChannel::~G4EvaporationChannel ( )

virtual

Definition at line 69 of file G4EvaporationChannel.cc.

{}

Member Function Documentation

G4Fragment * G4EvaporationChannel::EmittedFragment ( G4Fragment *  theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 130 of file G4EvaporationChannel.cc.

{
  G4Fragment* evFragment = nullptr;
  G4double ekin = 0.0;
  if(ResA <= 4 && 
    ((ResA == 4 && ResZ == 2) || (ResA == 3 && ResZ == 2) ||
     (ResA == 3 && ResZ == 1) || (ResA == 2 && ResZ == 1) ||
     (ResA == 1 && ResZ == 1) || (ResA == 1 && ResZ == 0) )) {
    G4double mres = G4NucleiProperties::GetNuclearMass(ResA, ResZ);
    ekin = 0.5*(Mass*Mass - mres*mres + EvapMass*EvapMass)/Mass - EvapMass;
  } else {
    ekin = theProbability->SampleKineticEnergy(MinKinEnergy, MaxKinEnergy,
                           CoulombBarrier);
  }
  G4LorentzVector lv0 = theNucleus->GetMomentum();
  G4LorentzVector lv(std::sqrt(ekin*(ekin + 2.0*EvapMass))*G4RandomDirection(), 
             ekin + EvapMass);
  lv.boost(lv0.boostVector());

  evFragment = new G4Fragment(theA, theZ, lv);
  lv0 -= lv;
  theNucleus->SetZandA_asInt(ResZ, ResA);
  theNucleus->SetMomentum(lv0);

  return evFragment; 
} 

Here is the call graph for this function:

G4double G4EvaporationChannel::GetEmissionProbability ( G4Fragment *  fragment )

virtual

Implements G4VEvaporationChannel.

Definition at line 78 of file G4EvaporationChannel.cc.

{
  G4int FragA = fragment->GetA_asInt();
  G4int FragZ = fragment->GetZ_asInt();
  ResA = FragA - theA;
  ResZ = FragZ - theZ;

  G4double FragmentMass = fragment->GetGroundStateMass();
  G4double ExEnergy = fragment->GetExcitationEnergy();
  Mass = FragmentMass + ExEnergy;
  //G4cout << "G4EvaporationChannel::Initialize Z= " << theZ << " A= " << theA 
  //     << " FragZ= " << FragZ << " FragA= " << FragA << G4endl;
  EmissionProbability = 0.0;

  // Only channels which are physically allowed are taken into account 
  if (ResA >= ResZ && ResZ > 0 && ResA >= theA) {
  
    //Effective excitation energy
    G4double ResMass = G4NucleiProperties::GetNuclearMass(ResA, ResZ);

    CoulombBarrier = (0 == theZ) ? 0.0 : 
      theCoulombBarrier->GetCoulombBarrier(ResA,ResZ,ExEnergy);

    G4double delta0 = 
      std::max(0.0,pairingCorrection->GetPairingCorrection(FragA,FragZ));
    G4double delta1 = 
      std::max(0.0,pairingCorrection->GetPairingCorrection(ResA,ResZ));
    ResMass += delta1;  
    /*
    G4cout << "ExEnergy= " << ExEnergy << " Ec= " << CoulombBarrier
       << " delta0= " << delta0 << " delta1= " << delta1
       << " Free= " << Mass - ResMass - EvapMass 
       << G4endl;
    */
    // for OPTxs >0 penetration under the barrier is taken into account
    G4double elim = (0 == OPTxs) ? CoulombBarrier : CoulombBarrier*0.7;
    if(ExEnergy >= delta0 && Mass >= ResMass + EvapMass + elim) {
      G4double xm2 = (Mass - EvapMass)*(Mass - EvapMass);
      G4double xm  = Mass - EvapMass - elim;
      MinKinEnergy = (0.0 == elim) ? 0.0 : std::max(0.5*(xm2 - xm*xm)/Mass, 0.0);
      MaxKinEnergy = std::max(0.5*(xm2 - ResMass*ResMass)/Mass, 0.0);
      //G4cout << "Emin= " << MinKinEnergy << " Emax= " << MaxKinEnergy 
      //     << "  xm= " << xm  << G4endl;
      EmissionProbability = theProbability->
    TotalProbability(*fragment, MinKinEnergy, MaxKinEnergy, CoulombBarrier);
    }
  }
  //G4cout << "G4EvaporationChannel:: probability= " 
  //    << EmissionProbability << G4endl;   
  return EmissionProbability;
}

Here is the call graph for this function:

void G4EvaporationChannel::Initialise ( )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 72 of file G4EvaporationChannel.cc.

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

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/hadronic/models/de_excitation/evaporation/include/G4EvaporationChannel.hh
• geant4.10.03.p01/source/processes/hadronic/models/de_excitation/evaporation/src/G4EvaporationChannel.cc