fixed GetSibyllMass, do clang

Processes/Sibyll/ParticleConversion.cc

@@ -20,7 +20,10 @@ corsika::units::si::HEPMassType corsika::process::sibyll::GetSibyllMass(
   if (pCode == corsika::particles::Code::Nucleus)
     throw std::runtime_error("Cannot GetMass() of particle::Nucleus -> unspecified");
   auto sCode = ConvertToSibyllRaw(pCode);
-  return sqrt(get_sibyll_mass2(sCode)) * 1_GeV;
+  if (sCode == 0)
+    return std::numeric_limits<double>::quiet_NaN() * 1_GeV;
+  else
+    return sqrt(get_sibyll_mass2(sCode)) * 1_GeV;
 }
 
 // const std::map<sibyll::PID, ParticleProperties::InternalParticleCode>

Processes/Sibyll/sibyll2.3c.cc

@@ -14,7 +14,7 @@
 #include <random>
 
 int get_nwounded() { return s_chist_.nwd; }
-double get_sibyll_mass2(int& id) { return s_mass1_.am2[id]; }
+double get_sibyll_mass2(int& id) { return s_mass1_.am2[abs(id) - 1]; }
 
 double s_rndm_(int&) {
   static corsika::random::RNG& rng =

Processes/Sibyll/testSibyll.cc

@@ -40,7 +40,6 @@ TEST_CASE("Sibyll", "[processes]") {
     REQUIRE(process::sibyll::ConvertToSibyllRaw(particles::Proton::GetCode()) == 13);
     REQUIRE(process::sibyll::ConvertToSibyll(particles::XiStarC0::GetCode()) ==
             process::sibyll::SibyllCode::XiStarC0);
-
   }
 
   SECTION("canInteractInSibyll") {