G4NeutronHPNBodyPhaseSpace Class Reference

#include <G4NeutronHPNBodyPhaseSpace.hh>

Inheritance diagram for G4NeutronHPNBodyPhaseSpace:

Public Member Functions
	G4NeutronHPNBodyPhaseSpace ()
	~G4NeutronHPNBodyPhaseSpace ()
void	Init (G4double aMass, G4int aCount)
void	Init (std::ifstream &aDataFile)
G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)
Public Member Functions inherited from G4VNeutronHPEnergyAngular
	G4VNeutronHPEnergyAngular ()
virtual	~G4VNeutronHPEnergyAngular ()
virtual void	Init (std::ifstream &aDataFile)=0
virtual G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)=0
virtual G4double	MeanEnergyOfThisInteraction ()=0
void	SetNeutron (G4ReactionProduct *aNeutron)
void	SetTarget (G4ReactionProduct *aTarget)
G4ReactionProduct *	GetTarget ()
G4ReactionProduct *	GetNeutron ()
G4ReactionProduct *	GetCMS ()
void	SetQValue (G4double aValue)
virtual void	ClearHistories ()

Additional Inherited Members
Protected Member Functions inherited from G4VNeutronHPEnergyAngular
G4double	GetQValue ()

Detailed Description

Definition at line 41 of file G4NeutronHPNBodyPhaseSpace.hh.

Constructor & Destructor Documentation

G4NeutronHPNBodyPhaseSpace()

G4NeutronHPNBodyPhaseSpace::G4NeutronHPNBodyPhaseSpace ( )

inline

Definition at line 45 of file G4NeutronHPNBodyPhaseSpace.hh.

{}

~G4NeutronHPNBodyPhaseSpace()

G4NeutronHPNBodyPhaseSpace::~G4NeutronHPNBodyPhaseSpace ( )

inline

Definition at line 46 of file G4NeutronHPNBodyPhaseSpace.hh.

{}

Member Function Documentation

Init() [1/2]

void G4NeutronHPNBodyPhaseSpace::Init ( G4double  aMass, G4int  aCount  )

inline

Definition at line 50 of file G4NeutronHPNBodyPhaseSpace.hh.

  {
    theTotalMass=aMass;
    theTotalCount=aCount;
  }

Referenced by G4NeutronHPInelasticBaseFS::BaseApply().

Init() [2/2]

void G4NeutronHPNBodyPhaseSpace::Init ( std::ifstream &  aDataFile )

inlinevirtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 56 of file G4NeutronHPNBodyPhaseSpace.hh.

  {
    aDataFile >> theTotalMass >> theTotalCount;
    theTotalMass *= G4Neutron::Neutron()->GetPDGMass();
  }

Sample()

G4ReactionProduct * G4NeutronHPNBodyPhaseSpace::Sample ( G4double  anEnergy, G4double  massCode, G4double  mass  )

virtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 43 of file G4NeutronHPNBodyPhaseSpace.cc.

{
   G4ReactionProduct * result = new G4ReactionProduct;
   G4int Z = static_cast<G4int>(massCode/1000);
   G4int A = static_cast<G4int>(massCode-1000*Z);
 
   if(massCode==0)
   {
     result->SetDefinition(G4Gamma::Gamma());
   }
   else if(A==0)
   {
     result->SetDefinition(G4Electron::Electron());     
     if(Z==1) result->SetDefinition(G4Positron::Positron());
   }
   else if(A==1)
   {
     result->SetDefinition(G4Neutron::Neutron());
     if(Z==1) result->SetDefinition(G4Proton::Proton());
   }
   else if(A==2)
   {
     result->SetDefinition(G4Deuteron::Deuteron());      
   }
   else if(A==3)
   {
     result->SetDefinition(G4Triton::Triton());  
     if(Z==2) result->SetDefinition(G4He3::He3());
   }
   else if(A==4)
   {
     result->SetDefinition(G4Alpha::Alpha());
     if(Z!=2) throw G4HadronicException(__FILE__, __LINE__, "Unknown ion case 1");    
   }
   else
   {
     throw G4HadronicException(__FILE__, __LINE__, "G4NeutronHPNBodyPhaseSpace: Unknown ion case 2");
   }
 
// Get the energy from phase-space distribution
   // in CMS
   // P = Cn*std::sqrt(E')*(Emax-E')**(3*n/2-4)
   G4double maxE = GetEmax(anEnergy, result->GetMass());
   G4double energy;
   G4double max(0);
   if(theTotalCount<=3)
   {
     max = maxE/2.;
   }
   else if(theTotalCount==4)
   {
     max = maxE/5.;
   }
   else if(theTotalCount==5)
   {
     max = maxE/8.;
   }
   else
   {
     throw G4HadronicException(__FILE__, __LINE__, "NeutronHP Phase-space distribution cannot cope with this number of particles");
   }
   G4double testit;
   G4double rand0 = Prob(max, maxE, theTotalCount);
   G4double rand;
   
   do
   {
     rand = rand0*G4UniformRand();
     energy = maxE*G4UniformRand();
     testit = Prob(energy, maxE, theTotalCount);
   }
   while(rand > testit);
   result->SetKineticEnergy(energy);
   
// now do random direction
   G4double cosTh = 2.*G4UniformRand()-1.;
   G4double phi = twopi*G4UniformRand();
   G4double theta = std::acos(cosTh);
   G4double sinth = std::sin(theta);
   G4double mtot = result->GetTotalMomentum(); 
   G4ThreeVector tempVector(mtot*sinth*std::cos(phi), mtot*sinth*std::sin(phi), mtot*std::cos(theta) );
   result->SetMomentum(tempVector);
   G4ReactionProduct aCMS = *GetTarget()+*GetNeutron();
   result->Lorentz(*result, -1.*aCMS);
   return result;
}

Referenced by G4NeutronHPInelasticBaseFS::BaseApply().

---

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

• G4NeutronHPNBodyPhaseSpace.hh
• G4NeutronHPNBodyPhaseSpace.cc