G4ParticleHPNBodyPhaseSpace Class Reference

#include <G4ParticleHPNBodyPhaseSpace.hh>

Inheritance diagram for G4ParticleHPNBodyPhaseSpace:

Collaboration diagram for G4ParticleHPNBodyPhaseSpace:

Public Member Functions
	G4ParticleHPNBodyPhaseSpace ()
	~G4ParticleHPNBodyPhaseSpace ()
void	Init (G4double aMass, G4int aCount)
void	Init (std::istream &aDataFile)
G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)
Public Member Functions inherited from G4VParticleHPEnergyAngular
	G4VParticleHPEnergyAngular ()
virtual	~G4VParticleHPEnergyAngular ()
void	SetProjectileRP (G4ReactionProduct *aIncidentParticleRP)
void	SetTarget (G4ReactionProduct *aTarget)
G4ReactionProduct *	GetTarget ()
G4ReactionProduct *	GetProjectileRP ()
G4ReactionProduct *	GetCMS ()
void	SetQValue (G4double aValue)
virtual void	ClearHistories ()

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

Detailed Description

Definition at line 42 of file G4ParticleHPNBodyPhaseSpace.hh.

Constructor & Destructor Documentation

G4ParticleHPNBodyPhaseSpace::G4ParticleHPNBodyPhaseSpace ( )

inline

Definition at line 46 of file G4ParticleHPNBodyPhaseSpace.hh.

                               {
   theTotalMass = 0.0;
   theTotalCount = 0;
  }

G4ParticleHPNBodyPhaseSpace::~G4ParticleHPNBodyPhaseSpace ( )

inline

Definition at line 50 of file G4ParticleHPNBodyPhaseSpace.hh.

{}

Member Function Documentation

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

inline

Definition at line 54 of file G4ParticleHPNBodyPhaseSpace.hh.

  {
    theTotalMass=aMass;
    theTotalCount=aCount;
  }

Here is the caller graph for this function:

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

inlinevirtual

Implements G4VParticleHPEnergyAngular.

Definition at line 60 of file G4ParticleHPNBodyPhaseSpace.hh.

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

Here is the call graph for this function:

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

virtual

Implements G4VParticleHPEnergyAngular.

Definition at line 44 of file G4ParticleHPNBodyPhaseSpace.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__, "G4ParticleHPNBodyPhaseSpace: 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;
   
   G4int icounter=0;
   G4int icounter_max=1024;
   do
   {
      icounter++;
      if ( icounter > icounter_max ) {
         G4cout << "Loop-counter exceeded the threshold value at " << __LINE__ << "th line of " << __FILE__ << "." << G4endl;
         break;
      }
     rand = rand0*G4UniformRand();
     energy = maxE*G4UniformRand();
     testit = Prob(energy, maxE, theTotalCount);
   }
   while(rand > testit); // Loop checking, 11.05.2015, T. Koi
   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()+*GetProjectileRP();
   result->Lorentz(*result, -1.*aCMS);
   return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/particle_hp/include/G4ParticleHPNBodyPhaseSpace.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/particle_hp/src/G4ParticleHPNBodyPhaseSpace.cc