G4PreCompoundEmission Class Reference

#include <G4PreCompoundEmission.hh>

Public Member Functions
	G4PreCompoundEmission ()
	~G4PreCompoundEmission ()
void	SetDefaultModel ()
void	SetHETCModel ()
G4ReactionProduct *	PerformEmission (G4Fragment &aFragment)
G4double	GetTotalProbability (const G4Fragment &aFragment)
void	Initialize (const G4Fragment &aFragment)
void	SetOPTxs (G4int)
void	UseSICB (G4bool)

Detailed Description

Definition at line 50 of file G4PreCompoundEmission.hh.

Constructor & Destructor Documentation

G4PreCompoundEmission::G4PreCompoundEmission

(

)

Definition at line 57 of file G4PreCompoundEmission.cc.

References G4PreCompoundParameters::GetFermiEnergy(), G4VPreCompoundEmissionFactory::GetFragmentVector(), G4Pow::GetInstance(), and G4PreCompoundParameters::GetLevelDensity().

{
  theFragmentsFactory = new G4PreCompoundEmissionFactory();
  theFragmentsVector = 
    new G4PreCompoundFragmentVector(theFragmentsFactory->GetFragmentVector());
  g4pow = G4Pow::GetInstance();
  G4PreCompoundParameters param;
  fLevelDensity = param.GetLevelDensity();
  fFermiEnergy  = param.GetFermiEnergy();
}

G4PreCompoundEmission::~G4PreCompoundEmission

(

)

Definition at line 68 of file G4PreCompoundEmission.cc.

{
  delete theFragmentsFactory; 
  delete theFragmentsVector; 
}

Member Function Documentation

G4double G4PreCompoundEmission::GetTotalProbability ( const G4Fragment &  aFragment )

inline

Definition at line 103 of file G4PreCompoundEmission.hh.

References G4PreCompoundFragmentVector::CalculateProbabilities().

Referenced by G4PreCompoundModel::DeExcite().

{
  return theFragmentsVector->CalculateProbabilities(aFragment);
}

void G4PreCompoundEmission::Initialize ( const G4Fragment &  aFragment )

inline

Definition at line 109 of file G4PreCompoundEmission.hh.

References G4PreCompoundFragmentVector::Initialize().

Referenced by G4PreCompoundModel::DeExcite().

{
  theFragmentsVector->Initialize(aFragment);
}

G4ReactionProduct * G4PreCompoundEmission::PerformEmission

(

G4Fragment &

aFragment

)

Definition at line 99 of file G4PreCompoundEmission.cc.

References CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4PreCompoundFragmentVector::ChooseFragment(), G4cout, G4endl, G4VPreCompoundFragment::GetA(), G4VPreCompoundFragment::GetKineticEnergy(), G4Fragment::GetMomentum(), G4VPreCompoundFragment::GetNuclearMass(), G4Fragment::GetNumberOfCharged(), G4Fragment::GetNumberOfParticles(), G4VPreCompoundFragment::GetReactionProduct(), G4VPreCompoundFragment::GetRestA(), G4VPreCompoundFragment::GetRestZ(), G4VPreCompoundFragment::GetZ(), G4Fragment::SetMomentum(), G4VPreCompoundFragment::SetMomentum(), G4Fragment::SetNumberOfCharged(), G4Fragment::SetNumberOfParticles(), and G4Fragment::SetZandA_asInt().

Referenced by G4PreCompoundModel::DeExcite().

{
  // Choose a Fragment for emission
  G4VPreCompoundFragment * thePreFragment = 
    theFragmentsVector->ChooseFragment();
  if (thePreFragment == 0)
    {
      G4cout << "G4PreCompoundEmission::PerformEmission : "
         << "I couldn't choose a fragment\n"
         << "while trying to de-excite\n" 
         << aFragment << G4endl;
      throw G4HadronicException(__FILE__, __LINE__, "");
    }

  //G4cout << "Chosen fragment: " << G4endl;
  //G4cout << *thePreFragment << G4endl;

  // Kinetic Energy of emitted fragment
  G4double kinEnergyOfEmittedFragment = 
    thePreFragment->GetKineticEnergy(aFragment);
  //  if(kinEnergyOfEmittedFragment < MeV) {
  //  G4cout << "Chosen fragment: " << G4endl;
  //  G4cout << *thePreFragment << G4endl;
  //  G4cout << "Ekin= " << kinEnergyOfEmittedFragment << G4endl;
    // }
  if(kinEnergyOfEmittedFragment < 0.0) { kinEnergyOfEmittedFragment = 0.0; }
  
  // Calculate the fragment momentum (three vector)
  AngularDistribution(thePreFragment,aFragment,kinEnergyOfEmittedFragment);
  
  // Mass of emittef fragment
  G4double EmittedMass = thePreFragment->GetNuclearMass();
  // Now we can calculate the four momentum 
  // both options are valid and give the same result but 2nd one is faster
  G4LorentzVector Emitted4Momentum(theFinalMomentum,
                   EmittedMass + kinEnergyOfEmittedFragment);
    
  // Perform Lorentz boost
  G4LorentzVector Rest4Momentum = aFragment.GetMomentum();
  Emitted4Momentum.boost(Rest4Momentum.boostVector());  

  // Set emitted fragment momentum
  thePreFragment->SetMomentum(Emitted4Momentum);    

  // NOW THE RESIDUAL NUCLEUS
  // ------------------------

  Rest4Momentum -= Emitted4Momentum;
    
  // Update nucleus parameters:
  // --------------------------

  // Z and A
  aFragment.SetZandA_asInt(thePreFragment->GetRestZ(),
               thePreFragment->GetRestA());
    
  // Number of excitons
  aFragment.SetNumberOfParticles(aFragment.GetNumberOfParticles()-
                 thePreFragment->GetA());
  // Number of charges
  aFragment.SetNumberOfCharged(aFragment.GetNumberOfCharged()-
                   thePreFragment->GetZ());
    
  // Update nucleus momentum 
  // A check on consistence of Z, A, and mass will be performed
  aFragment.SetMomentum(Rest4Momentum);
    
  // Create a G4ReactionProduct 
  G4ReactionProduct * MyRP = thePreFragment->GetReactionProduct();

  //  if(kinEnergyOfEmittedFragment < MeV) {
  //  G4cout << "G4PreCompoundEmission::Fragment emitted" << G4endl;
  //  G4cout << thePreFragment << G4endl;
    // }
  return MyRP;
}

void G4PreCompoundEmission::SetDefaultModel

(

)

Definition at line 74 of file G4PreCompoundEmission.cc.

References G4VPreCompoundEmissionFactory::GetFragmentVector(), and G4PreCompoundFragmentVector::SetVector().

Referenced by G4PreCompoundModel::G4PreCompoundModel(), and G4PreCompoundModel::UseDefaultEmission().

{
  if (theFragmentsFactory) { delete theFragmentsFactory; }
  theFragmentsFactory = new G4PreCompoundEmissionFactory();
  if (theFragmentsVector) {
    theFragmentsVector->SetVector(theFragmentsFactory->GetFragmentVector());
  } else {
    theFragmentsVector = 
      new G4PreCompoundFragmentVector(theFragmentsFactory->GetFragmentVector());
  }
}

void G4PreCompoundEmission::SetHETCModel

(

)

Definition at line 86 of file G4PreCompoundEmission.cc.

References G4VPreCompoundEmissionFactory::GetFragmentVector(), and G4PreCompoundFragmentVector::SetVector().

Referenced by G4PreCompoundModel::G4PreCompoundModel(), and G4PreCompoundModel::UseHETCEmission().

{
  if (theFragmentsFactory) delete theFragmentsFactory;
  theFragmentsFactory = new G4HETCEmissionFactory();
  if (theFragmentsVector) {
    theFragmentsVector->SetVector(theFragmentsFactory->GetFragmentVector());
  } else {
    theFragmentsVector = 
      new G4PreCompoundFragmentVector(theFragmentsFactory->GetFragmentVector());
  }
}

void G4PreCompoundEmission::SetOPTxs ( G4int  opt )

inline

Definition at line 114 of file G4PreCompoundEmission.hh.

References G4PreCompoundFragmentVector::SetOPTxs().

Referenced by G4PreCompoundModel::G4PreCompoundModel(), and G4PreCompoundModel::SetOPTxs().

{
  theFragmentsVector->SetOPTxs(opt);
}

void G4PreCompoundEmission::UseSICB ( G4bool  use )

inline

Definition at line 119 of file G4PreCompoundEmission.hh.

References G4PreCompoundFragmentVector::UseSICB().

Referenced by G4PreCompoundModel::G4PreCompoundModel(), and G4PreCompoundModel::UseSICB().

{
  theFragmentsVector->UseSICB(use);
}

---

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

• G4PreCompoundEmission.hh
• G4PreCompoundEmission.cc