G4FermiBreakUpVI Class Reference

#include <G4FermiBreakUpVI.hh>

Inheritance diagram for G4FermiBreakUpVI:

Collaboration diagram for G4FermiBreakUpVI:

Public Member Functions
	G4FermiBreakUpVI ()
virtual	~G4FermiBreakUpVI ()
virtual void	Initialise () final
virtual G4bool	IsApplicable (G4int ZZ, G4int AA, G4double etot) const final
virtual void	BreakFragment (G4FragmentVector *, G4Fragment *theNucleus) final
Public Member Functions inherited from G4VFermiBreakUp
	G4VFermiBreakUp ()
virtual	~G4VFermiBreakUp ()

Detailed Description

Definition at line 47 of file G4FermiBreakUpVI.hh.

Constructor & Destructor Documentation

G4FermiBreakUpVI::G4FermiBreakUpVI ( )

explicit

Definition at line 47 of file G4FermiBreakUpVI.cc.

  : theDecay(nullptr), rndmEngine(nullptr), verbose(0), maxZ(9), maxA(17)
{
  prob.reserve(10);
  frag.reserve(10);
  lvect.reserve(10);
  Z = A = spin = 0;
  mass = elim = excitation = tolerance = 0.0;  
  frag1 = frag2 = nullptr;
  Initialise();
}

Here is the call graph for this function:

G4FermiBreakUpVI::~G4FermiBreakUpVI ( )

virtual

Definition at line 59 of file G4FermiBreakUpVI.cc.

{
  if(G4Threading::IsMasterThread()) { 
    delete thePool;
    thePool = nullptr;
  }
}

Here is the call graph for this function:

Member Function Documentation

void G4FermiBreakUpVI::BreakFragment ( G4FragmentVector *  theResult, G4Fragment *  theNucleus  )

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 93 of file G4FermiBreakUpVI.cc.

{
  if(verbose > 0) {
    G4cout << "### G4FermiBreakUpVI::BreakFragment start new fragment " << G4endl;
    G4cout << *theNucleus << G4endl;
  }
  frag.clear();
  lvect.clear();

  // initial fragment
  Z = theNucleus->GetZ_asInt();
  A = theNucleus->GetA_asInt();
  excitation = theNucleus->GetExcitationEnergy();
  mass = theNucleus->GetGroundStateMass() + excitation;
  spin = -1;
  G4double time = theNucleus->GetCreationTime();

  lv0 = theNucleus->GetMomentum();
  rndmEngine = G4Random::getTheEngine();

  // sample first decay of an initial state
  if(!SampleDecay()) {
    theResult->push_back(theNucleus);
    return; 
  }
  delete theNucleus;

  static const G4int imax = 100; 

  // loop over vector of Fermi fragments
  // vector may grouw at each iteraction
  for(size_t i=0; i<frag.size(); ++i) {
    Z = frag[i]->GetZ();
    A = frag[i]->GetA();
    spin = frag[i]->GetSpin();
    mass = frag[i]->GetTotalEnergy();
    excitation = frag[i]->GetExcitationEnergy();
    lv0 = lvect[i];
    if(verbose > 0) {
      G4cout << "# FermiFrag " << i << ".  Z= " << Z << " A= " << A 
         << " mass= " << mass << " exc= " << excitation << G4endl;
    }
    // stable fragment
    if(!SampleDecay()) {
      if(verbose > 0) { G4cout << "   New G4Fragment" << G4endl; }
      G4Fragment* f = new G4Fragment(A, Z, lv0);
      f->SetSpin(0.5*spin);
      f->SetCreationTime(time);
      theResult->push_back(f);
    }
    // limit the loop
    if(i == imax) {
      break;
    }
  }
}

Here is the call graph for this function:

void G4FermiBreakUpVI::Initialise ( )

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 67 of file G4FermiBreakUpVI.cc.

{
  if(thePool == nullptr) { InitialisePool(); }
  theDecay = thePool->FermiDecayProbability();
  elim = thePool->GetEnergyLimit();
  tolerance = thePool->GetTolerance();
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4FermiBreakUpVI::IsApplicable ( G4int  ZZ, G4int  AA, G4double  etot  ) const

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 88 of file G4FermiBreakUpVI.cc.

{
  return (ZZ < maxZ && AA < maxA && AA > 0 && eexc <= elim) ? true : false;
}

---

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

• geant4.10.03.p01/source/processes/hadronic/models/de_excitation/fermi_breakup/include/G4FermiBreakUpVI.hh
• geant4.10.03.p01/source/processes/hadronic/models/de_excitation/fermi_breakup/src/G4FermiBreakUpVI.cc