Merge branch 'sibyll' of gitlab.ikp.kit.edu:AirShowerPhysics/corsika into sibyll

Documentation/Examples/cascade_example.cc

@@ -51,6 +51,8 @@ public:
     auto test = static_cast<corsika::process::sibyll::SibyllCodeIntType>(s_id);
     std::cout << "MinStepLength: particle input " << "sibyll id: |" << (int)test << "|" <<std::endl;
 											  // std::cout << "MinStepLength: particle input " << "sibyll id: " << process::sibyll::ConvertToSibyll( p.GetPID() ) << std::endl;
+    cout << p.GetPID() << " --> " << process::sibyll::ConvertToSibyllRaw( p.GetPID() ) << endl;
+    
     int kBeam   = 1;
 
     /* 

Processes/Sibyll/ParticleConversion.h

@@ -42,6 +42,10 @@ namespace corsika::process::sibyll {
     return sibyll2corsika[static_cast<SibyllCodeIntType>(pCode) - minSibyll];
   }
 
+  int ConvertToSibyllRaw(corsika::particles::Code pCode){
+    return  (int)static_cast<corsika::process::sibyll::SibyllCodeIntType>( corsika::process::sibyll::ConvertToSibyll( pCode ) );
+  }
+  
 } // namespace corsika::process::sibyll
 
 #endif

Processes/Sibyll/testSibyll.cc

@@ -29,6 +29,8 @@ TEST_CASE("Sibyll", "[processes]") {
   SECTION("Corsika -> Sibyll") {
     REQUIRE(process::sibyll::ConvertToSibyll(corsika::particles::Electron::GetCode()) ==
             process::sibyll::Code::Electron);
+    REQUIRE(process::sibyll::ConvertToSibyllRaw(corsika::particles::Proton::GetCode()) ==
+            13 );
   }
 
   SECTION("handledBySibyll") {