G4INCL::NuclearDensity Class Reference

#include <G4INCLNuclearDensity.hh>

Collaboration diagram for G4INCL::NuclearDensity:

Public Member Functions
	NuclearDensity (const G4int A, const G4int Z, InterpolationTable const *const rpCorrelationTableProton, InterpolationTable const *const rpCorrelationTableNeutron)
	~NuclearDensity ()
	NuclearDensity (const NuclearDensity &rhs)
	Copy constructor. More...
NuclearDensity &	operator= (const NuclearDensity &rhs)
	Assignment operator. More...
void	swap (NuclearDensity &rhs)
	Helper method for the assignment operator. More...
G4double	getMaxRFromP (const ParticleType t, const G4double p) const
	Get the maximum allowed radius for a given momentum. More...
G4double	getMinPFromR (const ParticleType t, const G4double r) const
G4double	getMaximumRadius () const
G4double	getTransmissionRadius (Particle const *const p) const
	The radius used for calculating the transmission coefficient. More...
G4double	getTransmissionRadius (ParticleType type) const
	The radius used for calculating the transmission coefficient. More...
G4int	getA () const
	Get the mass number. More...
G4int	getZ () const
	Get the charge number. More...
G4double	getProtonNuclearRadius () const
void	setProtonNuclearRadius (const G4double r)

Private Member Functions
void	initializeTransmissionRadii ()
	Initialize the transmission radius. More...

Private Attributes
G4int	theA
G4int	theZ
G4double	theMaximumRadius
G4double	theProtonNuclearRadius
	Represents INCL4.5's R0 variable. More...
G4double	transmissionRadius [UnknownParticle]
InterpolationTable const *	rFromP [UnknownParticle]
InterpolationTable const *	pFromR [UnknownParticle]

Detailed Description

Definition at line 54 of file G4INCLNuclearDensity.hh.

Constructor & Destructor Documentation

NuclearDensity() [1/2]

G4INCL::NuclearDensity::NuclearDensity	(	const G4int	A,
		const G4int	Z,
		InterpolationTable const *const	rpCorrelationTableProton,
		InterpolationTable const *const	rpCorrelationTableNeutron
	)

Definition at line 45 of file G4INCLNuclearDensity.cc.

                                                                                                                                                                                    :
    theA(A),
    theZ(Z),
    theMaximumRadius(std::min((*rpCorrelationTableProton)(1.), (*rpCorrelationTableNeutron)(1.))),
    theProtonNuclearRadius(ParticleTable::getNuclearRadius(Proton,theA,theZ))
  {
    std::fill(rFromP, rFromP + UnknownParticle, static_cast<InterpolationTable*>(NULL));
    rFromP[Proton] = rpCorrelationTableProton;
    rFromP[Neutron] = rpCorrelationTableNeutron;
    rFromP[DeltaPlusPlus] = rpCorrelationTableProton;
    rFromP[DeltaPlus] = rpCorrelationTableProton;
    rFromP[DeltaZero] = rpCorrelationTableNeutron;
    rFromP[DeltaMinus] = rpCorrelationTableNeutron;
    // The interpolation table for local-energy look-ups is simply obtained by
    // inverting the r-p correlation table.
    std::fill(pFromR, pFromR + UnknownParticle, static_cast<InterpolationTable*>(NULL));
    pFromR[Proton] = new InterpolationTable(rFromP[Proton]->getNodeValues(), rFromP[Proton]->getNodeAbscissae());
    pFromR[Neutron] = new InterpolationTable(rFromP[Neutron]->getNodeValues(), rFromP[Neutron]->getNodeAbscissae());
    pFromR[DeltaPlusPlus] = new InterpolationTable(rFromP[DeltaPlusPlus]->getNodeValues(), rFromP[DeltaPlusPlus]->getNodeAbscissae());
    pFromR[DeltaPlus] = new InterpolationTable(rFromP[DeltaPlus]->getNodeValues(), rFromP[DeltaPlus]->getNodeAbscissae());
    pFromR[DeltaZero] = new InterpolationTable(rFromP[DeltaZero]->getNodeValues(), rFromP[DeltaZero]->getNodeAbscissae());
    pFromR[DeltaMinus] = new InterpolationTable(rFromP[DeltaMinus]->getNodeValues(), rFromP[DeltaMinus]->getNodeAbscissae());
    INCL_DEBUG("Interpolation table for proton local energy (A=" << theA << ", Z=" << theZ << ") initialised:"
          << '\n'
          << pFromR[Proton]->print()
          << '\n'
          << "Interpolation table for neutron local energy (A=" << theA << ", Z=" << theZ << ") initialised:"
          << '\n'
          << pFromR[Neutron]->print()
          << '\n'
          << "Interpolation table for delta++ local energy (A=" << theA << ", Z=" << theZ << ") initialised:"
          << '\n'
          << pFromR[DeltaPlusPlus]->print()
          << '\n'
          << "Interpolation table for delta+ local energy (A=" << theA << ", Z=" << theZ << ") initialised:"
          << '\n'
          << pFromR[DeltaPlus]->print()
          << '\n'
          << "Interpolation table for delta0 local energy (A=" << theA << ", Z=" << theZ << ") initialised:"
          << '\n'
          << pFromR[DeltaZero]->print()
          << '\n'
          << "Interpolation table for delta- local energy (A=" << theA << ", Z=" << theZ << ") initialised:"
          << '\n'
          << pFromR[DeltaMinus]->print()
          << '\n');
    initializeTransmissionRadii();
  }

Here is the call graph for this function:

~NuclearDensity()

G4INCL::NuclearDensity::~NuclearDensity

(

)

Definition at line 94 of file G4INCLNuclearDensity.cc.

                                  {
    // We don't delete the rFromP tables, which are cached in the
    // NuclearDensityFactory
    delete pFromR[Proton];
    delete pFromR[Neutron];
    delete pFromR[DeltaPlusPlus];
    delete pFromR[DeltaPlus];
    delete pFromR[DeltaZero];
    delete pFromR[DeltaMinus];
  }

NuclearDensity() [2/2]

G4INCL::NuclearDensity::NuclearDensity

(

const NuclearDensity &

rhs

)

Copy constructor.

Definition at line 105 of file G4INCLNuclearDensity.cc.

                                                          :
    theA(rhs.theA),
    theZ(rhs.theZ),
    theMaximumRadius(rhs.theMaximumRadius),
    theProtonNuclearRadius(rhs.theProtonNuclearRadius)
  {
    // rFromP is owned by NuclearDensityFactory, so shallow copy is sufficient
    std::fill(rFromP, rFromP + UnknownParticle, static_cast<InterpolationTable*>(NULL));
    rFromP[Proton] = rhs.rFromP[Proton];
    rFromP[Neutron] = rhs.rFromP[Neutron];
    rFromP[DeltaPlusPlus] = rhs.rFromP[DeltaPlusPlus];
    rFromP[DeltaPlus] = rhs.rFromP[DeltaPlus];
    rFromP[DeltaZero] = rhs.rFromP[DeltaZero];
    rFromP[DeltaMinus] = rhs.rFromP[DeltaMinus];
    // deep copy for pFromR
    std::fill(pFromR, pFromR + UnknownParticle, static_cast<InterpolationTable*>(NULL));
    pFromR[Proton] = new InterpolationTable(*(rhs.pFromR[Proton]));
    pFromR[Neutron] = new InterpolationTable(*(rhs.pFromR[Neutron]));
    pFromR[DeltaPlusPlus] = new InterpolationTable(*(rhs.pFromR[DeltaPlusPlus]));
    pFromR[DeltaPlus] = new InterpolationTable(*(rhs.pFromR[DeltaPlus]));
    pFromR[DeltaZero] = new InterpolationTable(*(rhs.pFromR[DeltaZero]));
    pFromR[DeltaMinus] = new InterpolationTable(*(rhs.pFromR[DeltaMinus]));
    std::copy(rhs.transmissionRadius, rhs.transmissionRadius+UnknownParticle, transmissionRadius);
  }

Member Function Documentation

getA()

G4int G4INCL::NuclearDensity::getA ( ) const

inline

Get the mass number.

Definition at line 104 of file G4INCLNuclearDensity.hh.

{ return theA; }

getMaximumRadius()

G4double G4INCL::NuclearDensity::getMaximumRadius ( ) const

inline

Definition at line 78 of file G4INCLNuclearDensity.hh.

{ return theMaximumRadius; };

Here is the caller graph for this function:

getMaxRFromP()

G4double G4INCL::NuclearDensity::getMaxRFromP	(	const ParticleType	t,
		const G4double	p
	)		const

Get the maximum allowed radius for a given momentum.

Parameters

t	type of the particle
p	absolute value of the particle momentum, divided by the relevant Fermi momentum.

Returns

maximum allowed radius.

Definition at line 170 of file G4INCLNuclearDensity.cc.

                                                                                    {
// assert(t==Proton || t==Neutron || t==DeltaPlusPlus || t==DeltaPlus || t==DeltaZero || t==DeltaMinus);
    return (*(rFromP[t]))(p);
  }

Here is the caller graph for this function:

getMinPFromR()

G4double G4INCL::NuclearDensity::getMinPFromR	(	const ParticleType	t,
		const G4double	r
	)		const

Definition at line 175 of file G4INCLNuclearDensity.cc.

                                                                                    {
// assert(t==Proton || t==Neutron || t==DeltaPlusPlus || t==DeltaPlus || t==DeltaZero || t==DeltaMinus);
    return (*(pFromR[t]))(r);
  }

Here is the caller graph for this function:

getProtonNuclearRadius()

G4double G4INCL::NuclearDensity::getProtonNuclearRadius ( ) const

inline

Definition at line 109 of file G4INCLNuclearDensity.hh.

{ return theProtonNuclearRadius; }

Here is the caller graph for this function:

getTransmissionRadius() [1/2]

G4double G4INCL::NuclearDensity::getTransmissionRadius ( Particle const *const  p ) const

inline

The radius used for calculating the transmission coefficient.

Returns

the radius

Definition at line 84 of file G4INCLNuclearDensity.hh.

                                                                   {
      const ParticleType t = p->getType();
// assert(t!=Neutron && t!=PiZero && t!=DeltaZero); // no neutral particles here
      if(t==Composite) {
        return transmissionRadius[t] +
          ParticleTable::getNuclearRadius(t, p->getA(), p->getZ());
      } else
        return transmissionRadius[t];
    };

Here is the call graph for this function:

Here is the caller graph for this function:

getTransmissionRadius() [2/2]

G4double G4INCL::NuclearDensity::getTransmissionRadius ( ParticleType  type ) const

inline

The radius used for calculating the transmission coefficient.

Returns

the radius

Definition at line 98 of file G4INCLNuclearDensity.hh.

                                                            {
// assert(type!=Composite);
      return transmissionRadius[type];
    };

getZ()

G4int G4INCL::NuclearDensity::getZ ( ) const

inline

Get the charge number.

Definition at line 107 of file G4INCLNuclearDensity.hh.

{ return theZ; }

initializeTransmissionRadii()

void G4INCL::NuclearDensity::initializeTransmissionRadii ( )

private

Initialize the transmission radius.

Definition at line 156 of file G4INCLNuclearDensity.cc.

                                                   {
    const G4double theProtonRadius = 0.88; // fm
    const G4double theProtonTransmissionRadius = theProtonNuclearRadius + theProtonRadius;

    transmissionRadius[Proton] = theProtonTransmissionRadius;
    transmissionRadius[PiPlus] = theProtonNuclearRadius;
    transmissionRadius[PiMinus] = theProtonNuclearRadius;
    transmissionRadius[DeltaPlusPlus] = theProtonTransmissionRadius;
    transmissionRadius[DeltaPlus] = theProtonTransmissionRadius;
    transmissionRadius[DeltaMinus] = theProtonTransmissionRadius;
    transmissionRadius[Composite] = theProtonNuclearRadius;
    // transmission radii for neutral particles intentionally left uninitialised
  }

Here is the caller graph for this function:

operator=()

NuclearDensity & G4INCL::NuclearDensity::operator=

(

const NuclearDensity &

rhs

)

Assignment operator.

Definition at line 130 of file G4INCLNuclearDensity.cc.

                                                                     {
    NuclearDensity temporaryDensity(rhs);
    swap(temporaryDensity);
    return *this;
  }

Here is the call graph for this function:

setProtonNuclearRadius()

void G4INCL::NuclearDensity::setProtonNuclearRadius ( const G4double  r )

inline

Definition at line 110 of file G4INCLNuclearDensity.hh.

{ theProtonNuclearRadius = r; }

Here is the call graph for this function:

swap()

void G4INCL::NuclearDensity::swap

(

NuclearDensity &

rhs

)

Helper method for the assignment operator.

Definition at line 136 of file G4INCLNuclearDensity.cc.

                                               {
    std::swap(theA, rhs.theA);
    std::swap(theZ, rhs.theZ);
    std::swap(theMaximumRadius, rhs.theMaximumRadius);
    std::swap(theProtonNuclearRadius, rhs.theProtonNuclearRadius);
    std::swap_ranges(transmissionRadius, transmissionRadius+UnknownParticle, rhs.transmissionRadius);
    std::swap(rFromP[Proton], rhs.rFromP[Proton]);
    std::swap(rFromP[Neutron], rhs.rFromP[Neutron]);
    std::swap(rFromP[DeltaPlusPlus], rhs.rFromP[DeltaPlusPlus]);
    std::swap(rFromP[DeltaPlus], rhs.rFromP[DeltaPlus]);
    std::swap(rFromP[DeltaZero], rhs.rFromP[DeltaZero]);
    std::swap(rFromP[DeltaMinus], rhs.rFromP[DeltaMinus]);
    std::swap(pFromR[Proton], rhs.pFromR[Proton]);
    std::swap(pFromR[Neutron], rhs.pFromR[Neutron]);
    std::swap(pFromR[DeltaPlusPlus], rhs.pFromR[DeltaPlusPlus]);
    std::swap(pFromR[DeltaPlus], rhs.pFromR[DeltaPlus]);
    std::swap(pFromR[DeltaZero], rhs.pFromR[DeltaZero]);
    std::swap(pFromR[DeltaMinus], rhs.pFromR[DeltaMinus]);
 }

Here is the caller graph for this function:

Member Data Documentation

pFromR

InterpolationTable const* G4INCL::NuclearDensity::pFromR[UnknownParticle]

private

Definition at line 126 of file G4INCLNuclearDensity.hh.

rFromP

InterpolationTable const* G4INCL::NuclearDensity::rFromP[UnknownParticle]

private

Definition at line 125 of file G4INCLNuclearDensity.hh.

theA

G4int G4INCL::NuclearDensity::theA

private

Definition at line 117 of file G4INCLNuclearDensity.hh.

theMaximumRadius

G4double G4INCL::NuclearDensity::theMaximumRadius

private

Definition at line 118 of file G4INCLNuclearDensity.hh.

theProtonNuclearRadius

G4double G4INCL::NuclearDensity::theProtonNuclearRadius

private

Represents INCL4.5's R0 variable.

Definition at line 120 of file G4INCLNuclearDensity.hh.

theZ

G4int G4INCL::NuclearDensity::theZ

private

Definition at line 117 of file G4INCLNuclearDensity.hh.

transmissionRadius

G4double G4INCL::NuclearDensity::transmissionRadius[UnknownParticle]

private

Definition at line 123 of file G4INCLNuclearDensity.hh.

---

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

• G4INCLNuclearDensity.hh
• G4INCLNuclearDensity.cc