diff --git a/Processes/CONEXSourceCut/CONEXSourceCut.cc b/Processes/CONEXSourceCut/CONEXSourceCut.cc index a312349085fe00e08c2df5123c8a8e8d65b0d208..e1455a3b2469de7c49dfcf1f26841671f0552d9d 100644 --- a/Processes/CONEXSourceCut/CONEXSourceCut.cc +++ b/Processes/CONEXSourceCut/CONEXSourceCut.cc @@ -89,7 +89,7 @@ void CONEXSourceCut::SolveCE() { conex::conexcascade_(); } CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis showerAxis, units::si::LengthType groundDist, - units::si::GrammageType Xcut, + //units::si::GrammageType Xcut, units::si::HEPEnergyType primaryEnergy, particles::PDGCode primaryID) : center_{center} @@ -125,7 +125,7 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s showerAxis_.GetDirection().GetComponents(conexObservationCS_); double phi = 180 * std::atan2(-showerAxisConex.GetY().magnitude(), showerAxisConex.GetX().magnitude()); - double XmaxP_ = Xcut / (1_g / 1_cm / 1_cm); + //double XmaxP_ = Xcut / (1_g / 1_cm / 1_cm); int ipart = static_cast<int>(primaryID); auto rng = corsika::random::RNGManager::GetInstance().GetRandomStream("cascade"); diff --git a/Processes/CONEXSourceCut/CONEXSourceCut.h b/Processes/CONEXSourceCut/CONEXSourceCut.h index 918f6b1f3156c5fb260cde9aef1373d39a29f824..477261c585c24f01002aceebaf8e82b824beaf72 100644 --- a/Processes/CONEXSourceCut/CONEXSourceCut.h +++ b/Processes/CONEXSourceCut/CONEXSourceCut.h @@ -11,6 +11,8 @@ #ifndef _corsika_process_particle_cut_CONEXSourceCut_h_ #define _corsika_process_particle_cut_CONEXSourceCut_h_ +#include <ConexDynamicInterface.h> + #include <corsika/environment/ShowerAxis.h> #include <corsika/geometry/Point.h> #include <corsika/geometry/Vector.h> @@ -70,7 +72,7 @@ namespace corsika::process { public: CONEXSourceCut(geometry::Point center, environment::ShowerAxis showerAxis, - units::si::LengthType groundDist, units::si::GrammageType Xcut, + units::si::LengthType groundDist, /*units::si::GrammageType Xcut,*/ units::si::HEPEnergyType primaryEnergy, particles::PDGCode primaryID); corsika::process::EProcessReturn DoSecondaries(corsika::setup::StackView&);