IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 897bdc97 authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by ralfulrich
Browse files

updated modules usage of get_mass()

parent 7a3f337a
No related branches found
No related tags found
No related merge requests found
...@@ -100,7 +100,7 @@ namespace corsika::pythia8 { ...@@ -100,7 +100,7 @@ namespace corsika::pythia8 {
double py = pcomp[1] / 1_GeV; double py = pcomp[1] / 1_GeV;
double pz = pcomp[2] / 1_GeV; double pz = pcomp[2] / 1_GeV;
double en = vP.GetEnergy() / 1_GeV; double en = vP.GetEnergy() / 1_GeV;
double m = corsika::mass(vP.GetPID()) / 1_GeV; double m = corsika::get_mass(vP.GetPID()) / 1_GeV;
// add particle to pythia stack // add particle to pythia stack
event.append(pdgCode, 1, 0, 0, px, py, pz, en, m); event.append(pdgCode, 1, 0, 0, px, py, pz, en, m);
......
...@@ -295,7 +295,7 @@ namespace corsika::qgsjetII { ...@@ -295,7 +295,7 @@ namespace corsika::qgsjetII {
(projectileEnergyLab - corsika::Proton::mass))}); (projectileEnergyLab - corsika::Proton::mass))});
auto const energy = auto const energy =
sqrt(momentum.squaredNorm() + square(corsika::mass(idFragm))); sqrt(momentum.squaredNorm() + square(corsika::get_mass(idFragm)));
momentum.rebase(originalCS); // transform back into standard lab frame momentum.rebase(originalCS); // transform back into standard lab frame
std::cout << "secondary fragment> id=" << idFragm std::cout << "secondary fragment> id=" << idFragm
<< " p=" << momentum.GetComponents() << std::endl; << " p=" << momentum.GetComponents() << std::endl;
......
...@@ -123,7 +123,7 @@ namespace corsika::sibyll { ...@@ -123,7 +123,7 @@ namespace corsika::sibyll {
vP.GetMomentum(), vP.GetMomentum(),
// setting particle mass with Corsika values, may be inconsistent // setting particle mass with Corsika values, may be inconsistent
// with sibyll internal values // with sibyll internal values
corsika::mass(pCode)); corsika::get_mass(pCode));
// remember position // remember position
Point const decayPoint = vP.GetPosition(); Point const decayPoint = vP.GetPosition();
TimeType const t0 = vP.GetTime(); TimeType const t0 = vP.GetTime();
......
...@@ -560,7 +560,7 @@ namespace corsika::sibyll { ...@@ -560,7 +560,7 @@ namespace corsika::sibyll {
// CORSIKA 7 way // CORSIKA 7 way
// elastic nucleons inherit momentum from original projectile // elastic nucleons inherit momentum from original projectile
// neglecting momentum transfer in interaction // neglecting momentum transfer in interaction
const double mass_ratio = corsika::mass(elaNucCode) / ProjMass; const double mass_ratio = corsika::get_mass(elaNucCode) / ProjMass;
auto const Plab = PprojLab * mass_ratio; auto const Plab = PprojLab * mass_ratio;
vP.AddSecondary(std::tuple<corsika::Code, si::HEPEnergyType, vP.AddSecondary(std::tuple<corsika::Code, si::HEPEnergyType,
......
...@@ -39,8 +39,8 @@ namespace corsika::urqmd { ...@@ -39,8 +39,8 @@ namespace corsika::urqmd {
// the following is a translation of ptsigtot() into C++ // the following is a translation of ptsigtot() into C++
if (vProjectileCode != corsika::Code::Nucleus && if (vProjectileCode != corsika::Code::Nucleus &&
!corsika::is_nucleus(vTargetCode)) { // both particles are "special" !corsika::is_nucleus(vTargetCode)) { // both particles are "special"
auto const mProj = corsika::mass(vProjectileCode); auto const mProj = corsika::get_mass(vProjectileCode);
auto const mTar = corsika::mass(vTargetCode); auto const mTar = corsika::get_mass(vTargetCode);
double sqrtS = double sqrtS =
sqrt(static_pow<2>(mProj) + static_pow<2>(mTar) + 2 * vLabEnergy * mTar) * sqrt(static_pow<2>(mProj) + static_pow<2>(mTar) + 2 * vLabEnergy * mTar) *
(1 / 1_GeV); (1 / 1_GeV);
...@@ -226,7 +226,7 @@ namespace corsika::urqmd { ...@@ -226,7 +226,7 @@ namespace corsika::urqmd {
::urqmd::coor_.px[i], ::urqmd::coor_.py[i], ::urqmd::coor_.pz[i]} * ::urqmd::coor_.px[i], ::urqmd::coor_.py[i], ::urqmd::coor_.pz[i]} *
1_GeV); 1_GeV);
auto const energy = sqrt(momentum.squaredNorm() + square(corsika::mass(code))); auto const energy = sqrt(momentum.squaredNorm() + square(corsika::get_mass(code)));
momentum.rebase(originalCS); // transform back into standard lab frame momentum.rebase(originalCS); // transform back into standard lab frame
std::cout << i << " " << code << " " << momentum.GetComponents() << std::endl; std::cout << i << " " << code << " " << momentum.GetComponents() << std::endl;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment