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 = nullptr;
      theNumberOfPartials = 0;
      theRepresentationType = 0;
    }

~G4ParticleHPEnergyDistribution()

G4ParticleHPEnergyDistribution::~G4ParticleHPEnergyDistribution ( )

inline

Definition at line 56 of file G4ParticleHPEnergyDistribution.hh.

    {
      if (theEnergyDistribution != nullptr) {
        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 66 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;
    default:
            theEnergyDistribution[i] = new G4ParticleHPArbitaryTab;
        }
        theEnergyDistribution[i]->Init(theData);
      }
    }

Sample()

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

inline

Definition at line 99 of file G4ParticleHPEnergyDistribution.hh.

    {
      G4double result = 0;
      it = 0;
      if (theNumberOfPartials != 0) {
        G4double sum = 0;
        auto 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:

• G4ParticleHPEnergyDistribution.hh