diff --git a/Processes/CONEX/testConex.cc b/Processes/CONEX/testConex.cc
index 070b460dbf6bd991602abfe953a75dd591543a69..830bc3f204d6e03ef7440ffcdac5b3a547dae464 100644
--- a/Processes/CONEX/testConex.cc
+++ b/Processes/CONEX/testConex.cc
@@ -27,8 +27,8 @@ TEST_CASE("CONEX", "[processes]") {
using std::endl;
std::string parameterPathName = "";
- //auto cxModel = eSibyll23;
- //ConexDynamicInterface cx(cxModel);
+ // auto cxModel = eSibyll23;
+ // ConexDynamicInterface cx(cxModel);
int randomSeeds[3];
randomSeeds[0] = 1234;
@@ -38,7 +38,7 @@ TEST_CASE("CONEX", "[processes]") {
int nShower = 1; // large to avoid final stats.
int maxDetail = 0;
int particleListMode = 0;
- //cx.Init(nShower, randomSeeds, maxDetail, particleListMode, parameterPathName);
+ // cx.Init(nShower, randomSeeds, maxDetail, particleListMode, parameterPathName);
double energyInGeV = 100.;
double zenith = 60;
@@ -46,6 +46,7 @@ TEST_CASE("CONEX", "[processes]") {
double impactParameter = 0;
int particleType = 100;
- //cx.RunConex(randomSeeds, energyInGeV, zenith, azimuth, impactParameter, particleType);
+ // cx.RunConex(randomSeeds, energyInGeV, zenith, azimuth, impactParameter,
+ // particleType);
}
}
diff --git a/Processes/CONEXSourceCut/CONEXSourceCut.cc b/Processes/CONEXSourceCut/CONEXSourceCut.cc
index 1f84fbf257dd969f5452601aad06c13e12d7eee1..0808151768f8767a86e2355a151fbf9eaa82df98 100644
--- a/Processes/CONEXSourceCut/CONEXSourceCut.cc
+++ b/Processes/CONEXSourceCut/CONEXSourceCut.cc
@@ -76,10 +76,11 @@ corsika::process::EProcessReturn CONEXSourceCut::DoSecondaries(
<< std::endl;
p.Delete();
- //conex_.Shower(egs_pid, E, x, y, altitude, slantDistance, lateralX, lateralY, slantX,
+ // conex_.Shower(egs_pid, E, x, y, altitude, slantDistance, lateralX, lateralY,
+ // slantX,
// time, u, v, w, iri, weight, latchin);
conex::show_(egs_pid, E, x, y, altitude, slantDistance, lateralX, lateralY, slantX,
- time, u, v, w, iri, weight, latchin);
+ time, u, v, w, iri, weight, latchin);
}
return corsika::process::EProcessReturn::eOk;
@@ -92,14 +93,14 @@ void CONEXSourceCut::SolveCE() {
int iCEmode = 1;
int id = 0; // RU: max, fix this
int nshtot_ = 0; // RU: max, fix this
- //conex_.HadronCascade(id, nshtot_, zero, iCEmode);
- //conex_.SolveMomentEquations(zero);
+ // conex_.HadronCascade(id, nshtot_, zero, iCEmode);
+ // conex_.SolveMomentEquations(zero);
conex::hadroncascade_(id, nshtot_, zero, iCEmode);
conex::solvemomentequations_(zero);
// RU: this here is from cxroot,
- //int nX = conex_.GetNumberOfDepthBins(); // make sure this works!
+ // int nX = conex_.GetNumberOfDepthBins(); // make sure this works!
int nX = conex::get_number_of_depth_bins_(); // make sure this works!
int icut = 1;
@@ -130,12 +131,12 @@ void CONEXSourceCut::SolveCE() {
// conex_.GetGammaProfile(icutg, nX, Gamma[0]);
// conex_.GetElectronProfile(icute, nX, Electrons[0]);
// conex_.GetHadronProfile(icuth, nX, Hadrons[0]);
- conex::get_shower_data_(icut, iSec, nX, X[0], N[0], fitpars[0], H[0], D[0]);
- conex::get_shower_edep_(icut, nX, dEdX[0], EGround[0]);
- conex::get_shower_muon_(icutm, nX, Mu[0], dMu[0]);
- conex::get_shower_gamma_(icutg, nX, Gamma[0]);
- conex::get_shower_electron_(icute, nX, Electrons[0]);
- conex::get_shower_hadron_(icuth, nX, Hadrons[0]);
+ conex::get_shower_data_(icut, iSec, nX, X[0], N[0], fitpars[0], H[0], D[0]);
+ conex::get_shower_edep_(icut, nX, dEdX[0], EGround[0]);
+ conex::get_shower_muon_(icutm, nX, Mu[0], dMu[0]);
+ conex::get_shower_gamma_(icutg, nX, Gamma[0]);
+ conex::get_shower_electron_(icute, nX, Electrons[0]);
+ conex::get_shower_hadron_(icuth, nX, Hadrons[0]);
}
// RU: move all the non-C8 code from the following c++ function into a new file. Here we
@@ -146,8 +147,8 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s
// units::si::GrammageType Xcut,
units::si::HEPEnergyType primaryEnergy,
particles::PDGCode primaryID)
- : //conex_{ConexDynamicInterface(eSibyll23)}
- center_{center}
+ : // conex_{ConexDynamicInterface(eSibyll23)}
+ center_{center}
, showerAxis_{showerAxis}
, groundDist_{groundDist}
, conexObservationCS_{std::invoke([&]() {
@@ -171,21 +172,18 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s
int randomSeeds[3] = {1234, 0, 0}; // will be overwritten later??
int heModel = eSibyll23;
-
+
int nShower = 1; // large to avoid final stats.
int maxDetail = 0;
int particleListMode = 0;
- //conex_.Init(nShower, randomSeeds, maxDetail, particleListMode, parameterPathName);
+ // conex_.Init(nShower, randomSeeds, maxDetail, particleListMode, parameterPathName);
std::string configPath = CONEX_CONFIG_PATH;
- conex::initconex_(nShower, randomSeeds,
- heModel,
- maxDetail,
+ conex::initconex_(nShower, randomSeeds, heModel, maxDetail,
#ifdef CONEX_EXTENSIONS
- particleListMode,
+ particleListMode,
#endif
- configPath.c_str(),
- configPath.size());
+ configPath.c_str(), configPath.size());
double eprima = primaryEnergy / 1_GeV;
@@ -209,6 +207,6 @@ CONEXSourceCut::CONEXSourceCut(geometry::Point center, environment::ShowerAxis s
std::array<int, 3> ioseed{static_cast<int>(rng()), static_cast<int>(rng()),
static_cast<int>(rng())};
- //conex_.ConexRun(ipart, eprima, theta, phi, dimpact, ioseed.data());
+ // conex_.ConexRun(ipart, eprima, theta, phi, dimpact, ioseed.data());
conex::conexrun_(ipart, eprima, theta, phi, dimpact, ioseed.data());
}
diff --git a/Processes/CONEXSourceCut/CONEXSourceCut.h b/Processes/CONEXSourceCut/CONEXSourceCut.h
index 1686897a507e2b50396e46fa9df07e45072f97ad..1258be3954d70b5df339651caf7040401035010e 100644
--- a/Processes/CONEXSourceCut/CONEXSourceCut.h
+++ b/Processes/CONEXSourceCut/CONEXSourceCut.h
@@ -26,7 +26,7 @@
namespace conex {
corsika::units::si::LengthType constexpr earthRadius{6371315 *
- corsika::units::si::meter};
+ corsika::units::si::meter};
} // namespace conex
namespace corsika::process {
@@ -45,7 +45,7 @@ namespace corsika::process {
void SolveCE();
private:
- //ConexDynamicInterface conex_;
+ // ConexDynamicInterface conex_;
//! CONEX e.m. particle codes
static std::array<std::pair<particles::Code, int>, 3> constexpr egs_em_codes_{
diff --git a/Processes/CONEXSourceCut/CONEX_f.h b/Processes/CONEXSourceCut/CONEX_f.h
index a00bc7e35dac6351060f7f6db2dd3f509b093105..53719575142dd44917b5f600b98ae73470b5faca 100644
--- a/Processes/CONEXSourceCut/CONEX_f.h
+++ b/Processes/CONEXSourceCut/CONEX_f.h
@@ -19,33 +19,32 @@
#include <conexHEModels.h>
namespace conex {
-
- extern "C" {
- // ipart,energy,theta,phi,dimpact,ioseed
- void initconex_(int&, int*, int&, int&,
+
+ extern "C" {
+ // ipart,energy,theta,phi,dimpact,ioseed
+ void initconex_(int&, int*, int&, int&,
#ifdef CONEX_EXTENSIONS
- int&,
+ int&,
#endif
- const char*, int);
- void conexrun_(int& ipart, double& energy, double& theta, double& phi, double&
- dimpact, int ioseed[3]);
- void conexcascade_();
- void hadroncascade_(int&, int&, int&,
- int&);
- void solvemomentequations_(int&);
- void show_(int& iqi, double& ei, double& xmi,
- double& ymi, double& zmi, double& dmi, double& xi, double& yi, double& zi, double&
- tmi, double& ui, double& vi, double& wi, int& iri, double& wti, int& latchi);
-
- int get_number_of_depth_bins_();
-
- void get_shower_data_(const int&, const int&, const int&, float&, float&,
- float&, float&, float&);
- void get_shower_edep_(const int&, const int&, float&, float&);
- void get_shower_muon_(const int&, const int&, float&, float&);
- void get_shower_gamma_(const int&, const int&, float&);
- void get_shower_electron_(const int&, const int&, float&);
- void get_shower_hadron_(const int&, const int&, float&);
- }
-
-}
+ const char*, int);
+ void conexrun_(int& ipart, double& energy, double& theta, double& phi, double& dimpact,
+ int ioseed[3]);
+ void conexcascade_();
+ void hadroncascade_(int&, int&, int&, int&);
+ void solvemomentequations_(int&);
+ void show_(int& iqi, double& ei, double& xmi, double& ymi, double& zmi, double& dmi,
+ double& xi, double& yi, double& zi, double& tmi, double& ui, double& vi,
+ double& wi, int& iri, double& wti, int& latchi);
+
+ int get_number_of_depth_bins_();
+
+ void get_shower_data_(const int&, const int&, const int&, float&, float&, float&,
+ float&, float&);
+ void get_shower_edep_(const int&, const int&, float&, float&);
+ void get_shower_muon_(const int&, const int&, float&, float&);
+ void get_shower_gamma_(const int&, const int&, float&);
+ void get_shower_electron_(const int&, const int&, float&);
+ void get_shower_hadron_(const int&, const int&, float&);
+ }
+
+} // namespace conex