G4FermiConfigurationList Class Reference

#include <G4FermiConfigurationList.hh>

Public Member Functions
	G4FermiConfigurationList ()
	~G4FermiConfigurationList ()
G4FragmentVector *	GetFragments (const G4Fragment &theNucleus)

Detailed Description

Definition at line 47 of file G4FermiConfigurationList.hh.

Constructor & Destructor Documentation

G4FermiConfigurationList::G4FermiConfigurationList

(

)

Definition at line 47 of file G4FermiConfigurationList.cc.

References G4Pow::GetInstance(), python.hepunit::hbarc, G4FermiFragmentsPool::Instance(), python.hepunit::pi, G4Pow::powN(), and G4Pow::Z13().

{
  thePool = G4FermiFragmentsPool::Instance(); 
  g4pow = G4Pow::GetInstance();
  Coef = 0.6*(CLHEP::elm_coupling/r0)/g4pow->Z13(1+G4int(Kappa));
  ConstCoeff = g4pow->powN(r0/hbarc,3)*Kappa*std::sqrt(2.0/pi)/3.0;

  // 16 is the max number
  nmax = 50;
  NormalizedWeights.resize(nmax,0.0);
}

G4FermiConfigurationList::~G4FermiConfigurationList

(

)

Definition at line 59 of file G4FermiConfigurationList.cc.

{}

Member Function Documentation

G4FragmentVector * G4FermiConfigurationList::GetFragments

(

const G4Fragment &

theNucleus

)

Definition at line 175 of file G4FermiConfigurationList.cc.

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4FermiPhaseSpaceDecay::Decay(), G4Fragment::GetA_asInt(), G4Fragment::GetMomentum(), G4Fragment::GetZ_asInt(), CLHEP::HepLorentzVector::m(), and G4InuclParticleNames::nn.

Referenced by G4FermiBreakUp::BreakItUp().

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


  G4FragmentVector* theResult = new G4FragmentVector();
  size_t nn = conf->size();
  if(1 >= nn) {
    theResult->push_back(new G4Fragment(theNucleus));
    delete conf;
    return theResult; 
  }

  G4ThreeVector boostVector = theNucleus.GetMomentum().boostVector();  
  std::vector<G4double> mr;
  mr.reserve(nn);
  for(size_t i=0; i<nn; ++i) {
    mr.push_back( (*conf)[i]->GetTotalEnergy() );
  }
  std::vector<G4LorentzVector*>* mom = thePhaseSpace.Decay(M,mr);
  if(!mom) { 
    delete conf;
    return theResult; 
  }

  size_t nmom = mom->size();  

  // Go back to the Lab Frame
  if(0 < nmom) { 
    for (size_t j=0; j<nmom; ++j) {
      G4LorentzVector* FourMomentum = (*mom)[j];
    
      // Lorentz boost
      FourMomentum->boost(boostVector);
      
      G4FragmentVector* fragment = (*conf)[j]->GetFragment(*FourMomentum);
 
      size_t nfrag = fragment->size();
      for (size_t k=0; k<nfrag; ++k) { theResult->push_back((*fragment)[k]); }
      delete fragment;
      delete (*mom)[j];
    }
  }
  
  delete mom;
  delete conf;
  return theResult;
}

---

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

• G4FermiConfigurationList.hh
• G4FermiConfigurationList.cc