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 ()
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 40 of file G4GEMChannelVI.hh.

Constructor & Destructor Documentation

G4GEMChannelVI()

G4GEMChannelVI::G4GEMChannelVI ( G4int  theA, G4int  theZ  )

explicit

Definition at line 41 of file G4GEMChannelVI.cc.

  : A(theA), Z(theZ)
{ 
  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; 
}

~G4GEMChannelVI()

G4GEMChannelVI::~G4GEMChannelVI ( )

final

Definition at line 58 of file G4GEMChannelVI.cc.

{
  delete cBarrier;
  delete fProbability;
}

Member Function Documentation

Dump()

void G4GEMChannelVI::Dump ( ) const

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 122 of file G4GEMChannelVI.cc.

{}

EmittedFragment()

G4Fragment * G4GEMChannelVI::EmittedFragment ( G4Fragment *  theNucleus )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 94 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; 
  }
  theNucleus->SetZandA_asInt(resZ, resA);
  theNucleus->SetMomentum(lv0);
 
  return evFragment;  
} 

GetEmissionProbability()

G4double G4GEMChannelVI::GetEmissionProbability ( G4Fragment *  theNucleus )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 64 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