IAP GITLAB

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

some processes

parent 0a262e3d
No related branches found
No related tags found
1 merge request!280Refactory 2020
...@@ -35,7 +35,7 @@ namespace corsika::observation_plane { ...@@ -35,7 +35,7 @@ namespace corsika::observation_plane {
return corsika::EProcessReturn::eOk; return corsika::EProcessReturn::eOk;
} }
outputStream_ << static_cast<int>(corsika::GetPDG(particle.GetPID())) << ' ' outputStream_ << static_cast<int>(corsika::PDG(particle.GetPID())) << ' '
<< particle.GetEnergy() * (1 / 1_eV) << ' ' << particle.GetEnergy() * (1 / 1_eV) << ' '
<< (trajectory.GetPosition(1) - plane_.GetCenter()).norm() / 1_m << (trajectory.GetPosition(1) - plane_.GetCenter()).norm() / 1_m
<< std::endl; << std::endl;
...@@ -63,4 +63,4 @@ namespace corsika::observation_plane { ...@@ -63,4 +63,4 @@ namespace corsika::observation_plane {
return (trajectory.GetR0() - pointOfIntersection).norm() * 1.0001; return (trajectory.GetR0() - pointOfIntersection).norm() * 1.0001;
} }
} // namespace corsika::observation_plane } // namespace corsika::observation_plane
\ No newline at end of file
...@@ -35,7 +35,7 @@ namespace corsika::track_writer { ...@@ -35,7 +35,7 @@ namespace corsika::track_writer {
const TTrack& vT) { const TTrack& vT) {
auto const start = vT.GetPosition(0).GetCoordinates(); auto const start = vT.GetPosition(0).GetCoordinates();
auto const delta = vT.GetPosition(1).GetCoordinates() - start; auto const delta = vT.GetPosition(1).GetCoordinates() - start;
auto const pdg = static_cast<int>(corsika::GetPDG(vP.GetPID())); auto const pdg = static_cast<int>(corsika::PDG(vP.GetPID()));
// clang-format off // clang-format off
fFile << std::setw(7) << pdg fFile << std::setw(7) << pdg
......
...@@ -70,7 +70,7 @@ namespace corsika::energy_loss { ...@@ -70,7 +70,7 @@ namespace corsika::energy_loss {
double const gamma = E / m; double const gamma = E / m;
int const Z = p.GetChargeNumber(); int const Z = p.GetChargeNumber();
int const Z2 = Z * Z; int const Z2 = Z * Z;
HEPMassType constexpr me = corsika::Electron::GetMass(); HEPMassType constexpr me = corsika::Electron::mass();
auto const m2 = m * m; auto const m2 = m * m;
auto constexpr me2 = me * me; auto constexpr me2 = me * me;
double const gamma2 = gamma * gamma; double const gamma2 = gamma * gamma;
......
...@@ -47,7 +47,7 @@ TEST_CASE("ContinuousProcess interface", "[proccesses][observation_plane]") { ...@@ -47,7 +47,7 @@ TEST_CASE("ContinuousProcess interface", "[proccesses][observation_plane]") {
std::tuple<Code, HEPEnergyType, corsika::MomentumVector, Point, TimeType>{ std::tuple<Code, HEPEnergyType, corsika::MomentumVector, Point, TimeType>{
Code::NuMu, 1_GeV, Code::NuMu, 1_GeV,
corsika::MomentumVector(rootCS, corsika::MomentumVector(rootCS,
{0_GeV, 0_GeV, -elab2plab(1_GeV, NuMu::GetMass())}), {0_GeV, 0_GeV, -elab2plab(1_GeV, NuMu::mass())}),
Point(rootCS, {1_m, 1_m, 10_m}), 0_ns}); Point(rootCS, {1_m, 1_m, 10_m}), 0_ns});
} }
auto particle = stack.GetNextParticle(); auto particle = stack.GetNextParticle();
......
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