Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Functions
G4INCL::KinematicsUtils Namespace Reference

Functions

G4double fiveParFit (const G4double a, const G4double b, const G4double c, const G4double d, const G4double e, const G4double x)
 
G4double compute_xs (const std::vector< G4double > coefficients, const G4double pLab)
 
void transformToLocalEnergyFrame (Nucleus const *const n, Particle *const p)
 
G4double getLocalEnergy (Nucleus const *const n, Particle *const p)
 
ThreeVector makeBoostVector (Particle const *const p1, Particle const *const p2)
 
G4double totalEnergyInCM (Particle const *const p1, Particle const *const p2)
 
G4double squareTotalEnergyInCM (Particle const *const p1, Particle const *const p2)
 
G4double momentumInCM (Particle const *const p1, Particle const *const p2)
 gives the momentum in the CM frame of two particles.
 
G4double momentumInCM (const G4double E, const G4double M1, const G4double M2)
 
G4double momentumInLab (Particle const *const p1, Particle const *const p2)
 gives the momentum in the lab frame of two particles.
 
G4double momentumInLab (const G4double s, const G4double m1, const G4double m2)
 
G4double sumTotalEnergies (const ParticleList &)
 
ThreeVector sumMomenta (const ParticleList &)
 
G4double energy (const ThreeVector &p, const G4double m)
 
G4double invariantMass (const G4double E, const ThreeVector &p)
 
G4double squareInvariantMass (const G4double E, const ThreeVector &p)
 
G4double gammaFromKineticEnergy (const ParticleSpecies &p, const G4double EKin)
 

Function Documentation

◆ compute_xs()

G4double G4INCL::KinematicsUtils::compute_xs ( const std::vector< G4double > coefficients,
const G4double pLab )

Definition at line 49 of file G4INCLKinematicsUtils.cc.

49 {
50 G4double sigma = 0.;
51 G4double Ethreshold = 0.0;
52 if(coefficients.size() == 6){
53 Ethreshold = coefficients[5];
54 if(Ethreshold >= 5){ //there are no Ethreshold even close to 5 GeV.
55 if(pLab > Ethreshold){ // E is E cutoff, not threshold, we use it when sigma should be zero.
56 return 0.;
57 }
58 }
59 else{
60 if(pLab < Ethreshold){
61 return 0.;
62 }
63 }
64 }
65
66 sigma = fiveParFit(coefficients[0],coefficients[1],coefficients[2],coefficients[3],coefficients[4], pLab);
67 if(sigma < 0.){
68 return 0.;
69 };
70 return sigma;
71 }
G4double
double G4double
Definition G4Types.hh:83
G4INCL::KinematicsUtils::fiveParFit
G4double fiveParFit(const G4double a, const G4double b, const G4double c, const G4double d, const G4double e, const G4double x)
Definition G4INCLKinematicsUtils.cc:45

Referenced by G4INCL::NNbarToAnnihilationChannel::fillFinalState(), G4INCL::NNbarToLLbarChannel::fillFinalState(), G4INCL::NNbarToNNbar2piChannel::fillFinalState(), G4INCL::NNbarToNNbar3piChannel::fillFinalState(), G4INCL::NNbarToNNbarpiChannel::fillFinalState(), G4INCL::CrossSectionsAntiparticles::NNbarCEX(), G4INCL::CrossSectionsAntiparticles::NNbarElastic(), G4INCL::CrossSectionsAntiparticles::NNbarToAnnihilation(), G4INCL::CrossSectionsAntiparticles::NNbarToLLbar(), G4INCL::CrossSectionsAntiparticles::NNbarToNNbar2pi(), G4INCL::CrossSectionsAntiparticles::NNbarToNNbar3pi(), and G4INCL::CrossSectionsAntiparticles::NNbarToNNbarpi().

◆ energy()

G4double G4INCL::KinematicsUtils::energy ( const ThreeVector & p,
const G4double m )

Definition at line 185 of file G4INCLKinematicsUtils.cc.

185 {
186 return std::sqrt(p.mag2() + m*m);
187 }
G4INCL::ThreeVector::mag2
G4double mag2() const
Definition G4INCLThreeVector.hh:78

◆ fiveParFit()

G4double G4INCL::KinematicsUtils::fiveParFit ( const G4double a,
const G4double b,
const G4double c,
const G4double d,
const G4double e,
const G4double x )

Definition at line 45 of file G4INCLKinematicsUtils.cc.

45 {
46 return a+b*std::pow(x, c)+d*std::log(x)+e*std::log(x)*std::log(x);
47 }

Referenced by compute_xs().

◆ gammaFromKineticEnergy()

G4double G4INCL::KinematicsUtils::gammaFromKineticEnergy ( const ParticleSpecies & p,
const G4double EKin )

Definition at line 197 of file G4INCLKinematicsUtils.cc.

197 {
198 G4double mass;
199 if(p.theType==Composite)
200 mass = ParticleTable::getTableMass(p.theA, p.theZ, p.theS);
201 else
202 mass = ParticleTable::getTableParticleMass(p.theType);
203 return (1.+EKin/mass);
204 }

◆ getLocalEnergy()

G4double G4INCL::KinematicsUtils::getLocalEnergy ( Nucleus const *const n,
Particle *const p )

Definition at line 81 of file G4INCLKinematicsUtils.cc.

81 {
82// assert(!p->isMeson() && !p->isPhoton() && !p->isAntiNucleon()); // No local energy for mesons //D photons are bad too!
83 G4double vloc = 0.0;
84 const G4double r = p->getPosition().mag();
85 const G4double mass = p->getMass();
86
87 // Local energy is constant outside the surface
88 if(r > n->getUniverseRadius()) {
89 INCL_WARN("Tried to evaluate local energy for a particle outside the maximum radius."
90 << '\n' << p->print() << '\n'
91 << "Maximum radius = " << n->getDensity()->getMaximumRadius() << '\n'
92 << "Universe radius = " << n->getUniverseRadius() << '\n');
93 return 0.0;
94 }
95
96 G4double pfl0 = 0.0;
97 const ParticleType t = p->getType();
98 const G4double kinE = p->getKineticEnergy();
99 if(kinE <= n->getPotential()->getFermiEnergy(t)) {
100 pfl0 = n->getPotential()->getFermiMomentum(p);
101 } else {
102 const G4double tf0 = p->getPotentialEnergy() - n->getPotential()->getSeparationEnergy(p);
103 if(tf0<0.0) return 0.0;
104 pfl0 = std::sqrt(tf0*(tf0 + 2.0*mass));
105 }
106 const G4double pReflection = p->getReflectionMomentum()/pfl0;
107 const G4double reflectionRadius = n->getDensity()->getMaxRFromP(p->getType(), pReflection);
108 const G4double pNominal = p->getMomentum().mag()/pfl0;
109 const G4double nominalReflectionRadius = n->getDensity()->getMaxRFromP(p->getType(), pNominal);
110 const G4double pl = pfl0*n->getDensity()->getMinPFromR(t, r*nominalReflectionRadius/reflectionRadius);
111 vloc = std::sqrt(pl*pl + mass*mass) - mass;
112
113 return vloc;
114 }
INCL_WARN
#define INCL_WARN(x)
Definition G4INCLLogger.hh:232
G4INCL::Particle::getPotentialEnergy
G4double getPotentialEnergy() const
Get the particle potential energy.
Definition G4INCLParticle.hh:902
G4INCL::Particle::getPosition
const G4INCL::ThreeVector & getPosition() const
Definition G4INCLParticle.hh:956
G4INCL::Particle::getKineticEnergy
G4double getKineticEnergy() const
Get the particle kinetic energy.
Definition G4INCLParticle.hh:899
G4INCL::Particle::getReflectionMomentum
G4double getReflectionMomentum() const
Return the reflection momentum.
Definition G4INCLParticle.hh:1137
G4INCL::Particle::getMomentum
const G4INCL::ThreeVector & getMomentum() const
Definition G4INCLParticle.hh:934
G4INCL::Particle::getType
G4INCL::ParticleType getType() const
Definition G4INCLParticle.hh:187
G4INCL::Particle::print
std::string print() const
Definition G4INCLParticle.hh:1093
G4INCL::Particle::getMass
G4double getMass() const
Get the cached particle mass.
Definition G4INCLParticle.hh:534

Referenced by transformToLocalEnergyFrame().

◆ invariantMass()

G4double G4INCL::KinematicsUtils::invariantMass ( const G4double E,
const ThreeVector & p )

Definition at line 189 of file G4INCLKinematicsUtils.cc.

189 {
190 return std::sqrt(squareInvariantMass(E, p));
191 }
G4INCL::KinematicsUtils::squareInvariantMass
G4double squareInvariantMass(const G4double E, const ThreeVector &p)
Definition G4INCLKinematicsUtils.cc:193

◆ makeBoostVector()

ThreeVector G4INCL::KinematicsUtils::makeBoostVector ( Particle const *const p1,
Particle const *const p2 )

Definition at line 116 of file G4INCLKinematicsUtils.cc.

116 {
117 const G4double totalEnergy = p1->getEnergy() + p2->getEnergy();
118 return ((p1->getMomentum() + p2->getMomentum())/totalEnergy);
119 }

Referenced by G4INCL::InteractionAvatar::preInteraction(), and squareTotalEnergyInCM().

◆ momentumInCM() [1/2]

G4double G4INCL::KinematicsUtils::momentumInCM ( const G4double E,
const G4double M1,
const G4double M2 )

Definition at line 146 of file G4INCLKinematicsUtils.cc.

146 {
147 return 0.5*std::sqrt((E*E - std::pow(M1 + M2, 2))
148 *(E*E - std::pow(M1 - M2, 2)))/E;
149 }

◆ momentumInCM() [2/2]

G4double G4INCL::KinematicsUtils::momentumInCM ( Particle const *const p1,
Particle const *const p2 )

gives the momentum in the CM frame of two particles.

The formula is the following:

\[ p_{CM}^2 = \frac{z^2 - m_1^2 m_2^2}{2 z + m_1^2 + m_2^2} \]

where $z$ is the scalar product of the momentum four-vectors:

\[ z = E_1 E_2 - \vec{p}_1\cdot\vec{p}_2 \]

Parameters
p1pointer to particle 1
p2pointer to particle 2
Returns
the absolute value of the momentum of any of the two particles in the CM frame, in MeV/c.

Definition at line 134 of file G4INCLKinematicsUtils.cc.

134 {
135 const G4double m1sq = std::pow(p1->getMass(),2);
136 const G4double m2sq = std::pow(p2->getMass(),2);
137 const G4double z = p1->getEnergy()*p2->getEnergy() - p1->getMomentum().dot(p2->getMomentum());
138 G4double pcm2 = (z*z-m1sq*m2sq)/(2*z+m1sq+m2sq);
139 if(pcm2 < 0.0) {
140 INCL_ERROR("momentumInCM: pcm2 == " << pcm2 << " < 0.0" << '\n');
141 pcm2 = 0.0;
142 }
143 return std::sqrt(pcm2);
144 }
INCL_ERROR
#define INCL_ERROR(x)
Definition G4INCLLogger.hh:225

Referenced by G4INCL::DeltaDecayChannel::computeDecayTime(), G4INCL::Nucleus::decayOutgoingDeltas(), G4INCL::Nucleus::decayOutgoingPionResonances(), G4INCL::Nucleus::decayOutgoingSigmaZero(), G4INCL::CrossSectionsMultiPionsAndResonances::etaNToPiN(), G4INCL::DeltaDecayChannel::fillFinalState(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::NKbToLpiChannel::fillFinalState(), G4INCL::NKbToNKbChannel::fillFinalState(), G4INCL::NKbToSpiChannel::fillFinalState(), G4INCL::NKElasticChannel::fillFinalState(), G4INCL::NKToNKChannel::fillFinalState(), G4INCL::NpiToLKChannel::fillFinalState(), G4INCL::NpiToSKChannel::fillFinalState(), G4INCL::NYElasticChannel::fillFinalState(), G4INCL::PionResonanceDecayChannel::fillFinalState(), G4INCL::RecombinationChannel::fillFinalState(), G4INCL::SigmaZeroDecayChannel::fillFinalState(), G4INCL::StrangeAbsorbtionChannel::fillFinalState(), G4INCL::PhaseSpaceKopylov::generate(), G4INCL::NKbElasticChannel::KaonMomentum(), and G4INCL::CrossSectionsMultiPionsAndResonances::omegaNToPiN().

◆ momentumInLab() [1/2]

G4double G4INCL::KinematicsUtils::momentumInLab ( const G4double s,
const G4double m1,
const G4double m2 )

Definition at line 151 of file G4INCLKinematicsUtils.cc.

151 {
152 const G4double m1sq = m1*m1;
153 const G4double m2sq = m2*m2;
154 G4double plab2 = (s*s-2*s*(m1sq+m2sq)+(m1sq-m2sq)*(m1sq-m2sq))/(4*m2sq);
155 if(plab2 < 0.0) {
156 INCL_ERROR("momentumInLab: plab2 == " << plab2 << " < 0.0; m1sq == " << m1sq << "; m2sq == " << m2sq << "; s == " << s << '\n');
157 plab2 = 0.0;
158 }
159 return std::sqrt(plab2);
160 }

◆ momentumInLab() [2/2]

G4double G4INCL::KinematicsUtils::momentumInLab ( Particle const *const p1,
Particle const *const p2 )

gives the momentum in the lab frame of two particles.

Assumes particle 1 carries all the momentum and particle 2 is at rest.

The formula is the following:

\[ p_{lab}^2 = \frac{s^2 - 2 s (m_1^2 + m_2^2) + {(m_1^2 - m_2^2)}^2}{4 m_2^2} \]

Parameters
p1pointer to particle 1
p2pointer to particle 2
Returns
the absolute value of the momentum of particle 1 in the lab frame, in MeV/c

Definition at line 162 of file G4INCLKinematicsUtils.cc.

162 {
163 const G4double m1 = p1->getMass();
164 const G4double m2 = p2->getMass();
165 const G4double s = squareTotalEnergyInCM(p1, p2);
166 return momentumInLab(s, m1, m2);
167 }
G4INCL::KinematicsUtils::squareTotalEnergyInCM
G4double squareTotalEnergyInCM(Particle const *const p1, Particle const *const p2)
Definition G4INCLKinematicsUtils.cc:125
G4INCL::KinematicsUtils::momentumInLab
G4double momentumInLab(Particle const *const p1, Particle const *const p2)
gives the momentum in the lab frame of two particles.
Definition G4INCLKinematicsUtils.cc:162

Referenced by G4INCL::CrossSectionsINCL46::elasticNNLegacy(), G4INCL::CrossSectionsMultiPionsAndResonances::etaNElastic(), G4INCL::CrossSectionsMultiPionsAndResonances::etaNToPiN(), G4INCL::CrossSectionsMultiPionsAndResonances::etaNToPiPiN(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::ElasticChannel::fillFinalState(), G4INCL::EtaNElasticChannel::fillFinalState(), G4INCL::EtaNToPiNChannel::fillFinalState(), G4INCL::NNbarToAnnihilationChannel::fillFinalState(), G4INCL::NNbarToLLbarChannel::fillFinalState(), G4INCL::NNbarToNNbar2piChannel::fillFinalState(), G4INCL::NNbarToNNbar3piChannel::fillFinalState(), G4INCL::NNbarToNNbarpiChannel::fillFinalState(), G4INCL::NNToMissingStrangenessChannel::fillFinalState(), G4INCL::NpiToMissingStrangenessChannel::fillFinalState(), G4INCL::NKbElasticChannel::KaonMomentum(), G4INCL::NKbToLpiChannel::KaonMomentum(), G4INCL::NKbToNKbChannel::KaonMomentum(), G4INCL::NKbToSpiChannel::KaonMomentum(), G4INCL::NpiToLKChannel::KaonMomentum(), G4INCL::NpiToSKChannel::KaonMomentum(), momentumInLab(), G4INCL::CrossSectionsINCL46::NDeltaToNN(), G4INCL::CrossSectionsStrangeness::NKbelastic(), G4INCL::CrossSectionsStrangeness::NKbToNKb(), G4INCL::CrossSectionsStrangeness::NKbToNKb2pi(), G4INCL::CrossSectionsStrangeness::NKbToNKbpi(), G4INCL::CrossSectionsStrangeness::NKbToS2pi(), G4INCL::CrossSectionsStrangeness::NKbToSpi(), G4INCL::CrossSectionsStrangeness::NKelastic(), G4INCL::CrossSectionsStrangeness::NKToNK(), G4INCL::CrossSectionsStrangeness::NKToNK2pi(), G4INCL::CrossSectionsStrangeness::NKToNKpi(), G4INCL::CrossSectionsStrangeness::NLToNS(), G4INCL::CrossSectionsAntiparticles::NNbarCEX(), G4INCL::CrossSectionsAntiparticles::NNbarElastic(), G4INCL::CrossSectionsAntiparticles::NNbarToAnnihilation(), G4INCL::CrossSectionsAntiparticles::NNbarToLLbar(), G4INCL::CrossSectionsAntiparticles::NNbarToNNbar2pi(), G4INCL::CrossSectionsAntiparticles::NNbarToNNbar3pi(), G4INCL::CrossSectionsAntiparticles::NNbarToNNbarpi(), G4INCL::CrossSectionsMultiPions::NNElastic(), G4INCL::CrossSectionsMultiPions::NNElasticFixed(), G4INCL::CrossSectionsMultiPions::NNOnePiOrDelta(), G4INCL::CrossSectionsMultiPions::NNThreePi(), G4INCL::CrossSectionsStrangeness::NNToMissingStrangeness(), G4INCL::CrossSectionsINCL46::NNToNDelta(), G4INCL::CrossSectionsStrangeness::NNToNLK(), G4INCL::CrossSectionsStrangeness::NNToNSK(), G4INCL::CrossSectionsMultiPions::NNTotFixed(), G4INCL::CrossSectionsMultiPions::NNTwoPi(), G4INCL::CrossSectionsStrangeness::NpiToLK2pi(), G4INCL::CrossSectionsStrangeness::NpiToLKpi(), G4INCL::CrossSectionsStrangeness::NpiToMissingStrangeness(), G4INCL::CrossSectionsStrangeness::NpiToNKKb(), G4INCL::CrossSectionsStrangeness::NpiToSK2pi(), G4INCL::CrossSectionsStrangeness::NpiToSKpi(), G4INCL::CrossSectionsStrangeness::NSToNL(), G4INCL::CrossSectionsStrangeness::NSToNS(), G4INCL::CrossSectionsStrangeness::NYelastic(), G4INCL::CrossSectionsMultiPionsAndResonances::omegaNElastic(), G4INCL::CrossSectionsMultiPionsAndResonances::omegaNInelastic(), G4INCL::CrossSectionsMultiPionsAndResonances::omegaNToPiN(), G4INCL::CrossSectionsStrangeness::p_kmToL_pp_pm(), G4INCL::CrossSectionsStrangeness::p_kmToL_pz(), G4INCL::CrossSectionsStrangeness::p_pimToLK0(), G4INCL::CrossSectionsStrangeness::p_pimToSmKp(), G4INCL::CrossSectionsStrangeness::p_pimToSzKz(), G4INCL::CrossSectionsStrangeness::p_pipToSpKp(), G4INCL::CrossSectionsStrangeness::p_pizToSzKp(), G4INCL::CrossSectionsMultiPions::piMinuspIne(), G4INCL::CrossSectionsMultiPions::piMinuspOnePi(), G4INCL::CrossSectionsMultiPionsAndResonances::piMinuspToEtaN(), G4INCL::CrossSectionsMultiPionsAndResonances::piMinuspToEtaN(), G4INCL::CrossSectionsMultiPionsAndResonances::piMinuspToOmegaN(), G4INCL::CrossSectionsMultiPionsAndResonances::piMinuspToOmegaN(), G4INCL::CrossSectionsMultiPions::piMinuspTwoPi(), G4INCL::CrossSectionsMultiPions::piNIne(), G4INCL::CrossSectionsMultiPions::piNOnePi(), G4INCL::CrossSectionsMultiPions::piNToxPiN(), G4INCL::CrossSectionsMultiPions::piNTwoPi(), G4INCL::CrossSectionsMultiPions::piPluspIne(), G4INCL::CrossSectionsMultiPions::piPluspOnePi(), and G4INCL::CrossSectionsMultiPions::piPluspTwoPi().

◆ squareInvariantMass()

G4double G4INCL::KinematicsUtils::squareInvariantMass ( const G4double E,
const ThreeVector & p )

Definition at line 193 of file G4INCLKinematicsUtils.cc.

193 {
194 return E*E - p.mag2();
195 }

Referenced by invariantMass().

◆ squareTotalEnergyInCM()

G4double G4INCL::KinematicsUtils::squareTotalEnergyInCM ( Particle const *const p1,
Particle const *const p2 )

Definition at line 125 of file G4INCLKinematicsUtils.cc.

125 {
126 G4double beta2 = makeBoostVector(p1, p2).mag2();
127 if(beta2 > 1.0) {
128 INCL_ERROR("squareTotalEnergyInCM: beta2 == " << beta2 << " > 1.0" << '\n');
129 beta2 = 0.0;
130 }
131 return (1.0 - beta2)*std::pow(p1->getEnergy() + p2->getEnergy(), 2);
132 }
G4INCL::KinematicsUtils::makeBoostVector
ThreeVector makeBoostVector(Particle const *const p1, Particle const *const p2)
Definition G4INCLKinematicsUtils.cc:116

Referenced by G4INCL::CrossSectionsINCL46::elasticNNLegacy(), G4INCL::ElasticChannel::fillFinalState(), G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::BinaryCollisionAvatar::getChannel(), momentumInLab(), G4INCL::CrossSectionsStrangeness::NDeltaToDeltaLK(), G4INCL::CrossSectionsStrangeness::NDeltaToDeltaSK(), G4INCL::CrossSectionsStrangeness::NDeltaToNLK(), G4INCL::CrossSectionsINCL46::NDeltaToNN(), G4INCL::CrossSectionsMultiPions::NDeltaToNN(), G4INCL::CrossSectionsStrangeness::NDeltaToNSK(), G4INCL::CrossSectionsMultiPions::NNElastic(), G4INCL::CrossSectionsMultiPions::NNFourPi(), G4INCL::CrossSectionsMultiPions::NNTot(), and totalEnergyInCM().

◆ sumMomenta()

ThreeVector G4INCL::KinematicsUtils::sumMomenta ( const ParticleList & pl)

Definition at line 177 of file G4INCLKinematicsUtils.cc.

177 {
178 ThreeVector p(0.0, 0.0, 0.0);
179 for(ParticleIter i=pl.begin(), e=pl.end(); i!=e; ++i) {
180 p += (*i)->getMomentum();
181 }
182 return p;
183 }
G4INCL::ParticleIter
ParticleList::const_iterator ParticleIter
Definition G4INCLParticle.hh:72

◆ sumTotalEnergies()

G4double G4INCL::KinematicsUtils::sumTotalEnergies ( const ParticleList & pl)

Definition at line 169 of file G4INCLKinematicsUtils.cc.

169 {
170 G4double E = 0.0;
171 for(ParticleIter i=pl.begin(), e=pl.end(); i!=e; ++i) {
172 E += (*i)->getEnergy();
173 }
174 return E;
175 }

◆ totalEnergyInCM()

G4double G4INCL::KinematicsUtils::totalEnergyInCM ( Particle const *const p1,
Particle const *const p2 )

Definition at line 121 of file G4INCLKinematicsUtils.cc.

121 {
122 return std::sqrt(squareTotalEnergyInCM(p1,p2));
123 }

Referenced by G4INCL::CrossSectionsMultiPionsAndResonances::etaNToPiN(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::EtaNToPiPiNChannel::fillFinalState(), G4INCL::NDeltaEtaProductionChannel::fillFinalState(), G4INCL::NDeltaOmegaProductionChannel::fillFinalState(), G4INCL::NDeltaToDeltaLKChannel::fillFinalState(), G4INCL::NDeltaToDeltaSKChannel::fillFinalState(), G4INCL::NDeltaToNLKChannel::fillFinalState(), G4INCL::NDeltaToNNKKbChannel::fillFinalState(), G4INCL::NDeltaToNSKChannel::fillFinalState(), G4INCL::NKbToL2piChannel::fillFinalState(), G4INCL::NKbToNKb2piChannel::fillFinalState(), G4INCL::NKbToNKbpiChannel::fillFinalState(), G4INCL::NKbToS2piChannel::fillFinalState(), G4INCL::NKToNK2piChannel::fillFinalState(), G4INCL::NKToNKpiChannel::fillFinalState(), G4INCL::NLToNSChannel::fillFinalState(), G4INCL::NNbarCEXChannel::fillFinalState(), G4INCL::NNbarElasticChannel::fillFinalState(), G4INCL::NNbarToAnnihilationChannel::fillFinalState(), G4INCL::NNbarToLLbarChannel::fillFinalState(), G4INCL::NNbarToNNbar2piChannel::fillFinalState(), G4INCL::NNbarToNNbar3piChannel::fillFinalState(), G4INCL::NNbarToNNbarpiChannel::fillFinalState(), G4INCL::NNEtaToMultiPionsChannel::fillFinalState(), G4INCL::NNOmegaToMultiPionsChannel::fillFinalState(), G4INCL::NNToMissingStrangenessChannel::fillFinalState(), G4INCL::NNToMultiPionsChannel::fillFinalState(), G4INCL::NNToNLK2piChannel::fillFinalState(), G4INCL::NNToNLKChannel::fillFinalState(), G4INCL::NNToNLKpiChannel::fillFinalState(), G4INCL::NNToNNEtaChannel::fillFinalState(), G4INCL::NNToNNKKbChannel::fillFinalState(), G4INCL::NNToNNOmegaChannel::fillFinalState(), G4INCL::NNToNSK2piChannel::fillFinalState(), G4INCL::NNToNSKChannel::fillFinalState(), G4INCL::NNToNSKpiChannel::fillFinalState(), G4INCL::NpiToLK2piChannel::fillFinalState(), G4INCL::NpiToLKpiChannel::fillFinalState(), G4INCL::NpiToMissingStrangenessChannel::fillFinalState(), G4INCL::NpiToNKKbChannel::fillFinalState(), G4INCL::NpiToSK2piChannel::fillFinalState(), G4INCL::NpiToSKpiChannel::fillFinalState(), G4INCL::NSToNLChannel::fillFinalState(), G4INCL::NSToNSChannel::fillFinalState(), G4INCL::OmegaNToPiPiNChannel::fillFinalState(), G4INCL::PiNToEtaChannel::fillFinalState(), G4INCL::PiNToMultiPionsChannel::fillFinalState(), G4INCL::RecombinationChannel::fillFinalState(), G4INCL::StrangeAbsorbtionChannel::fillFinalState(), G4INCL::CrossSectionsStrangeness::NDeltaToNNKKb(), G4INCL::CrossSectionsMultiPions::NNOnePi(), G4INCL::CrossSectionsMultiPions::NNOnePiOrDelta(), G4INCL::CrossSectionsMultiPions::NNThreePi(), G4INCL::CrossSectionsINCL46::NNToNDelta(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNDeltaEta(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNDeltaOmega(), G4INCL::CrossSectionsStrangeness::NNToNLK2pi(), G4INCL::CrossSectionsStrangeness::NNToNLKpi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEta(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEtaExclu(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEtaFourPi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEtaOnePi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEtaOnePiOrDelta(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEtaThreePi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEtaTwoPi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNEtaxPi(), G4INCL::CrossSectionsStrangeness::NNToNNKKb(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmega(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmegaExclu(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmegaFourPi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmegaOnePi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmegaOnePiOrDelta(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmegaThreePi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmegaTwoPi(), G4INCL::CrossSectionsMultiPionsAndResonances::NNToNNOmegaxPi(), G4INCL::CrossSectionsStrangeness::NNToNSK2pi(), G4INCL::CrossSectionsStrangeness::NNToNSKpi(), G4INCL::CrossSectionsMultiPions::NNTwoPi(), G4INCL::CrossSectionsMultiPionsAndResonances::omegaNToPiN(), G4INCL::CrossSectionsMultiPionsAndResonances::piMinuspToEtaN(), G4INCL::CrossSectionsMultiPionsAndResonances::piMinuspToOmegaN(), G4INCL::CrossSectionsINCL46::piNToDelta(), G4INCL::CrossSectionsMultiPions::piNToDelta(), and G4INCL::CrossSectionsMultiPions::piNTot().

◆ transformToLocalEnergyFrame()

void G4INCL::KinematicsUtils::transformToLocalEnergyFrame ( Nucleus const *const n,
Particle *const p )

Definition at line 73 of file G4INCLKinematicsUtils.cc.

73 {
74// assert(!p->isMeson() && !p->isPhoton() && !p->isAntiNucleon()); // No local energy for mesons //D nor for photons!
75 const G4double localEnergy = getLocalEnergy(n, p);
76 const G4double localTotalEnergy = p->getEnergy() - localEnergy;
77 p->setEnergy(localTotalEnergy);
78 p->adjustMomentumFromEnergy();
79 }
G4INCL::Particle::getEnergy
G4double getEnergy() const
Definition G4INCLParticle.hh:910
G4INCL::Particle::adjustMomentumFromEnergy
const ThreeVector & adjustMomentumFromEnergy()
Rescale the momentum to match the total energy.
Definition G4INCLParticle.cc:154
G4INCL::Particle::setEnergy
void setEnergy(G4double energy)
Definition G4INCLParticle.hh:926
G4INCL::KinematicsUtils::getLocalEnergy
G4double getLocalEnergy(Nucleus const *const n, Particle *const p)
Definition G4INCLKinematicsUtils.cc:81

Referenced by G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), and G4INCL::InteractionAvatar::preInteractionLocalEnergy().