Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_noPU_HFnose_v0.tcl@ 861ad5a

ImprovedOutputFile Timing
Last change on this file since 861ad5a was 190cfa0, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago

first attempt at describing HFnose in Delphes
  • Property mode set to 100644
File size: 57.6 KB
Line 
1
2###########################################
3#
4# Main authors: Michele Selvaggi (CERN)
5#
6# Released on: May 2018
7#
8# Based on CMS_PhaseII_200PU_v03
9#
10# Test card for CMS PhaseII detector without pile-up (use for testing only)
11#
12############################################
13
14
15#######################################
16# Order of execution of various modules
17#######################################
18
19set ExecutionPath {
20 ParticlePropagator
21
22 ChargedHadronTrackingEfficiency
23 ElectronTrackingEfficiency
24 MuonTrackingEfficiency
25
26 ChargedHadronMomentumSmearing
27 ElectronEnergySmearing
28 MuonMomentumSmearing
29
30 TrackMerger
31
32 ECal
33 HCal
34
35 PhotonEnergySmearing
36 ElectronFilter
37
38 Calorimeter
39 EFlowMerger
40 EFlowFilter
41
42 PhotonEfficiency
43 PhotonIsolation
44
45 ElectronFilter
46 ElectronEfficiency
47 ElectronIsolation
48
49 ChargedHadronFilter
50
51 MuonEfficiency
52 MuonIsolation
53
54 MissingET
55
56 NeutrinoFilter
57 GenJetFinder
58 GenMissingET
59
60 FastJetFinder
61 FatJetFinder
62
63 JetEnergyScale
64
65 JetFlavorAssociation
66
67 BTagging
68 TauTagging
69
70 UniqueObjectFinder
71
72 ScalarHT
73
74 TreeWriter
75}
76
77#################################
78# Propagate particles in cylinder
79#################################
80
81module ParticlePropagator ParticlePropagator {
82 set InputArray Delphes/stableParticles
83
84 set OutputArray stableParticles
85 set ChargedHadronOutputArray chargedHadrons
86 set ElectronOutputArray electrons
87 set MuonOutputArray muons
88
89 # radius of the magnetic field coverage, in m
90 set Radius 1.29
91 # half-length of the magnetic field coverage, in m
92 set HalfLength 3.00
93
94 # magnetic field
95 set Bz 3.8
96}
97
98
99####################################
100# Charged hadron tracking efficiency
101####################################
102
103module Efficiency ChargedHadronTrackingEfficiency {
104 ## particles after propagation
105 set InputArray ParticlePropagator/chargedHadrons
106 set OutputArray chargedHadrons
107 # tracking efficiency formula for charged hadrons
108 set EfficiencyFormula {
109 (pt <= 0.2) * (0.00) + \
110 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
111 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
112 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
113 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
114 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
115 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
116 (abs(eta) > 4.0) * (0.00)
117 }
118}
119
120
121#####################################
122# Electron tracking efficiency - ID
123####################################
124
125module Efficiency ElectronTrackingEfficiency {
126 set InputArray ParticlePropagator/electrons
127 set OutputArray electrons
128 # tracking efficiency formula for electrons
129 set EfficiencyFormula {
130 (pt <= 0.2) * (0.00) + \
131 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
132 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
133 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
134 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
135 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
136 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
137 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
138 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
139 (abs(eta) > 4.0) * (0.00)
140
141 }
142}
143
144##########################
145# Muon tracking efficiency
146##########################
147
148module Efficiency MuonTrackingEfficiency {
149 set InputArray ParticlePropagator/muons
150 set OutputArray muons
151 # tracking efficiency formula for muons
152 set EfficiencyFormula {
153 (pt <= 0.2) * (0.00) + \
154 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
155 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
156 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
157 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
158 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
159 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
160 (abs(eta) > 4.0) * (0.00)
161
162 }
163}
164
165
166########################################
167# Momentum resolution for charged tracks
168########################################
169
170module MomentumSmearing ChargedHadronMomentumSmearing {
171 ## hadrons after having applied the tracking efficiency
172 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
173 set OutputArray chargedHadrons
174 # resolution formula for charged hadrons ,
175
176 #
177 # Automatically generated tracker resolution formula for layout: OT612IT4025
178 #
179 # By Unknown author on: 2017-06-30.17:03:00
180 #
181 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
182 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
183 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
184 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
185 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
186 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
187 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
188 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
189 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
190 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
191 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
192 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
193 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
194 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
195 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
196 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
197 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
198 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
199 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
200 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
201 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
202 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
203 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
204 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
205 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
206 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
207 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
208 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
209 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
210 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
211 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
212 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
213 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
214 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
215 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
216 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
217 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
218 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
219 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
220 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
221 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
222 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
223 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
224 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
225 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
226 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
227 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
228 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
229 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
230 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
231 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
232 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
233 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
234 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
235 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
236 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
237 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
238 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
239 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
240 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
241 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
242 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
243 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
244 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
245 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
246 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
247 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
248 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
249 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
250 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
251 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
252 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
253 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
254 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
255 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
256 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
257 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
258 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
259 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
260 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000)
261 }
262
263
264}
265
266
267#################################
268# Energy resolution for electrons
269#################################
270
271module EnergySmearing ElectronEnergySmearing {
272 set InputArray ElectronTrackingEfficiency/electrons
273 set OutputArray electrons
274
275 # set ResolutionFormula {resolution formula as a function of eta and energy}
276
277 # resolution formula for electrons
278
279 # taking something flat in energy for now, ECAL will take over at high energy anyway.
280 # inferred from hep-ex/1306.2016 and 1502.02701
281 set ResolutionFormula {
282
283 (abs(eta) <= 1.5) * (energy*0.028) +
284 (abs(eta) > 1.5 && abs(eta) <= 1.75) * (energy*0.037) +
285 (abs(eta) > 1.75 && abs(eta) <= 2.15) * (energy*0.038) +
286 (abs(eta) > 2.15 && abs(eta) <= 3.00) * (energy*0.044) +
287 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)}
288
289}
290
291###############################
292# Momentum resolution for muons
293###############################
294
295module MomentumSmearing MuonMomentumSmearing {
296 set InputArray MuonTrackingEfficiency/muons
297 set OutputArray muons
298 # resolution formula for muons
299
300 # up to |eta| < 2.8 take measurement from tracking + muon chambers
301 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
302 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
303 source muonMomentumResolution.tcl
304}
305
306
307##############
308# Track merger
309##############
310
311module Merger TrackMerger {
312# add InputArray InputArray
313 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
314 add InputArray ElectronEnergySmearing/electrons
315 add InputArray MuonMomentumSmearing/muons
316 set OutputArray tracks
317}
318
319
320#############
321# ECAL
322#############
323
324module SimpleCalorimeter ECal {
325 set ParticleInputArray ParticlePropagator/stableParticles
326 set TrackInputArray TrackMerger/tracks
327
328 set TowerOutputArray ecalTowers
329 set EFlowTrackOutputArray eflowTracks
330 set EFlowTowerOutputArray eflowPhotons
331
332 set IsEcal true
333
334 set EnergyMin 0.5
335 set EnergySignificanceMin 1.0
336
337 set SmearTowerCenter true
338
339 set pi [expr {acos(-1)}]
340
341 # lists of the edges of each tower in eta and phi
342 # each list starts with the lower edge of the first tower
343 # the list ends with the higher edged of the last tower
344
345 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
346
347 set PhiBins {}
348 for {set i -180} {$i <= 180} {incr i} {
349 add PhiBins [expr {$i * $pi/180.0}]
350 }
351
352 # 0.02 unit in eta up to eta = 1.5 (barrel)
353 for {set i -85} {$i <= 86} {incr i} {
354 set eta [expr {$i * 0.0174}]
355 add EtaPhiBins $eta $PhiBins
356 }
357
358 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 4.175 (HGCAL- ECAL)
359
360 set PhiBins {}
361 for {set i -180} {$i <= 180} {incr i} {
362 add PhiBins [expr {$i * $pi/180.0}]
363 }
364
365 # 0.02 unit in eta up to eta = 3
366 for {set i 1} {$i <= 154} {incr i} {
367 set eta [expr { -4.176 + $i * 0.0174}]
368 add EtaPhiBins $eta $PhiBins
369 }
370
371 for {set i 1} {$i <= 154} {incr i} {
372 set eta [expr { 1.4964 + $i * 0.0174}]
373 add EtaPhiBins $eta $PhiBins
374 }
375
376 # remaining HF part
377 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 4.2 < |eta| < 5.0
378
379 set PhiBins {}
380 for {set i -18} {$i <= 18} {incr i} {
381 add PhiBins [expr {$i * $pi/18.0}]
382 }
383
384 foreach eta {-5 -4.7 -4.525 -4.35 -4.176 4.35 4.525 4.7 5} {
385 add EtaPhiBins $eta $PhiBins
386 }
387
388 add EnergyFraction {0} {0.0}
389 # energy fractions for e, gamma and pi0
390 add EnergyFraction {11} {1.0}
391 add EnergyFraction {22} {1.0}
392 add EnergyFraction {111} {1.0}
393 # energy fractions for muon, neutrinos and neutralinos
394 add EnergyFraction {12} {0.0}
395 add EnergyFraction {13} {0.0}
396 add EnergyFraction {14} {0.0}
397 add EnergyFraction {16} {0.0}
398 add EnergyFraction {1000022} {0.0}
399 add EnergyFraction {1000023} {0.0}
400 add EnergyFraction {1000025} {0.0}
401 add EnergyFraction {1000035} {0.0}
402 add EnergyFraction {1000045} {0.0}
403 # energy fractions for K0short and Lambda
404 add EnergyFraction {310} {0.3}
405 add EnergyFraction {3122} {0.3}
406
407 # set ResolutionFormula {resolution formula as a function of eta and energy}
408
409 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
410 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
411
412 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
413 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
414 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
415 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
416 (abs(eta) >= 3.0 && abs(eta) <= 4.176) * sqrt(energy^2*0.01^2 + energy*0.67^2) + \
417 (abs(eta) >= 4.176 && abs(eta) <= 5.0) * sqrt(energy^2*0.10^2 + energy*1.82^2)}
418
419}
420
421#############
422# HCAL
423#############
424
425module SimpleCalorimeter HCal {
426 set ParticleInputArray ParticlePropagator/stableParticles
427 set TrackInputArray ECal/eflowTracks
428
429 set TowerOutputArray hcalTowers
430 set EFlowTrackOutputArray eflowTracks
431 set EFlowTowerOutputArray eflowNeutralHadrons
432
433 set IsEcal false
434
435 set EnergyMin 1.0
436 set EnergySignificanceMin 1.0
437
438 set SmearTowerCenter true
439
440 set pi [expr {acos(-1)}]
441
442 # lists of the edges of each tower in eta and phi
443 # each list starts with the lower edge of the first tower
444 # the list ends with the higher edged of the last tower
445
446 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
447
448 set PhiBins {}
449 for {set i -36} {$i <= 36} {incr i} {
450 add PhiBins [expr {$i * $pi/36.0}]
451 }
452 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.65} {
453 add EtaPhiBins $eta $PhiBins
454 }
455
456 # assume 0.07 x 0.07 resolution in eta,phi in the endcaps 1.5 < |eta| < 4.2 (HGCAL- HCAL)
457
458 set PhiBins {}
459 for {set i -45} {$i <= 45} {incr i} {
460 add PhiBins [expr {$i * $pi/45.0}]
461 }
462
463 # 0.07 unit in eta up to eta = 3
464 for {set i 1} {$i <= 39} {incr i} {
465 set eta [expr { -4.2108 + $i * 0.0696}]
466 add EtaPhiBins $eta $PhiBins
467 }
468
469 for {set i 1} {$i <= 39} {incr i} {
470 set eta [expr { 1.4964 + $i * 0.0696}]
471 add EtaPhiBins $eta $PhiBins
472 }
473
474 # take present CMS granularity for HF
475
476 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
477 set PhiBins {}
478 for {set i -18} {$i <= 18} {incr i} {
479 add PhiBins [expr {$i * $pi/18.0}]
480 }
481
482 foreach eta {-5 -4.7 -4.525 -4.35 -4.2108 4.35 4.525 4.7 5} {
483 add EtaPhiBins $eta $PhiBins
484 }
485
486
487 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
488 add EnergyFraction {0} {1.0}
489 # energy fractions for e, gamma and pi0
490 add EnergyFraction {11} {0.0}
491 add EnergyFraction {22} {0.0}
492 add EnergyFraction {111} {0.0}
493 # energy fractions for muon, neutrinos and neutralinos
494 add EnergyFraction {12} {0.0}
495 add EnergyFraction {13} {0.0}
496 add EnergyFraction {14} {0.0}
497 add EnergyFraction {16} {0.0}
498 add EnergyFraction {1000022} {0.0}
499 add EnergyFraction {1000023} {0.0}
500 add EnergyFraction {1000025} {0.0}
501 add EnergyFraction {1000035} {0.0}
502 add EnergyFraction {1000045} {0.0}
503 # energy fractions for K0short and Lambda
504 add EnergyFraction {310} {0.7}
505 add EnergyFraction {3122} {0.7}
506
507# set ResolutionFormula {resolution formula as a function of eta and energy}
508 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
509 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
510 (abs(eta) > 3.0 && abs(eta) <= 4.2108) * sqrt(energy^2*0.05^2 + energy*1.96^2) + \
511 (abs(eta) > 4.2108 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.40^2)}
512
513}
514
515#################################
516# Energy resolution for electrons
517#################################
518
519module EnergySmearing PhotonEnergySmearing {
520 set InputArray ECal/eflowPhotons
521 set OutputArray eflowPhotons
522
523 # adding 1% extra photon smearing
524 set ResolutionFormula {energy*0.01}
525
526}
527
528
529
530#################
531# Electron filter
532#################
533
534module PdgCodeFilter ElectronFilter {
535 set InputArray HCal/eflowTracks
536 set OutputArray electrons
537 set Invert true
538 add PdgCode {11}
539 add PdgCode {-11}
540}
541
542
543######################
544# ChargedHadronFilter
545######################
546
547module PdgCodeFilter ChargedHadronFilter {
548 set InputArray HCal/eflowTracks
549 set OutputArray chargedHadrons
550
551 add PdgCode {11}
552 add PdgCode {-11}
553 add PdgCode {13}
554 add PdgCode {-13}
555}
556
557
558###################################################
559# Tower Merger (in case not using e-flow algorithm)
560###################################################
561
562module Merger Calorimeter {
563# add InputArray InputArray
564 add InputArray ECal/ecalTowers
565 add InputArray HCal/hcalTowers
566 set OutputArray towers
567}
568
569
570
571####################
572# Energy flow merger
573####################
574
575module Merger EFlowMerger {
576# add InputArray InputArray
577 add InputArray HCal/eflowTracks
578 add InputArray PhotonEnergySmearing/eflowPhotons
579 add InputArray HCal/eflowNeutralHadrons
580 set OutputArray eflow
581}
582
583######################
584# EFlowFilter
585######################
586
587module PdgCodeFilter EFlowFilter {
588 set InputArray EFlowMerger/eflow
589 set OutputArray eflow
590
591 add PdgCode {11}
592 add PdgCode {-11}
593 add PdgCode {13}
594 add PdgCode {-13}
595}
596
597
598###################
599# Photon efficiency
600###################
601
602module Efficiency PhotonEfficiency {
603 set InputArray PhotonEnergySmearing/eflowPhotons
604 set OutputArray photons
605
606 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
607
608 # efficiency formula for photons
609 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
610 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
611 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
612 (abs(eta) > 4.0) * (0.00)}
613}
614
615
616##################
617# Photon isolation
618##################
619
620module Isolation PhotonIsolation {
621 set CandidateInputArray PhotonEfficiency/photons
622 set IsolationInputArray EFlowFilter/eflow
623
624 set OutputArray photons
625
626 set DeltaRMax 0.5
627
628 set PTMin 0.5
629
630 set PTRatioMax 0.12
631}
632
633
634#####################
635# Electron efficiency
636#####################
637
638module Efficiency ElectronEfficiency {
639 set InputArray ElectronFilter/electrons
640 set OutputArray electrons
641
642 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
643 set EfficiencyFormula {
644 (pt <= 4.0) * (0.00) +
645 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.018) +
646 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.252) +
647 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.480) +
648 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.681) +
649 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.792) +
650 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.862) +
651 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.859) +
652 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.016) +
653 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.198) +
654 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.446) +
655 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.598) +
656 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.759) +
657 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.847) +
658 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.872) +
659 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 4.0 && pt <= 6.0) * (0.005) +
660 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 6.0 && pt <= 8.0) * (0.029) +
661 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 8.0 && pt <= 10.0) * (0.108) +
662 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 10.0 && pt <= 20.0) * (0.289) +
663 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 20.0 && pt <= 35.0) * (0.570) +
664 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 35.0 && pt <= 50.0) * (0.743) +
665 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 50.0 && pt <= 14000.0) * (0.828) +
666 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 6.0) * (0.026) +
667 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 6.0 && pt <= 8.0) * (0.045) +
668 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 8.0 && pt <= 10.0) * (0.133) +
669 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 20.0) * (0.411) +
670 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 20.0 && pt <= 35.0) * (0.629) +
671 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 35.0 && pt <= 50.0) * (0.761) +
672 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 50.0 && pt <= 14000.0) * (0.752) +
673 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.061) +
674 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.191) +
675 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.337) +
676 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.475) +
677 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.605) +
678 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.713) +
679 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.794) +
680 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.100) +
681 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.223) +
682 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.427) +
683 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.590) +
684 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.720) +
685 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.800) +
686 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.840) +
687 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 4.0 && pt <= 6.0) * (0.049) +
688 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 6.0 && pt <= 8.0) * (0.152) +
689 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 8.0 && pt <= 10.0) * (0.436) +
690 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.679) +
691 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 35.0) * (0.778) +
692 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 35.0 && pt <= 50.0) * (0.830) +
693 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 14000.0) * (0.919) +
694
695 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 4.0 && pt <= 6.0) * (0.049) +
696 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 6.0 && pt <= 8.0) * (0.152) +
697 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 8.0 && pt <= 10.0) * (0.436) +
698 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 20.0) * (0.679) +
699 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 20.0 && pt <= 35.0) * (0.778) +
700 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 35.0 && pt <= 50.0) * (0.830) +
701 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 50.0 && pt <= 14000.0) * (0.919)}
702}
703
704####################
705# Electron isolation
706####################
707
708module Isolation ElectronIsolation {
709 set CandidateInputArray ElectronEfficiency/electrons
710 set IsolationInputArray EFlowFilter/eflow
711
712 set OutputArray electrons
713
714 set DeltaRMax 0.5
715
716 set PTMin 0.5
717
718 set PTRatioMax 0.12
719}
720
721#################
722# Muon efficiency
723#################
724
725module Efficiency MuonEfficiency {
726 set InputArray MuonMomentumSmearing/muons
727 set OutputArray muons
728
729 # set EfficiencyFormula {efficiency as a function of eta and pt}
730
731 # efficiency formula for muons
732 set EfficiencyFormula {
733 (pt <= 2.0) * (0.00) +
734 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
735 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
736 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
737 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.67) +
738 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
739 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.90) +
740 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.92) +
741 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.90) +
742 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
743 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.46) +
744 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
745 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.65) +
746 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
747 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
748 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.93) +
749 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.92) +
750 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.15) +
751 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.47) +
752 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.55) +
753 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.64) +
754 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.78) +
755 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
756 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
757 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
758 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.23) +
759 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
760 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.53) +
761 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.68) +
762 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.78) +
763 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.89) +
764 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.95) +
765 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.88) +
766 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 2.0 && pt <= 4.0) * (0.22) +
767 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 4.0 && pt <= 6.0) * (0.36) +
768 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 6.0 && pt <= 8.0) * (0.44) +
769 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 8.0 && pt <= 10.0) * (0.57) +
770 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.63) +
771 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 35.0) * (0.71) +
772 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 35.0 && pt <= 50.0) * (0.76) +
773 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 50.0 && pt <= 14000.0) * (0.82) +
774 (abs(eta) > 2.8) * (0.00)
775 }}
776
777################
778# Muon isolation
779################
780
781module Isolation MuonIsolation {
782 set CandidateInputArray MuonEfficiency/muons
783 set IsolationInputArray EFlowFilter/eflow
784
785 set OutputArray muons
786
787 set DeltaRMax 0.5
788
789 set PTMin 0.5
790
791 set PTRatioMax 0.25
792}
793
794###################
795# Missing ET merger
796###################
797
798module Merger MissingET {
799# add InputArray InputArray
800 add InputArray EFlowMerger/eflow
801 set MomentumOutputArray momentum
802}
803
804##################
805# Scalar HT merger
806##################
807
808module Merger ScalarHT {
809# add InputArray InputArray
810 add InputArray UniqueObjectFinder/jets
811 add InputArray UniqueObjectFinder/electrons
812 add InputArray UniqueObjectFinder/photons
813 add InputArray UniqueObjectFinder/muons
814 set EnergyOutputArray energy
815}
816
817
818#####################
819# Neutrino Filter
820#####################
821
822module PdgCodeFilter NeutrinoFilter {
823
824 set InputArray Delphes/stableParticles
825 set OutputArray filteredParticles
826
827 set PTMin 0.0
828
829 add PdgCode {12}
830 add PdgCode {14}
831 add PdgCode {16}
832 add PdgCode {-12}
833 add PdgCode {-14}
834 add PdgCode {-16}
835
836}
837
838
839#####################
840# MC truth jet finder
841#####################
842
843module FastJetFinder GenJetFinder {
844 set InputArray NeutrinoFilter/filteredParticles
845
846 set OutputArray jets
847
848 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
849 set JetAlgorithm 6
850 set ParameterR 0.5
851
852 set JetPTMin 15.0
853}
854
855#########################
856# Gen Missing ET merger
857########################
858
859module Merger GenMissingET {
860# add InputArray InputArray
861 add InputArray NeutrinoFilter/filteredParticles
862 set MomentumOutputArray momentum
863}
864
865
866
867############
868# Jet finder
869############
870
871module FastJetFinder FastJetFinder {
872# set InputArray Calorimeter/towers
873 set InputArray EFlowMerger/eflow
874
875 set OutputArray jets
876
877 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
878 set JetAlgorithm 6
879 set ParameterR 0.5
880
881 set JetPTMin 20.0
882}
883
884##################
885# Fat Jet finder
886##################
887
888module FastJetFinder FatJetFinder {
889 set InputArray EFlowMerger/eflow
890
891 set OutputArray jets
892
893 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
894 set JetAlgorithm 6
895 set ParameterR 0.8
896
897 set ComputeNsubjettiness 1
898 set Beta 1.0
899 set AxisMode 4
900
901 set ComputeTrimming 1
902 set RTrim 0.2
903 set PtFracTrim 0.05
904
905 set ComputePruning 1
906 set ZcutPrun 0.1
907 set RcutPrun 0.5
908 set RPrun 0.8
909
910 set ComputeSoftDrop 1
911 set BetaSoftDrop 0.0
912 set SymmetryCutSoftDrop 0.1
913 set R0SoftDrop 0.8
914
915 set JetPTMin 200.0
916}
917
918
919
920
921##################
922# Jet Energy Scale
923##################
924
925module EnergyScale JetEnergyScale {
926 set InputArray FastJetFinder/jets
927 set OutputArray jets
928
929 # scale formula for jets
930 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
931}
932
933########################
934# Jet Flavor Association
935########################
936
937module JetFlavorAssociation JetFlavorAssociation {
938
939 set PartonInputArray Delphes/partons
940 set ParticleInputArray Delphes/allParticles
941 set ParticleLHEFInputArray Delphes/allParticlesLHEF
942 set JetInputArray JetEnergyScale/jets
943
944 set DeltaR 0.5
945 set PartonPTMin 1.0
946 set PartonEtaMax 2.5
947
948}
949
950###########
951# b-tagging
952###########
953
954module BTagging BTagging {
955 set JetInputArray JetEnergyScale/jets
956
957 set BitNumber 0
958
959 add EfficiencyFormula {0} {0.001}
960
961 add EfficiencyFormula {5} {
962 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.527) +
963 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.598) +
964 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.632) +
965 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.647) +
966 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.652) +
967 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.653) +
968 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.653) +
969 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.655) +
970 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.635) +
971 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.626) +
972 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.614) +
973 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.591) +
974 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.576) +
975 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.550) +
976 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.489) +
977 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.444) +
978 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.408) +
979 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.369) +
980 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.334) +
981 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.269) +
982 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.253) +
983 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.247) +
984 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.230) +
985 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.209) +
986 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.209) +
987 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.209) +
988 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.363) +
989 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.465) +
990 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.504) +
991 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.525) +
992 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.531) +
993 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.526) +
994 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.511) +
995 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.506) +
996 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.492) +
997 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.474) +
998 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.452) +
999 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.443) +
1000 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.433) +
1001 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.385) +
1002 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.349) +
1003 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.325) +
1004 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.316) +
1005 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.280) +
1006 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.206) +
1007 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.206) +
1008 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.206) +
1009 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.206) +
1010 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.206) +
1011 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.206) +
1012 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.206) +
1013 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.206) +
1014 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.211) +
1015 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.278) +
1016 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.310) +
1017 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.319) +
1018 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.327) +
1019 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.320) +
1020 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.318) +
1021 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.314) +
1022 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.304) +
1023 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.296) +
1024 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.271) +
1025 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.278) +
1026 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.277) +
1027 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.255) +
1028 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.218) +
1029 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.218) +
1030 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.218) +
1031 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.218) +
1032 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.218) +
1033 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.218) +
1034 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.218) +
1035 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.218) +
1036 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.218) +
1037 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.218) +
1038 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.218) +
1039 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.218)
1040 }
1041
1042 add EfficiencyFormula {4} {
1043 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.023) +
1044 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.025) +
1045 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.029) +
1046 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.031) +
1047 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.034) +
1048 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.035) +
1049 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.036) +
1050 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.038) +
1051 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.039) +
1052 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.039) +
1053 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.042) +
1054 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.042) +
1055 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.045) +
1056 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.040) +
1057 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.038) +
1058 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.034) +
1059 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.031) +
1060 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.026) +
1061 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.025) +
1062 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.019) +
1063 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.017) +
1064 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.018) +
1065 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.016) +
1066 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.015) +
1067 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.015) +
1068 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.015) +
1069 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.017) +
1070 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.018) +
1071 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.023) +
1072 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.026) +
1073 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.029) +
1074 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.028) +
1075 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.026) +
1076 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.032) +
1077 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.032) +
1078 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.033) +
1079 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.034) +
1080 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.036) +
1081 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.035) +
1082 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.035) +
1083 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.032) +
1084 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.025) +
1085 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.027) +
1086 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.020) +
1087 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.017) +
1088 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.017) +
1089 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.017) +
1090 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.017) +
1091 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.017) +
1092 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.017) +
1093 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.017) +
1094 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.017) +
1095 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.018) +
1096 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.017) +
1097 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.018) +
1098 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.018) +
1099 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.020) +
1100 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.020) +
1101 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.019) +
1102 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.020) +
1103 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.019) +
1104 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.018) +
1105 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.021) +
1106 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.016) +
1107 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.022) +
1108 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.023) +
1109 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.014) +
1110 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.014) +
1111 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.014) +
1112 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.014) +
1113 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.014) +
1114 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.014) +
1115 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.014) +
1116 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.014) +
1117 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.014) +
1118 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.014) +
1119 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.014) +
1120 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.014)
1121 }
1122}
1123
1124#############
1125# tau-tagging
1126#############
1127
1128module TauTagging TauTagging {
1129 set ParticleInputArray Delphes/allParticles
1130 set PartonInputArray Delphes/partons
1131 set JetInputArray JetEnergyScale/jets
1132
1133 set DeltaR 0.5
1134
1135 set TauPTMin 20.0
1136
1137 set TauEtaMax 2.3
1138
1139 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1140
1141 add EfficiencyFormula {0} { (abs(eta) < 2.3) * ((( -0.00621816+0.00130097*pt-2.19642e-5*pt^2+1.49393e-7*pt^3-4.58972e-10*pt^4+5.27983e-13*pt^5 )) * (pt<250) + 0.0032*(pt>250)) + \
1142 (abs(eta) > 2.3) * (0.000)
1143 }
1144 add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.97*0.77*( (0.32 + 0.01*pt - 0.000054*pt*pt )*(pt<100)+0.78*(pt>100) ) + \
1145 (abs(eta) > 2.3) * (0.000)
1146 }
1147}
1148
1149#####################################################
1150# Find uniquely identified photons/electrons/tau/jets
1151#####################################################
1152
1153module UniqueObjectFinder UniqueObjectFinder {
1154# earlier arrays take precedence over later ones
1155# add InputArray InputArray OutputArray
1156 add InputArray PhotonIsolation/photons photons
1157 add InputArray ElectronIsolation/electrons electrons
1158 add InputArray MuonIsolation/muons muons
1159 add InputArray JetEnergyScale/jets jets
1160}
1161
1162##################
1163# ROOT tree writer
1164##################
1165
1166# tracks, towers and eflow objects are not stored by default in the output.
1167# if needed (for jet constituent or other studies), uncomment the relevant
1168# "add Branch ..." lines.
1169
1170module TreeWriter TreeWriter {
1171# add Branch InputArray BranchName BranchClass
1172 add Branch Delphes/allParticles Particle GenParticle
1173
1174 add Branch TrackMerger/tracks Track Track
1175 add Branch Calorimeter/towers Tower Tower
1176
1177 add Branch HCal/eflowTracks EFlowTrack Track
1178 add Branch PhotonEnergySmearing/eflowPhotons EFlowPhoton Tower
1179 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1180
1181 add Branch GenJetFinder/jets GenJet Jet
1182 add Branch GenMissingET/momentum GenMissingET MissingET
1183
1184 add Branch UniqueObjectFinder/jets Jet Jet
1185 add Branch UniqueObjectFinder/electrons Electron Electron
1186 add Branch UniqueObjectFinder/photons Photon Photon
1187 add Branch UniqueObjectFinder/muons Muon Muon
1188
1189 add Branch FatJetFinder/jets FatJet Jet
1190
1191 add Branch MissingET/momentum MissingET MissingET
1192 add Branch ScalarHT/energy ScalarHT ScalarHT
1193}
Note: See TracBrowser for help on using the repository browser.