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	pionNucleon (Particle const *const p1, Particle const *const p2)
G4double	recombination (Particle const *const p1, Particle const *const p2)
G4double	deltaProduction (Particle const *const p1, Particle const *const p2)
G4double	calculateNNAngularSlope (G4double energyCM, G4int iso)
	Calculate the slope of the NN DDXS. More...
G4double	interactionDistanceNN (const ParticleSpecies &aSpecies, const G4double kineticEnergy)
	Compute the "interaction distance". More...
G4double	interactionDistancePiN (const G4double projectileKineticEnergy)
	Compute the "interaction distance". More...
void	setCrossSections (ICrossSections *c)
void	deleteCrossSections ()

Function Documentation

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

Calculate the slope of the NN DDXS.

Parameters

energyCM	energy in the CM frame, in MeV
iso	total isospin of the system

Returns

the slope of the angular distribution

Definition at line 70 of file G4INCLCrossSections.cc.

Referenced by G4INCL::ElasticChannel::getFinalState().

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

void G4INCL::CrossSections::deleteCrossSections

(

)

Definition at line 142 of file G4INCLCrossSections.cc.

Referenced by G4INCL::INCL::~INCL().

                               {
      delete theCrossSections;
      theCrossSections = NULL;
    }

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

Definition at line 66 of file G4INCLCrossSections.cc.

Referenced by G4INCL::BinaryCollisionAvatar::getChannel().

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

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

Definition at line 50 of file G4INCLCrossSections.cc.

Referenced by G4QMDReaction::ApplyYourself(), G4QMDCollision::CalFinalStateOfTheBinaryCollisionJQMD(), and G4INCL::BinaryCollisionAvatar::getChannel().

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

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

Compute the "interaction distance".

Defined on the basis of the average value of the N-N cross sections at the given kinetic energy.

Returns

the interaction distance

Definition at line 107 of file G4INCLCrossSections.cc.

References G4INCL::Particle::adjustMomentumFromEnergy(), G4INCL::Particle::getMass(), G4INCL::Math::max(), G4INCL::Neutron, G4INCL::Proton, G4INCL::Particle::setEnergy(), G4INCL::Math::tenPi, G4INCL::ParticleSpecies::theA, and total().

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

G4double G4INCL::CrossSections::interactionDistancePiN

(

const G4double

projectileKineticEnergy

)

Compute the "interaction distance".

Defined on the basis of the average value of the pi-N cross sections at the given kinetic energy.

Returns

the interaction distance

Definition at line 74 of file G4INCLCrossSections.cc.

References G4INCL::Particle::adjustMomentumFromEnergy(), G4INCL::Particle::getMass(), G4INCL::Math::max(), G4INCL::Neutron, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, G4INCL::Particle::setEnergy(), G4INCL::Math::tenPi, and total().

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

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

Definition at line 58 of file G4INCLCrossSections.cc.

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

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

Definition at line 62 of file G4INCLCrossSections.cc.

Referenced by G4INCL::BinaryCollisionAvatar::getChannel().

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

void G4INCL::CrossSections::setCrossSections

(

ICrossSections *

c

)

Definition at line 138 of file G4INCLCrossSections.cc.

References test::c.

Referenced by G4INCL::INCL::INCL().

                                             {
      theCrossSections = c;
    }

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

Definition at line 54 of file G4INCLCrossSections.cc.

Referenced by G4ScreeningMottCrossSection::AngleDistribution(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4Volant::getTotalMass(), G4NeutronHPFissionBaseFS::Init(), G4NeutronHPList::Init(), G4NeutronHPInelasticBaseFS::Init(), G4NeutronHPInelasticCompFS::Init(), G4NeutronHPVector::Init(), G4NeutronHPVector::IntegrateAndNormalise(), interactionDistanceNN(), interactionDistancePiN(), G4SPBaryon::MatchDiQuarkAndGetQuark(), G4INCL::InteractionAvatar::preInteractionBlocking(), and CexmcRunAction::PrintResults().

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