Merge pull request #10549 from heppye/patch-5

add DQM Validation plots for MSSM Hbb Trigger Paths.
cms-sw · Sep 1, 2015 · 55ac7a3 · 55ac7a3
2 parents de24a50 + a6f14f0
commit 55ac7a3
Show file tree

Hide file tree

Showing 2 changed files with 42 additions and 7 deletions.
diff --git a/HLTriggerOffline/Higgs/python/hltHiggsPostProcessors_cff.py b/HLTriggerOffline/Higgs/python/hltHiggsPostProcessors_cff.py
@@ -182,7 +182,7 @@ def get_reco_strings(strings):
 #dEtaqq, mqq, dPhibb, CVS1, maxCSV_jets, maxCSV_E, MET, pt1, pt2, pt3, pt4
 NminOneCutNames = ("EffdEtaqq", "Effmqq", "EffdPhibb", "EffCSV1", "EffCSV2", "EffCSV3",  "EffmaxCSV", "", "", "TurnOn1", "TurnOn2", "TurnOn3", "TurnOn4")
 plot_types = []
-NminOneCuts = (_config.__getattribute__("VBFHbb_2btag")).__getattribute__("NminOneCuts")
+NminOneCuts = _config.VBFHbb_2btag.NminOneCuts
 if NminOneCuts: 
     for iCut in range(0,len(NminOneCuts)):
         if( NminOneCuts[iCut] and NminOneCutNames[iCut] ):
@@ -206,7 +206,7 @@ def get_reco_strings(strings):
 #Specific plots for VBFHbb_1btag  
 #dEtaqq, mqq, dPhibb, CVS1, maxCSV_jets, maxCSV_E, MET, pt1, pt2, pt3, pt4
 plot_types = []
-NminOneCuts = (_config.__getattribute__("VBFHbb_1btag")).__getattribute__("NminOneCuts")
+NminOneCuts = _config.VBFHbb_1btag.NminOneCuts
 if NminOneCuts: 
     for iCut in range(0,len(NminOneCuts)):
         if( NminOneCuts[iCut] and NminOneCutNames[iCut] ):
@@ -230,7 +230,7 @@ def get_reco_strings(strings):
 #Specific plots for VBFHbb_0btag  
 #dEtaqq, mqq, dPhibb, CVS1, maxCSV_jets, maxCSV_E, MET, pt1, pt2, pt3, pt4
 plot_types = []
-NminOneCuts = (_config.__getattribute__("VBFHbb_0btag")).__getattribute__("NminOneCuts")
+NminOneCuts = _config.VBFHbb_0btag.NminOneCuts
 if NminOneCuts: 
     for iCut in range(0,len(NminOneCuts)):
         if( NminOneCuts[iCut] and NminOneCutNames[iCut] ):
@@ -255,7 +255,7 @@ def get_reco_strings(strings):
 #Specific plots for ZnnHbb
 #Jet plots
 plot_types = ["EffEta", "EffPhi"]
-NminOneCuts = (_config.__getattribute__("ZnnHbb")).__getattribute__("NminOneCuts")
+NminOneCuts = _config.ZnnHbb.NminOneCuts
 if NminOneCuts: 
     for iCut in range(0,len(NminOneCuts)):
         if( NminOneCuts[iCut] and NminOneCutNames[iCut] ):
@@ -288,7 +288,7 @@ def get_reco_strings(strings):
 
 #Specific plots for X4b
 #Jet plots
-NminOneCuts = (_config.__getattribute__("X4b")).__getattribute__("NminOneCuts")
+NminOneCuts = _config.X4b.NminOneCuts
 if NminOneCuts: 
     for iCut in range(0,len(NminOneCuts)):
         if( NminOneCuts[iCut] and NminOneCutNames[iCut] ):
@@ -320,6 +320,27 @@ def get_reco_strings(strings):
 hltHiggsPostWHToENuBB.subDirs = ['HLT/Higgs/WHToENuBB']
 hltHiggsPostWHToENuBB.efficiencyProfile = efficiency_strings_WHToENuBB
 
+#Specific plots for MSSMHbb
+#Jet plots
+NMinOneCuts = _config.MSSMHbb.NminOneCuts
+if NminOneCuts:
+    for iCut in range(0,len(NminOneCuts)):
+        if( NminOneCuts[iCut] and NminOneCutNames[iCut] ):
+            plot_types.append(NminOneCutNames[iCut])
+
+efficiency_strings = []
+for type in plot_types:
+    for obj in ["Jet"]:
+        for trig in triggers:
+            efficiency_strings.append(efficiency_string(obj,type,trig))
+
+efficiency_strings = get_reco_strings(efficiency_strings)
+
+hltHiggsPostMSSMHbb = hltHiggsPostProcessor.clone()
+hltHiggsPostMSSMHbb.subDirs = ['HLT/Higgs/MSSMHbb']
+hltHiggsPostMSSMHbb.efficiencyProfile = efficiency_strings
+
+
 hltHiggsPostProcessors = cms.Sequence(
         hltHiggsPostHWW+
         hltHiggsPostHZZ+
@@ -336,7 +357,8 @@ def get_reco_strings(strings):
         hltHiggsPostDoubleHinTaus+
         hltHiggsPostHiggsDalitz+
         hltHiggsPostX4b+
-        hltHiggsPostWHToENuBB
+        hltHiggsPostWHToENuBB+
+        hltHiggsPostMSSMHbb
 
 )
 

diff --git a/HLTriggerOffline/Higgs/python/hltHiggsValidator_cfi.py b/HLTriggerOffline/Higgs/python/hltHiggsValidator_cfi.py
@@ -4,7 +4,7 @@
 hltHiggsValidator = cms.EDAnalyzer("HLTHiggsValidator",
 
     hltProcessName = cms.string("HLT"),
-    analysis       = cms.vstring("HWW", "HZZ", "Hgg", "HggControlPaths", "Htaunu", "H2tau", "VBFHbb_0btag", "VBFHbb_1btag", "VBFHbb_2btag",  "ZnnHbb","DoubleHinTaus","HiggsDalitz","X4b","TTHbbej","AHttH","WHToENuBB"),
+    analysis       = cms.vstring("HWW", "HZZ", "Hgg", "HggControlPaths", "Htaunu", "H2tau", "VBFHbb_0btag", "VBFHbb_1btag", "VBFHbb_2btag",  "ZnnHbb","DoubleHinTaus","HiggsDalitz","X4b","TTHbbej","AHttH","WHToENuBB","MSSMHbb"),
     histDirectory  = cms.string("HLT/Higgs"),
 
     # -- The instance name of the reco::GenParticles collection 
@@ -358,4 +358,17 @@
         recJetLabel  = cms.string("ak4PFJetsCHS"),
         minCandidates = cms.uint32(1),
         ),
+    MSSMHbb  = cms.PSet(
+        hltPathsToCheck = cms.vstring(
+            "HLT_DoubleJetsC100_DoubleBTagCSV0p85_DoublePFJetsC160_v",
+            "HLT_DoubleJetsC100_DoubleBTagCSV0p9_DoublePFJetsC100MaxDeta1p6_v",
+            "HLT_DoubleJetsC112_DoubleBTagCSV0p85_DoublePFJetsC172_v",
+            "HLT_DoubleJetsC112_DoubleBTagCSV0p9_DoublePFJetsC112MaxDeta1p6_v",
+            ),
+        recJetLabel  = cms.string("ak4PFJetsCHS"),
+        jetTagLabel  = cms.string("combinedInclusiveSecondaryVertexV2BJetTags"),
+        # -- Analysis specific cuts
+        minCandidates = cms.uint32(3),
+        NminOneCuts = cms.untracked.vdouble(0, 0, 0, 0.941, 0.941 , 0.00, 0, 0, 0, 100, 100, 0.0, 0.0),
+        ),
 )