G4NeutronHPProduct Class Reference

#include <G4NeutronHPProduct.hh>

Public Member Functions
	G4NeutronHPProduct ()
	~G4NeutronHPProduct ()
G4double	GetMassCode ()
G4double	GetMass ()
void	Init (std::ifstream &aDataFile)
G4ReactionProductVector *	Sample (G4double anEnergy)
G4double	GetMeanYield (G4double anEnergy)
void	SetNeutron (G4ReactionProduct *aNeutron)
void	SetTarget (G4ReactionProduct *aTarget)
G4ReactionProduct *	GetTarget ()
G4ReactionProduct *	GetNeutron ()
G4double	MeanEnergyOfThisInteraction ()
G4double	GetQValue ()

Detailed Description

Definition at line 49 of file G4NeutronHPProduct.hh.

Constructor & Destructor Documentation

G4NeutronHPProduct()

G4NeutronHPProduct::G4NeutronHPProduct ( )

inline

Definition at line 52 of file G4NeutronHPProduct.hh.

  {
    theDist = 0;
  }

~G4NeutronHPProduct()

G4NeutronHPProduct::~G4NeutronHPProduct ( )

inline

Definition at line 56 of file G4NeutronHPProduct.hh.

  {
    if(theDist != 0) delete theDist;
  }

Member Function Documentation

GetMass()

G4double G4NeutronHPProduct::GetMass ( )

inline

Definition at line 65 of file G4NeutronHPProduct.hh.

{return theMass;};

GetMassCode()

G4double G4NeutronHPProduct::GetMassCode ( )

inline

Definition at line 64 of file G4NeutronHPProduct.hh.

{return theMassCode;};

GetMeanYield()

G4double G4NeutronHPProduct::GetMeanYield ( G4double  anEnergy )

inline

Definition at line 132 of file G4NeutronHPProduct.hh.

  {
    return theYield.GetY(anEnergy);
  }

GetNeutron()

G4ReactionProduct * G4NeutronHPProduct::GetNeutron ( )

inline

Definition at line 149 of file G4NeutronHPProduct.hh.

{ return theNeutron; }

GetQValue()

G4double G4NeutronHPProduct::GetQValue ( )

inline

Definition at line 166 of file G4NeutronHPProduct.hh.

{ return theActualStateQValue; }

Referenced by G4NeutronHPEnAngCorrelation::Sample().

GetTarget()

G4ReactionProduct * G4NeutronHPProduct::GetTarget ( )

inline

Definition at line 147 of file G4NeutronHPProduct.hh.

{ return theTarget; }

Init()

void G4NeutronHPProduct::Init ( std::ifstream &  aDataFile )

inline

Definition at line 67 of file G4NeutronHPProduct.hh.

  {
    aDataFile >> theMassCode>>theMass>>theIsomerFlag>>theDistLaw
              >> theGroundStateQValue>>theActualStateQValue;
    theGroundStateQValue*= CLHEP::eV;
    theActualStateQValue*= CLHEP::eV;
    theYield.Init(aDataFile, CLHEP::eV);
    if(theDistLaw==0)
    {
      // distribution not known, use E-independent, isotropic angular distribution
      theDist = new G4NeutronHPIsotropic;
    }
    else if(theDistLaw == 1)
    {
      // Continuum energy-angular distribution
      theDist = new G4NeutronHPContEnergyAngular;
    }
    else if(theDistLaw == 2)
    {
      // Discrete 2-body scattering
      theDist = new G4NeutronHPDiscreteTwoBody;
    }
    else if(theDistLaw == 3)
    {
      // Isotropic emission
      theDist = new G4NeutronHPIsotropic;
    }
    else if(theDistLaw == 4)
    {
      // Discrete 2-body recoil modification
      // not used for now. @@@@
      theDist = new G4NeutronHPDiscreteTwoBody; 
      // the above is only temporary;
      // recoils need to be addressed
      // properly
      delete theDist;
      theDist = 0;
    }
    else if(theDistLaw == 5)
    {
      // charged particles only, to be used in a later stage. @@@@
    }
    else if(theDistLaw == 6)
    {
      // N-Body phase space
      theDist = new G4NeutronHPNBodyPhaseSpace;
    }
    else if(theDistLaw == 7)
    {
      // Laboratory angular energy paraetrisation
      theDist = new G4NeutronHPLabAngularEnergy;
    }
    else
    {
      throw G4HadronicException(__FILE__, __LINE__, "distribution law unknown to G4NeutronHPProduct");
    }
    if(theDist!=0)
    {
      theDist->SetQValue(theActualStateQValue);      
      theDist->Init(aDataFile);
    }
  }

Referenced by G4NeutronHPEnAngCorrelation::Init().

MeanEnergyOfThisInteraction()

G4double G4NeutronHPProduct::MeanEnergyOfThisInteraction ( )

inline

Definition at line 151 of file G4NeutronHPProduct.hh.

  { 
    G4double result;
    if(theDist == 0)
    {
      result = 0;
    }
    else
    {
      result=theDist->MeanEnergyOfThisInteraction();
      result *= theCurrentMultiplicity;
    }
    return result;
  }

Referenced by G4NeutronHPEnAngCorrelation::Sample().

Sample()

G4ReactionProductVector * G4NeutronHPProduct::Sample

(

G4double

anEnergy

)

Definition at line 42 of file G4NeutronHPProduct.cc.

{
  if(theDist == 0) { return 0; }
  G4ReactionProductVector * result = new G4ReactionProductVector;
  G4double mean = theYield.GetY(anEnergy);
  G4int multi;
  multi = G4int(mean+0.0001);
  //if(theMassCode==0) multi = G4Poisson(mean); // @@@@gammas. please X-check this
  //080718
  if ( theMassCode == 0 )
  { 
     if ( G4int ( mean ) == mean )
     {
        multi = (G4int) mean;
     }
     else
     {
        multi = G4Poisson ( mean ); 
     }
  }
  theDist->SetTarget(theTarget);
  theDist->SetNeutron(theNeutron);
  G4int i;
//  G4double eMax = GetTarget()->GetMass()+GetNeutron()->GetMass()
//                  - theActualStateQValue;
  theCurrentMultiplicity = static_cast<G4int>(mean);
  G4ReactionProduct * tmp;
  theDist->ClearHistories();
  for(i=0;i<multi;i++)
  {
    tmp = theDist->Sample(anEnergy, theMassCode, theMass);
    if(tmp != 0) { result->push_back(tmp); }
  }
  if(multi == 0) 
  {
    tmp = theDist->Sample(anEnergy, theMassCode, theMass);
    delete  tmp;
  }
/*
//080901 TK Comment out, too many secondaries are produced in deuteron reactions
  if(theTarget->GetMass()<2*GeV) // @@@ take care of residuals in all cases
  {
    tmp = theDist->Sample(anEnergy, theMassCode, theMass);
    tmp->SetDefinition(G4Proton::Proton());
    if(tmp != 0) { result->push_back(tmp); }
  }
*/
  return result;
}

Referenced by G4NeutronHPEnAngCorrelation::Sample(), and G4NeutronHPEnAngCorrelation::SampleOne().

SetNeutron()

void G4NeutronHPProduct::SetNeutron ( G4ReactionProduct *  aNeutron )

inline

Definition at line 137 of file G4NeutronHPProduct.hh.

  { 
    theNeutron = aNeutron; 
  }

SetTarget()

void G4NeutronHPProduct::SetTarget ( G4ReactionProduct *  aTarget )

inline

Definition at line 142 of file G4NeutronHPProduct.hh.

  { 
    theTarget = aTarget; 
  }

---

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

• G4NeutronHPProduct.hh
• G4NeutronHPProduct.cc