Generate momenta using the RauboldLynch method. More...

#include <G4INCLPhaseSpaceRauboldLynch.hh>

Inheritance diagram for G4INCL::PhaseSpaceRauboldLynch:

Public Member Functions
	PhaseSpaceRauboldLynch ()

virtual	~PhaseSpaceRauboldLynch ()

	PhaseSpaceRauboldLynch (PhaseSpaceRauboldLynch const &other)
	Dummy copy constructor to silence Coverity warning.

PhaseSpaceRauboldLynch &	operator= (PhaseSpaceRauboldLynch const &rhs)
	Dummy assignment operator to silence Coverity warning.

void	generate (const G4double sqrtS, ParticleList &particles)
	Generate momenta according to a uniform, Lorentz-invariant phase-space model.

G4double	getMaxGeneratedWeight () const
	Return the largest generated weight.

Public Member Functions inherited from G4INCL::IPhaseSpaceGenerator
	IPhaseSpaceGenerator ()

virtual	~IPhaseSpaceGenerator ()

virtual void	generate (const G4double sqrtS, ParticleList &particles)=0
	Generate an event in the CM frame.

Detailed Description

Generate momenta using the RauboldLynch method.

Definition at line 49 of file G4INCLPhaseSpaceRauboldLynch.hh.

Constructor & Destructor Documentation

◆ PhaseSpaceRauboldLynch() [1/2]

G4INCL::PhaseSpaceRauboldLynch::PhaseSpaceRauboldLynch ( )

Definition at line 180 of file G4INCLPhaseSpaceRauboldLynch.cc.

                                                 :
    nParticles(0),
    sqrtS(0.),
    availableEnergy(0.),
    maxGeneratedWeight(0.)
  {
    std::vector<G4double> wMaxMasslessXV(wMaxMasslessX, wMaxMasslessX + wMaxNE);
    std::vector<G4double> wMaxMasslessYV(wMaxMasslessY, wMaxMasslessY + wMaxNE);
    wMaxMassless = new InterpolationTable(wMaxMasslessXV, wMaxMasslessYV);
    std::vector<G4double> wMaxCorrectionXV(wMaxCorrectionX, wMaxCorrectionX + wMaxNE);
    std::vector<G4double> wMaxCorrectionYV(wMaxCorrectionY, wMaxCorrectionY + wMaxNE);
    wMaxCorrection = new InterpolationTable(wMaxCorrectionXV, wMaxCorrectionYV);
 
    // initialize the precalculated coefficients
    prelog[0] = 0.;
    for(size_t i=1; i<wMaxNP; ++i) {
      prelog[i] = -std::log(G4double(i));
    }
  }

◆ ~PhaseSpaceRauboldLynch()

G4INCL::PhaseSpaceRauboldLynch::~PhaseSpaceRauboldLynch ( )

virtual

Definition at line 200 of file G4INCLPhaseSpaceRauboldLynch.cc.

                                                  {
    delete wMaxMassless;
    delete wMaxCorrection;
  }

◆ PhaseSpaceRauboldLynch() [2/2]

G4INCL::PhaseSpaceRauboldLynch::PhaseSpaceRauboldLynch ( PhaseSpaceRauboldLynch const & other )

Dummy copy constructor to silence Coverity warning.

Member Function Documentation

◆ generate()

void G4INCL::PhaseSpaceRauboldLynch::generate	(	const G4double	sqrtS,
		ParticleList &	particles
	)

virtual

Generate momenta according to a uniform, Lorentz-invariant phase-space model.

This function will assign momenta to the particles in the list that is passed as an argument. The event is generated in the CM frame.

Parameters

sqrtS	total centre-of-mass energy of the system
particles	list of particles

Implements G4INCL::IPhaseSpaceGenerator.

Definition at line 205 of file G4INCLPhaseSpaceRauboldLynch.cc.

                                                                                   {
    maxGeneratedWeight = 0.;
 
    sqrtS = sqs;
 
    // initialize the structures containing the particle masses
    initialize(particles);
 
    // initialize the maximum weight
    const G4double weightMax = computeMaximumWeightParam();
 
    const G4int maxIter = 500;
    G4int iter = 0;
    G4double weight, r;
    do {
      weight = computeWeight();
      maxGeneratedWeight = std::max(weight, maxGeneratedWeight);
// assert(weight<=weightMax);
      r = Random::shoot();
    } while(++iter<maxIter && r*weightMax>weight); /* Loop checking, 10.07.2015, D.Mancusi */
 
#ifndef INCLXX_IN_GEANT4_MODE
    if(iter>=maxIter) {
      INCL_WARN("Number of tries exceeded in PhaseSpaceRauboldLynch::generate()\n"
                << "nParticles=" << nParticles << ", sqrtS=" << sqrtS << '\n');
    }
#endif
 
    generateEvent(particles);
  }

◆ getMaxGeneratedWeight()

G4double G4INCL::PhaseSpaceRauboldLynch::getMaxGeneratedWeight ( ) const

Return the largest generated weight.

Definition at line 345 of file G4INCLPhaseSpaceRauboldLynch.cc.

                                                               {
    return maxGeneratedWeight;
  }

◆ operator=()

PhaseSpaceRauboldLynch & G4INCL::PhaseSpaceRauboldLynch::operator= ( PhaseSpaceRauboldLynch const & rhs )

Dummy assignment operator to silence Coverity warning.

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

Public Member Functions