#include <G4INCLParticleSampler.hh>

Public Member Functions
	ParticleSampler (const G4int A, const G4int Z, InverseInterpolationTable const const rCDFTable, InverseInterpolationTable const const pCDFTable)
	Constructor.

	~ParticleSampler ()
	Destructor.

NuclearDensity const *	getDensity () const
	Getter for theDensity.

NuclearPotential::INuclearPotential const *	getPotential () const
	Getter for thePotential.

void	setDensity (NuclearDensity const *const d)
	Setter for theDensity.

void	setPotential (NuclearPotential::INuclearPotential const *const p)
	Setter for thePotential.

ParticleList	sampleParticles (ThreeVector const &position) const

Detailed Description

Definition at line 54 of file G4INCLParticleSampler.hh.

Constructor & Destructor Documentation

◆ ParticleSampler()

G4INCL::ParticleSampler::ParticleSampler	(	const G4int	A,
		const G4int	Z,
		InverseInterpolationTable const *const	rCDFTable,
		InverseInterpolationTable const *const	pCDFTable
	)

Constructor.

Parameters

A	the mass number
Z	the charge number
rCDFTable	the interpolation table for the inverse CDF in r-space
pCDFTable	the interpolation table for the inverse CDF in p-space

Definition at line 50 of file G4INCLParticleSampler.cc.

                                                                                                                                                                     :
    sampleOneParticle(&ParticleSampler::sampleOneParticleWithoutRPCorrelation),
    theA(A),
    theZ(Z),
    theRCDFTable(rCDFTable),
    thePCDFTable(pCDFTable),
    theDensity(NULL),
    thePotential(NULL)
  {}

◆ ~ParticleSampler()

G4INCL::ParticleSampler::~ParticleSampler ( )

Destructor.

Definition at line 60 of file G4INCLParticleSampler.cc.

60 {

61 }

Member Function Documentation

◆ getDensity()

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

inline

Getter for theDensity.

Definition at line 70 of file G4INCLParticleSampler.hh.

70{ return theDensity; }

◆ getPotential()

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

inline

Getter for thePotential.

Definition at line 73 of file G4INCLParticleSampler.hh.

73{ return thePotential; }

◆ sampleParticles()

ParticleList G4INCL::ParticleSampler::sampleParticles ( ThreeVector const & position ) const

Definition at line 80 of file G4INCLParticleSampler.cc.

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

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

◆ setDensity()

void G4INCL::ParticleSampler::setDensity ( NuclearDensity const *const d )

Setter for theDensity.

Definition at line 63 of file G4INCLParticleSampler.cc.

                                                                 {
    theDensity = d;
    updateSampleOneParticleMethod();
  }

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

◆ setPotential()

void G4INCL::ParticleSampler::setPotential ( NuclearPotential::INuclearPotential const *const p )

Setter for thePotential.

Definition at line 68 of file G4INCLParticleSampler.cc.

                                                                                      {
    thePotential = p;
    updateSampleOneParticleMethod();
  }