G4DistributionGenerator Class Reference

#include <G4DistributionGenerator.hh>

Public Member Functions
	G4DistributionGenerator (std::vector< G4double > &x, std::vector< G4double > &values)
	G4DistributionGenerator ()
	~G4DistributionGenerator ()
G4double	Generate (G4double ranflat)

Detailed Description

Definition at line 41 of file G4DistributionGenerator.hh.

Constructor & Destructor Documentation

G4DistributionGenerator() [1/2]

G4DistributionGenerator::G4DistributionGenerator	(	std::vector< G4double > &	x,
		std::vector< G4double > &	values
	)

Definition at line 42 of file G4DistributionGenerator.cc.

{
  G4HadronicDeprecate("G4DistributionGenerator");
  _x = x;
 
  // Check boundaries: must be size(x) = size(values) + 1
  if (x.size() != (values.size() + 1))
    { G4cout << " Inconsistent parameters in G4DistributionGenerator "
       << G4endl;
    }
  assert (x.size() == (values.size() + 1));
 
  G4double tot = 0.;
  unsigned int i;
  for (i=0; i<values.size(); i++) { tot += values[i]; }
  assert (tot > 0.);
  
  _cumProb.push_back(0.);    
  //  _cumProb.push_back(values[0] / tot);
  G4double sum = 0.;
  for (i=0; i<values.size(); i++) 
    { 
      sum += values[i];
      _cumProb.push_back(sum / tot); }
 
  // Debugging
  /*
  for (i=0; i<values.size(); i++)
    { G4cout << values[i] << "  " ; }
  G4cout << "  Integral = " << tot << G4endl;
  for (i=0; i<_cumProb.size(); i++)
    { 
      G4cout << "Variable " << _x[i]  
       << " --- cumProb = " << _cumProb[i] << G4endl;
    }
  */
  // End of debugging
 
}

G4DistributionGenerator() [2/2]

G4DistributionGenerator::G4DistributionGenerator

(

)

~G4DistributionGenerator()

G4DistributionGenerator::~G4DistributionGenerator

(

)

Definition at line 86 of file G4DistributionGenerator.cc.

{
}

Member Function Documentation

Generate()

G4double G4DistributionGenerator::Generate

(

G4double

ranflat

)

Definition at line 91 of file G4DistributionGenerator.cc.

{
  G4double xRandom = _x[0];
 
  G4int bin = _cumProb.size() - 1;
  unsigned int i;
  for (i=1; i<_cumProb.size(); i++)
    {
      if (ranflat >= _cumProb[i-1] && ranflat < _cumProb[i])
    {
      bin = i - 1;
    }
    }
 
  if (bin >= 0 && bin < static_cast<G4int>(_cumProb.size()-1) && bin < static_cast<G4int>(_x.size()-1)) 
    {
      G4double coeff = (ranflat - _cumProb[bin]) *  (_x[bin+1] - _x[bin]) / 
                   (_cumProb[bin+1] - _cumProb[bin]);
      xRandom = _x[bin] + coeff;
 
      // Deugging
      /*
      G4cout << "Random = " << ranflat << " Generated " << xRandom << G4endl;
      */
      // Endo of Debugging
 
    }
  else
    {   
      // Debugging
      /*
      G4cout << "Bin " << bin << " "
       << _cumProb.size() << " " 
       << _x.size()
       << G4endl;
      */
      // End of debugging
    }
 
  return xRandom;
 
}

Referenced by G4PiMinusStopMaterial::P4Vector().

---

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

• G4DistributionGenerator.hh
• G4DistributionGenerator.cc