G4NeutronHPNBodyPhaseSpace Class Reference

#include <G4NeutronHPNBodyPhaseSpace.hh>

Inheritance diagram for G4NeutronHPNBodyPhaseSpace:

Public Member Functions
	G4NeutronHPNBodyPhaseSpace ()
	~G4NeutronHPNBodyPhaseSpace ()
void	Init (G4double aMass, G4int aCount)
void	Init (std::istream &aDataFile)
G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)
Public Member Functions inherited from G4VNeutronHPEnergyAngular
	G4VNeutronHPEnergyAngular ()
virtual	~G4VNeutronHPEnergyAngular ()
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 40 of file G4NeutronHPNBodyPhaseSpace.hh.

Constructor & Destructor Documentation

G4NeutronHPNBodyPhaseSpace::G4NeutronHPNBodyPhaseSpace ( )

inline

Definition at line 44 of file G4NeutronHPNBodyPhaseSpace.hh.

{}

G4NeutronHPNBodyPhaseSpace::~G4NeutronHPNBodyPhaseSpace ( )

inline

Definition at line 45 of file G4NeutronHPNBodyPhaseSpace.hh.

{}

Member Function Documentation

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

inline

Definition at line 49 of file G4NeutronHPNBodyPhaseSpace.hh.

Referenced by G4NeutronHPInelasticBaseFS::BaseApply().

  {
    theTotalMass=aMass;
    theTotalCount=aCount;
  }

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

inlinevirtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 55 of file G4NeutronHPNBodyPhaseSpace.hh.

References G4ParticleDefinition::GetPDGMass(), and G4Neutron::Neutron().

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

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

virtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 42 of file G4NeutronHPNBodyPhaseSpace.cc.

References G4Alpha::Alpha(), G4Deuteron::Deuteron(), G4Electron::Electron(), energy(), G4UniformRand, G4Gamma::Gamma(), G4ReactionProduct::GetMass(), G4VNeutronHPEnergyAngular::GetNeutron(), G4VNeutronHPEnergyAngular::GetTarget(), G4ReactionProduct::GetTotalMomentum(), G4He3::He3(), G4ReactionProduct::Lorentz(), G4INCL::Math::max(), G4Neutron::Neutron(), G4Positron::Positron(), G4Proton::Proton(), G4ReactionProduct::SetDefinition(), G4ReactionProduct::SetKineticEnergy(), G4ReactionProduct::SetMomentum(), G4Triton::Triton(), and python.hepunit::twopi.

Referenced by G4NeutronHPInelasticBaseFS::BaseApply().

{
   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;
}

---

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

• G4NeutronHPNBodyPhaseSpace.hh
• G4NeutronHPNBodyPhaseSpace.cc