G4NeutronHPIsotropic Class Reference

#include <G4NeutronHPIsotropic.hh>

Inheritance diagram for G4NeutronHPIsotropic:

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

Constructor & Destructor Documentation

G4NeutronHPIsotropic::G4NeutronHPIsotropic ( )

inline

Definition at line 41 of file G4NeutronHPIsotropic.hh.

{}

G4NeutronHPIsotropic::~G4NeutronHPIsotropic ( )

inline

Definition at line 43 of file G4NeutronHPIsotropic.hh.

{}

Member Function Documentation

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

virtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 46 of file G4NeutronHPIsotropic.cc.

{
}

G4double G4NeutronHPIsotropic::MeanEnergyOfThisInteraction ( )

inlinevirtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 49 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(), G4IonTable::GetIonTable(), G4VNeutronHPEnergyAngular::GetQValue(), G4ReactionProduct::GetTotalMomentum(), G4He3::He3(), G4INCL::Math::max(), python.hepunit::MeV, G4Neutron::Neutron(), G4Positron::Positron(), G4Proton::Proton(), G4ReactionProduct::SetDefinition(), G4ReactionProduct::SetKineticEnergy(), G4ReactionProduct::SetMomentum(), G4Triton::Triton(), and python.hepunit::twopi.

{
   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__, "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