IAP GITLAB

Skip to content
Snippets Groups Projects
Commit ccb14a00 authored by ralfulrich's avatar ralfulrich
Browse files

fixed some warnings

parent 9216c9b3
No related branches found
No related tags found
2 merge requests!234WIP: Initial example of python as script language from C++,!232Conex EM
Pipeline #1830 passed with warnings
......@@ -67,4 +67,3 @@ namespace corsika::process {
} // namespace conex_source_cut
} // namespace corsika::process
#endif
......@@ -51,23 +51,10 @@ TEST_CASE("CONEXSourceCut") {
builder.assemble(env);
const HEPEnergyType mass = particles::GetMass(particles::Code::Proton);
const HEPEnergyType E0 = 1_PeV;
double thetaDeg = 60.;
auto const thetaRad = thetaDeg / 180. * M_PI;
auto elab2plab = [](HEPEnergyType Elab, HEPMassType m) {
return sqrt((Elab - m) * (Elab + m));
};
HEPMomentumType P0 = elab2plab(E0, mass);
auto momentumComponents = [](double thetaRad, HEPMomentumType ptot) {
return std::make_tuple(ptot * sin(thetaRad), 0_eV, -ptot * cos(thetaRad));
};
auto const [px, py, pz] = momentumComponents(thetaRad, P0);
auto plab = corsika::stack::MomentumVector(rootCS, {px, py, pz});
auto const observationHeight = 1.4_km + conex::earthRadius;
auto const injectionHeight = 112.75_km + conex::earthRadius;
auto const t =
......
......@@ -119,6 +119,8 @@ TEST_CASE("InteractionCounter") {
}
auto [env, csPtr, nodePtr] = setupEnvironment(particles::Code::Oxygen);
[[maybe_unused]] auto& env_dummy = env;
SECTION("DoInteraction nucleus") {
unsigned short constexpr A = 14, Z = 7;
auto [stackPtr, secViewPtr] = setupStack(A, Z, 105_TeV, nodePtr, *csPtr);
......
......@@ -315,8 +315,6 @@ namespace corsika::process::sibyll {
auto const pCoM = Vector<hepmomentum_d>(csPrime, tmp);
HEPEnergyType const eCoM = psib.GetEnergy();
auto const Plab = boost.fromCoM(FourVector(eCoM, pCoM));
auto const p3lab = Plab.GetSpaceLikeComponents();
assert(p3lab.GetCoordinateSystem() == originalCS); // just to be sure!
auto const pid = process::sibyll::ConvertFromSibyll(psib.GetPID());
......@@ -324,8 +322,8 @@ namespace corsika::process::sibyll {
auto pnew = vP.AddSecondary(
tuple<particles::Code, units::si::HEPEnergyType, stack::MomentumVector,
geometry::Point, units::si::TimeType>{
process::sibyll::ConvertFromSibyll(psib.GetPID()),
Plab.GetTimeLikeComponent(), p3lab, pOrig, tOrig});
pid, Plab.GetTimeLikeComponent(), Plab.GetSpaceLikeComponents(), pOrig,
tOrig});
Plab_final += pnew.GetMomentum();
Elab_final += pnew.GetEnergy();
......
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