G4NeutronHPEnergyDistribution Class Reference

#include <G4NeutronHPEnergyDistribution.hh>

Public Member Functions
	G4NeutronHPEnergyDistribution ()
	~G4NeutronHPEnergyDistribution ()
void	Init (std::istream &theData)
G4double	Sample (G4double anEnergy, G4int &it)

Detailed Description

Definition at line 46 of file G4NeutronHPEnergyDistribution.hh.

Constructor & Destructor Documentation

G4NeutronHPEnergyDistribution::G4NeutronHPEnergyDistribution ( )

inline

Definition at line 49 of file G4NeutronHPEnergyDistribution.hh.

  {
    theEnergyDistribution = 0;
    theNumberOfPartials = 0;
    theRepresentationType = 0;
  }

G4NeutronHPEnergyDistribution::~G4NeutronHPEnergyDistribution ( )

inline

Definition at line 55 of file G4NeutronHPEnergyDistribution.hh.

  {
    if(theEnergyDistribution != 0)
    {
      for(G4int i=0; i<theNumberOfPartials; i++) 
      {
        delete theEnergyDistribution[i];
      }
      delete [] theEnergyDistribution;
    }
  }

Member Function Documentation

void G4NeutronHPEnergyDistribution::Init ( std::istream &  theData )

inline

Definition at line 67 of file G4NeutronHPEnergyDistribution.hh.

References G4VNeutronHPEDis::Init().

Referenced by G4NeutronHPFSFissionFS::Init(), G4NeutronHPFissionBaseFS::Init(), G4NeutronHPInelasticBaseFS::Init(), G4NeutronHPInelasticCompFS::Init(), and G4FissionLibrary::Init().

  {
    G4double dummy;
    theData >> dummy >> theNumberOfPartials;
    theEnergyDistribution = new G4VNeutronHPEDis * [theNumberOfPartials];
    for(G4int i=0; i<theNumberOfPartials; i++) 
    {
      theData >> theRepresentationType;
      switch(theRepresentationType)
      {
    case 1:
          theEnergyDistribution[i] = new G4NeutronHPArbitaryTab;
          break;
    case 5:        
          theEnergyDistribution[i] = new G4NeutronHPEvapSpectrum;
          break;
    case 7:
          theEnergyDistribution[i] = new G4NeutronHPFissionSpectrum;
          break;
    case 9:
          theEnergyDistribution[i] = new G4NeutronHPSimpleEvapSpectrum;
          break;
    case 11:
          theEnergyDistribution[i] = new G4NeutronHPWattSpectrum;
          break;
    case 12:
          theEnergyDistribution[i] = new G4NeutronHPMadlandNixSpectrum;
          break;
      }
      theEnergyDistribution[i]->Init(theData);
    }
  }

G4double G4NeutronHPEnergyDistribution::Sample ( G4double  anEnergy, G4int &  it  )

inline

Definition at line 100 of file G4NeutronHPEnergyDistribution.hh.

References G4UniformRand, G4VNeutronHPEDis::GetFractionalProbability(), and G4VNeutronHPEDis::Sample().

Referenced by G4NeutronHPFSFissionFS::ApplyYourself(), G4NeutronHPFissionBaseFS::ApplyYourself(), G4NeutronHPInelasticBaseFS::BaseApply(), and G4NeutronHPInelasticCompFS::CompositeApply().

  {
    G4double result = 0;
    it = 0;
    if (theNumberOfPartials != 0)
    {
      G4double sum=0;
      G4double * running = new G4double[theNumberOfPartials];
      running[0] = 0;
      G4int i;
      for (i=0; i<theNumberOfPartials; i++)
      { 
    if (i!=0) running[i]=running[i-1];
    running[i]+=theEnergyDistribution[i]->GetFractionalProbability(anEnergy);
      }
      sum = running[theNumberOfPartials-1];
      G4double random = G4UniformRand();
      for(i=0; i<theNumberOfPartials; i++)
      {
    it = i;
    if(running[i]/sum>random) break;
      }
      delete [] running;
      if(it==theNumberOfPartials) it--;
      result = theEnergyDistribution[it]->Sample(anEnergy);
    }
    return result;
  }

---

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

• G4NeutronHPEnergyDistribution.hh