G4NeutronHPContEnergyAngular Class Reference

#include <G4NeutronHPContEnergyAngular.hh>

Inheritance diagram for G4NeutronHPContEnergyAngular:

Public Member Functions
	G4NeutronHPContEnergyAngular ()
	~G4NeutronHPContEnergyAngular ()
void	Init (std::ifstream &aDataFile)
G4double	MeanEnergyOfThisInteraction ()
G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)
void	ClearHistories ()
Public Member Functions inherited from G4VNeutronHPEnergyAngular
	G4VNeutronHPEnergyAngular ()
virtual	~G4VNeutronHPEnergyAngular ()
virtual void	Init (std::ifstream &aDataFile)=0
virtual G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)=0
virtual G4double	MeanEnergyOfThisInteraction ()=0
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 43 of file G4NeutronHPContEnergyAngular.hh.

Constructor & Destructor Documentation

G4NeutronHPContEnergyAngular()

G4NeutronHPContEnergyAngular::G4NeutronHPContEnergyAngular ( )

inline

Definition at line 47 of file G4NeutronHPContEnergyAngular.hh.

  {
    theAngular = 0;
    currentMeanEnergy = -2;
  }

~G4NeutronHPContEnergyAngular()

G4NeutronHPContEnergyAngular::~G4NeutronHPContEnergyAngular ( )

inline

Definition at line 53 of file G4NeutronHPContEnergyAngular.hh.

  {
    if(theAngular!=0) delete [] theAngular;
  }

Member Function Documentation

ClearHistories()

void G4NeutronHPContEnergyAngular::ClearHistories ( )

virtual

Reimplemented from G4VNeutronHPEnergyAngular.

Definition at line 94 of file G4NeutronHPContEnergyAngular.cc.

{ 
   if ( theAngular!= NULL )
   { 
      for ( G4int i = 0 ; i< nEnergy ; i++ ) 
         theAngular[i].ClearHistories(); 
   } 
}

Referenced by ClearHistories().

Init()

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

inlinevirtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 60 of file G4NeutronHPContEnergyAngular.hh.

  {
    aDataFile >> theTargetCode >> theAngularRep >> theInterpolation >> nEnergy;
    theAngular = new G4NeutronHPContAngularPar[nEnergy];
    theManager.Init(aDataFile);
    for(G4int i=0; i<nEnergy; i++)
    {
      theAngular[i].Init(aDataFile);
      theAngular[i].SetInterpolation(theInterpolation);
    }
  }

MeanEnergyOfThisInteraction()

G4double G4NeutronHPContEnergyAngular::MeanEnergyOfThisInteraction ( )

virtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 76 of file G4NeutronHPContEnergyAngular.cc.

{
   G4double result(0);
   if(currentMeanEnergy<-1)
   {
     throw G4HadronicException(__FILE__, __LINE__, "G4NeutronHPContEnergyAngular: Logical error in Product class");
   }
   else
   {
     result = currentMeanEnergy;
   }
   currentMeanEnergy = -2;
   return result;
}

Sample()

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

virtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 35 of file G4NeutronHPContEnergyAngular.cc.

{
   G4ReactionProduct * result;
   G4int i(0);
   G4int it(0);
   for(i=0;i<nEnergy;i++)
   {
     it = i;
     if(theAngular[i].GetEnergy()>anEnergy) break;
   }
   G4double targetMass = GetTarget()->GetMass();
   if(it==0)
   {
     theAngular[0].SetTarget(GetTarget());
     theAngular[0].SetTargetCode(theTargetCode);
     theAngular[0].SetPrimary(GetNeutron());
     result = theAngular[0].Sample(anEnergy, massCode, targetMass, 
                                  theAngularRep, theInterpolation);
     currentMeanEnergy = theAngular[0].MeanEnergyOfThisInteraction();
   }
   else
   {
     // interpolation through alternating sampling. This needs improvement @@@
     // This is the cause of the He3 problem !!!!!!!!
     // See to it, if you can improve this.
     //080714 TK commnet Randomizing use angular distribution
     //080714 TK Always use the upper side distribution. enabling ClearHistories method.
     //G4double random = G4UniformRand();
     //G4double deltaE = theAngular[it].GetEnergy()-theAngular[it-1].GetEnergy();
     //G4double offset = theAngular[it].GetEnergy()-anEnergy;
     //if(random<offset/deltaE) it--;
     theAngular[it].SetTarget(GetTarget());
     theAngular[it].SetTargetCode(theTargetCode);
     theAngular[it].SetPrimary(GetNeutron());
     result = theAngular[it].Sample(anEnergy, massCode, targetMass, 
                                    theAngularRep, theInterpolation);
     currentMeanEnergy = theAngular[it].MeanEnergyOfThisInteraction();
   }
   return result;
}

---

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

• G4NeutronHPContEnergyAngular.hh
• G4NeutronHPContEnergyAngular.cc