#include <G4ParticleHPNBodyPhaseSpace.hh>

Inheritance diagram for G4ParticleHPNBodyPhaseSpace:

Public Member Functions
virtual void	ClearHistories ()

	G4ParticleHPNBodyPhaseSpace ()

G4ReactionProduct *	GetCMS ()

G4ReactionProduct *	GetProjectileRP ()

G4ReactionProduct *	GetTarget ()

void	Init (G4double aMass, G4int aCount)

void	Init (std::istream &aDataFile)

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

void	SetProjectileRP (G4ReactionProduct *aIncidentParticleRP)

void	SetQValue (G4double aValue)

void	SetTarget (G4ReactionProduct *aTarget)

	~G4ParticleHPNBodyPhaseSpace ()

Protected Member Functions
G4double	GetQValue ()

Private Member Functions
G4double	C (G4double anEnergy, G4double mass)

G4double	GetEmax (G4double anEnergy, G4double mass)

G4double	MeanEnergyOfThisInteraction ()

G4double	Prob (G4double anEnergy, G4double eMax, G4int n)

Private Attributes
G4Cache< toBeCached >	fCache

G4double	theQValue

G4int	theTotalCount

G4double	theTotalMass

Detailed Description

Definition at line 42 of file G4ParticleHPNBodyPhaseSpace.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPNBodyPhaseSpace()

G4ParticleHPNBodyPhaseSpace::G4ParticleHPNBodyPhaseSpace ( )

inline

Definition at line 46 of file G4ParticleHPNBodyPhaseSpace.hh.

                               {
   theTotalMass = 0.0;
   theTotalCount = 0;
  }

References theTotalCount, and theTotalMass.

◆ ~G4ParticleHPNBodyPhaseSpace()

G4ParticleHPNBodyPhaseSpace::~G4ParticleHPNBodyPhaseSpace ( )

inline

Definition at line 50 of file G4ParticleHPNBodyPhaseSpace.hh.

50{}

Member Function Documentation

◆ C()

G4double G4ParticleHPNBodyPhaseSpace::C	(	G4double	anEnergy,
		G4double	mass
	)

inlineprivate

Definition at line 77 of file G4ParticleHPNBodyPhaseSpace.hh.

  {
    G4double result(0);
    if(theTotalCount==3) result = 4./CLHEP::pi/G4Pow::GetInstance()->powN(GetEmax(anEnergy, mass),2);
    if(theTotalCount==4) result = 105./32./G4Pow::GetInstance()->powA(GetEmax(anEnergy, mass), 3.5);
    //if(theTotalCount==5) result = 256./14./CLHEP::pi/G4Pow::GetInstance()->powA(GetEmax(anEnergy, mass), 5.);
    if(theTotalCount==5) result = 256./14./CLHEP::pi/G4Pow::GetInstance()->powN(GetEmax(anEnergy, mass), 5);
    return result;
  }

References GetEmax(), G4Pow::GetInstance(), CLHEP::pi, G4Pow::powA(), G4Pow::powN(), and theTotalCount.

◆ ClearHistories()

virtual void G4VParticleHPEnergyAngular::ClearHistories ( )

inlinevirtualinherited

Reimplemented in G4ParticleHPContEnergyAngular.

Definition at line 106 of file G4VParticleHPEnergyAngular.hh.

107 {

108 }

Referenced by G4ParticleHPProduct::Sample().

◆ GetCMS()

G4ReactionProduct * G4VParticleHPEnergyAngular::GetCMS ( )

inlineinherited

Definition at line 94 of file G4VParticleHPEnergyAngular.hh.

  { 
     *fCache.Get().theCMS = *fCache.Get().theProjectileRP
                          + *fCache.Get().theTarget;
     return fCache.Get().theCMS;
  }

References G4VParticleHPEnergyAngular::fCache.

◆ GetEmax()

G4double G4ParticleHPNBodyPhaseSpace::GetEmax	(	G4double	anEnergy,
		G4double	mass
	)

inlineprivate

Definition at line 87 of file G4ParticleHPNBodyPhaseSpace.hh.

  {
    G4double result;
    G4double tMass = GetTarget()->GetMass();
    G4double pMass = GetProjectileRP()->GetMass();
    G4double availableEnergy = GetQValue() + anEnergy*tMass/(pMass+tMass);
    result = availableEnergy*(theTotalMass-mass)/theTotalMass;
    return result;
  }

References G4ReactionProduct::GetMass(), G4VParticleHPEnergyAngular::GetProjectileRP(), G4VParticleHPEnergyAngular::GetQValue(), G4VParticleHPEnergyAngular::GetTarget(), and theTotalMass.

Referenced by C(), and Sample().

◆ GetProjectileRP()

G4ReactionProduct * G4VParticleHPEnergyAngular::GetProjectileRP ( )

inlineinherited

Definition at line 89 of file G4VParticleHPEnergyAngular.hh.

  {
    return fCache.Get().theProjectileRP;
  }

References G4VParticleHPEnergyAngular::fCache.

Referenced by GetEmax(), G4ParticleHPContEnergyAngular::Sample(), G4ParticleHPDiscreteTwoBody::Sample(), and Sample().

◆ GetQValue()

G4double G4VParticleHPEnergyAngular::GetQValue ( )

inlineprotectedinherited

Definition at line 112 of file G4VParticleHPEnergyAngular.hh.

112{ return theQValue; }

G4VParticleHPEnergyAngular::theQValue

G4double theQValue

Definition: G4VParticleHPEnergyAngular.hh:116

References G4VParticleHPEnergyAngular::theQValue.

Referenced by GetEmax(), G4ParticleHPDiscreteTwoBody::Sample(), and G4ParticleHPIsotropic::Sample().

◆ GetTarget()

G4ReactionProduct * G4VParticleHPEnergyAngular::GetTarget ( void )

inlineinherited

Definition at line 84 of file G4VParticleHPEnergyAngular.hh.

  {
    return fCache.Get().theTarget;
  }

References G4VParticleHPEnergyAngular::fCache.

Referenced by GetEmax(), G4ParticleHPContEnergyAngular::Sample(), G4ParticleHPDiscreteTwoBody::Sample(), and Sample().

◆ Init() [1/2]

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

inline

Definition at line 54 of file G4ParticleHPNBodyPhaseSpace.hh.

  {
    theTotalMass=aMass;
    theTotalCount=aCount;
  }

References theTotalCount, and theTotalMass.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply().

◆ Init() [2/2]

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

inlinevirtual

Implements G4VParticleHPEnergyAngular.

Definition at line 60 of file G4ParticleHPNBodyPhaseSpace.hh.

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

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

◆ MeanEnergyOfThisInteraction()

G4double G4ParticleHPNBodyPhaseSpace::MeanEnergyOfThisInteraction ( )

inlineprivatevirtual

Implements G4VParticleHPEnergyAngular.

Definition at line 97 of file G4ParticleHPNBodyPhaseSpace.hh.

97{return -1; }

◆ Prob()

G4double G4ParticleHPNBodyPhaseSpace::Prob	(	G4double	anEnergy,
		G4double	eMax,
		G4int	n
	)

inlineprivate

Definition at line 70 of file G4ParticleHPNBodyPhaseSpace.hh.

  {
    G4double result;
    result = std::sqrt(anEnergy)*G4Pow::GetInstance()->powA(eMax-anEnergy, 3.*n/2.-4.);
    return result;
  }

References G4Pow::GetInstance(), CLHEP::detail::n, and G4Pow::powA().

Referenced by Sample().

◆ Sample()

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());
   if(maxE<=0){
     maxE=1.*CLHEP::eV;
   }
   G4double energy=0.;
   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;
}

References A, G4Alpha::Alpha(), G4Deuteron::Deuteron(), G4Electron::Electron(), G4INCL::KinematicsUtils::energy(), CLHEP::eV, G4cout, G4endl, G4UniformRand, G4Gamma::Gamma(), GetEmax(), G4ReactionProduct::GetMass(), G4VParticleHPEnergyAngular::GetProjectileRP(), G4VParticleHPEnergyAngular::GetTarget(), G4ReactionProduct::GetTotalMomentum(), G4He3::He3(), G4ReactionProduct::Lorentz(), G4INCL::Math::max(), G4Neutron::Neutron(), G4Positron::Positron(), Prob(), G4Proton::Proton(), G4ReactionProduct::SetDefinition(), G4ReactionProduct::SetKineticEnergy(), G4ReactionProduct::SetMomentum(), theTotalCount, G4Triton::Triton(), twopi, and Z.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply().

◆ SetProjectileRP()

void G4VParticleHPEnergyAngular::SetProjectileRP ( G4ReactionProduct * aIncidentParticleRP )

inlineinherited

Definition at line 74 of file G4VParticleHPEnergyAngular.hh.

  { 
    fCache.Get().theProjectileRP = aIncidentParticleRP;
  }

References G4VParticleHPEnergyAngular::fCache.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPProduct::Sample().

◆ SetQValue()

void G4VParticleHPEnergyAngular::SetQValue ( G4double aValue )

inlineinherited

Definition at line 101 of file G4VParticleHPEnergyAngular.hh.

  {
    theQValue = aValue;
  }

References G4VParticleHPEnergyAngular::theQValue.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPProduct::Init().

◆ SetTarget()

void G4VParticleHPEnergyAngular::SetTarget ( G4ReactionProduct * aTarget )

inlineinherited

Definition at line 79 of file G4VParticleHPEnergyAngular.hh.

  { 
    fCache.Get().theTarget = aTarget; 
  }

References G4VParticleHPEnergyAngular::fCache.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPProduct::Sample().

Field Documentation

◆ fCache

G4Cache<toBeCached> G4VParticleHPEnergyAngular::fCache

privateinherited

Definition at line 118 of file G4VParticleHPEnergyAngular.hh.

Referenced by G4VParticleHPEnergyAngular::G4VParticleHPEnergyAngular(), G4VParticleHPEnergyAngular::GetCMS(), G4VParticleHPEnergyAngular::GetProjectileRP(), G4VParticleHPEnergyAngular::GetTarget(), G4VParticleHPEnergyAngular::SetProjectileRP(), and G4VParticleHPEnergyAngular::SetTarget().

◆ theQValue

G4double G4VParticleHPEnergyAngular::theQValue

privateinherited

Definition at line 116 of file G4VParticleHPEnergyAngular.hh.

Referenced by G4VParticleHPEnergyAngular::G4VParticleHPEnergyAngular(), G4VParticleHPEnergyAngular::GetQValue(), and G4VParticleHPEnergyAngular::SetQValue().

◆ theTotalCount

G4int G4ParticleHPNBodyPhaseSpace::theTotalCount

private

Definition at line 102 of file G4ParticleHPNBodyPhaseSpace.hh.

Referenced by C(), G4ParticleHPNBodyPhaseSpace(), Init(), and Sample().

◆ theTotalMass

G4double G4ParticleHPNBodyPhaseSpace::theTotalMass

private

Definition at line 101 of file G4ParticleHPNBodyPhaseSpace.hh.

Referenced by G4ParticleHPNBodyPhaseSpace(), GetEmax(), and Init().

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

source/processes/hadronic/models/particle_hp/include/G4ParticleHPNBodyPhaseSpace.hh
source/processes/hadronic/models/particle_hp/src/G4ParticleHPNBodyPhaseSpace.cc

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4ParticleHPNBodyPhaseSpace()

◆ ~G4ParticleHPNBodyPhaseSpace()

Member Function Documentation

◆ C()

◆ ClearHistories()

◆ GetCMS()

◆ GetEmax()

◆ GetProjectileRP()

◆ GetQValue()

◆ GetTarget()

◆ Init() [1/2]

◆ Init() [2/2]

◆ MeanEnergyOfThisInteraction()

◆ Prob()

◆ Sample()

◆ SetProjectileRP()

◆ SetQValue()

◆ SetTarget()

Field Documentation

◆ fCache

◆ theQValue

◆ theTotalCount

◆ theTotalMass