G4FermiBreakUpVI Class Reference

#include <G4FermiBreakUpVI.hh>

Inheritance diagram for G4FermiBreakUpVI:

Public Member Functions
	G4FermiBreakUpVI ()
	~G4FermiBreakUpVI () override
void	Initialise () override
G4bool	IsApplicable (G4int ZZ, G4int AA, G4double eexc) const override
void	BreakFragment (G4FragmentVector *, G4Fragment *theNucleus) override
	G4FermiBreakUpVI (const G4FermiBreakUpVI &right)=delete
const G4FermiBreakUpVI &	operator= (const G4FermiBreakUpVI &right)=delete
G4bool	operator== (const G4FermiBreakUpVI &right) const =delete
G4bool	operator!= (const G4FermiBreakUpVI &right) const =delete
Public Member Functions inherited from G4VFermiBreakUp
	G4VFermiBreakUp ()
virtual	~G4VFermiBreakUp ()=default
	G4VFermiBreakUp (const G4VFermiBreakUp &right)=delete
const G4VFermiBreakUp &	operator= (const G4VFermiBreakUp &right)=delete
G4bool	operator== (const G4VFermiBreakUp &right) const =delete
G4bool	operator!= (const G4VFermiBreakUp &right) const =delete
void	SetVerbose (G4int val)

Additional Inherited Members
Protected Attributes inherited from G4VFermiBreakUp
G4int	verbose {0}

Detailed Description

Definition at line 43 of file G4FermiBreakUpVI.hh.

Constructor & Destructor Documentation

G4FermiBreakUpVI() [1/2]

G4FermiBreakUpVI::G4FermiBreakUpVI

(

)

Definition at line 45 of file G4FermiBreakUpVI.cc.

{
  frag.reserve(10);
  lvect.reserve(10);
  secID = G4PhysicsModelCatalog::GetModelID("model_G4FermiBreakUpVI");
  prob.resize(12,0.0);
  if (nullptr == fPool) {
    fPool = new G4FermiFragmentsPoolVI();
    fPool->Initialise();
    isFirst = true;
  }
}

~G4FermiBreakUpVI()

G4FermiBreakUpVI::~G4FermiBreakUpVI ( )

override

Definition at line 58 of file G4FermiBreakUpVI.cc.

{
  if (isFirst) { 
    delete fPool;
    fPool = nullptr;
  }
}

G4FermiBreakUpVI() [2/2]

G4FermiBreakUpVI::G4FermiBreakUpVI ( const G4FermiBreakUpVI & right )

delete

Member Function Documentation

BreakFragment()

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

overridevirtual

Implements G4VFermiBreakUp.

Definition at line 85 of file G4FermiBreakUpVI.cc.

{
  if (verbose > 1) {
    G4cout << "### G4FermiBreakUpVI::BreakFragment start new fragment " 
           << G4endl;
    G4cout << *theNucleus << G4endl;
  }
  if (!fPool->IsInitialized()) { fPool->Initialise(); } 
 
  // initial fragment
  G4int Z = theNucleus->GetZ_asInt();
  G4int A = theNucleus->GetA_asInt();
  G4double excitation = theNucleus->GetExcitationEnergy();
  if (!IsApplicable(Z, A, excitation)) { return; }
  G4double mass = theNucleus->GetGroundStateMass() + excitation;
  G4LorentzVector lv0 = theNucleus->GetMomentum();
 
  // sample first decay of an initial state
  // if not possible to decay - exit
  if (!SampleDecay(Z, A, mass, excitation, lv0)) { 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 (std::size_t i=0; i<frag.size(); ++i) {
    Z = frag[i]->GetZ();
    A = frag[i]->GetA();
    excitation = frag[i]->GetExcitationEnergy();
    lv0 = lvect[i];
    G4bool unstable = (IsApplicable(Z, A, excitation) && frag[i]->GetLifeTime() < fTimeLim);
    if (unstable) {
      mass = frag[i]->GetTotalEnergy();
      if (verbose > 1) {
    G4cout << "# FermiFrag " << i << ".  Z= " << Z << " A= " << A 
           << " mass= " << mass << " exc= " 
           << frag[i]->GetExcitationEnergy() << G4endl;
      }
      unstable = SampleDecay(Z, A, mass, excitation, lv0);
    }
    // stable fragment
    if (!unstable) {
      if(verbose > 1) { G4cout << "   New G4Fragment" << G4endl; }
      G4Fragment* f = new G4Fragment(A, Z, lv0);
      f->SetCreationTime(time);
      f->SetCreatorModelID(secID);
      theResult->push_back(f);
    }
    // limit the loop
    if (i == imax) { break; }
  }
  frag.clear();
  lvect.clear();
}

Initialise()

void G4FermiBreakUpVI::Initialise ( )

overridevirtual

Implements G4VFermiBreakUp.

Definition at line 66 of file G4FermiBreakUpVI.cc.

{
  G4DeexPrecoParameters* param = 
    G4NuclearLevelData::GetInstance()->GetParameters();
  fTolerance = param->GetMinExcitation();
  fElim = param->GetFBUEnergyLimit();
  fTimeLim = param->GetMaxLifeTime();
  if (verbose > 1) {
    G4cout << "### G4FermiBreakUpVI::Initialise(): the pool is initilized=" 
       << fPool->IsInitialized() << " fTolerance(eV)=" << fTolerance/CLHEP::eV
           << " Elim(MeV)=" << fElim/CLHEP::MeV << G4endl;
  }
}

IsApplicable()

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

overridevirtual

Implements G4VFermiBreakUp.

Definition at line 80 of file G4FermiBreakUpVI.cc.

{
  return (Z < maxZ && A < maxA && eexc <= fElim && fPool->HasDecay(Z, A, eexc));
}

Referenced by BreakFragment().

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

---

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

• G4FermiBreakUpVI.hh
• G4FermiBreakUpVI.cc