G4INCL::CrossSections Namespace Reference

Functions
G4double	elastic (Particle const *const p1, Particle const *const p2)
G4double	total (Particle const *const p1, Particle const *const p2)
G4double	NDeltaToNN (Particle const *const p1, Particle const *const p2)
G4double	NNToNDelta (Particle const *const p1, Particle const *const p2)
G4double	NNToxPiNN (const G4int xpi, Particle const *const p1, Particle const *const p2)
G4double	piNToDelta (Particle const *const p1, Particle const *const p2)
G4double	piNToxPiN (const G4int xpi, Particle const *const p1, Particle const *const p2)
G4double	calculateNNAngularSlope (G4double energyCM, G4int iso)
G4double	interactionDistancePiN (const G4double projectileKineticEnergy)
G4double	interactionDistanceNN (const ParticleSpecies &aSpecies, const G4double kineticEnergy)
void	setCrossSections (ICrossSections *c)
void	deleteCrossSections ()
void	initialize (Config const *const theConfig)
	Initialize the Coulomb-distortion algorithm. More...

Function Documentation

calculateNNAngularSlope()

G4double G4INCL::CrossSections::calculateNNAngularSlope	(	G4double	energyCM,
		G4int	iso
	)

Definition at line 82 of file G4INCLCrossSections.cc.

                                                                   {
      return theCrossSections->calculateNNAngularSlope(energyCM, iso);
    }

Here is the caller graph for this function:

deleteCrossSections()

void G4INCL::CrossSections::deleteCrossSections

(

)

Definition at line 154 of file G4INCLCrossSections.cc.

                               {
      delete theCrossSections;
      theCrossSections = NULL;
    }

Here is the caller graph for this function:

elastic()

G4double G4INCL::CrossSections::elastic	(	Particle const *const	p1,
		Particle const *const	p2
	)

Definition at line 54 of file G4INCLCrossSections.cc.

                                                                           {
      return theCrossSections->elastic(p1,p2);
    }

Here is the caller graph for this function:

initialize()

void G4INCL::CrossSections::initialize

(

Config const *const

theConfig

)

Initialize the Coulomb-distortion algorithm.

Initialize the Coulomb-distortion algorithm.

Definition at line 159 of file G4INCLCrossSections.cc.

                                                    {
      CrossSectionsType crossSections = theConfig->getCrossSectionsType();
      if(crossSections == INCL46CrossSections)
        setCrossSections(new CrossSectionsINCL46);
      else if(crossSections == MultiPionsCrossSections)
        setCrossSections(new CrossSectionsMultiPions);
      else if(crossSections == TruncatedMultiPionsCrossSections) {
        const G4int nMaxPi = theConfig->getMaxNumberMultipions();
        if(nMaxPi>0)
          setCrossSections(new CrossSectionsTruncatedMultiPions(nMaxPi));
        else {
          INCL_WARN("Truncated multipion cross sections were requested, but the specified maximum\n"
                    << "number of pions is <=0. Falling back to standard multipion cross-sections.\n");
          setCrossSections(new CrossSectionsMultiPions);
        }
      }
    }

Here is the call graph for this function:

Here is the caller graph for this function:

interactionDistanceNN()

G4double G4INCL::CrossSections::interactionDistanceNN	(	const ParticleSpecies &	aSpecies,
		const G4double	kineticEnergy
	)

Definition at line 119 of file G4INCLCrossSections.cc.

                                                                                                  {
// assert(aSpecies.theType==Proton || aSpecies.theType==Neutron || aSpecies.theType==Composite);
// assert(aSpecies.theA>0);
      ThreeVector nullVector;
      ThreeVector unitVector(0.,0.,1.);

      const G4double kineticEnergyPerNucleon = kineticEnergy / aSpecies.theA;

      Particle protonProjectile(Proton, unitVector, nullVector);
      protonProjectile.setEnergy(protonProjectile.getMass()+kineticEnergyPerNucleon);
      protonProjectile.adjustMomentumFromEnergy();
      Particle neutronProjectile(Neutron, unitVector, nullVector);
      neutronProjectile.setEnergy(neutronProjectile.getMass()+kineticEnergyPerNucleon);
      neutronProjectile.adjustMomentumFromEnergy();

      Particle protonTarget(Proton, nullVector, nullVector);
      Particle neutronTarget(Neutron, nullVector, nullVector);
      const G4double sigmapp = total(&protonProjectile, &protonTarget);
      const G4double sigmapn = total(&protonProjectile, &neutronTarget);
      const G4double sigmann = total(&neutronProjectile, &neutronTarget);
      /* We compute the interaction distance from the largest of the NN cross
       * sections. Note that this is different from INCL4.6, which just takes the
       * average of the four, and will in general lead to a different geometrical
       * cross section.
       */
      const G4double largestSigma = std::max(sigmapp, std::max(sigmapn, sigmann));
      const G4double interactionDistance = std::sqrt(largestSigma/Math::tenPi);

      return interactionDistance;
    }

Here is the call graph for this function:

Here is the caller graph for this function:

interactionDistancePiN()

G4double G4INCL::CrossSections::interactionDistancePiN

(

const G4double

projectileKineticEnergy

)

Definition at line 86 of file G4INCLCrossSections.cc.

                                                                            {
      ThreeVector nullVector;
      ThreeVector unitVector(0., 0., 1.);

      Particle piPlusProjectile(PiPlus, unitVector, nullVector);
      piPlusProjectile.setEnergy(piPlusProjectile.getMass()+projectileKineticEnergy);
      piPlusProjectile.adjustMomentumFromEnergy();
      Particle piZeroProjectile(PiZero, unitVector, nullVector);
      piZeroProjectile.setEnergy(piZeroProjectile.getMass()+projectileKineticEnergy);
      piZeroProjectile.adjustMomentumFromEnergy();
      Particle piMinusProjectile(PiMinus, unitVector, nullVector);
      piMinusProjectile.setEnergy(piMinusProjectile.getMass()+projectileKineticEnergy);
      piMinusProjectile.adjustMomentumFromEnergy();

      Particle protonTarget(Proton, nullVector, nullVector);
      Particle neutronTarget(Neutron, nullVector, nullVector);
      const G4double sigmapipp = total(&piPlusProjectile, &protonTarget);
      const G4double sigmapipn = total(&piPlusProjectile, &neutronTarget);
      const G4double sigmapi0p = total(&piZeroProjectile, &protonTarget);
      const G4double sigmapi0n = total(&piZeroProjectile, &neutronTarget);
      const G4double sigmapimp = total(&piMinusProjectile, &protonTarget);
      const G4double sigmapimn = total(&piMinusProjectile, &neutronTarget);
      /* We compute the interaction distance from the largest of the pi-N cross
       * sections. Note that this is different from INCL4.6, which just takes the
       * average of the six, and will in general lead to a different geometrical
       * cross section.
       */
      const G4double largestSigma = std::max(sigmapipp, std::max(sigmapipn, std::max(sigmapi0p, std::max(sigmapi0n, std::max(sigmapimp,sigmapimn)))));
      const G4double interactionDistance = std::sqrt(largestSigma/Math::tenPi);

      return interactionDistance;
    }

Here is the call graph for this function:

Here is the caller graph for this function:

NDeltaToNN()

G4double G4INCL::CrossSections::NDeltaToNN	(	Particle const *const	p1,
		Particle const *const	p2
	)

Definition at line 62 of file G4INCLCrossSections.cc.

                                                                              {
      return theCrossSections->NDeltaToNN(p1,p2);
    }

Here is the caller graph for this function:

NNToNDelta()

G4double G4INCL::CrossSections::NNToNDelta	(	Particle const *const	p1,
		Particle const *const	p2
	)

Definition at line 66 of file G4INCLCrossSections.cc.

                                                                              {
      return theCrossSections->NNToNDelta(p1,p2);
    }

Here is the caller graph for this function:

NNToxPiNN()

G4double G4INCL::CrossSections::NNToxPiNN	(	const G4int	xpi,
		Particle const *const	p1,
		Particle const *const	p2
	)

Definition at line 70 of file G4INCLCrossSections.cc.

                                                                                                {
          return theCrossSections->NNToxPiNN(xpi,p1,p2);
      }

Here is the caller graph for this function:

piNToDelta()

G4double G4INCL::CrossSections::piNToDelta	(	Particle const *const	p1,
		Particle const *const	p2
	)

Definition at line 74 of file G4INCLCrossSections.cc.

                                                                                {
          return theCrossSections->piNToDelta(p1,p2);
      }

Here is the caller graph for this function:

piNToxPiN()

G4double G4INCL::CrossSections::piNToxPiN	(	const G4int	xpi,
		Particle const *const	p1,
		Particle const *const	p2
	)

Definition at line 78 of file G4INCLCrossSections.cc.

                                                                                                {
          return theCrossSections->piNToxPiN(xpi,p1,p2);
      }

Here is the caller graph for this function:

setCrossSections()

void G4INCL::CrossSections::setCrossSections

(

ICrossSections *

c

)

Definition at line 150 of file G4INCLCrossSections.cc.

                                             {
      theCrossSections = c;
    }

Here is the caller graph for this function:

total()

G4double G4INCL::CrossSections::total	(	Particle const *const	p1,
		Particle const *const	p2
	)

Definition at line 58 of file G4INCLCrossSections.cc.

                                                                         {
      return theCrossSections->total(p1,p2);
    }

Here is the caller graph for this function: