G4INCL::ParticleSampler Class Reference

#include <G4INCLParticleSampler.hh>

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

Detailed Description

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

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

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

sampleParticles()

ParticleList G4INCL::ParticleSampler::sampleParticles

(

ThreeVector const &

position

)

Definition at line 99 of file G4INCLParticleSampler.cc.

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

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

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

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

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

---

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

• G4INCLParticleSampler.hh
• G4INCLParticleSampler.cc