#include <G4ParticleHPEnergyDistribution.hh>

Public Member Functions
	G4ParticleHPEnergyDistribution ()

void	Init (std::istream &theData)

G4double	Sample (G4double anEnergy, G4int &it)

	~G4ParticleHPEnergyDistribution ()

Private Attributes
G4VParticleHPEDis **	theEnergyDistribution

G4int	theNumberOfPartials

G4int	theRepresentationType

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

References theEnergyDistribution, theNumberOfPartials, and theRepresentationType.

◆ ~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;
    }
  }

References theEnergyDistribution, and theNumberOfPartials.

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

References G4VParticleHPEDis::Init(), theEnergyDistribution, theNumberOfPartials, and theRepresentationType.

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