G4INCL::ParticleSampler Class Reference

#include <G4INCLParticleSampler.hh>

Public Member Functions
NuclearDensity const *	getDensity () const
	Getter for theDensity. More...
NuclearPotential::INuclearPotential const *	getPotential () const
	Getter for thePotential. More...
G4double	getRPCorrelationCoefficient (const ParticleType t) const
	Getter for rpCorrelationCoefficient. More...
	ParticleSampler (const G4int A, const G4int Z, const G4int S)
	Constructor. More...
ParticleList	sampleParticles (ThreeVector const &position)
void	sampleParticlesIntoList (ThreeVector const &position, ParticleList &theList)
void	setDensity (NuclearDensity const *const d)
	Setter for theDensity. More...
void	setPotential (NuclearPotential::INuclearPotential const *const p)
	Setter for thePotential. More...
void	setRPCorrelationCoefficient (const ParticleType t, const G4double corrCoeff)
	Setter for rpCorrelationCoefficient. More...
	~ParticleSampler ()
	Destructor. More...

Private Types
typedef Particle *(ParticleSampler::*	ParticleSamplerMethod) (const ParticleType t) const

Private Member Functions
Particle *	sampleOneParticleWithFuzzyRPCorrelation (const ParticleType t) const
	Sample one particle with a fuzzy rp-correlation. More...
Particle *	sampleOneParticleWithoutRPCorrelation (const ParticleType t) const
	Sample one particle not taking into account the rp-correlation. More...
Particle *	sampleOneParticleWithRPCorrelation (const ParticleType t) const
	Sample one particle taking into account the rp-correlation. More...
void	updateSampleOneParticleMethods ()

Private Attributes
G4double	rpCorrelationCoefficient [UnknownParticle]
	Correlation coefficients for the r-p correlation. More...
ParticleSamplerMethod	sampleOneNeutron
	Sample a list of particles. More...
ParticleSamplerMethod	sampleOneProton
	Sample a list of particles. More...
const G4int	theA
	Mass number. More...
NuclearDensity const *	theDensity
	Pointer to the Cluster's NuclearDensity. More...
InterpolationTable const *	thePCDFTable [UnknownParticle]
	Array of pointers to the p-space CDF table. More...
NuclearPotential::INuclearPotential const *	thePotential
	Pointer to the Cluster's NuclearPotential. More...
InterpolationTable const *	theRCDFTable [UnknownParticle]
	Array of pointers to the r-space CDF table. More...
const G4int	theS
	Strangeness number. More...
const G4int	theZ
	Charge number. More...

Detailed Description

Definition at line 54 of file G4INCLParticleSampler.hh.

Member Typedef Documentation

ParticleSamplerMethod

typedef Particle *(ParticleSampler::* G4INCL::ParticleSampler::ParticleSamplerMethod) (const ParticleType t) const

private

Definition at line 98 of file G4INCLParticleSampler.hh.

Constructor & Destructor Documentation

ParticleSampler()

G4INCL::ParticleSampler::ParticleSampler	(	const G4int	A,
		const G4int	Z,
		const G4int	S
	)

Constructor.

Parameters

A	the mass number
Z	the charge number

Definition at line 51 of file G4INCLParticleSampler.cc.

                                                                              :
    sampleOneProton(&ParticleSampler::sampleOneParticleWithoutRPCorrelation),
    sampleOneNeutron(&ParticleSampler::sampleOneParticleWithoutRPCorrelation),
    theA(A),
    theZ(Z),
    theS(S),
    theDensity(NULL),
    thePotential(NULL)
  {
    std::fill(theRCDFTable, theRCDFTable + UnknownParticle, static_cast<InterpolationTable *>(NULL));
    std::fill(thePCDFTable, thePCDFTable + UnknownParticle, static_cast<InterpolationTable *>(NULL));
    std::fill(rpCorrelationCoefficient, rpCorrelationCoefficient + UnknownParticle, 1.);
    rpCorrelationCoefficient[Proton] = ParticleTable::getRPCorrelationCoefficient(Proton);
    rpCorrelationCoefficient[Neutron] = ParticleTable::getRPCorrelationCoefficient(Neutron);
    rpCorrelationCoefficient[Lambda] = ParticleTable::getRPCorrelationCoefficient(Lambda);
  }

References G4INCL::ParticleTable::getRPCorrelationCoefficient(), G4INCL::Lambda, G4INCL::Neutron, G4INCL::Proton, rpCorrelationCoefficient, thePCDFTable, theRCDFTable, and G4INCL::UnknownParticle.

~ParticleSampler()

G4INCL::ParticleSampler::~ParticleSampler

(

)

Destructor.

Definition at line 68 of file G4INCLParticleSampler.cc.

                                    {
  }

Member Function Documentation

getDensity()

NuclearDensity const * G4INCL::ParticleSampler::getDensity ( ) const

inline

Getter for theDensity.

Definition at line 68 of file G4INCLParticleSampler.hh.

{ return theDensity; }

References theDensity.

getPotential()

NuclearPotential::INuclearPotential const * G4INCL::ParticleSampler::getPotential ( ) const

inline

Getter for thePotential.

Definition at line 71 of file G4INCLParticleSampler.hh.

{ return thePotential; }

References thePotential.

getRPCorrelationCoefficient()

G4double G4INCL::ParticleSampler::getRPCorrelationCoefficient ( const ParticleType  t ) const

inline

Getter for rpCorrelationCoefficient.

Definition at line 74 of file G4INCLParticleSampler.hh.

                                                                       {
// assert(t==Proton || t==Neutron);
        return rpCorrelationCoefficient[t];
      }

References rpCorrelationCoefficient.

sampleOneParticleWithFuzzyRPCorrelation()

Particle * G4INCL::ParticleSampler::sampleOneParticleWithFuzzyRPCorrelation ( const ParticleType  t ) const

private

Sample one particle with a fuzzy rp-correlation.

Definition at line 165 of file G4INCLParticleSampler.cc.

                                                                                               {
// assert(theDensity && thePotential);
    std::pair<G4double,G4double> ranNumbers = Random::correlatedUniform(rpCorrelationCoefficient[t]);
    const G4double x = Math::pow13(ranNumbers.first);
    const G4double y = Math::pow13(ranNumbers.second);
    const G4double theFermiMomentum = thePotential->getFermiMomentum(t);
    const ThreeVector momentumVector = Random::normVector(y*theFermiMomentum);
    const G4double reflectionRadius = theDensity->getMaxRFromP(t, x);
    const ThreeVector positionVector = Random::sphereVector(reflectionRadius);
    Particle *aParticle = new Particle(t, momentumVector, positionVector);
    aParticle->setUncorrelatedMomentum(x*theFermiMomentum);
    return aParticle;
  }

References G4INCL::Random::correlatedUniform(), G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::NuclearDensity::getMaxRFromP(), G4INCL::Random::normVector(), G4INCL::Math::pow13(), rpCorrelationCoefficient, G4INCL::Particle::setUncorrelatedMomentum(), G4INCL::Random::sphereVector(), theDensity, and thePotential.

Referenced by updateSampleOneParticleMethods().

sampleOneParticleWithoutRPCorrelation()

Particle * G4INCL::ParticleSampler::sampleOneParticleWithoutRPCorrelation ( const ParticleType  t ) const

private

Sample one particle not taking into account the rp-correlation.

Definition at line 157 of file G4INCLParticleSampler.cc.

                                                                                             {
    const G4double position = (*(theRCDFTable[t]))(Random::shoot());
    const G4double momentum = (*(thePCDFTable[t]))(Random::shoot());
    ThreeVector positionVector = Random::normVector(position);
    ThreeVector momentumVector = Random::normVector(momentum);
    return new Particle(t, momentumVector, positionVector);
  }

References G4INCL::Random::normVector(), G4INCL::Random::shoot(), thePCDFTable, and theRCDFTable.

Referenced by sampleParticlesIntoList(), and updateSampleOneParticleMethods().

sampleOneParticleWithRPCorrelation()

Particle * G4INCL::ParticleSampler::sampleOneParticleWithRPCorrelation ( const ParticleType  t ) const

private

Sample one particle taking into account the rp-correlation.

Definition at line 144 of file G4INCLParticleSampler.cc.

                                                                                          {
// assert(theDensity && thePotential);
    const G4double theFermiMomentum = thePotential->getFermiMomentum(t);
    const ThreeVector momentumVector = Random::sphereVector(theFermiMomentum);
    const G4double momentumAbs = momentumVector.mag();
    const G4double momentumRatio = momentumAbs/theFermiMomentum;
    const G4double reflectionRadius = theDensity->getMaxRFromP(t, momentumRatio);
    const ThreeVector positionVector = Random::sphereVector(reflectionRadius);
    Particle *aParticle = new Particle(t, momentumVector, positionVector);
    aParticle->setUncorrelatedMomentum(momentumAbs);
    return aParticle;
  }

References G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::NuclearDensity::getMaxRFromP(), G4INCL::ThreeVector::mag(), G4INCL::Particle::setUncorrelatedMomentum(), G4INCL::Random::sphereVector(), theDensity, and thePotential.

Referenced by updateSampleOneParticleMethods().

sampleParticles()

ParticleList G4INCL::ParticleSampler::sampleParticles

(

ThreeVector const &

position

)

Definition at line 99 of file G4INCLParticleSampler.cc.

                                                                           {
    ParticleList aList;
    sampleParticlesIntoList(position, aList);
    return aList;
  }

References sampleParticlesIntoList().

sampleParticlesIntoList()

void G4INCL::ParticleSampler::sampleParticlesIntoList	(	ThreeVector const &	position,
		ParticleList &	theList
	)

Definition at line 105 of file G4INCLParticleSampler.cc.

                                                                                                  {
 
    if(sampleOneProton == &ParticleSampler::sampleOneParticleWithoutRPCorrelation) {
      // sampling without correlation, we need to initialize the CDF tables
      theRCDFTable[Proton] = NuclearDensityFactory::createRCDFTable(Proton, theA, theZ);
      thePCDFTable[Proton] = NuclearDensityFactory::createPCDFTable(Proton, theA, theZ);
      theRCDFTable[Neutron] = NuclearDensityFactory::createRCDFTable(Neutron, theA, theZ);
      thePCDFTable[Neutron] = NuclearDensityFactory::createPCDFTable(Neutron, theA, theZ);
      theRCDFTable[Lambda] = NuclearDensityFactory::createRCDFTable(Lambda, theA, theZ);
      thePCDFTable[Lambda] = NuclearDensityFactory::createPCDFTable(Lambda, theA, theZ);
    }
 
    theList.resize(theA);
    if(theA > 2) {
      ParticleType type = Proton;
      ParticleSamplerMethod sampleOneParticle = sampleOneProton;
      for(G4int i = 0; i < theA; ++i) {
        if(i == theZ) { // Nucleons [Z..A-1] are neutrons
          type = Lambda;
          sampleOneParticle = sampleOneNeutron; // hypothesis: Lambdas follow the same rules than neutrons
        }
        if(i == theZ - theS) type = Neutron;
        Particle *p = (this->*sampleOneParticle)(type);
        p->setPosition(position + p->getPosition());
        theList[i] = p;
      }
    } else {
      // For deuterons, only sample the proton position and momentum. The
      // neutron position and momenta are determined by the conditions of
      // vanishing CM position and total momentum.
// assert(theZ==1);
      Particle *aProton = (this->*(this->sampleOneProton))(Proton);
      Particle *aNeutron = new Particle(Neutron, -aProton->getMomentum(), position - aProton->getPosition());
      aProton->setPosition(position + aProton->getPosition());
      theList[0] = aProton;
      theList[1] = aNeutron;
    }
  }

References G4INCL::NuclearDensityFactory::createPCDFTable(), G4INCL::NuclearDensityFactory::createRCDFTable(), G4INCL::Particle::getMomentum(), G4INCL::Particle::getPosition(), G4INCL::Lambda, G4INCL::Neutron, G4INCL::Proton, sampleOneNeutron, sampleOneParticleWithoutRPCorrelation(), sampleOneProton, G4INCL::Particle::setPosition(), theA, thePCDFTable, theRCDFTable, theS, and theZ.

Referenced by G4INCL::Cluster::initializeParticles(), and sampleParticles().

setDensity()

void G4INCL::ParticleSampler::setDensity

(

NuclearDensity const *const

d

)

Setter for theDensity.

Definition at line 71 of file G4INCLParticleSampler.cc.

                                                                 {
    theDensity = d;
    updateSampleOneParticleMethods();
  }

References theDensity, and updateSampleOneParticleMethods().

Referenced by G4INCL::Nucleus::Nucleus(), and G4INCL::Nucleus::setDensity().

setPotential()

void G4INCL::ParticleSampler::setPotential

(

NuclearPotential::INuclearPotential const *const

p

)

Setter for thePotential.

Definition at line 76 of file G4INCLParticleSampler.cc.

                                                                                      {
    thePotential = p;
    updateSampleOneParticleMethods();
  }

References thePotential, and updateSampleOneParticleMethods().

Referenced by G4INCL::Nucleus::Nucleus().

setRPCorrelationCoefficient()

void G4INCL::ParticleSampler::setRPCorrelationCoefficient ( const ParticleType  t, const G4double  corrCoeff  )

inline

Setter for rpCorrelationCoefficient.

Definition at line 86 of file G4INCLParticleSampler.hh.

                                                                                       {
// assert(t==Proton || t==Neutron);
        rpCorrelationCoefficient[t] = corrCoeff;
      }

References rpCorrelationCoefficient.

updateSampleOneParticleMethods()

void G4INCL::ParticleSampler::updateSampleOneParticleMethods ( )

private

Definition at line 81 of file G4INCLParticleSampler.cc.

                                                       {
    if(theDensity && thePotential) {
      if(rpCorrelationCoefficient[Proton]>0.99999) {
        sampleOneProton = &ParticleSampler::sampleOneParticleWithRPCorrelation;
      } else {
        sampleOneProton = &ParticleSampler::sampleOneParticleWithFuzzyRPCorrelation;
      }
      if(rpCorrelationCoefficient[Neutron]>0.99999) {
        sampleOneNeutron = &ParticleSampler::sampleOneParticleWithRPCorrelation;
      } else {
        sampleOneNeutron = &ParticleSampler::sampleOneParticleWithFuzzyRPCorrelation;
      }
    } else {
      sampleOneProton = &ParticleSampler::sampleOneParticleWithoutRPCorrelation;
      sampleOneNeutron = &ParticleSampler::sampleOneParticleWithoutRPCorrelation;
    }
  }

References G4INCL::Neutron, G4INCL::Proton, rpCorrelationCoefficient, sampleOneNeutron, sampleOneParticleWithFuzzyRPCorrelation(), sampleOneParticleWithoutRPCorrelation(), sampleOneParticleWithRPCorrelation(), sampleOneProton, theDensity, and thePotential.

Referenced by setDensity(), and setPotential().

Field Documentation

rpCorrelationCoefficient

G4double G4INCL::ParticleSampler::rpCorrelationCoefficient[UnknownParticle]

private

Correlation coefficients for the r-p correlation.

Definition at line 143 of file G4INCLParticleSampler.hh.

Referenced by getRPCorrelationCoefficient(), ParticleSampler(), sampleOneParticleWithFuzzyRPCorrelation(), setRPCorrelationCoefficient(), and updateSampleOneParticleMethods().

sampleOneNeutron

ParticleSamplerMethod G4INCL::ParticleSampler::sampleOneNeutron

private

Sample a list of particles.

This method is a pointer to the method that does the real work for neutrons.

Definition at line 110 of file G4INCLParticleSampler.hh.

Referenced by sampleParticlesIntoList(), and updateSampleOneParticleMethods().

sampleOneProton

ParticleSamplerMethod G4INCL::ParticleSampler::sampleOneProton

private

Sample a list of particles.

This method is a pointer to the method that does the real work for protons.

Definition at line 104 of file G4INCLParticleSampler.hh.

Referenced by sampleParticlesIntoList(), and updateSampleOneParticleMethods().

theA

const G4int G4INCL::ParticleSampler::theA

private

Mass number.

Definition at line 122 of file G4INCLParticleSampler.hh.

Referenced by sampleParticlesIntoList().

theDensity

NuclearDensity const* G4INCL::ParticleSampler::theDensity

private

Pointer to the Cluster's NuclearDensity.

Definition at line 137 of file G4INCLParticleSampler.hh.

Referenced by getDensity(), sampleOneParticleWithFuzzyRPCorrelation(), sampleOneParticleWithRPCorrelation(), setDensity(), and updateSampleOneParticleMethods().

thePCDFTable

InterpolationTable const* G4INCL::ParticleSampler::thePCDFTable[UnknownParticle]

private

Array of pointers to the p-space CDF table.

Definition at line 134 of file G4INCLParticleSampler.hh.

Referenced by ParticleSampler(), sampleOneParticleWithoutRPCorrelation(), and sampleParticlesIntoList().

thePotential

NuclearPotential::INuclearPotential const* G4INCL::ParticleSampler::thePotential

private

Pointer to the Cluster's NuclearPotential.

Definition at line 140 of file G4INCLParticleSampler.hh.

Referenced by getPotential(), sampleOneParticleWithFuzzyRPCorrelation(), sampleOneParticleWithRPCorrelation(), setPotential(), and updateSampleOneParticleMethods().

theRCDFTable

InterpolationTable const* G4INCL::ParticleSampler::theRCDFTable[UnknownParticle]

private

Array of pointers to the r-space CDF table.

Definition at line 131 of file G4INCLParticleSampler.hh.

Referenced by ParticleSampler(), sampleOneParticleWithoutRPCorrelation(), and sampleParticlesIntoList().

theS

const G4int G4INCL::ParticleSampler::theS

private

Strangeness number.

Definition at line 128 of file G4INCLParticleSampler.hh.

Referenced by sampleParticlesIntoList().

theZ

const G4int G4INCL::ParticleSampler::theZ

private

Charge number.

Definition at line 125 of file G4INCLParticleSampler.hh.

Referenced by sampleParticlesIntoList().

---

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

• source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLParticleSampler.hh
• source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLParticleSampler.cc