G4GEMChannelVI Class Reference

#include <G4GEMChannelVI.hh>

Inheritance diagram for G4GEMChannelVI:

Collaboration diagram for G4GEMChannelVI:

Public Member Functions
	G4GEMChannelVI (G4int theA, G4int theZ)
virtual	~G4GEMChannelVI ()
virtual void	Initialise () final
virtual G4double	GetEmissionProbability (G4Fragment *theNucleus) final
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus) final
virtual void	Dump () const
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	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 49 of file G4GEMChannelVI.hh.

Constructor & Destructor Documentation

G4GEMChannelVI::G4GEMChannelVI ( G4int  theA, G4int  theZ  )

explicit

Definition at line 71 of file G4GEMChannelVI.cc.

  : A(theA), Z(theZ), levelDensity(0.1)
{ 
  fG4pow = G4Pow::GetInstance(); 
  Z13 = fG4pow->Z13(Z);
  A13 = fG4pow->Z13(A);

  cBarrier = new G4CoulombBarrier(A, Z);
  pairingCorrection = G4PairingCorrection::GetInstance();

  nData = G4NuclearLevelData::GetInstance();
  levelManager = nData->GetLevelManager(Z, A);
  maxLevelE = levelManager->MaxLevelEnergy();

  massGround = G4NucleiProperties::GetNuclearMass(A, Z);

  resA = resZ = fragZ = fragA = nWarn = 0;
  massGround = maxLevelE = Z13 = A13 = massFrag = eCBarrier 
    = resMassGround = resZ13 = resA13 = delta0
    = delta1 = maxExc = maxProb = alphaP = betaP = maxKinEnergy = 0.0;
  coeff = CLHEP::fermi*CLHEP::fermi/(CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc);
}

Here is the call graph for this function:

G4GEMChannelVI::~G4GEMChannelVI ( )

virtual

Definition at line 94 of file G4GEMChannelVI.cc.

{
  delete cBarrier;
}

Member Function Documentation

void G4GEMChannelVI::Dump ( ) const

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 368 of file G4GEMChannelVI.cc.

{
}

G4Fragment * G4GEMChannelVI::EmittedFragment ( G4Fragment *  theNucleus )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 180 of file G4GEMChannelVI.cc.

{
  G4double exc;
  G4double resExc;

  static G4double factor = 1.2;
  maxProb *= factor;
 
  CLHEP::HepRandomEngine* rndm = G4Random::getTheEngine();
  for(G4int i=0; i<100; ++i) {
    do {
      exc = maxExc*rndm->flat();
      resExc = maxExc*rndm->flat();
    } while (exc + resExc > maxExc);

    G4double prob = ProbabilityDistributionFunction(exc, resExc);
    if(prob > maxProb && nWarn < 10) {
      ++nWarn;
      G4cout << "### G4GEMChannelVI::EmittedFragment WARNING: majoranta "
         << maxProb << " is exceeded " << prob << "\n"
         << " fragZ= " << fragZ << " fragA= " << fragA
         << " Z= " << Z << " A= " << A 
         << " resZ= " << resZ << " resA= " << resA << "\n"
             << " exc(MeV)= " << exc << " resExc(MeV)= " << resExc 
         << " maxExc(MeV)= " << maxExc << G4endl;
    }
    if(maxProb*rndm->flat() <= prob) { break; }
  }
  if(exc <= maxLevelE) {
    exc = FindLevel(levelManager, exc, maxExc - resExc);
  }
  if(resA >= nData->GetMinA(resZ) && resA <= nData->GetMaxA(resZ) 
     && resExc < nData->GetMaxLevelEnergy(Z, A)) {
    const G4LevelManager* lman = nData->GetLevelManager(Z, A);
    if(lman) { resExc = FindLevel(lman, resExc, maxExc - exc); }
  }

  G4LorentzVector lv0 = theNucleus->GetMomentum();
  G4double mass1 = massGround + exc;
  G4double mass2 = resMassGround + resExc;

  G4double e1 = 0.5*((massFrag - mass2)*(massFrag + mass2) 
             + mass1*mass1)/massFrag;

  G4double p1(0.0);
  if(e1 > mass1) {
    p1 = std::sqrt((e1 - mass1)*(e1 + mass1));
  } else {
    e1 = mass1;
  }
  G4ThreeVector v = G4RandomDirection();
  G4LorentzVector lv1(p1*v, e1);

  G4ThreeVector boostVector = lv0.boostVector();
  lv1.boost(boostVector);
  
  G4Fragment* frag = new G4Fragment(A, Z, lv1);

  G4double e2 = massFrag - e1;
  if(e2 < mass2) {
    e2 = mass2;
    p1 = 0.0;
  }
  lv0.set(-v*p1, e2);
  lv0.boost(boostVector);

  theNucleus->SetZandA_asInt(resZ, resA);
  theNucleus->SetMomentum(lv0);

  return frag; 
} 

Here is the call graph for this function:

G4double G4GEMChannelVI::GetEmissionProbability ( G4Fragment *  theNucleus )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 106 of file G4GEMChannelVI.cc.

{
  fragZ = fragment->GetZ_asInt();
  fragA = fragment->GetA_asInt();
  resZ = fragZ - Z;
  resA = fragA - A;
  G4double prob = 0.0;
  if(resA < A || resA < resZ || resZ < 0 || (resA == A && resZ < Z)) { 
    return prob; 
  }
 
  resMassGround = G4NucleiProperties::GetNuclearMass(resA, resZ);
  G4double exc = fragment->GetExcitationEnergy();
  massFrag = fragment->GetGroundStateMass() + exc;
  delta0 = pairingCorrection->GetPairingCorrection(fragA, fragZ);
  eCBarrier = cBarrier->GetCoulombBarrier(resA, resZ, exc);
  
  maxExc = massFrag - massGround - resMassGround - eCBarrier - delta0;
  if(maxExc < 0.0) { return prob; }

  resZ13 = fG4pow->Z13(resZ);
  resA13 = fG4pow->Z13(resA);
  delta1 = pairingCorrection->GetPairingCorrection(resA, resZ);

  G4double C = 0.0;
  if(resA >= 50) {
    C = -0.10/G4double(A);
  } else if(resZ > 20) {
    C = (0.123482-0.00534691*Z-0.0000610624*(Z*Z)+5.93719*1e-7*(Z*Z*Z)+
     1.95687*1e-8*(Z*Z*Z*Z))/G4double(A);
  }
  if(0 == Z) {
    alphaP = 0.76+1.93/resA13;
    betaP = (1.66/(resA13*resA13)-0.05)*CLHEP::MeV/alphaP;
  } else {
    alphaP = 1.0 + C;
    betaP = - eCBarrier;
  }

  maxProb = 0.0;
  G4double e0 = maxExc*0.5;

  // e is an excitation of emitted fragment 
  for (G4int i=0; i<NPOINTSGEM; ++i) {
    prob += ws[i]*IntegratedProbability(e0*(xs[i] + 1.0));
  }
  prob *= coeff*e0*e0;

  /*
  G4cout << "G4GEMChannelVI: Z= " << Z << "  A= " << A 
     << " FragmentZ= " << aZ << " FragmentA= " << anA
     << " Zres= " << ResidualZ << " Ares= " << ResidualA 
     << G4endl; 
  */
  
  //G4cout << "Prob= " << prob << G4endl;
  return prob;
}

Here is the call graph for this function:

void G4GEMChannelVI::Initialise ( )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 99 of file G4GEMChannelVI.cc.

{
  G4DeexPrecoParameters* param = 
    G4NuclearLevelData::GetInstance()->GetParameters();
  levelDensity = param->GetLevelDensity();
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/gem_evaporation/include/G4GEMChannelVI.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/gem_evaporation/src/G4GEMChannelVI.cc