fixed physics processes

Processes/Pythia/Decay.cc

@@ -148,9 +148,6 @@ namespace corsika::process::pythia {
     // set particle stable
     Decay::SetStable(vP.GetPID());
 
-    // remove original particle from corsika stack
-    vP.Delete();
-    //    if (fCount>10) throw std::runtime_error("stop here");
   }
 
 } // namespace corsika::process::pythia

Processes/Pythia/Interaction.cc

@@ -411,8 +411,6 @@ namespace corsika::process::pythia {
              << "Elab_final=" << Elab_final / 1_GeV
              << ", Plab_final=" << (Plab_final / 1_GeV).GetComponents() << endl;
       }
-      // delete current particle
-      vP.Delete();
     }
     return process::EProcessReturn::eOk;
   }

Processes/Sibyll/Interaction.cc

@@ -366,8 +366,6 @@ namespace corsika::process::sibyll {
              << ", Plab_final=" << (Plab_final / 1_GeV).GetComponents() << endl;
       }
     }
-    // delete current particle
-    vP.Delete();
     return process::EProcessReturn::eOk;
   }