Neutrino interactions with pythia 8310

Never mind the branch name. This uses the latest pythia 8.310.

applications/c8_air_shower.cpp

@@ -242,7 +242,6 @@ int main(int argc, char** argv) {
       ->default_val(0)
       ->check(CLI::Range(0, 20))
       ->group("Radio");
-
   // parse the command line options into the variables
   CLI11_PARSE(app, argc, argv);
 
@@ -399,8 +398,8 @@ int main(int argc, char** argv) {
 
   corsika::pythia8::Decay decayPythia;
 
-  // neutrino interactions
-  corsika::pythia8::NeutrinoInteraction neutrinoPrimaryPythia(E0);
+  // neutrino interactions (options are: NC, CC)
+  corsika::pythia8::NeutrinoInteraction neutrinoPrimaryPythia(true, true);
 
   // hadronic photon interactions in resonance region
   corsika::sophia::InteractionModel sophia;
@@ -591,7+590,7 @@
     auto const primaryProperties = std::make_tuple(
         beamCode, calculate_kinetic_energy(plab.getNorm(), get_mass(beamCode)),
         plab.normalized(), injectionPos, 0_ns);
     stack.addParticle(primaryProperties);
 
     // if we want to fix the first location of the shower
     if (force_interaction) {