Count BB4L FSR vetos

cms-sw · May 16, 2017 · 417c42f · 417c42f
1 parent 0fdfecb
commit 417c42f
Show file tree

Hide file tree

Showing 2 changed files with 32 additions and 10 deletions.
diff --git a/GeneratorInterface/Pythia8Interface/plugins/PowhegHooksBB4L.h b/GeneratorInterface/Pythia8Interface/plugins/PowhegHooksBB4L.h
@@ -280,31 +280,45 @@ class PowhegHooksBB4L : public UserHooks {
 			}
 
 			if (iTopCharge > 0) {
-				if (onlyDistance1)
-					return (distance == 1) && scale > topresscale;
+				if (onlyDistance1) {
+					if ((distance == 1) && scale > topresscale) {
+					  nFSRveto++;
+					  return true;
+					}
+				}
 				else
-					return scale > topresscale;
+					if (scale > topresscale) {
+					  nFSRveto++;
+					  return true;
+					}
 			}
 			else {
-				if (onlyDistance1)
-					return (distance == 1) && scale > atopresscale;
+				if (onlyDistance1) {
+					if ((distance == 1) && scale > atopresscale) {
+					  nFSRveto++;
+					  return true;
+					}
+				}
 				else
-					return scale > atopresscale;
+					if (scale > atopresscale) {
+					  nFSRveto++;
+					  return true;
+					}
 			}
-		} else {
-			return false;
 		}
+		return false;
 	}
 
 //--------------------------------------------------------------------------
 
   // Functions to return information
 
-//  inline int    getNFSRveto() { return nFSRveto; }
+  inline int    getNFSRveto() { return nFSRveto; }
 
 //--------------------------------------------------------------------------
 
 private:
+	int nFSRveto = 0;
 	double topresscale, atopresscale;
 	bool onlyDistance1, useScaleResonanceInstead;
 };

diff --git a/GeneratorInterface/Pythia8Interface/plugins/Pythia8Hadronizer.cc b/GeneratorInterface/Pythia8Interface/plugins/Pythia8Hadronizer.cc
@@ -157,6 +157,7 @@ class Pythia8Hadronizer : public BaseHadronizer, public Py8InterfaceBase {
 
     int nISRveto;
     int nFSRveto;
+    int nFSRvetoBB4L;
 
     int NHooks;
 
@@ -171,7 +172,7 @@ Pythia8Hadronizer::Pythia8Hadronizer(const edm::ParameterSet &params) :
   fInitialState(PP),
   fMultiUserHook(new MultiUserHook), fReweightUserHook(0),fReweightRapUserHook(0),fReweightPtHatRapUserHook(0),
   fJetMatchingHook(0),fJetMatchingPy8InternalHook(0), fMergingHook(0),
-  fEmissionVetoHook(0), fEmissionVetoHook1(0), fResonanceDecayFilterHook(0), fPTFilterHook(0), nME(-1), nMEFiltered(-1), nISRveto(0), nFSRveto(0),
+  fEmissionVetoHook(0), fEmissionVetoHook1(0), fResonanceDecayFilterHook(0), fPTFilterHook(0), nME(-1), nMEFiltered(-1), nISRveto(0), nFSRveto(0), nFSRvetoBB4L(0),
   NHooks(0)
 {
 
@@ -627,6 +628,10 @@ void Pythia8Hadronizer::statistics()
     edm::LogPrint("Pythia8Interface")
       << "Number of FSR vetoed = " << nFSRveto;
   }
+  if(fPowhegHooksBB4L) {
+    edm::LogInfo("Pythia8Interface") << "\n"
+      << "BB4L: Number of FSR vetoed = " << nFSRvetoBB4L;
+  }
 
   double xsec = fMasterGen->info.sigmaGen(); // cross section in mb
   xsec *= 1.0e9; // translate to pb (CMS/Gen "convention" as of May 2009)
@@ -662,6 +667,9 @@ bool Pythia8Hadronizer::generatePartonsAndHadronize()
     nISRveto += fEmissionVetoHook->getNISRveto();
     nFSRveto += fEmissionVetoHook->getNFSRveto();  
   }
+  if (fPowhegHooksBB4L) {
+    nFSRvetoBB4L += fPowhegHooksBB4L->getNFSRveto();
+  }
 
   return true;