G4INCL::ParticleSampler Class Reference

#include <G4INCLParticleSampler.hh>

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)

Detailed Description

Definition at line 52 of file G4INCLParticleSampler.hh.

Constructor & Destructor Documentation

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

Constructor.

Parameters

A	the mass number
Z	the charge number

Definition at line 50 of file G4INCLParticleSampler.cc.

References G4INCL::ParticleTable::getRPCorrelationCoefficient(), G4INCL::Neutron, G4INCL::Proton, and G4INCL::UnknownParticle.

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

G4INCL::ParticleSampler::~ParticleSampler

(

)

Destructor.

Definition at line 65 of file G4INCLParticleSampler.cc.

                                    {
  }

Member Function Documentation

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

inline

Getter for theDensity.

Definition at line 66 of file G4INCLParticleSampler.hh.

{ return theDensity; }

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

inline

Getter for thePotential.

Definition at line 69 of file G4INCLParticleSampler.hh.

{ return thePotential; }

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

inline

Getter for rpCorrelationCoefficient.

Definition at line 72 of file G4INCLParticleSampler.hh.

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

ParticleList G4INCL::ParticleSampler::sampleParticles

(

ThreeVector const &

position

)

Definition at line 96 of file G4INCLParticleSampler.cc.

References G4INCL::NuclearDensityFactory::createPCDFTable(), G4INCL::NuclearDensityFactory::createRCDFTable(), G4INCL::Particle::getMomentum(), G4INCL::Particle::getPosition(), G4INCL::Neutron, G4INCL::Proton, and G4INCL::Particle::setPosition().

Referenced by G4INCL::Cluster::initializeParticles().

                                                                           {

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

    ParticleList theList;
    if(theA > 2) {
      ParticleType type = Proton;
      ParticleSamplerMethod sampleOneParticle = sampleOneProton;
      for(G4int i = 1; i <= theA; ++i) {
        if(i == (theZ + 1)) { // Nucleons [Z+1..A] are neutrons
          type = Neutron;
          sampleOneParticle = sampleOneNeutron;
        }
        Particle *p = (this->*sampleOneParticle)(type);
        p->setPosition(position + p->getPosition());
        theList.push_back(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.push_back(aProton);
      theList.push_back(aNeutron);
    }

    return theList;
  }

void G4INCL::ParticleSampler::setDensity

(

NuclearDensity const *const

d

)

Setter for theDensity.

Definition at line 68 of file G4INCLParticleSampler.cc.

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

                                                                 {
    theDensity = d;
    updateSampleOneParticleMethods();
  }

void G4INCL::ParticleSampler::setPotential

(

NuclearPotential::INuclearPotential const *const

p

)

Setter for thePotential.

Definition at line 73 of file G4INCLParticleSampler.cc.

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

                                                                                      {
    thePotential = p;
    updateSampleOneParticleMethods();
  }

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

inline

Setter for rpCorrelationCoefficient.

Definition at line 84 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