Manually port Markus' commit, avoiding rivet version changes

copy of commit 1ec4f08 Clang format Fix typo Fix merge mistake
cms-sw · Nov 29, 2021 · ca42d89 · ca42d89
1 parent e9addd3
commit ca42d89
Show file tree

Hide file tree

Showing 4 changed files with 10 additions and 6 deletions.
diff --git a/GeneratorInterface/RivetInterface/interface/RivetAnalysis.h b/GeneratorInterface/RivetInterface/interface/RivetAnalysis.h
@@ -34,6 +34,7 @@ namespace Rivet {
       bool _usePromptFinalStates;
       bool _excludePromptLeptonsFromJetClustering;
       bool _excludeNeutrinosFromJetClustering;
+      bool _doJetClustering;
 
       double _particleMinPt, _particleMaxEta;
       double _lepConeSize, _lepMinPt, _lepMaxEta;
@@ -51,6 +52,7 @@ namespace Rivet {
       _usePromptFinalStates(pset.getParameter<bool>("usePromptFinalStates")),
       _excludePromptLeptonsFromJetClustering(pset.getParameter<bool>("excludePromptLeptonsFromJetClustering")),
       _excludeNeutrinosFromJetClustering(pset.getParameter<bool>("excludeNeutrinosFromJetClustering")),
+      _doJetClustering(pset.getParameter<bool>("doJetClustering")),
 
       _particleMinPt  (pset.getParameter<double>("particleMinPt")),
       _particleMaxEta (pset.getParameter<double>("particleMaxEta")),
@@ -201,8 +203,10 @@ namespace Rivet {
 
         }
 
-        _jets      = applyProjection<FastJets>(event, "Jets").jetsByPt(jet_cut);
-        _fatjets   = applyProjection<FastJets>(event, "FatJets").jetsByPt(fatjet_cut);
+        if (_doJetClustering) {
+          _jets = apply<FastJets>(event, "Jets").jetsByPt(jet_cut);
+          _fatjets = apply<FastJets>(event, "FatJets").jetsByPt(fatjet_cut);
+        }
         _neutrinos = applyProjection<FinalState>(event, "Neutrinos").particlesByPt();
         _met       = applyProjection<MissingMomentum>(event, "MET").missingMomentum().p3();
 

diff --git a/GeneratorInterface/RivetInterface/python/particleLevel_cfi.py b/GeneratorInterface/RivetInterface/python/particleLevel_cfi.py
@@ -6,6 +6,7 @@
     usePromptFinalStates = cms.bool(True), # for leptons, photons, neutrinos
     excludePromptLeptonsFromJetClustering = cms.bool(True),
     excludeNeutrinosFromJetClustering = cms.bool(True),
+    doJetClustering = cms.bool(True),
 
     particleMinPt  = cms.double(0.),
     particleMaxEta = cms.double(5.), # HF range. Maximum 6.0 on MiniAOD

diff --git a/PhysicsTools/NanoAOD/plugins/GenJetGenPartMerger.cc b/PhysicsTools/NanoAOD/plugins/GenJetGenPartMerger.cc
@@ -84,9 +84,9 @@ void GenJetGenPartMerger::produce(edm::Event& iEvent, const edm::EventSetup& iSe
   for (unsigned int ijet = 0; ijet < jetHandle->size(); ++ijet) {
     auto jet = jetHandle->at(ijet);
     if (cut_(jet)) {
-        merged->push_back(reco::GenJet(jet));
-        reco::GenJetRef jetRef(jetHandle, ijet);
-        hasTauAncValues.push_back((*tauAncHandle)[jetRef]);
+      merged->push_back(reco::GenJet(jet));
+      reco::GenJetRef jetRef(jetHandle, ijet);
+      hasTauAncValues.push_back((*tauAncHandle)[jetRef]);
     }
   }
 

diff --git a/PhysicsTools/NanoAOD/python/electrons_cff.py b/PhysicsTools/NanoAOD/python/electrons_cff.py
@@ -532,7 +532,6 @@ def _get_bitmapVIDForEle_docstring(modules,WorkingPoints):
 
 electronSequence = cms.Sequence(bitmapVIDForEle + bitmapVIDForEleHEEP + isoForEle + ptRatioRelForEle + seedGainEle + slimmedElectronsWithUserData + finalElectrons)
 electronTables = cms.Sequence (electronMVATTH + electronTable)
-<<<<<<< HEAD
 electronMCold = cms.Sequence(electronsMCMatchForTable + electronMCTable)
 electronMC = cms.Sequence(tautaggerForMatching + matchingElecPhoton + electronsMCMatchForTable + electronsMCMatchForTableAlt + electronMCTable)
 ( run2_nanoAOD_106Xv1 & ~run2_nanoAOD_devel).toModify( electronMCTable,