G4FermiBreakUpVI Class Reference

#include <G4FermiBreakUpVI.hh>

Inheritance diagram for G4FermiBreakUpVI:

Public Member Functions
	G4FermiBreakUpVI ()
	~G4FermiBreakUpVI () final
void	Initialise () final
G4bool	IsApplicable (G4int ZZ, G4int AA, G4double etot) const final
void	BreakFragment (G4FragmentVector *, G4Fragment *theNucleus) final
Public Member Functions inherited from G4VFermiBreakUp
	G4VFermiBreakUp ()
virtual	~G4VFermiBreakUp ()
virtual void	Initialise ()=0
virtual G4bool	IsApplicable (G4int Z, G4int A, G4double mass) const =0
virtual void	BreakFragment (G4FragmentVector *results, G4Fragment *theNucleus)=0
void	SetVerbose (G4int val)

Additional Inherited Members
Protected Attributes inherited from G4VFermiBreakUp
G4int	verbose

Detailed Description

Definition at line 46 of file G4FermiBreakUpVI.hh.

Constructor & Destructor Documentation

G4FermiBreakUpVI()

G4FermiBreakUpVI::G4FermiBreakUpVI ( )

explicit

Definition at line 46 of file G4FermiBreakUpVI.cc.

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

~G4FermiBreakUpVI()

G4FermiBreakUpVI::~G4FermiBreakUpVI ( )

final

Definition at line 59 of file G4FermiBreakUpVI.cc.

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

Member Function Documentation

BreakFragment()

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

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 96 of file G4FermiBreakUpVI.cc.

{
  if(verbose > 1) {
    G4cout << "### G4FermiBreakUpVI::BreakFragment start new fragment " 
           << G4endl;
    G4cout << *theNucleus << G4endl;
  }
 
  // initial fragment
  Z = theNucleus->GetZ_asInt();
  A = theNucleus->GetA_asInt();
  excitation = theNucleus->GetExcitationEnergy();
  mass = theNucleus->GetGroundStateMass() + excitation;
  spin = -1;
 
  lv0 = theNucleus->GetMomentum();
  rndmEngine = G4Random::getTheEngine();
 
  // sample first decay of an initial state
  // if not possible to decay - exit
  if(!SampleDecay()) {
    return; 
  }
 
  G4double time = theNucleus->GetCreationTime();
  delete theNucleus;
 
  static const G4int imax = 100; 
 
  // loop over vector of Fermi fragments
  // vector may grow 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();
    lv0 = lvect[i];
    if(verbose > 1) {
      G4cout << "# FermiFrag " << i << ".  Z= " << Z << " A= " << A 
         << " mass= " << mass << " exc= " 
         << frag[i]->GetExcitationEnergy() << G4endl;
    }
    // stable fragment
    if(!SampleDecay()) {
      if(verbose > 1) { 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;
    }
  }
  frag.clear();
  lvect.clear();
}

Initialise()

void G4FermiBreakUpVI::Initialise ( )

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 67 of file G4FermiBreakUpVI.cc.

{
  if(verbose > 1) {
    G4cout << "### G4FermiBreakUpVI::Initialise(): " << thePool << G4endl;
  }
  if(thePool == nullptr) { InitialisePool(); }
  theDecay = thePool->FermiDecayProbability();
  elim = thePool->GetEnergyLimit();
}

Referenced by G4FermiBreakUpVI().

IsApplicable()

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

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 90 of file G4FermiBreakUpVI.cc.

{
  return (ZZ < maxZ && AA < maxA && AA > 0 && eexc <= elim 
      && thePool->HasChannels(ZZ, AA, eexc));
}

---

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

• G4FermiBreakUpVI.hh
• G4FermiBreakUpVI.cc