G4ParticleHPIsotropic Class Reference

#include <G4ParticleHPIsotropic.hh>

Inheritance diagram for G4ParticleHPIsotropic:

Collaboration diagram for G4ParticleHPIsotropic:

Public Member Functions
	G4ParticleHPIsotropic ()
	~G4ParticleHPIsotropic ()
void	Init (std::istream &aDataFile)
G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)
G4double	MeanEnergyOfThisInteraction ()
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 38 of file G4ParticleHPIsotropic.hh.

Constructor & Destructor Documentation

G4ParticleHPIsotropic::G4ParticleHPIsotropic ( )

inline

Definition at line 42 of file G4ParticleHPIsotropic.hh.

{}

G4ParticleHPIsotropic::~G4ParticleHPIsotropic ( )

inline

Definition at line 44 of file G4ParticleHPIsotropic.hh.

{}

Member Function Documentation

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

virtual

Implements G4VParticleHPEnergyAngular.

Definition at line 48 of file G4ParticleHPIsotropic.cc.

{
}

G4double G4ParticleHPIsotropic::MeanEnergyOfThisInteraction ( )

inlinevirtual

Implements G4VParticleHPEnergyAngular.

Definition at line 50 of file G4ParticleHPIsotropic.hh.

  {
    return -1.;
  }

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

virtual

Implements G4VParticleHPEnergyAngular.

Definition at line 52 of file G4ParticleHPIsotropic.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());
     //110607 TK modified following parts for migration to G4NDL3.15 (ENDF VII.r0)
     //if(Z!=2) throw G4HadronicException(__FILE__, __LINE__, "Unknown ion case 1");    
     if(Z!=2) 
     {
        result->SetDefinition( G4IonTable::GetIonTable()->GetIon ( Z ,  A , 0.0 ) );
     }
   }
   else
   {
     //110607 TK modified following parts for migration to G4NDL3.15 (ENDF VII.r0)
     result->SetDefinition( G4IonTable::GetIonTable()->GetIon ( Z ,  A , 0.0 ) );
     //throw G4HadronicException(__FILE__, __LINE__, "G4ParticleHPIsotropic: Unknown ion case 2");
   }

   G4double cosTh = G4UniformRand();
   G4double phi = twopi*G4UniformRand();
   G4double theta = std::acos(cosTh);
   G4double sinth = std::sin(theta);
   
// we need the the Q value of the reaction
   result->SetKineticEnergy(std::max(0.001*MeV, anEnergy+GetQValue()));
   G4double mtot = result->GetTotalMomentum(); 
   G4ThreeVector tempVector(mtot*sinth*std::cos(phi), mtot*sinth*std::sin(phi), mtot*std::cos(theta) );
   result->SetMomentum(tempVector);

   return result;
}

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/processes/hadronic/models/particle_hp/include/G4ParticleHPIsotropic.hh
• geant4.10.03.p01/source/processes/hadronic/models/particle_hp/src/G4ParticleHPIsotropic.cc