G4FermiBreakUp Class Reference

#include <G4FermiBreakUp.hh>

Inheritance diagram for G4FermiBreakUp:

Collaboration diagram for G4FermiBreakUp:

Public Member Functions
	G4FermiBreakUp ()
virtual	~G4FermiBreakUp ()
virtual void	Initialise () final
virtual G4bool	IsApplicable (G4int Z, G4int A, G4double mass) const final
virtual void	BreakFragment (G4FragmentVector *, G4Fragment *theNucleus) final
Public Member Functions inherited from G4VFermiBreakUp
	G4VFermiBreakUp ()
virtual	~G4VFermiBreakUp ()

Detailed Description

Definition at line 44 of file G4FermiBreakUp.hh.

Constructor & Destructor Documentation

G4FermiBreakUp::G4FermiBreakUp ( )

explicit

Definition at line 41 of file G4FermiBreakUp.cc.

                               : thePhaseSpace(nullptr)
{
  g4calc = G4Pow::GetInstance();
  Coef = ConstCoeff = 0.0;

  nmax = 16;
  NormalizedWeights.resize(nmax,0.0);
  massRes.reserve(4);
  frag.resize(4, 0);
  thePool = G4FermiFragmentsPool::Instance(); 
  Initialise();
}

Here is the call graph for this function:

G4FermiBreakUp::~G4FermiBreakUp ( )

virtual

Definition at line 54 of file G4FermiBreakUp.cc.

{}

Member Function Documentation

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

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 95 of file G4FermiBreakUp.cc.

{
  if(thePool == nullptr) { Initialise(); }
  // Calculate Momenta of K fragments
  G4double M = theNucleus->GetMomentum().m();
  const std::vector<const G4VFermiFragment*>* conf = 
    SelectConfiguration(theNucleus->GetZ_asInt(), 
            theNucleus->GetA_asInt(), M);

  // should never happen
  if(!conf) { 
    theResult->push_back(theNucleus);
    return; 
  }

  size_t nn = conf->size();

  // should never happen
  if(0 == nn) {
    theResult->push_back(theNucleus);
    return; 
  }

  G4LorentzVector fourMomentum = theNucleus->GetMomentum();

  // one unstable fragment
  if(1 == nn) {
    (*conf)[0]->FillFragment(theResult, fourMomentum);

    // normal case
  } else {
    G4ThreeVector boostVector = fourMomentum.boostVector();  
    massRes.clear();
    for(size_t i=0; i<nn; ++i) {
      massRes.push_back( (*conf)[i]->GetTotalEnergy() );
    }
    std::vector<G4LorentzVector*>* mom = thePhaseSpace->Decay(M, massRes);

    //  size_t nmom = mom->size();  
    // G4cout << "nmom= " << nmom << G4endl;

    // Go back to the Lab Frame
    for (size_t j=0; j<nn; ++j) {    
      (*mom)[j]->boost(boostVector); 
      (*conf)[j]->FillFragment(theResult, *((*mom)[j]));
      delete (*mom)[j];
    }
    delete mom;
  }
  delete theNucleus;
}

Here is the call graph for this function:

void G4FermiBreakUp::Initialise ( )

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 57 of file G4FermiBreakUp.cc.

{
  if(thePhaseSpace != nullptr) { return; }
  // Kappa = V/V_0 it is used in calculation of Coulomb energy
  // Nuclear radius r0 (is a model parameter)
  G4double Kappa = 6.0;
  G4double r0 = 1.3*CLHEP::fermi;

  Coef = 0.6*(CLHEP::elm_coupling/r0)/g4calc->Z13(1+G4int(Kappa));
  ConstCoeff = g4calc->powN(r0/hbarc,3)*Kappa/(6.0*pi*pi);

  thePhaseSpace = thePool->GetFermiPhaseSpaceDecay();
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4FermiBreakUp::IsApplicable ( G4int  Z, G4int  A, G4double  mass  ) const

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 71 of file G4FermiBreakUp.cc.

{
  return thePool->IsApplicable(Z, A, mass);
}

Here is the call graph for this function:

---

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

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