Skip to content

Commit

Permalink
Setup Run2, Run3 era dependent JMENanoDQM
Browse files Browse the repository at this point in the history
  • Loading branch information
nurfikri89 committed Dec 2, 2022
1 parent 635bebb commit 4aa7802
Showing 1 changed file with 59 additions and 6 deletions.
65 changes: 59 additions & 6 deletions PhysicsTools/NanoAOD/python/nanojmeDQM_cff.py
Expand Up @@ -12,7 +12,14 @@
# Add more variables for AK4 Puppi jets
#
#============================================
nanojmeDQM.vplots.Jet.plots.extend([
_ak4puppiplots = cms.VPSet(
Count1D('_size', 20, -0.5, 19.5, 'AK4 PF Puppi jets with JECs applied.')
)
for plot in nanojmeDQM.vplots.Jet.plots:
if plot.name.value()=="_size": continue
_ak4puppiplots.append(plot)

_ak4puppiplots.extend([
Plot1D('nConstChHads', 'nConstChHads', 20, 14.5, 34.5, 'number of charged hadrons in the jet'),
Plot1D('nConstElecs', 'nConstElecs', 3, -0.5, 2.5, 'number of electrons in the jet'),
Plot1D('nConstHFEMs', 'nConstHFEMs', 1, -0.5, 0.5, 'number of HF EMs in the jet'),
Expand Down Expand Up @@ -57,7 +64,7 @@
_ak4chsplots = cms.VPSet(
Count1D('_size', 20, -0.5, 19.5, 'AK4 PF CHS jets with JECs applied.')
)
for plot in nanojmeDQM.vplots.Jet.plots:
for plot in _ak4puppiplots:
if plot.name.value()=="_size": continue
_ak4chsplots.append(plot)
_ak4chsplots.extend([
Expand All @@ -66,11 +73,24 @@
Plot1D('chFPV3EF', 'chFPV3EF', 20, 0, 2, 'charged fromPV==3 Energy Fraction (component of the total charged Energy Fraction).'),
])

nanojmeDQM.vplots.JetCHS = cms.PSet(
#============================================
#
# Setup all extra AK4 collections. Will remove
# collection depending on era.
#
#============================================
nanojmeDQM.vplots.Jet.plots = _ak4puppiplots #Puppi is default "Jet collection" for Run-3
nanojmeDQM.vplots.JetPuppi = cms.PSet( # This is for the Run-2 extra "JetPuppi" collection
sels = nanojmeDQM.vplots.Jet.sels,
plots = _ak4puppiplots
)
nanojmeDQM.vplots.JetCHS = cms.PSet( # This is for the Run-3 extra "JetCHS" collection
sels = nanojmeDQM.vplots.Jet.sels,
plots = _ak4chsplots
)



#============================================
#
# Add more variables for AK8 Puppi jets
Expand Down Expand Up @@ -172,7 +192,6 @@
),
)


##MC
nanojmeDQMMC = nanojmeDQM.clone()
#nanojmeDQMMC.vplots.Electron.sels.Prompt = cms.string("genPartFlav == 1")
Expand All @@ -182,8 +201,42 @@
nanojmeDQMMC.vplots.Tau.sels.Prompt = cms.string("genPartFlav == 5")
nanojmeDQMMC.vplots.Jet.sels.Prompt = cms.string("genJetIdx != 1")
nanojmeDQMMC.vplots.Jet.sels.PromptB = cms.string("genJetIdx != 1 && hadronFlavour == 5")
nanojmeDQMMC.vplots.JetCHS.sels.Prompt = cms.string("genJetIdx != 1")
nanojmeDQMMC.vplots.JetCHS.sels.PromptB = cms.string("genJetIdx != 1 && hadronFlavour == 5")

#============================================
#
# Era dependent customization
#
#============================================
#
# Run 3
#
(~run2_nanoAOD_ANY).toModify(
nanojmeDQM.vplots.Jet,
plots = _ak4puppiplots,
).toModify(
nanojmeDQM.vplots,
JetPuppi = None # Remove "JetPuppi" from DQM
)
(~run2_nanoAOD_ANY).toModify(
nanojmeDQMMC.vplots.JetCHS.sels,
Prompt = nanojmeDQMMC.vplots.Jet.sels.Prompt,
PromptB = nanojmeDQMMC.vplots.Jet.sels.PromptB
)
#
# Run 2
#
run2_nanoAOD_ANY.toModify(
nanojmeDQM.vplots.Jet,
plots = _ak4chsplots, #
).toModify(
nanojmeDQM.vplots,
JetCHS = None # Remove "JetCHS" from DQM
)
run2_nanoAOD_ANY.toModify(
nanojmeDQMMC.vplots.JetPuppi.sels,
Prompt = nanojmeDQMMC.vplots.Jet.sels.Prompt,
PromptB = nanojmeDQMMC.vplots.Jet.sels.PromptB
)

from DQMServices.Core.DQMQualityTester import DQMQualityTester
nanoDQMQTester = DQMQualityTester(
Expand Down

0 comments on commit 4aa7802

Please sign in to comment.