G4GEMChannelVI Class Reference

#include <G4GEMChannelVI.hh>

Inheritance diagram for G4GEMChannelVI:

Public Member Functions
	G4GEMChannelVI (G4int theA, G4int theZ)
	~G4GEMChannelVI () final
G4double	GetEmissionProbability (G4Fragment *theNucleus) final
G4Fragment *	EmittedFragment (G4Fragment *theNucleus) final
void	Dump () const final
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()=default
virtual void	Initialise ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus)
G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)
virtual G4double	ComputeInverseXSection (G4Fragment *theNucleus, G4double kinEnergy)
virtual G4double	ComputeProbability (G4Fragment *theNucleus, G4double kinEnergy)
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 40 of file G4GEMChannelVI.hh.

Constructor & Destructor Documentation

G4GEMChannelVI()

G4GEMChannelVI::G4GEMChannelVI ( G4int theA, G4int theZ )

explicit

Definition at line 42 of file G4GEMChannelVI.cc.

  : A(theA), Z(theZ), secID(-1)
{ 
  G4NuclearLevelData* nData = G4NuclearLevelData::GetInstance();
  pairingCorrection = nData->GetPairingCorrection();
  const G4LevelManager* lManager = nullptr;
  if(A > 4) { lManager = nData->GetLevelManager(Z, A); }
  evapMass = G4NucleiProperties::GetNuclearMass(A, Z);
  evapMass2 = evapMass*evapMass;
 
  cBarrier = new G4CoulombBarrier(A, Z);
  fProbability = new G4GEMProbabilityVI(A, Z, lManager);
 
  resA = resZ = fragZ = fragA = 0;
  mass = resMass = 0.0;
  secID = G4PhysicsModelCatalog::GetModelID("model_G4GEMChannelVI");
}

~G4GEMChannelVI()

G4GEMChannelVI::~G4GEMChannelVI ( )

final

Definition at line 60 of file G4GEMChannelVI.cc.

{
  delete cBarrier;
  delete fProbability;
}

Member Function Documentation

Dump()

void G4GEMChannelVI::Dump ( ) const

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 126 of file G4GEMChannelVI.cc.

{}

EmittedFragment()

G4Fragment * G4GEMChannelVI::EmittedFragment ( G4Fragment * theNucleus )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 96 of file G4GEMChannelVI.cc.

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

GetEmissionProbability()

G4double G4GEMChannelVI::GetEmissionProbability ( G4Fragment * theNucleus )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 66 of file G4GEMChannelVI.cc.

{
  fProbability->ResetProbability();
  fragZ = fragment->GetZ_asInt();
  fragA = fragment->GetA_asInt();
  resZ = fragZ - Z;
  resA = fragA - A;
  if(resA < A || resA < resZ || resZ < 0 || (resA == A && resZ < Z)) { 
    return 0.0; 
  }
 
  const G4double exc = fragment->GetExcitationEnergy();
  const G4double delta0 = 
    std::max(pairingCorrection->GetPairingCorrection(fragA, fragZ),0.0);
  if(exc < delta0) { return 0.0; }
 
  resMass = G4NucleiProperties::GetNuclearMass(resA, resZ);
  const G4double fragM = fragment->GetGroundStateMass() + exc;
  const G4double CB = cBarrier->GetCoulombBarrier(resA, resZ, exc);
  
  const G4double delta1 = 
    std::max(0.0,pairingCorrection->GetPairingCorrection(resA,resZ));
  if(fragM <= resMass + CB + delta1) { return 0.0; }
 
  fProbability->SetDecayKinematics(resZ, resA, resMass, fragM);
  G4double prob = fProbability->ComputeTotalProbability(*fragment, CB);
  //G4cout<<"G4EvaporationChannel: probability= "<< prob <<G4endl;
  return prob;
}

---

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

• G4GEMChannelVI.hh
• G4GEMChannelVI.cc