IAP GITLAB

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

comment

parent 8e630180
No related branches found
No related tags found
2 merge requests!234WIP: Initial example of python as script language from C++,!232Conex EM
...@@ -89,7 +89,7 @@ void CONEXSourceCut::SolveCE() { conex::conexcascade_(); } ...@@ -89,7 +89,7 @@ void CONEXSourceCut::SolveCE() { conex::conexcascade_(); }
CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis showerAxis, CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis showerAxis,
units::si::LengthType groundDist, units::si::LengthType groundDist,
//units::si::GrammageType Xcut, // units::si::GrammageType Xcut,
units::si::HEPEnergyType primaryEnergy, units::si::HEPEnergyType primaryEnergy,
particles::PDGCode primaryID) particles::PDGCode primaryID)
: center_{center} : center_{center}
...@@ -114,6 +114,8 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s ...@@ -114,6 +114,8 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s
})} })}
, y_sf_{showerAxis_.GetDirection().cross(x_sf_)} { , y_sf_{showerAxis_.GetDirection().cross(x_sf_)} {
// TODO: call initconex
double eprima = primaryEnergy / 1_GeV; double eprima = primaryEnergy / 1_GeV;
// set phi, theta // set phi, theta
...@@ -125,7 +127,7 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s ...@@ -125,7 +127,7 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s
showerAxis_.GetDirection().GetComponents(conexObservationCS_); showerAxis_.GetDirection().GetComponents(conexObservationCS_);
double phi = 180 * std::atan2(-showerAxisConex.GetY().magnitude(), double phi = 180 * std::atan2(-showerAxisConex.GetY().magnitude(),
showerAxisConex.GetX().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); int ipart = static_cast<int>(primaryID);
auto rng = corsika::random::RNGManager::GetInstance().GetRandomStream("cascade"); auto rng = corsika::random::RNGManager::GetInstance().GetRandomStream("cascade");
......
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