Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156
import FWCore.ParameterSet.Config as cms

from DQMOffline.Trigger.JetMonitor_cfi import hltJetMETmonitoring
from DQMOffline.Trigger.HTMonitor_cfi import hltHTmonitoring
from DQMOffline.Trigger.MjjMonitor_cfi import hltMjjmonitoring
from DQMOffline.Trigger.SoftdropMonitor_cfi import hltSoftdropmonitoring
from DQMOffline.Trigger.B2GTnPMonitor_cfi import B2GegmGsfElectronIDsForDQM,B2GegHLTDQMOfflineTnPSource
from DQMOffline.Trigger.TopMonitor_cfi import hltTOPmonitoring

### B2G triggers:
# HLT_AK8PFJet*_SoftDropMass*
# HLT_AK8DiPFJet*_*_SoftDropMass*
# HLT_Mu37_Ele27_CaloIdL_MW
# HLT_Mu27_Ele37_CaloIdL_MW
# HLT_Mu37_TkMu27
#
# Additionally, we monitor mjj and mSD for PFHT1050 and AK8PFJet500

# HT and AK8jet monitoring

PFHT1050_Mjjmonitoring = hltMjjmonitoring.clone(
    FolderName = 'HLT/B2G/PFHT1050',
    jets = "ak8PFJetsPuppi",
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet= dict(hltPaths = ["HLT_PFHT1050_v*"])     
)

PFHT1050_Softdropmonitoring = hltSoftdropmonitoring.clone(
    FolderName = 'HLT/B2G/PFHT1050',
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet = dict(hltPaths = ["HLT_PFHT1050_v*"])
)

AK8PFJet500_Mjjmonitoring = hltMjjmonitoring.clone(
    FolderName = 'HLT/B2G/AK8PFJet500',
    jets = "ak8PFJetsPuppi",
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet = dict(hltPaths = ["HLT_AK8PFJet500_v*"])
)

AK8PFJet380_SoftDropMass30_Mjjmonitoring = hltMjjmonitoring.clone(
    FolderName = 'HLT/B2G/AK8PFJet380_SoftDropMass30',
    jets = "ak8PFJetsPuppi",
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet= dict(hltPaths = ["HLT_AK8PFJet380_SoftDropMass30_v*"])
)

AK8DiPFJet260_260_SoftDropMass30_Mjjmonitoring = hltMjjmonitoring.clone(
    FolderName = 'HLT/B2G/AK8DiPFJet260_260_SoftDropMass30',
    jets = "ak8PFJetsPuppi",
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet= dict(hltPaths = ["HLT_AK8DiPFJet260_260_SoftDropMass30_v*"])
)

AK8PFJet500_Softdropmonitoring = hltSoftdropmonitoring.clone(
    FolderName = 'HLT/B2G/AK8PFJet500',
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet= dict(hltPaths = ["HLT_AK8PFJet500_v*"]),
    histoPSet = dict(
        htBinning = [0., 10., 20., 30., 40., 50., 60., 70., 80., 90., 100., 110., 120., 130., 140., 150., 160., 170., 180., 190., 200., 210., 220., 230., 240., 250., 260., 270., 280., 290., 300., 310., 320., 330., 340., 350.],
        htPSet = dict(nbins = 200, xmin = -0.5, xmax = 19999.5)
    )
)

# AK8PFJet380_SoftDropMass30 monitoring

AK8PFJet380_SoftDropMass30_PromptMonitoring = hltJetMETmonitoring.clone(
    FolderName = 'HLT/B2G/AK8PFJet380_SoftDropMass30',
    ptcut = 200,
    numGenericTriggerEventPSet = dict(hltPaths = ["HLT_AK8PFJet380_SoftDropMass30_v*"])
)


AK8PFJet380_SoftDropMass30_Softdropmonitoring = hltSoftdropmonitoring.clone(
    FolderName = 'HLT/B2G/AK8PFJet380_SoftDropMass30',
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet = dict(hltPaths = ["HLT_AK8PFJet380_SoftDropMass30_v*"]),
    histoPSet = dict(
        htBinning = [0., 10., 20., 30., 40., 50., 60., 70., 80., 90., 100., 110., 120., 130., 140., 150., 160., 170., 180., 190., 200., 210., 220., 230., 240., 250., 260., 270., 280., 290., 300., 310., 320., 330., 340., 350.],
        htPSet = dict(nbins = 200, xmin = -0.5, xmax = 19999.5)
    )
)

# AK8DiPFJet260_260_SoftDropMass30 monitoring

AK8DiPFJet260_260_SoftDropMass30_PromptMonitoring = hltJetMETmonitoring.clone(
    FolderName = 'HLT/B2G/AK8DiPFJet260_260_SoftDropMass30',
    ptcut = 200,
    numGenericTriggerEventPSet = dict(hltPaths = ["HLT_AK8DiPFJet260_260_SoftDropMass30_v*"])
)


AK8DiPFJet260_260_SoftDropMass30_Softdropmonitoring = hltSoftdropmonitoring.clone(
    FolderName = 'HLT/B2G/AK8DiPFJet260_260_SoftDropMass30',
    jetSelection = "pt > 200 && eta < 2.4",
    numGenericTriggerEventPSet = dict(hltPaths = ["HLT_AK8DiPFJet260_260_SoftDropMass30_v*"]),
    histoPSet = dict(
        htBinning = [0., 10., 20., 30., 40., 50., 60., 70., 80., 90., 100., 110., 120., 130., 140., 150., 160., 170., 180., 190., 200., 210., 220., 230., 240., 250., 260., 270., 280., 290., 300., 310., 320., 330., 340., 350.],
        htPSet = dict(nbins = 200, xmin = -0.5, xmax = 19999.5)
    )
)

# Lepton cross trigger monitoring

hltDQMonitorB2G_MuEle = hltTOPmonitoring.clone(
    FolderName = 'HLT/B2G/Dileptonic/HLT_MuXX_EleXX_CaloIdL_MW',
    nelectrons = 1,
    eleSelection = 'pt>20 & abs(eta)<2.4',
    nmuons = 1,
    muoSelection = 'pt>20 & abs(eta)<2.4 & ((pfIsolationR04.sumChargedHadronPt + max(pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - (pfIsolationR04.sumPUPt)/2.,0.))/pt < 0.25)  & isPFMuon & (isTrackerMuon || isGlobalMuon)',
    numGenericTriggerEventPSet = dict(hltPaths = ['HLT_Mu37_Ele27_CaloIdL_MW_v*', 'HLT_Mu27_Ele37_CaloIdL_MW_v*'])
)

hltDQMonitorB2G_MuTkMu = hltTOPmonitoring.clone(
    FolderName = 'HLT/B2G/Dileptonic/HLT_Mu37_TkMu27',
    nmuons = 2,
    muoSelection = 'pt>20 & abs(eta)<2.4 & ((pfIsolationR04.sumChargedHadronPt + max(pfIsolationR04.sumNeutralHadronEt + pfIsolationR04.sumPhotonEt - (pfIsolationR04.sumPUPt)/2.,0.))/pt < 0.25)  & isPFMuon & (isTrackerMuon || isGlobalMuon)',
    numGenericTriggerEventPSet = dict(hltPaths = ['HLT_Mu37_TkMu27_v*'])
)


# the sequence

b2gMonitorHLT = cms.Sequence(

    PFHT1050_Mjjmonitoring +

    AK8PFJet500_Mjjmonitoring +

    AK8PFJet380_SoftDropMass30_Mjjmonitoring +
    AK8DiPFJet260_260_SoftDropMass30_Mjjmonitoring +

    AK8PFJet380_SoftDropMass30_PromptMonitoring +
    AK8DiPFJet260_260_SoftDropMass30_PromptMonitoring +

    B2GegHLTDQMOfflineTnPSource

  * hltDQMonitorB2G_MuEle
  * hltDQMonitorB2G_MuTkMu

  , cms.Task(B2GegmGsfElectronIDsForDQM) ## unschedule execution [Use of electron VID requires this module being executed first]
)

## as reported in https://github.com/cms-sw/cmssw/issues/24444
## it turned out that all softdrop modules rely on a jet collection which is available only if the miniAOD step is run @Tier0
## ==> it is fine in the PromptReco workflow, but this collection is not available in the Express reconstruction
## in addition, it is not available in the AOD (!!!!) ==> these modules needs to be run *WithRECO* step workflow (actually w/ the miniAOD step ....)
b2gHLTDQMSourceWithRECO = cms.Sequence(
    PFHT1050_Softdropmonitoring +
    AK8PFJet500_Softdropmonitoring +
    AK8PFJet380_SoftDropMass30_Softdropmonitoring +
    AK8DiPFJet260_260_SoftDropMass30_Softdropmonitoring
)

b2gHLTDQMSourceExtra = cms.Sequence(
)