G4INCL::ParticleSampler Class Reference

#include <G4INCLParticleSampler.hh>

Collaboration diagram for G4INCL::ParticleSampler:

Public Member Functions
	ParticleSampler (const G4int A, const G4int Z)
	Constructor. More...
	~ParticleSampler ()
	Destructor. More...
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...
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...
ParticleList	sampleParticles (ThreeVector const &position)
void	sampleParticlesIntoList (ThreeVector const &position, ParticleList &theList)

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

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

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

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),
    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);
  }

Here is the call graph for this function:

~ParticleSampler()

G4INCL::ParticleSampler::~ParticleSampler

(

)

Destructor.

Definition at line 66 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; }

getPotential()

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

inline

Getter for thePotential.

Definition at line 71 of file G4INCLParticleSampler.hh.

{ return 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];
      }

Here is the call graph for this function:

sampleOneParticleWithFuzzyRPCorrelation()

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

private

Sample one particle with a fuzzy rp-correlation.

Definition at line 160 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;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

sampleOneParticleWithoutRPCorrelation()

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

private

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

Definition at line 152 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);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

sampleOneParticleWithRPCorrelation()

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

private

Sample one particle taking into account the rp-correlation.

Definition at line 139 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;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

sampleParticles()

ParticleList G4INCL::ParticleSampler::sampleParticles

(

ThreeVector const &

position

)

Definition at line 97 of file G4INCLParticleSampler.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

sampleParticlesIntoList()

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

Definition at line 103 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);
    }

    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 = Neutron;
          sampleOneParticle = sampleOneNeutron;
        }
        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;
    }
  }

Here is the call graph for this function:

Here is the caller graph for this function:

setDensity()

void G4INCL::ParticleSampler::setDensity

(

NuclearDensity const *const

d

)

Setter for theDensity.

Definition at line 69 of file G4INCLParticleSampler.cc.

                                                                 {
    theDensity = d;
    updateSampleOneParticleMethods();
  }

Here is the call graph for this function:

Here is the caller graph for this function:

setPotential()

void G4INCL::ParticleSampler::setPotential

(

NuclearPotential::INuclearPotential const *const

p

)

Setter for thePotential.

Definition at line 74 of file G4INCLParticleSampler.cc.

                                                                                      {
    thePotential = p;
    updateSampleOneParticleMethods();
  }

Here is the call graph for this function:

Here is the caller graph for this function:

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;
      }

Here is the call graph for this function:

updateSampleOneParticleMethods()

void G4INCL::ParticleSampler::updateSampleOneParticleMethods ( )

private

Definition at line 79 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;
    }
  }

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

rpCorrelationCoefficient

G4double G4INCL::ParticleSampler::rpCorrelationCoefficient[UnknownParticle]

private

Correlation coefficients for the r-p correlation.

Definition at line 140 of file G4INCLParticleSampler.hh.

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.

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.

theA

const G4int G4INCL::ParticleSampler::theA

private

Mass number.

Definition at line 122 of file G4INCLParticleSampler.hh.

theDensity

NuclearDensity const* G4INCL::ParticleSampler::theDensity

private

Pointer to the Cluster's NuclearDensity.

Definition at line 134 of file G4INCLParticleSampler.hh.

thePCDFTable

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

private

Array of pointers to the p-space CDF table.

Definition at line 131 of file G4INCLParticleSampler.hh.

thePotential

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

private

Pointer to the Cluster's NuclearPotential.

Definition at line 137 of file G4INCLParticleSampler.hh.

theRCDFTable

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

private

Array of pointers to the r-space CDF table.

Definition at line 128 of file G4INCLParticleSampler.hh.

theZ

const G4int G4INCL::ParticleSampler::theZ

private

Charge number.

Definition at line 125 of file G4INCLParticleSampler.hh.

---

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

• G4INCLParticleSampler.hh
• G4INCLParticleSampler.cc