G4GEMChannelVI Class Reference

#include <G4GEMChannelVI.hh>

Inheritance diagram for G4GEMChannelVI:

Public Member Functions
	G4GEMChannelVI (G4int theA, G4int theZ)
	~G4GEMChannelVI () override
void	Initialise () override
G4double	GetEmissionProbability (G4Fragment *theNucleus) override
G4Fragment *	EmittedFragment (G4Fragment *theNucleus) override
void	Dump () const override
	G4GEMChannelVI (const G4GEMChannelVI &right)=delete
const G4GEMChannelVI &	operator= (const G4GEMChannelVI &right)=delete
G4bool	operator== (const G4GEMChannelVI &right) const =delete
G4bool	operator!= (const G4GEMChannelVI &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 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() [1/2]

G4GEMChannelVI::G4GEMChannelVI ( G4int theA, G4int theZ )

explicit

Definition at line 49 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); }
  fEvapMass = G4NucleiProperties::GetNuclearMass(A, Z);
  fEvapMass2 = fEvapMass*fEvapMass;
 
  cBarrier = new G4CoulombBarrier(A, Z);
  fProbability = new G4GEMProbabilityVI(A, Z, lManager);
 
  fCoeff = CLHEP::millibarn/((CLHEP::pi*CLHEP::hbarc)*(CLHEP::pi*CLHEP::hbarc)); 
 
  secID = G4PhysicsModelCatalog::GetModelID("model_G4GEMChannelVI");
  if (Z == 0 && A == 1) {
    indexC = 0;
    fCoeff *= 2.0;
  } else if (Z == 1 && A == 1) {
    indexC = 1;
    fCoeff *= 2.0;
  } else if (Z == 1 && A == 2) {
    indexC = 2;
    fCoeff *= 3.0;
  } else if (Z == 1 && A == 3) {
    indexC = 3;
    fCoeff *= 2.0;
  } else if (Z == 2 && A == 3) {
    indexC = 4;
    fCoeff *= 2.0;
  } else if (Z == 2 && A == 4) {
    indexC = 5;
  } else {
    indexC = 6;
  }
}

Referenced by G4GEMChannelVI(), operator!=(), operator=(), and operator==().

~G4GEMChannelVI()

G4GEMChannelVI::~G4GEMChannelVI ( )

override

Definition at line 87 of file G4GEMChannelVI.cc.

{
  delete cBarrier;
  delete fProbability;
}

G4GEMChannelVI() [2/2]

G4GEMChannelVI::G4GEMChannelVI ( const G4GEMChannelVI & right )

delete

Member Function Documentation

Dump()

void G4GEMChannelVI::Dump ( ) const

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 202 of file G4GEMChannelVI.cc.

{}

EmittedFragment()

G4Fragment * G4GEMChannelVI::EmittedFragment ( G4Fragment * theNucleus )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 159 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();
  G4double ekin;
  G4double exc = 0.0;
  G4double probMax = std::max(fEData[nProb - 1].prob, 0.0); 
  if (0.0 >= probMax) {
    ekin = std::max(0.5*(fMass*fMass - fResMass*fResMass + fEvapMass2)
            /fMass - fEvapMass, 0.0);
  } else if (1 == nProb) {
    ekin = fProbability->SampleEnergy(fEData[0].ekin1, fEData[0].ekin2,
                      bCoulomb, fExc, 0.0);
  } else {
    G4double p = G4UniformRand()*probMax;
    G4int i{1};
    for (; i<nProb; ++i) {
      if (p <= fEData[i].prob) { break; }
    }
    G4double e1 = fEData[i - 1].exc;
    G4double e2 = fEData[i].exc;
    G4double p1 = fEData[i - 1].prob;
    G4double p2 = fEData[i].prob;
    exc = e1 + (e2 - e1)*(p - p1)/(p2 - p1);
    ekin = fProbability->SampleEnergy(fEData[i].ekin1, fEData[i].ekin2,
                                      bCoulomb, fExc, exc);
  }
  G4double m1 = fEvapMass + exc;
  G4LorentzVector lv(std::sqrt(ekin*(ekin + 2.0*m1))
             *G4RandomDirection(), ekin + m1);
  lv.boost(lv0.boostVector());
  evFragment = new G4Fragment(A, Z, lv);
  lv0 -= lv;
  evFragment->SetCreatorModelID(secID);
  theNucleus->SetZandA_asInt(resZ, resA);
  theNucleus->SetMomentum(lv0);
  theNucleus->SetCreatorModelID(secID);
  
  return evFragment;  
} 

GetEmissionProbability()

G4double G4GEMChannelVI::GetEmissionProbability ( G4Fragment * theNucleus )

overridevirtual

Implements G4VEvaporationChannel.

Definition at line 99 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; 
  }
 
  fExc = fragment->GetExcitationEnergy();
  fMass = fragment->GetGroundStateMass() + fExc;
  fResMass = G4NucleiProperties::GetNuclearMass(resA, resZ);
 
  // limit for the case when both evaporation and residual 
  // fragments are in ground states
  if (fMass <= fEvapMass + fResMass) { return 0.0; } 
 
  if (Z > 0) {
    bCoulomb = cBarrier->GetCoulombBarrier(resA, resZ, 0.0);
  }
  G4double de = fMass - fEvapMass - fResMass - bCoulomb;
  nProb = (G4int)(de/minExc);
  if (nProb <= 1 || indexC < 6 || resA <= 4) {
    nProb = 1;
  } else {
    nProb = std::min(nProb, nProbMax);
  }
  if (2 < fVerbose) {
    G4cout << "## G4GEMChannelVI::GetEmissionProbability fragZ="
       << fragZ << " fragA=" << fragA << " Z=" << Z << " A=" << A
       << " Eex(MeV)=" << fExc << " nProb=" << nProb 
       << G4endl;
  }
  fProbability->SetDecayKinematics(resZ, resA, fResMass, fMass);
  G4double sump = 0.0; 
  for (G4int i=0; i<nProb; ++i) {
    G4double exc = std::min(minExc*i, de);
    G4double m1 = fEvapMass + exc;
    G4double e2 = 0.5*((fMass-fResMass)*(fMass+fResMass) + m1*m1)/fMass - m1;
    G4double m2 = fMass - m1 - 0.5*bCoulomb;
    if (m2 < fResMass) {
      nProb = i;
      break;
    }
    G4double e1 = std::max(0.5*((fMass-m2)*(fMass+m2) + m1*m1)/fMass - m1, 0.0);
    if (e1 >= e2) {
      nProb = i;
      break;
    }
    sump += fProbability->TotalProbability(*fragment, e1, e2, bCoulomb, fExc, exc);
    fEData[i].exc = exc;
    fEData[i].ekin1 = e1;
    fEData[i].ekin2 = e2;
    fEData[i].prob = sump;
  }
  return sump;
}

Initialise()

void G4GEMChannelVI::Initialise ( )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 93 of file G4GEMChannelVI.cc.

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

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

---

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

• G4GEMChannelVI.hh
• G4GEMChannelVI.cc