IAP GITLAB

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

update modules usage of get_PDG()

parent 0609f9e2
No related branches found
No related tags found
No related merge requests found
......@@ -35,7 +35,7 @@ namespace corsika::observation_plane {
return corsika::EProcessReturn::eOk;
}
outputStream_ << static_cast<int>(corsika::PDG(particle.GetPID())) << ' '
outputStream_ << static_cast<int>(corsika::get_PDG(particle.GetPID())) << ' '
<< particle.GetEnergy() * (1 / 1_eV) << ' '
<< (trajectory.GetPosition(1) - plane_.GetCenter()).norm() / 1_m
<< std::endl;
......
......@@ -35,7 +35,7 @@ namespace corsika::track_writer {
const TTrack& vT) {
auto const start = vT.GetPosition(0).GetCoordinates();
auto const delta = vT.GetPosition(1).GetCoordinates() - start;
auto const pdg = static_cast<int>(corsika::PDG(vP.GetPID()));
auto const pdg = static_cast<int>(corsika::get_PDG(vP.GetPID()));
// clang-format off
fFile << std::setw(7) << pdg
......
......@@ -50,12 +50,12 @@ namespace corsika::pythia8 {
void Decay::SetUnstable(const corsika::Code pCode) {
std::cout << "Pythia::Decay: setting " << pCode << " unstable.." << std::endl;
fPythia.particleData.mayDecay(static_cast<int>(corsika::PDG(pCode)), true);
fPythia.particleData.mayDecay(static_cast<int>(corsika::get_PDG(pCode)), true);
}
void Decay::SetStable(const corsika::Code pCode) {
std::cout << "Pythia::Decay: setting " << pCode << " stable.." << std::endl;
fPythia.particleData.mayDecay(static_cast<int>(corsika::PDG(pCode)), false);
fPythia.particleData.mayDecay(static_cast<int>(corsika::get_PDG(pCode)), false);
}
template <typename TParticle>
......@@ -93,7 +93,7 @@ namespace corsika::pythia8 {
Decay::SetUnstable(vP.GetPID());
// input particle PDG
auto const pdgCode = static_cast<int>(corsika::PDG(vP.GetPID()));
auto const pdgCode = static_cast<int>(corsika::get_PDG(vP.GetPID()));
auto const pcomp = vP.GetMomentum().GetComponents();
double px = pcomp[0] / 1_GeV;
......
......@@ -77,12 +77,12 @@ namespace corsika::pythia8 {
void Interaction::SetUnstable(const corsika::Code pCode) {
std::cout << "Pythia::Interaction: setting " << pCode << " unstable.." << std::endl;
fPythia.particleData.mayDecay(static_cast<int>(corsika::PDG(pCode)), true);
fPythia.particleData.mayDecay(static_cast<int>(corsika::get_PDG(pCode)), true);
}
void Interaction::SetStable(const corsika::Code pCode) {
std::cout << "Pythia::Interaction: setting " << pCode << " stable.." << std::endl;
fPythia.particleData.mayDecay(static_cast<int>(corsika::PDG(pCode)), false);
fPythia.particleData.mayDecay(static_cast<int>(corsika::get_PDG(pCode)), false);
}
void Interaction::ConfigureLabFrameCollision(const corsika::Code BeamId,
......@@ -93,15 +93,15 @@ namespace corsika::pythia8 {
// set beam
// beam id for pythia
auto const pdgBeam = static_cast<int>(corsika::PDG(BeamId));
auto const pdgBeam = static_cast<int>(corsika::get_PDG(BeamId));
std::stringstream stBeam;
stBeam << "Beams:idA = " << pdgBeam;
fPythia.readString(stBeam.str());
// set target
auto pdgTarget = static_cast<int>(corsika::PDG(TargetId));
auto pdgTarget = static_cast<int>(corsika::get_PDG(TargetId));
// replace hydrogen with proton, otherwise pythia goes into heavy ion mode!
if (TargetId == corsika::Code::Hydrogen)
pdgTarget = static_cast<int>(corsika::PDG(corsika::Code::Proton));
pdgTarget = static_cast<int>(corsika::get_PDG(corsika::Code::Proton));
std::stringstream stTarget;
stTarget << "Beams:idB = " << pdgTarget;
fPythia.readString(stTarget.str());
......@@ -131,8 +131,8 @@ namespace corsika::pythia8 {
if (TargetId == corsika::Code::Proton || TargetId == corsika::Code::Hydrogen) {
if (CanInteract(BeamId) && ValidCoMEnergy(CoMenergy)) {
// input particle PDG
auto const pdgCodeBeam = static_cast<int>(corsika::PDG(BeamId));
auto const pdgCodeTarget = static_cast<int>(corsika::PDG(TargetId));
auto const pdgCodeBeam = static_cast<int>(corsika::get_PDG(BeamId));
auto const pdgCodeTarget = static_cast<int>(corsika::get_PDG(TargetId));
const double ecm = CoMenergy / 1_GeV;
// calculate cross section
......
......@@ -310,7 +310,7 @@ namespace corsika::urqmd {
{443, {135, 0}}, // jpsi
};
return mapPDGToUrQMD.at(static_cast<int>(PDG(code)));
return mapPDGToUrQMD.at(static_cast<int>(get_PDG(code)));
}
} // namespace corsika::urqmd
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