G4NeutronHPIsotropic Class Reference

#include <G4NeutronHPIsotropic.hh>

Inheritance diagram for G4NeutronHPIsotropic:

Public Member Functions
	G4NeutronHPIsotropic ()
	~G4NeutronHPIsotropic ()
void	Init (std::ifstream &aDataFile)
G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)
G4double	MeanEnergyOfThisInteraction ()

---

Detailed Description

Definition at line 38 of file G4NeutronHPIsotropic.hh.

---

Constructor & Destructor Documentation

G4NeutronHPIsotropic::G4NeutronHPIsotropic

(

)

[inline]

Definition at line 42 of file G4NeutronHPIsotropic.hh.

{}

G4NeutronHPIsotropic::~G4NeutronHPIsotropic

(

)

[inline]

Definition at line 44 of file G4NeutronHPIsotropic.hh.

{}

---

Member Function Documentation

void G4NeutronHPIsotropic::Init

(

std::ifstream &

aDataFile

)

[virtual]

Implements G4VNeutronHPEnergyAngular.

Definition at line 46 of file G4NeutronHPIsotropic.cc.

{
}

G4double G4NeutronHPIsotropic::MeanEnergyOfThisInteraction

(

)

[inline, virtual]

Implements G4VNeutronHPEnergyAngular.

Definition at line 50 of file G4NeutronHPIsotropic.hh.

  {
    return -1.;
  }

G4ReactionProduct * G4NeutronHPIsotropic::Sample	(	G4double	anEnergy,
		G4double	massCode,
		G4double	mass
	)			[virtual]

Implements G4VNeutronHPEnergyAngular.

Definition at line 50 of file G4NeutronHPIsotropic.cc.

References G4Alpha::Alpha(), G4Deuteron::Deuteron(), G4Electron::Electron(), G4UniformRand, G4Gamma::Gamma(), G4ParticleTable::GetParticleTable(), G4VNeutronHPEnergyAngular::GetQValue(), G4ReactionProduct::GetTotalMomentum(), G4He3::He3(), G4Neutron::Neutron(), G4Positron::Positron(), G4Proton::Proton(), G4ReactionProduct::SetDefinition(), G4ReactionProduct::SetKineticEnergy(), G4ReactionProduct::SetMomentum(), and G4Triton::Triton().

{
   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( G4ParticleTable::GetParticleTable()->GetIon ( Z ,  A , 0.0 ) );
     }
   }
   else
   {
     //110607 TK modified following parts for migration to G4NDL3.15 (ENDF VII.r0)
     result->SetDefinition( G4ParticleTable::GetParticleTable()->GetIon ( Z ,  A , 0.0 ) );
     //throw G4HadronicException(__FILE__, __LINE__, "G4NeutronHPIsotropic: 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;
}

---

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

• G4NeutronHPIsotropic.hh
• G4NeutronHPIsotropic.cc

---