Changing TauDecay to not use intermediate resonance for pythia8

cms-sw · Sep 14, 2013 · 4281919 · 4281919
1 parent f1cc058
commit 4281919
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/Validation/EventGenerator/src/TauDecay.cc b/Validation/EventGenerator/src/TauDecay.cc
@@ -195,7 +195,7 @@ void TauDecay::ClassifyDecayMode(unsigned int &JAK_ID,unsigned int &TauBitMask){
     JAK_ID=JAK_PION;
     return;
   }
-  if(n_e==0 && n_mu==0 && n_pi==1 && n_pi0==1 && n_K==0 && n_K0L==0 && n_K0S==0 && n_nu==1 && n_rho==1){
+  if(n_e==0 && n_mu==0 && n_pi==1 && n_pi0==1 && n_K==0 && n_K0L==0 && n_K0S==0 && n_nu==1){// && n_rho==1){ removing intermediate resoance to be compatible with pythia8
     JAK_ID=JAK_RHO_PIPI0;
     return;
   }
@@ -247,10 +247,11 @@ void TauDecay::ClassifyDecayMode(unsigned int &JAK_ID,unsigned int &TauBitMask){
     JAK_ID=JAK_KPIPI;
     return;
   }
-  if(n_e==0 && n_mu==0 && n_pi==1 && n_pi0==1 && n_K==0 && n_K0L==0 && n_K0S==0 && n_nu==1 && n_gamma>=1 && n_rho==0){
+  // removing JAKID 21 to allow for compatibility with Pythia8 
+  /*  if(n_e==0 && n_mu==0 && n_pi==1 && n_pi0==1 && n_K==0 && n_K0L==0 && n_K0S==0 && n_nu==1 && n_gamma>=1 && n_rho==0){
     JAK_ID=JAK_PIPI0GAM;
     return;
-  }
+    }*/
   std::cout << "Tau Mode not found: n_e " <<  n_e << " n_mu " << n_mu << " n_pi " << n_pi << " n_pi0 " << n_pi0 << " n_K " << n_K << "  n_K0L " << n_K0L << "  n_K0S " << n_K0S << " n_nu  " << n_nu << " n_gamma " << n_gamma << std::endl;
   JAK_ID=JAK_UNKNOWN;
 }