Merge pull request #36 from JETSCAPE/fix_HH_bug_decays

Fix bug when not decaying hadrons Hybrid Hadronization (Pythia)
JETSCAPE · Sep 19, 2024 · 369570b · 369570b
2 parents 8637288 + 3d02b5b
commit 369570b
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/hadronization/HybridHadronization.cc b/src/hadronization/HybridHadronization.cc
@@ -709,6 +709,10 @@ void HybridHadronization::DoHadronization(vector<vector<shared_ptr<Parton>>>& sh
       } //do this only if there are negative partons, otherwise the baryon recombination would be switched off after the first event without negative partons
     }
     if(HH_shower.num() == 0){
+      th_recofactor = tmp_threco;
+      maxB_level = tmp_maxB_level;
+      reco_hadrons_pythia = tmp_reco_hadrons_pythia;
+      if(pythia_decays == "off"){pythia.readString("HadronLevel:Decay = off"); pythia.init();}
       continue;
     } //attempting to handle events/configurations with 0 partons will result in a crash