diff --git a/Processes/CONEXSourceCut/CONEXSourceCut.cc b/Processes/CONEXSourceCut/CONEXSourceCut.cc
index e1455a3b2469de7c49dfcf1f26841671f0552d9d..094f993d3b7a46fd2192dd17076715c060180f59 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}
@@ -114,6 +114,8 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s
     })}
     , y_sf_{showerAxis_.GetDirection().cross(x_sf_)} {
 
+  // TODO: call initconex
+
   double eprima = primaryEnergy / 1_GeV;
 
   // set phi, theta
@@ -125,7 +127,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");