CMS 3D CMS Logo

triggerObjects_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
4 from PhysicsTools.NanoAOD.triggerObjectTableProducer_cfi import triggerObjectTableProducer
5 from PhysicsTools.NanoAOD.globalVariablesTableProducer_cfi import globalVariablesTableProducer
6 import copy
7 
8 unpackedPatTrigger = cms.EDProducer("PATTriggerObjectStandAloneUnpacker",
9  patTriggerObjectsStandAlone = cms.InputTag('slimmedPatTrigger'),
10  triggerResults = cms.InputTag('TriggerResults::HLT'),
11  unpackFilterLabels = cms.bool(True)
12 )
13 
14 def mksel( selection, doc=None, bit=None):
15  ddoc=""
16  if 'OR' in selection:
17  selection,ddoc=OR(selection)
18  if 'AND' in selection:
19  selection,ddoc=AND(selection)
20  if type(selection)==list:
21  selection,ddoc=OR(selection)
22  if type(selection)==tuple:
23  selection,ddoc=AND(selection)
24  if doc==None:
25  doc=ddoc
26  #print("creating an entry with",selection,doc)
27  return cms.PSet(selection=cms.string(selection),doc=cms.string(doc),bit=cms.uint32(bit)) if bit is not None else cms.PSet(selection=cms.string(selection),doc=cms.string(doc))
28 
29 def chaintoken_(tokens,f,OR_or_AND=None):
30  if (not type(tokens) in [list,tuple]):
31  doc=tokens
32  tokens= [t.strip().rstrip() for t in tokens.split(OR_or_AND)]
33  else:
34  doc=f' {OR_or_AND} '.join(tokens)
35  selection = f"filter('{tokens[0]}')"
36  for token in tokens[1:]:
37  filter = f"filter('{token}')"
38  selection = f"{f}({filter},{selection})"
39  return selection,doc
40 def OR(tokens):
41  return chaintoken_(tokens,"max","OR")
42 def AND(tokens):
43  return chaintoken_(tokens,"min","AND")
44 
45 triggerObjectTable = triggerObjectTableProducer.clone(
46  name= cms.string("TrigObj"),
47  src = cms.InputTag("unpackedPatTrigger"),
48  l1EG = cms.InputTag("caloStage2Digis","EGamma"),
49  l1Sum = cms.InputTag("caloStage2Digis","EtSum"),
50  l1Jet = cms.InputTag("caloStage2Digis","Jet"),
51  l1Muon = cms.InputTag("gmtStage2Digis","Muon"),
52  l1Tau = cms.InputTag("caloStage2Digis","Tau"),
53  selections = cms.PSet(
54  Electron = cms.PSet(
55  doc = cms.string("PixelMatched e/gamma"), # this may also select photons!
56  id = cms.int32(11),
57  sel = cms.string("type(92) && pt > 7 && coll('hltEgammaCandidates') && filter('*PixelMatchFilter')"),
58  l1seed = cms.string("type(-98)"), l1deltaR = cms.double(0.3),
59  #l2seed = cms.string("type(92) && coll('')"), l2deltaR = cms.double(0.5),
60  skipObjectsNotPassingQualityBits = cms.bool(True),
61  qualityBits = cms.VPSet(
62  mksel("filter('*CaloIdLTrackIdLIsoVL*TrackIso*Filter')","CaloIdL_TrackIdL_IsoVL"),
63  mksel("filter('hltEle*WPTight*TrackIsoFilter*')","1e (WPTight)"),
64  mksel("filter('hltEle*WPLoose*TrackIsoFilter')","1e (WPLoose)"),
65  mksel("filter('*OverlapFilter*IsoEle*PFTau*')","OverlapFilter PFTau"),
66  mksel("filter('hltEle*Ele*CaloIdLTrackIdLIsoVL*Filter')","2e"),
67  mksel("filter('hltMu*TrkIsoVVL*Ele*CaloIdLTrackIdLIsoVL*Filter*')","1e-1mu"),
68  mksel("filter('hlt*OverlapFilterIsoEle*PFTau*')","1e-1tau"),
69  mksel("filter('hltEle*Ele*Ele*CaloIdLTrackIdLDphiLeg*Filter')","3e"),
70  mksel(["hltL3fL1Mu*DoubleEG*Filtered*","hltMu*DiEle*CaloIdLTrackIdLElectronleg*Filter"],"2e-1mu"),
71  mksel(["hltL3fL1DoubleMu*EG*Filter*","hltDiMu*Ele*CaloIdLTrackIdLElectronleg*Filter"],"1e-2mu"),
72  mksel(("hltEle32L1DoubleEGWPTightGsfTrackIsoFilter","hltEGL1SingleEGOrFilter"),"1e (32_L1DoubleEG_AND_L1SingleEGOr)"),
73  mksel("filter('hltEle*CaloIdVTGsfTrkIdTGsfDphiFilter')","1e (CaloIdVT_GsfTrkIdT)"),
74  mksel("path('HLT_Ele*PFJet*')","1e (PFJet)"),
75  mksel(["hltEG175HEFilter","hltEG200HEFilter"],"1e (Photon175_OR_Photon200)"),
76  )
77  ),
78  Photon = cms.PSet(
79  id = cms.int32(22),
80  sel = cms.string("type(92) && pt > 15 && coll('hltEgammaCandidates')"),
81  l1seed = cms.string("type(-98)"), l1deltaR = cms.double(0.3),
82  #l2seed = cms.string("type(92) && coll('')"), l2deltaR = cms.double(0.5),
83  skipObjectsNotPassingQualityBits = cms.bool(True),
84  qualityBits = cms.VPSet(
85  mksel("filter('hltEG33L1EG26HEFilter')","hltEG33L1EG26HEFilter"),
86  mksel("filter('hltEG50HEFilter')","hltEG50HEFilter"),
87  mksel("filter('hltEG75HEFilter')","hltEG75HEFilter"),
88  mksel("filter('hltEG90HEFilter')","hltEG90HEFilter"),
89  mksel("filter('hltEG120HEFilter')","hltEG120HEFilter"),
90  mksel("filter('hltEG150HEFilter')","hltEG150HEFilter"),
91  mksel("filter('hltEG175HEFilter')","hltEG150HEFilter"),
92  mksel("filter('hltEG200HEFilter')","hltEG200HEFilter"),
93  mksel("filter('hltHtEcal800')","hltHtEcal800"),
94  mksel("filter('hltEG110EBTightIDTightIsoTrackIsoFilter')","hltEG110EBTightIDTightIsoTrackIsoFilter"),
95  mksel("filter('hltEG120EBTightIDTightIsoTrackIsoFilter')","hltEG120EBTightIDTightIsoTrackIsoFilter"),
96  mksel("filter('hltMu17Photon30IsoCaloIdPhotonlegTrackIsoFilter')","1mu-1photon"),
97  mksel("filter('hltEG30LR9Id85b90eHE12R9Id50b80eR9IdLastFilter')","hltEG30LR9Id85b90eHE12R9Id50b80eR9IdLastFilter"),
98  mksel("filter('hltEG30LIso60CaloId15b35eHE12R9Id50b80eEcalIsoLastFilter')","hltEG30LIso60CaloId15b35eHE12R9Id50b80eEcalIsoLastFilter"),
99  mksel("filter('hltEG22Iso60CaloId15b35eHE12R9Id50b80eTrackIsoUnseededLastFilter')","hltEG22Iso60CaloId15b35eHE12R9Id50b80eTrackIsoUnseededLastFilter"),
100  mksel("filter('hltEG22R9Id85b90eHE12R9Id50b80eR9UnseededLastFilter')","hltEG22R9Id85b90eHE12R9Id50b80eR9UnseededLastFilter"),
101  mksel("filter('hltEG30Iso60CaloId15b35eR9Id50b90eHE12b10eR9Id50b80eEcalIsoFilter')","hltEG30Iso60CaloId15b35eR9Id50b90eHE12b10eR9Id50b80eEcalIsoFilter"),
102  mksel("filter('hltEG18TrackIso60Iso60CaloId15b35eR9Id50b90eHE12b10eR9Id50b80eTrackIsoUnseededFilter')","hltEG18TrackIso60Iso60CaloId15b35eR9Id50b90eHE12b10eR9Id50b80eTrackIsoUnseededFilter")
103  )
104  ),
105  Muon = cms.PSet(
106  id = cms.int32(13),
107  sel = cms.string("type(83) && pt > 5 && (coll('hltIterL3MuonCandidates') || (pt > 45 && coll('hltHighPtTkMuonCands')) || (pt > 95 && coll('hltOldL3MuonCandidates')))"),
108  l1seed = cms.string("type(-81)"), l1deltaR = cms.double(0.5),
109  l2seed = cms.string("type(83) && coll('hltL2MuonCandidates')"), l2deltaR = cms.double(0.3),
110  skipObjectsNotPassingQualityBits = cms.bool(True),
111  qualityBits = cms.VPSet(
112  mksel(["*RelTrkIsoVVLFiltered0p4","*RelTrkIsoVVLFiltered"],"TrkIsoVVL"),
113  mksel(["hltL3crIso*IsoFiltered0p07","hltL3crIso*IsoFiltered0p08","hltL3crIso*IsoFiltered"],"Iso"),
114  mksel("filter('*OverlapFilterIsoMu*PFTau*')","OverlapFilter PFTau"),
115  mksel(["hltL3crIsoL1*SingleMu*IsoFiltered0p07","hltL3crIsoL1sMu*IsoFiltered0p07","hltL3crIsoL1*SingleMu*IsoFiltered0p08","hltL3crIsoL1sMu*IsoFiltered0p08","hltL3crIsoL1*SingleMu*IsoFiltered","hltL3crIsoL1sMu*IsoFiltered"],"1mu"),
116  mksel("filter('hltDiMuon*Filtered*')","2mu"),
117  mksel("filter('hltMu*TrkIsoVVL*Ele*CaloIdLTrackIdLIsoVL*Filter*')","1mu-1e"),
118  mksel("filter('hlt*OverlapFilterIsoMu*PFTau*')","1mu-1tau"),
119  mksel("filter('hltL3fL1TripleMu*')","3mu"),
120  mksel(["hltL3fL1DoubleMu*EG*Filtered*","hltDiMu*Ele*CaloIdLTrackIdLElectronleg*Filter"],"2mu-1e"),
121  mksel(["hltL3fL1Mu*DoubleEG*Filtered*","hltMu*DiEle*CaloIdLTrackIdLElectronleg*Filter"],"1mu-2e"),
122  mksel(["hltL3fL1sMu*L3Filtered50*","hltL3fL1sMu*TkFiltered50*"],"1mu (Mu50)"),
123  mksel(["hltL3fL1sMu*L3Filtered100*","hltL3fL1sMu*TkFiltered100*"],"1mu (Mu100)"),
124  mksel("filter('hltMu17Photon30IsoCaloIdMuonlegL3Filtered17Q')","1mu-1photon")
125  )
126  ),
127  Tau = cms.PSet(
128  id = cms.int32(15),
129  sel = cms.string("type(84) && pt > 5 && coll('*Tau*') && ( filter('*LooseChargedIso*') || filter('*MediumChargedIso*') || filter('*DeepTau*') || filter('*TightChargedIso*') || filter('*TightOOSCPhotons*') || filter('hltL2TauIsoFilter') || filter('*OverlapFilterIsoMu*') || filter('*OverlapFilterIsoEle*') || filter('*L1HLTMatched*') || filter('*Dz02*') || filter('*DoublePFTau*') || filter('*SinglePFTau*') || filter('hlt*SelectedPFTau') || filter('*DisplPFTau*') )"), #All trigger objects from a Tau collection + passing at least one filter
130  l1seed = cms.string("type(-100)"), l1deltaR = cms.double(0.3),
131  l2seed = cms.string("type(84) && coll('hltL2TauJetsL1IsoTauSeeded')"), l2deltaR = cms.double(0.3),
132  skipObjectsNotPassingQualityBits = cms.bool(True),
133  qualityBits = cms.VPSet(
134  mksel("filter('*LooseChargedIso*')","LooseChargedIso"),
135  mksel("filter('*MediumChargedIso*')","MediumChargedIso"),
136  mksel("filter('*TightChargedIso*')","TightChargedIso"),
137  mksel("filter('*DeepTau*')","DeepTau"),
138  mksel("filter('*TightOOSCPhotons*')","TightID OOSC photons"),
139  mksel("filter('*Hps*')","HPS"),
140  mksel("filter('hlt*DoublePFTau*TrackPt1*ChargedIsolation*Dz02*')","charged iso di-tau"),
141  mksel("filter('hlt*DoublePFTau*DeepTau*L1HLTMatched')","deeptau di-tau"),
142  mksel("filter('hlt*OverlapFilterIsoEle*WPTightGsf*PFTau*')","e-tau"),
143  mksel("filter('hlt*OverlapFilterIsoMu*PFTau*')","mu-tau"),
144  mksel("filter('hlt*SelectedPFTau*L1HLTMatched')","single-tau/tau+MET"),
145  mksel("filter('hlt*DoublePFTau*TrackPt1*ChargedIso*')","run 2 VBF+ditau"),
146  mksel("filter('hlt*DoublePFTau*Track*ChargedIso*AgainstMuon')","run 3 VBF+ditau"),
147  mksel("filter('hltHpsSinglePFTau*HLTMatched')","run 3 double PF jets + ditau"),
148  mksel("filter('hltHpsOverlapFilterDeepTauDoublePFTau*PFJet*')","di-tau + PFJet"),
149  mksel("filter('hlt*Double*ChargedIsoDisplPFTau*Dxy*')","Displaced Tau"),
150  mksel("filter('*Monitoring')","Monitoring"),
151  mksel("filter('*Reg')","regional paths"),
152  mksel("filter('*L1Seeded')","L1 seeded paths"),
153  mksel("filter('*1Prong')","1 prong tau paths")
154  )
155  ),
156  BoostedTau = cms.PSet(
157  id = cms.int32(1515),
158  sel = cms.string("type(85) && pt > 120 && coll('hltAK8PFJetsCorrected') && filter('hltAK8SinglePFJets*SoftDropMass40*ParticleNetTauTau')"),
159  l1seed = cms.string("type(-99)"), l1deltaR = cms.double(0.3),
160  l2seed = cms.string("type(85) && coll('hltAK8CaloJetsCorrectedIDPassed')"), l2deltaR = cms.double(0.3),
161  skipObjectsNotPassingQualityBits = cms.bool(True),
162  qualityBits = cms.VPSet(
163  mksel("filter('hltAK8SinglePFJets*SoftDropMass40*ParticleNetTauTau')","HLT_AK8PFJetX_SoftDropMass40_PFAK8ParticleNetTauTau0p30")
164  )
165  ),
166  Jet = cms.PSet(
167  id = cms.int32(1),
168  sel = cms.string("( type(0) || type(85) || type(86) || type(-99) )"),
169  l1seed = cms.string("type(-99)"), l1deltaR = cms.double(0.3),
170  l2seed = cms.string("type(85) || type(86) || type(-99)"), l2deltaR = cms.double(0.3),
171  skipObjectsNotPassingQualityBits = cms.bool(True),
172  qualityBits = cms.VPSet(
173  mksel(["hlt4PixelOnlyPFCentralJetTightIDPt20"]), # 0
174  mksel(["hlt3PixelOnlyPFCentralJetTightIDPt30"]), # 1
175  mksel(["hltPFJetFilterTwoC30"]), # 2
176  mksel(["hlt4PFCentralJetTightIDPt30"]), # 3
177  mksel(["hlt4PFCentralJetTightIDPt35"]), # 4
178  mksel(["hltQuadCentralJet30"]), # 5
179  mksel(["hlt2PixelOnlyPFCentralJetTightIDPt40"]), # 6
180  mksel(["hltL1sTripleJet1008572VBFIorHTTIorDoubleJetCIorSingleJet","hltL1sTripleJet1058576VBFIorHTTIorDoubleJetCIorSingleJet","hltL1sTripleJetVBFIorHTTIorSingleJet"]), # 7
181  mksel(["hlt3PFCentralJetTightIDPt40"]), # 8
182  mksel(["hlt3PFCentralJetTightIDPt45"]), # 9
183  mksel(["hltL1sQuadJetC60IorHTT380IorHTT280QuadJetIorHTT300QuadJet","hltL1sQuadJetC50to60IorHTT280to500IorHTT250to340QuadJet"]), # 10
184  mksel(["hltBTagCaloDeepCSVp17Double"]), # 11
185  mksel(["hltPFCentralJetLooseIDQuad30"]), # 12
186  mksel(["hlt1PFCentralJetLooseID75"]), # 13
187  mksel(["hlt2PFCentralJetLooseID60"]), # 14
188  mksel(["hlt3PFCentralJetLooseID45"]), # 15
189  mksel(["hlt4PFCentralJetLooseID40"]), # 16
190  mksel(["hltBTagPFDeepCSV4p5Triple"]), # 17
191  mksel("filter('hltHpsOverlapFilterDeepTauDoublePFTau*PFJet*')","(Double tau + jet) hltHpsOverlapFilterDeepTauDoublePFTau*PFJet*"), # 18
192  mksel("filter('*CrossCleaned*MediumDeepTauDitauWPPFTau*')","(VBF cross-cleaned from medium deeptau PFTau) *CrossCleaned*MediumDeepTauDitauWPPFTau*"), # 19
193  mksel("filter('*CrossCleanedUsingDiJetCorrChecker*')","(VBF cross-cleaned using dijet correlation checker) *CrossCleanedUsingDiJetCorrChecker*"), # 20
194  mksel("filter('hltHpsOverlapFilterDeepTauPFTau*PFJet*')","(monitoring muon + tau + jet) hltHpsOverlapFilterDeepTauPFTau*PFJet*"), # 21
195  mksel(["hlt2PFCentralJetTightIDPt50"]), # 22
196  mksel(["hlt1PixelOnlyPFCentralJetTightIDPt60"]), # 23
197  mksel(["hlt1PFCentralJetTightIDPt70"]), # 24
198  mksel(["hltBTagPFDeepJet1p5Single"]), # 25
199  mksel(["hltBTagPFDeepJet4p5Triple"]), # 26
200  mksel(["hltBTagCentralJetPt35PFParticleNet2BTagSum0p65","hltBTagCentralJetPt30PFParticleNet2BTagSum0p65","hltPFJetTwoC30PFBTagParticleNet2BTagSum0p65"]), # 27
201  mksel(["hltBTagPFDeepCSV1p5Single"]) # 28
202  ),
203  ),
204  FatJet = cms.PSet(
205  id = cms.int32(6),
206  sel = cms.string("type(85) && pt > 120 && coll('hltAK8PFJetsCorrected')"),
207  l1seed = cms.string("type(-99)"), l1deltaR = cms.double(0.3),
208  l2seed = cms.string("type(85) && coll('hltAK8CaloJetsCorrectedIDPassed')"), l2deltaR = cms.double(0.3),
209  skipObjectsNotPassingQualityBits = cms.bool(True),
210  qualityBits = cms.VPSet()
211  ),
212  MET = cms.PSet(
213  id = cms.int32(2),
214  sel = cms.string("type(87) && pt > 30 && coll('hltPFMETProducer')"),
215  l1seed = cms.string("type(-87) && coll('L1ETM')"), l1deltaR = cms.double(9999),
216  l1seed_2 = cms.string("type(-87) && coll('L1ETMHF')"), l1deltaR_2 = cms.double(9999),
217  l2seed = cms.string("type( 87) && coll('hltMetClean')"), l2deltaR = cms.double(9999),
218  skipObjectsNotPassingQualityBits = cms.bool(True),
219  qualityBits = cms.VPSet()
220  ),
221  HT = cms.PSet(
222  id = cms.int32(3),
223  sel = cms.string("type(89) || type(-89)"),
224  l1seed = cms.string("type(-89) && coll('L1HTT')"), l1deltaR = cms.double(9999),
225  l1seed_2 = cms.string("type(-89) && coll('L1HTTHF')"), l1deltaR_2 = cms.double(9999),
226  l2seed = cms.string("type(89) && coll('hltHtMhtJet30')"), l2deltaR = cms.double(9999),
227  skipObjectsNotPassingQualityBits = cms.bool(True),
228  qualityBits = cms.VPSet(
229  mksel(["hltL1sTripleJetVBFIorHTTIorDoubleJetCIorSingleJet"]),
230  mksel(["hltL1sQuadJetC50IorQuadJetC60IorHTT280IorHTT300IorHTT320IorTripleJet846848VBFIorTripleJet887256VBFIorTripleJet927664VBF","hltL1sQuadJetCIorTripleJetVBFIorHTT"]),
231  mksel(["hltL1sQuadJetC60IorHTT380IorHTT280QuadJetIorHTT300QuadJet","hltL1sQuadJetC50to60IorHTT280to500IorHTT250to340QuadJet"]),
232  mksel(["hltCaloQuadJet30HT300","hltCaloQuadJet30HT320"]),
233  mksel(["hltPFCentralJetsLooseIDQuad30HT300","hltPFCentralJetsLooseIDQuad30HT330"])
234  ),
235  ),
236  MHT = cms.PSet(
237  id = cms.int32(4),
238  sel = cms.string("type(90)"),
239  l1seed = cms.string("type(-90) && coll('L1HTM')"), l1deltaR = cms.double(9999),
240  l1seed_2 = cms.string("type(-90) && coll('L1HTMHF')"), l1deltaR_2 = cms.double(9999),
241  l2seed = cms.string("type(90) && coll('hltHtMhtJet30')"), l2deltaR = cms.double(9999),
242  skipObjectsNotPassingQualityBits = cms.bool(True),
243  qualityBits = cms.VPSet(
244  mksel(["hltCaloQuadJet30HT300","hltCaloQuadJet30HT320"]),
245  mksel(["hltPFCentralJetsLooseIDQuad30HT300","hltPFCentralJetsLooseIDQuad30HT330"])
246  ),
247  ),
248  ),
249 )
250 
251 # ERA-dependent configuration
252 # Tune filter and collection names to 2016 HLT menus
253 # FIXME: check non-lepton objects and cross check leptons
254 run2_HLTconditions_2016.toModify(
255  triggerObjectTable.selections.Muon,
256  sel = "type(83) && pt > 5 && (coll('hlt*L3MuonCandidates') || coll('hlt*TkMuonCands') || coll('hlt*TrkMuonCands'))",
257  qualityBits = cms.VPSet(
258  mksel("filter('*RelTrkIso*Filtered0p4')","TrkIsoVVL"),
259  mksel("filter('hltL3cr*IsoFiltered0p09')","Iso"),
260  mksel("filter('*OverlapFilter*IsoMu*PFTau*')","OverlapFilter PFTau"),
261  mksel("filter('hltL3f*IsoFiltered0p09')","IsoTkMu"),
262  mksel(["hltL3fL1sMu*L3Filtered50*","hltL3fL1sMu*TkFiltered50*"],"1mu (Mu50)", bit=10)
263  )
264 ).toModify(
265  triggerObjectTable.selections.Tau,
266  sel = "type(84) && pt > 5 && coll('*Tau*') && ( filter('*LooseIso*') || filter('*MediumIso*') || filter('*MediumComb*Iso*') || filter('hltL2TauIsoFilter') || filter('*OverlapFilter*IsoMu*') || filter('*OverlapFilter*IsoEle*') || filter('*L1HLTMatched*') || filter('*Dz02*') )",
267  qualityBits = cms.VPSet(
268  mksel("(filter('*LooseIso*')-filter('*VLooseIso*'))","LooseIso"),
269  mksel("filter('*Medium*Iso*')","Medium(Comb)Iso"),
270  mksel("filter('*VLooseIso*')","VLooseIso"),
271  mksel("0","None"),
272  mksel("filter('hltL2TauIsoFilter')","L2p5 pixel iso"),
273  mksel("filter('*OverlapFilter*IsoMu*')","OverlapFilter IsoMu"),
274  mksel("filter('*OverlapFilter*IsoEle*')","OverlapFilter IsoEle"),
275  mksel("filter('*L1HLTMatched*')","L1-HLT matched"),
276  mksel("filter('*Dz02*')","Dz")
277  )
278 )
279 
280 _run2_HLTconditions = run2_HLTconditions_2016 | run2_HLTconditions_2017 | run2_HLTconditions_2018
281 
282 _run2_2016_jet_filters = [
283  mksel(["*CrossCleaned*LooseChargedIsoPFTau*"], "VBF cross-cleaned from loose iso PFTau"), # 0
284  mksel(["hltBTagCaloCSVp087Triple"], "hltBTagCaloCSVp087Triple"), # 1
285  mksel(["hltDoubleCentralJet90"], "hltDoubleCentralJet90"), # 2
286  mksel(["hltDoublePFCentralJetLooseID90"], "hltDoublePFCentralJetLooseID90"), # 3
287  mksel(["hltL1sTripleJetVBFIorHTTIorDoubleJetCIorSingleJet"], "hltL1sTripleJetVBFIorHTTIorDoubleJetCIorSingleJet"), # 4
288  mksel(["hltQuadCentralJet30"], "hltQuadCentralJet30"), # 5
289  mksel(["hltQuadPFCentralJetLooseID30"], "hltQuadPFCentralJetLooseID30"), # 6
290  mksel(["hltL1sQuadJetC50IorQuadJetC60IorHTT280IorHTT300IorHTT320IorTripleJet846848VBFIorTripleJet887256VBFIorTripleJet927664VBF", "hltL1sQuadJetCIorTripleJetVBFIorHTT"], "hltL1sQuadJetC50IorQuadJetC60IorHTT280IorHTT300IorHTT320IorTripleJet846848VBFIorTripleJet887256VBFIorTripleJet927664VBF or hltL1sQuadJetCIorTripleJetVBFIorHTT"), # 7
291  mksel(["hltQuadCentralJet45"], "hltQuadCentralJet45"), # 8
292  mksel(["hltQuadPFCentralJetLooseID45"], "hltQuadPFCentralJetLooseID45"), # 9
293  mksel(["hltL1sQuadJetC60IorHTT380IorHTT280QuadJetIorHTT300QuadJet", "hltL1sQuadJetC50to60IorHTT280to500IorHTT250to340QuadJet"], "hltL1sQuadJetC60IorHTT380IorHTT280QuadJetIorHTT300QuadJet or hltL1sQuadJetC50to60IorHTT280to500IorHTT250to340QuadJet"), # 10
294  mksel(["hltBTagCaloCSVp05Double", "hltBTagCaloDeepCSVp17Double"], "hltBTagCaloCSVp05Double or hltBTagCaloDeepCSVp17Double"), # 11
295  mksel(["hltPFCentralJetLooseIDQuad30"], "hltPFCentralJetLooseIDQuad30"), # 12
296  mksel(["hlt1PFCentralJetLooseID75"], "hlt1PFCentralJetLooseID75"), # 13
297  mksel(["hlt2PFCentralJetLooseID60"], "hlt2PFCentralJetLooseID60"), # 14
298  mksel(["hlt3PFCentralJetLooseID45"], "hlt3PFCentralJetLooseID45"), # 15
299  mksel(["hlt4PFCentralJetLooseID40"], "hlt4PFCentralJetLooseID40"), # 16
300  mksel(["hltBTagPFCSVp070Triple", "hltBTagPFDeepCSVp24Triple", "hltBTagPFDeepCSV4p5Triple"], "hltBTagPFCSVp070Triple or hltBTagPFDeepCSVp24Triple or hltBTagPFDeepCSV4p5Triple"), # 17
301  mksel(["hltHpsOverlapFilterDeepTauDoublePFTau*PFJet*"], "Double tau + jet"), # 18
302  mksel(["*CrossCleaned*MediumDeepTauDitauWPPFTau*"], "VBF cross-cleaned from medium deeptau PFTau"), # 19
303  mksel(["*CrossCleanedUsingDiJetCorrChecker*"], "VBF cross-cleaned using dijet correlation checker"), # 20
304  mksel(["hltHpsOverlapFilterDeepTauPFTau*PFJet*"], "monitoring muon + tau + jet"), # 21
305 ]
306 run2_HLTconditions_2016.toModify(triggerObjectTable.selections.Jet, qualityBits = cms.VPSet(_run2_2016_jet_filters))
307 
308 _run2_2017_jet_filters = copy.deepcopy(_run2_2016_jet_filters)
309 _run2_2017_jet_filters[7] = mksel(["hltL1sTripleJet1008572VBFIorHTTIorDoubleJetCIorSingleJet","hltL1sTripleJet1058576VBFIorHTTIorDoubleJetCIorSingleJet","hltL1sTripleJetVBFIorHTTIorSingleJet"])
310 run2_HLTconditions_2017.toModify(triggerObjectTable.selections.Jet, qualityBits = cms.VPSet(_run2_2017_jet_filters))
311 
312 _run2_2018_jet_filters = copy.deepcopy(_run2_2017_jet_filters)
313 _run2_2018_jet_filters[2] = mksel(["hltPFJetFilterTwoC30"])
314 _run2_2018_jet_filters.append(mksel(["hltBTagPFDeepCSV1p5Single"])) # 22
315 run2_HLTconditions_2018.toModify(triggerObjectTable.selections.Jet, qualityBits = cms.VPSet(_run2_2018_jet_filters))
316 
317 
318 from PhysicsTools.PatUtils.L1PrefiringWeightProducer_cff import prefiringweight
319 #Next lines are for UL2016 maps
320 (run2_muon_2016 & tracker_apv_vfp30_2016).toModify(
321  prefiringweight,
322  DataEraECAL = cms.string("UL2016preVFP"),
323  DataEraMuon = cms.string("2016preVFP")
324 )
325 (run2_muon_2016 & ~tracker_apv_vfp30_2016).toModify(
326  prefiringweight,
327  DataEraECAL = cms.string("UL2016postVFP"),
328  DataEraMuon = cms.string("2016postVFP")
329 )
330 #Next line is for UL2017 maps
331 run2_jme_2017.toModify(
332  prefiringweight,
333  DataEraECAL = cms.string("UL2017BtoF"),
334  DataEraMuon = cms.string("20172018")
335 )
336 #Next line is for UL2018 maps
337 run2_muon_2018.toModify(
338  prefiringweight,
339  DataEraECAL = cms.string("None"),
340  DataEraMuon = cms.string("20172018")
341 )
342 
343 l1PreFiringEventWeightTable = globalVariablesTableProducer.clone(
344  name = cms.string("L1PreFiringWeight"),
345  variables = cms.PSet(
346  Nom = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProb"), "float", doc = "L1 pre-firing event correction weight (1-probability)", precision=8),
347  Up = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbUp"), "float", doc = "L1 pre-firing event correction weight (1-probability), up var.", precision=8),
348  Dn = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbDown"), "float", doc = "L1 pre-firing event correction weight (1-probability), down var.", precision=8),
349  Muon_Nom = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbMuon"), "float", doc = "Muon L1 pre-firing event correction weight (1-probability)", precision=8),
350  Muon_SystUp = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbMuonSystUp"), "float", doc = "Muon L1 pre-firing event correction weight (1-probability), up var. syst.", precision=8),
351  Muon_SystDn = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbMuonSystDown"), "float", doc = "Muon L1 pre-firing event correction weight (1-probability), down var. syst.", precision=8),
352  Muon_StatUp = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbMuonStatUp"), "float", doc = "Muon L1 pre-firing event correction weight (1-probability), up var. stat.", precision=8),
353  Muon_StatDn = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbMuonStatDown"), "float", doc = "Muon L1 pre-firing event correction weight (1-probability), down var. stat.", precision=8),
354  ECAL_Nom = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbECAL"), "float", doc = "ECAL L1 pre-firing event correction weight (1-probability)", precision=8),
355  ECAL_Up = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbECALUp"), "float", doc = "ECAL L1 pre-firing event correction weight (1-probability), up var.", precision=8),
356  ECAL_Dn = ExtVar(cms.InputTag("prefiringweight:nonPrefiringProbECALDown"), "float", doc = "ECAL L1 pre-firing event correction weight (1-probability), down var.", precision=8),
357  )
358 )
359 
360 l1bits=cms.EDProducer("L1TriggerResultsConverter",
361  src=cms.InputTag("gtStage2Digis"),
362  legacyL1=cms.bool(False),
363  storeUnprefireableBit=cms.bool(True),
364  src_ext=cms.InputTag("simGtExtUnprefireable"))
365 
366 triggerObjectTablesTask = cms.Task( unpackedPatTrigger,triggerObjectTable,l1bits)
367 
368 _run2_HLTconditions.toReplaceWith(
369  triggerObjectTablesTask, triggerObjectTablesTask.copyAndAdd(prefiringweight,l1PreFiringEventWeightTable)
370 )
def ExtVar(tag, valtype, doc=None, precision=-1)
Definition: common_cff.py:27
def chaintoken_(tokens, f, OR_or_AND=None)
def mksel(selection, doc=None, bit=None)
static std::string join(char **cmd)
Definition: RemoteFile.cc:19