G4INCL::PhaseSpaceGenerator Namespace Reference

Functions
void	generate (const G4double sqrtS, ParticleList &particles)
	Generate an event in the CM system.
void	generateBiased (const G4double sqrtS, ParticleList &particles, const size_t index, const G4double slope)
	Generate a biased event in the CM system.
void	setPhaseSpaceGenerator (IPhaseSpaceGenerator *g)
IPhaseSpaceGenerator *	getPhaseSpaceGenerator ()
void	deletePhaseSpaceGenerator ()
void	initialize (Config const *const theConfig)

Function Documentation

deletePhaseSpaceGenerator()

void G4INCL::PhaseSpaceGenerator::deletePhaseSpaceGenerator

(

)

Definition at line 117 of file G4INCLPhaseSpaceGenerator.cc.

                                     {
      delete thePhaseSpaceGenerator;
      thePhaseSpaceGenerator = NULL;
    }

Referenced by G4INCL::INCL::~INCL().

generate()

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

Generate an event in the CM system.

Definition at line 94 of file G4INCLPhaseSpaceGenerator.cc.

                                                                 {
      return thePhaseSpaceGenerator->generate(sqrtS, particles);
    }

Referenced by G4INCL::NNbarToAnnihilationChannel::fillFinalState(), G4INCL::NNbarToLLbarChannel::fillFinalState(), G4INCL::NNbarToNNbar2piChannel::fillFinalState(), G4INCL::NNbarToNNbar3piChannel::fillFinalState(), G4INCL::NNbarToNNbarpiChannel::fillFinalState(), G4INCL::PionResonanceDecayChannel::fillFinalState(), generateBiased(), G4INCL::PbarAtrestEntryChannel::makeMesonStar(), and G4INCL::INCL::processEvent().

generateBiased()

void G4INCL::PhaseSpaceGenerator::generateBiased	(	const G4double	sqrtS,
		ParticleList &	particles,
		const size_t	index,
		const G4double	slope )

Generate a biased event in the CM system.

This method first generates a "flat" event by calling generate(). The particles are subsequently rotated in such a way that one of them (identified by the parameter index) is biased towards the collisionAxis with an exponential distribution of the form

where is the usual Mandelstam variable. The incoming momentum is taken to be the momentum of particles[index] at the moment of the call.

Parameters

sqrtS	total energy in the centre of mass, in MeV
particles	list of particles for which the event will be generated (modified on exit)
index	index of the particle to be biased; all the other particles will follow
slope	slope of the angular distribution: , in (GeV/c)^(-2)

Definition at line 98 of file G4INCLPhaseSpaceGenerator.cc.

                                                                                                                 {
// assert(index<particles.size());
      // store the incoming momentum of particle[index]; it will be used to
      // compute t when biasing
      biasMe = particles[index];
      const ThreeVector pInVec = biasMe->getMomentum();
      generate(sqrtS, particles);
      // Extremely rare event try to bias with vector null
      if(pInVec.mag() != 0.) bias(particles, pInVec, slope);
    }

Referenced by G4INCL::EtaNToPiPiNChannel::fillFinalState(), G4INCL::NDeltaEtaProductionChannel::fillFinalState(), G4INCL::NDeltaOmegaProductionChannel::fillFinalState(), G4INCL::NDeltaToDeltaLKChannel::fillFinalState(), G4INCL::NDeltaToDeltaSKChannel::fillFinalState(), G4INCL::NDeltaToNLKChannel::fillFinalState(), G4INCL::NDeltaToNNKKbChannel::fillFinalState(), G4INCL::NDeltaToNSKChannel::fillFinalState(), G4INCL::NKbToL2piChannel::fillFinalState(), G4INCL::NKbToNKb2piChannel::fillFinalState(), G4INCL::NKbToNKbpiChannel::fillFinalState(), G4INCL::NKbToS2piChannel::fillFinalState(), G4INCL::NKToNK2piChannel::fillFinalState(), G4INCL::NKToNKpiChannel::fillFinalState(), G4INCL::NNEtaToMultiPionsChannel::fillFinalState(), G4INCL::NNOmegaToMultiPionsChannel::fillFinalState(), G4INCL::NNToMissingStrangenessChannel::fillFinalState(), G4INCL::NNToMultiPionsChannel::fillFinalState(), G4INCL::NNToNLK2piChannel::fillFinalState(), G4INCL::NNToNLKChannel::fillFinalState(), G4INCL::NNToNLKpiChannel::fillFinalState(), G4INCL::NNToNNEtaChannel::fillFinalState(), G4INCL::NNToNNKKbChannel::fillFinalState(), G4INCL::NNToNNOmegaChannel::fillFinalState(), G4INCL::NNToNSK2piChannel::fillFinalState(), G4INCL::NNToNSKChannel::fillFinalState(), G4INCL::NNToNSKpiChannel::fillFinalState(), G4INCL::NpiToLK2piChannel::fillFinalState(), G4INCL::NpiToLKpiChannel::fillFinalState(), G4INCL::NpiToMissingStrangenessChannel::fillFinalState(), G4INCL::NpiToNKKbChannel::fillFinalState(), G4INCL::NpiToSK2piChannel::fillFinalState(), G4INCL::NpiToSKpiChannel::fillFinalState(), G4INCL::OmegaNToPiPiNChannel::fillFinalState(), and G4INCL::PiNToMultiPionsChannel::fillFinalState().

getPhaseSpaceGenerator()

IPhaseSpaceGenerator * G4INCL::PhaseSpaceGenerator::getPhaseSpaceGenerator

(

)

Definition at line 113 of file G4INCLPhaseSpaceGenerator.cc.

                                                   {
      return thePhaseSpaceGenerator;
    }

initialize()

void G4INCL::PhaseSpaceGenerator::initialize

(

Config const *const

theConfig

)

Definition at line 122 of file G4INCLPhaseSpaceGenerator.cc.

                                                    {
      PhaseSpaceGeneratorType psg = theConfig->getPhaseSpaceGeneratorType();
      if(psg==RauboldLynchType)
        setPhaseSpaceGenerator(new PhaseSpaceRauboldLynch);
      else if(psg==KopylovType)
        setPhaseSpaceGenerator(new PhaseSpaceKopylov);
      else
        setPhaseSpaceGenerator(NULL);
    }

Referenced by G4INCL::INCL::INCL().

setPhaseSpaceGenerator()

void G4INCL::PhaseSpaceGenerator::setPhaseSpaceGenerator

(

IPhaseSpaceGenerator *

g

)

Definition at line 109 of file G4INCLPhaseSpaceGenerator.cc.

                                                         {
      thePhaseSpaceGenerator = g;
    }

Referenced by initialize().