G4ParticleHPEnergyDistribution Class Reference

#include <G4ParticleHPEnergyDistribution.hh>

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

Detailed Description

Definition at line 47 of file G4ParticleHPEnergyDistribution.hh.

Constructor & Destructor Documentation

G4ParticleHPEnergyDistribution()

G4ParticleHPEnergyDistribution::G4ParticleHPEnergyDistribution ( )

inline

Definition at line 50 of file G4ParticleHPEnergyDistribution.hh.

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

~G4ParticleHPEnergyDistribution()

G4ParticleHPEnergyDistribution::~G4ParticleHPEnergyDistribution ( )

inline

Definition at line 56 of file G4ParticleHPEnergyDistribution.hh.

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

Member Function Documentation

Init()

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

inline

Definition at line 68 of file G4ParticleHPEnergyDistribution.hh.

  {
    G4double dummy;
    theData >> dummy >> theNumberOfPartials;
    theEnergyDistribution = new G4VParticleHPEDis * [theNumberOfPartials];
    for(G4int i=0; i<theNumberOfPartials; i++) 
    {
      theData >> theRepresentationType;
      switch(theRepresentationType)
      {
    case 1:
          theEnergyDistribution[i] = new G4ParticleHPArbitaryTab;
          break;
    case 5:        
          theEnergyDistribution[i] = new G4ParticleHPEvapSpectrum;
          break;
    case 7:
          theEnergyDistribution[i] = new G4ParticleHPFissionSpectrum;
          break;
    case 9:
          theEnergyDistribution[i] = new G4ParticleHPSimpleEvapSpectrum;
          break;
    case 11:
          theEnergyDistribution[i] = new G4ParticleHPWattSpectrum;
          break;
    case 12:
          theEnergyDistribution[i] = new G4ParticleHPMadlandNixSpectrum;
          break;
      }
      theEnergyDistribution[i]->Init(theData);
    }
  }

Referenced by G4FissionLibrary::Init(), G4ParticleHPFSFissionFS::Init(), G4ParticleHPInelasticCompFS::Init(), G4ParticleHPFissionBaseFS::Init(), and G4ParticleHPInelasticBaseFS::Init().

Sample()

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

inline

Definition at line 101 of file G4ParticleHPEnergyDistribution.hh.

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

Referenced by G4ParticleHPFissionBaseFS::ApplyYourself(), G4ParticleHPFSFissionFS::ApplyYourself(), G4ParticleHPInelasticBaseFS::BaseApply(), and G4ParticleHPInelasticCompFS::CompositeApply().

---

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

• G4ParticleHPEnergyDistribution.hh