G4NeutronHPContEnergyAngular Class Reference

#include <G4NeutronHPContEnergyAngular.hh>

Inheritance diagram for G4NeutronHPContEnergyAngular:

Public Member Functions
	G4NeutronHPContEnergyAngular ()
	~G4NeutronHPContEnergyAngular ()
void	Init (std::istream &aDataFile)
G4double	MeanEnergyOfThisInteraction ()
G4ReactionProduct *	Sample (G4double anEnergy, G4double massCode, G4double mass)
void	ClearHistories ()
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)

Additional Inherited Members
Protected Member Functions inherited from G4VNeutronHPEnergyAngular
G4double	GetQValue ()

Detailed Description

Definition at line 42 of file G4NeutronHPContEnergyAngular.hh.

Constructor & Destructor Documentation

G4NeutronHPContEnergyAngular::G4NeutronHPContEnergyAngular ( )

inline

Definition at line 46 of file G4NeutronHPContEnergyAngular.hh.

  {
    theAngular = 0;
    currentMeanEnergy = -2;
  }

G4NeutronHPContEnergyAngular::~G4NeutronHPContEnergyAngular ( )

inline

Definition at line 52 of file G4NeutronHPContEnergyAngular.hh.

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

Member Function Documentation

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(); 
   } 
}

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

inlinevirtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 59 of file G4NeutronHPContEnergyAngular.hh.

References G4NeutronHPContAngularPar::Init(), G4InterpolationManager::Init(), and G4NeutronHPContAngularPar::SetInterpolation().

  {
    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);
    }
  }

G4double G4NeutronHPContEnergyAngular::MeanEnergyOfThisInteraction ( )

virtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 77 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;
}

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

virtual

Implements G4VNeutronHPEnergyAngular.

Definition at line 35 of file G4NeutronHPContEnergyAngular.cc.

References G4ReactionProduct::GetMass(), G4VNeutronHPEnergyAngular::GetNeutron(), G4VNeutronHPEnergyAngular::GetTarget(), G4NeutronHPContAngularPar::MeanEnergyOfThisInteraction(), G4NeutronHPContAngularPar::Sample(), G4NeutronHPContAngularPar::SetPrimary(), G4NeutronHPContAngularPar::SetTarget(), and G4NeutronHPContAngularPar::SetTargetCode().

{
   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