Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_noPU_v0.tcl@ 6182516

ImprovedOutputFile Timing
Last change on this file since 6182516 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: 56.8 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| < 3.0 (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 <= 84} {incr i} {
367 set eta [expr { -2.958 + $i * 0.0174}]
368 add EtaPhiBins $eta $PhiBins
369 }
370
371 for {set i 1} {$i <= 84} {incr i} {
372 set eta [expr { 1.4964 + $i * 0.0174}]
373 add EtaPhiBins $eta $PhiBins
374 }
375
376 # take present CMS granularity for HF
377
378 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
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.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
385 add EtaPhiBins $eta $PhiBins
386 }
387
388
389 add EnergyFraction {0} {0.0}
390 # energy fractions for e, gamma and pi0
391 add EnergyFraction {11} {1.0}
392 add EnergyFraction {22} {1.0}
393 add EnergyFraction {111} {1.0}
394 # energy fractions for muon, neutrinos and neutralinos
395 add EnergyFraction {12} {0.0}
396 add EnergyFraction {13} {0.0}
397 add EnergyFraction {14} {0.0}
398 add EnergyFraction {16} {0.0}
399 add EnergyFraction {1000022} {0.0}
400 add EnergyFraction {1000023} {0.0}
401 add EnergyFraction {1000025} {0.0}
402 add EnergyFraction {1000035} {0.0}
403 add EnergyFraction {1000045} {0.0}
404 # energy fractions for K0short and Lambda
405 add EnergyFraction {310} {0.3}
406 add EnergyFraction {3122} {0.3}
407
408 # set ResolutionFormula {resolution formula as a function of eta and energy}
409
410 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
411 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
412
413 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
414 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
415 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
416 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
417 (abs(eta) >= 3.0 && 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| < 3.0 (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 <= 21} {incr i} {
465 set eta [expr { -2.958 + $i * 0.0696}]
466 add EtaPhiBins $eta $PhiBins
467 }
468
469 for {set i 1} {$i <= 21} {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.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 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) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.40^2)}
511
512}
513
514#################################
515# Energy resolution for electrons
516#################################
517
518module EnergySmearing PhotonEnergySmearing {
519 set InputArray ECal/eflowPhotons
520 set OutputArray eflowPhotons
521
522 # adding 1% extra photon smearing
523 set ResolutionFormula {energy*0.01}
524
525}
526
527
528
529#################
530# Electron filter
531#################
532
533module PdgCodeFilter ElectronFilter {
534 set InputArray HCal/eflowTracks
535 set OutputArray electrons
536 set Invert true
537 add PdgCode {11}
538 add PdgCode {-11}
539}
540
541
542######################
543# ChargedHadronFilter
544######################
545
546module PdgCodeFilter ChargedHadronFilter {
547 set InputArray HCal/eflowTracks
548 set OutputArray chargedHadrons
549
550 add PdgCode {11}
551 add PdgCode {-11}
552 add PdgCode {13}
553 add PdgCode {-13}
554}
555
556
557###################################################
558# Tower Merger (in case not using e-flow algorithm)
559###################################################
560
561module Merger Calorimeter {
562# add InputArray InputArray
563 add InputArray ECal/ecalTowers
564 add InputArray HCal/hcalTowers
565 set OutputArray towers
566}
567
568
569
570####################
571# Energy flow merger
572####################
573
574module Merger EFlowMerger {
575# add InputArray InputArray
576 add InputArray HCal/eflowTracks
577 add InputArray PhotonEnergySmearing/eflowPhotons
578 add InputArray HCal/eflowNeutralHadrons
579 set OutputArray eflow
580}
581
582######################
583# EFlowFilter
584######################
585
586module PdgCodeFilter EFlowFilter {
587 set InputArray EFlowMerger/eflow
588 set OutputArray eflow
589
590 add PdgCode {11}
591 add PdgCode {-11}
592 add PdgCode {13}
593 add PdgCode {-13}
594}
595
596
597###################
598# Photon efficiency
599###################
600
601module Efficiency PhotonEfficiency {
602 set InputArray PhotonEnergySmearing/eflowPhotons
603 set OutputArray photons
604
605 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
606
607 # efficiency formula for photons
608 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
609 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
610 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
611 (abs(eta) > 4.0) * (0.00)}
612}
613
614
615##################
616# Photon isolation
617##################
618
619module Isolation PhotonIsolation {
620 set CandidateInputArray PhotonEfficiency/photons
621 set IsolationInputArray EFlowFilter/eflow
622
623 set OutputArray photons
624
625 set DeltaRMax 0.5
626
627 set PTMin 0.5
628
629 set PTRatioMax 0.12
630}
631
632
633#####################
634# Electron efficiency
635#####################
636
637module Efficiency ElectronEfficiency {
638 set InputArray ElectronFilter/electrons
639 set OutputArray electrons
640
641 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
642 set EfficiencyFormula {
643 (pt <= 4.0) * (0.00) +
644 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.018) +
645 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.252) +
646 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.480) +
647 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.681) +
648 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.792) +
649 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.862) +
650 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.859) +
651 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.016) +
652 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.198) +
653 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.446) +
654 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.598) +
655 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.759) +
656 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.847) +
657 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.872) +
658 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 4.0 && pt <= 6.0) * (0.005) +
659 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 6.0 && pt <= 8.0) * (0.029) +
660 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 8.0 && pt <= 10.0) * (0.108) +
661 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 10.0 && pt <= 20.0) * (0.289) +
662 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 20.0 && pt <= 35.0) * (0.570) +
663 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 35.0 && pt <= 50.0) * (0.743) +
664 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 50.0 && pt <= 14000.0) * (0.828) +
665 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 6.0) * (0.026) +
666 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 6.0 && pt <= 8.0) * (0.045) +
667 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 8.0 && pt <= 10.0) * (0.133) +
668 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 20.0) * (0.411) +
669 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 20.0 && pt <= 35.0) * (0.629) +
670 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 35.0 && pt <= 50.0) * (0.761) +
671 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 50.0 && pt <= 14000.0) * (0.752) +
672 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.061) +
673 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.191) +
674 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.337) +
675 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.475) +
676 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.605) +
677 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.713) +
678 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.794) +
679 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.100) +
680 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.223) +
681 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.427) +
682 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.590) +
683 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.720) +
684 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.800) +
685 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.840) +
686 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 4.0 && pt <= 6.0) * (0.049) +
687 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 6.0 && pt <= 8.0) * (0.152) +
688 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 8.0 && pt <= 10.0) * (0.436) +
689 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.679) +
690 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 35.0) * (0.778) +
691 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 35.0 && pt <= 50.0) * (0.830) +
692 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 14000.0) * (0.919) +
693 (abs(eta) > 3.0) * (0.00)
694 }
695}
696
697####################
698# Electron isolation
699####################
700
701module Isolation ElectronIsolation {
702 set CandidateInputArray ElectronEfficiency/electrons
703 set IsolationInputArray EFlowFilter/eflow
704
705 set OutputArray electrons
706
707 set DeltaRMax 0.5
708
709 set PTMin 0.5
710
711 set PTRatioMax 0.12
712}
713
714#################
715# Muon efficiency
716#################
717
718module Efficiency MuonEfficiency {
719 set InputArray MuonMomentumSmearing/muons
720 set OutputArray muons
721
722 # set EfficiencyFormula {efficiency as a function of eta and pt}
723
724 # efficiency formula for muons
725 set EfficiencyFormula {
726 (pt <= 2.0) * (0.00) +
727 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
728 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
729 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
730 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.67) +
731 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
732 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.90) +
733 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.92) +
734 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.90) +
735 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
736 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.46) +
737 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
738 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.65) +
739 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
740 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
741 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.93) +
742 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.92) +
743 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.15) +
744 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.47) +
745 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.55) +
746 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.64) +
747 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.78) +
748 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
749 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
750 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
751 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.23) +
752 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
753 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.53) +
754 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.68) +
755 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.78) +
756 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.89) +
757 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.95) +
758 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.88) +
759 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 2.0 && pt <= 4.0) * (0.22) +
760 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 4.0 && pt <= 6.0) * (0.36) +
761 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 6.0 && pt <= 8.0) * (0.44) +
762 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 8.0 && pt <= 10.0) * (0.57) +
763 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 10.0 && pt <= 20.0) * (0.63) +
764 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 20.0 && pt <= 35.0) * (0.71) +
765 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 35.0 && pt <= 50.0) * (0.76) +
766 (abs(eta) > 2.0 && abs(eta) <= 2.8) * (pt > 50.0 && pt <= 14000.0) * (0.82) +
767 (abs(eta) > 2.8) * (0.00)
768 }}
769
770################
771# Muon isolation
772################
773
774module Isolation MuonIsolation {
775 set CandidateInputArray MuonEfficiency/muons
776 set IsolationInputArray EFlowFilter/eflow
777
778 set OutputArray muons
779
780 set DeltaRMax 0.5
781
782 set PTMin 0.5
783
784 set PTRatioMax 0.25
785}
786
787###################
788# Missing ET merger
789###################
790
791module Merger MissingET {
792# add InputArray InputArray
793 add InputArray EFlowMerger/eflow
794 set MomentumOutputArray momentum
795}
796
797##################
798# Scalar HT merger
799##################
800
801module Merger ScalarHT {
802# add InputArray InputArray
803 add InputArray UniqueObjectFinder/jets
804 add InputArray UniqueObjectFinder/electrons
805 add InputArray UniqueObjectFinder/photons
806 add InputArray UniqueObjectFinder/muons
807 set EnergyOutputArray energy
808}
809
810
811#####################
812# Neutrino Filter
813#####################
814
815module PdgCodeFilter NeutrinoFilter {
816
817 set InputArray Delphes/stableParticles
818 set OutputArray filteredParticles
819
820 set PTMin 0.0
821
822 add PdgCode {12}
823 add PdgCode {14}
824 add PdgCode {16}
825 add PdgCode {-12}
826 add PdgCode {-14}
827 add PdgCode {-16}
828
829}
830
831
832#####################
833# MC truth jet finder
834#####################
835
836module FastJetFinder GenJetFinder {
837 set InputArray NeutrinoFilter/filteredParticles
838
839 set OutputArray jets
840
841 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
842 set JetAlgorithm 6
843 set ParameterR 0.5
844
845 set JetPTMin 15.0
846}
847
848#########################
849# Gen Missing ET merger
850########################
851
852module Merger GenMissingET {
853# add InputArray InputArray
854 add InputArray NeutrinoFilter/filteredParticles
855 set MomentumOutputArray momentum
856}
857
858
859
860############
861# Jet finder
862############
863
864module FastJetFinder FastJetFinder {
865# set InputArray Calorimeter/towers
866 set InputArray EFlowMerger/eflow
867
868 set OutputArray jets
869
870 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
871 set JetAlgorithm 6
872 set ParameterR 0.5
873
874 set JetPTMin 20.0
875}
876
877##################
878# Fat Jet finder
879##################
880
881module FastJetFinder FatJetFinder {
882 set InputArray EFlowMerger/eflow
883
884 set OutputArray jets
885
886 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
887 set JetAlgorithm 6
888 set ParameterR 0.8
889
890 set ComputeNsubjettiness 1
891 set Beta 1.0
892 set AxisMode 4
893
894 set ComputeTrimming 1
895 set RTrim 0.2
896 set PtFracTrim 0.05
897
898 set ComputePruning 1
899 set ZcutPrun 0.1
900 set RcutPrun 0.5
901 set RPrun 0.8
902
903 set ComputeSoftDrop 1
904 set BetaSoftDrop 0.0
905 set SymmetryCutSoftDrop 0.1
906 set R0SoftDrop 0.8
907
908 set JetPTMin 200.0
909}
910
911
912
913
914##################
915# Jet Energy Scale
916##################
917
918module EnergyScale JetEnergyScale {
919 set InputArray FastJetFinder/jets
920 set OutputArray jets
921
922 # scale formula for jets
923 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
924}
925
926########################
927# Jet Flavor Association
928########################
929
930module JetFlavorAssociation JetFlavorAssociation {
931
932 set PartonInputArray Delphes/partons
933 set ParticleInputArray Delphes/allParticles
934 set ParticleLHEFInputArray Delphes/allParticlesLHEF
935 set JetInputArray JetEnergyScale/jets
936
937 set DeltaR 0.5
938 set PartonPTMin 1.0
939 set PartonEtaMax 2.5
940
941}
942
943###########
944# b-tagging
945###########
946
947module BTagging BTagging {
948 set JetInputArray JetEnergyScale/jets
949
950 set BitNumber 0
951
952 add EfficiencyFormula {0} {0.001}
953
954 add EfficiencyFormula {5} {
955 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.527) +
956 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.598) +
957 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.632) +
958 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.647) +
959 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.652) +
960 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.653) +
961 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.653) +
962 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.655) +
963 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.635) +
964 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.626) +
965 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.614) +
966 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.591) +
967 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.576) +
968 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.550) +
969 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.489) +
970 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.444) +
971 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.408) +
972 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.369) +
973 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.334) +
974 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.269) +
975 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.253) +
976 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.247) +
977 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.230) +
978 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.209) +
979 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.209) +
980 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.209) +
981 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.363) +
982 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.465) +
983 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.504) +
984 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.525) +
985 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.531) +
986 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.526) +
987 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.511) +
988 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.506) +
989 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.492) +
990 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.474) +
991 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.452) +
992 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.443) +
993 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.433) +
994 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.385) +
995 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.349) +
996 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.325) +
997 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.316) +
998 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.280) +
999 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.206) +
1000 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.206) +
1001 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.206) +
1002 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.206) +
1003 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.206) +
1004 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.206) +
1005 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.206) +
1006 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.206) +
1007 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.211) +
1008 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.278) +
1009 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.310) +
1010 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.319) +
1011 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.327) +
1012 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.320) +
1013 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.318) +
1014 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.314) +
1015 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.304) +
1016 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.296) +
1017 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.271) +
1018 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.278) +
1019 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.277) +
1020 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.255) +
1021 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.218) +
1022 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.218) +
1023 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.218) +
1024 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.218) +
1025 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.218) +
1026 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.218) +
1027 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.218) +
1028 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.218) +
1029 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.218) +
1030 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.218) +
1031 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.218) +
1032 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.218)
1033 }
1034
1035 add EfficiencyFormula {4} {
1036 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.023) +
1037 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.025) +
1038 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.029) +
1039 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.031) +
1040 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.034) +
1041 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.035) +
1042 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.036) +
1043 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.038) +
1044 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.039) +
1045 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.039) +
1046 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.042) +
1047 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.042) +
1048 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.045) +
1049 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.040) +
1050 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.038) +
1051 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.034) +
1052 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.031) +
1053 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.026) +
1054 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.025) +
1055 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.019) +
1056 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.017) +
1057 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.018) +
1058 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.016) +
1059 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.015) +
1060 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.015) +
1061 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.015) +
1062 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.017) +
1063 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.018) +
1064 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.023) +
1065 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.026) +
1066 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.029) +
1067 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.028) +
1068 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.026) +
1069 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.032) +
1070 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.032) +
1071 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.033) +
1072 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.034) +
1073 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.036) +
1074 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.035) +
1075 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.035) +
1076 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.032) +
1077 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.025) +
1078 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.027) +
1079 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.020) +
1080 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.017) +
1081 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.017) +
1082 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.017) +
1083 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.017) +
1084 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.017) +
1085 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.017) +
1086 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.017) +
1087 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.017) +
1088 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.018) +
1089 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.017) +
1090 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.018) +
1091 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.018) +
1092 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.020) +
1093 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.020) +
1094 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.019) +
1095 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.020) +
1096 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.019) +
1097 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.018) +
1098 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.021) +
1099 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.016) +
1100 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.022) +
1101 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.023) +
1102 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.014) +
1103 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.014) +
1104 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.014) +
1105 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.014) +
1106 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.014) +
1107 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.014) +
1108 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.014) +
1109 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.014) +
1110 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.014) +
1111 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.014) +
1112 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.014) +
1113 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.014)
1114 }
1115
1116}
1117
1118#############
1119# tau-tagging
1120#############
1121
1122module TauTagging TauTagging {
1123 set ParticleInputArray Delphes/allParticles
1124 set PartonInputArray Delphes/partons
1125 set JetInputArray JetEnergyScale/jets
1126
1127 set DeltaR 0.5
1128
1129 set TauPTMin 20.0
1130
1131 set TauEtaMax 2.3
1132
1133 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1134
1135 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)) + \
1136 (abs(eta) > 2.3) * (0.000)
1137 }
1138 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) ) + \
1139 (abs(eta) > 2.3) * (0.000)
1140 }
1141}
1142
1143#####################################################
1144# Find uniquely identified photons/electrons/tau/jets
1145#####################################################
1146
1147module UniqueObjectFinder UniqueObjectFinder {
1148# earlier arrays take precedence over later ones
1149# add InputArray InputArray OutputArray
1150 add InputArray PhotonIsolation/photons photons
1151 add InputArray ElectronIsolation/electrons electrons
1152 add InputArray MuonIsolation/muons muons
1153 add InputArray JetEnergyScale/jets jets
1154}
1155
1156##################
1157# ROOT tree writer
1158##################
1159
1160# tracks, towers and eflow objects are not stored by default in the output.
1161# if needed (for jet constituent or other studies), uncomment the relevant
1162# "add Branch ..." lines.
1163
1164module TreeWriter TreeWriter {
1165# add Branch InputArray BranchName BranchClass
1166 add Branch Delphes/allParticles Particle GenParticle
1167
1168 add Branch TrackMerger/tracks Track Track
1169 add Branch Calorimeter/towers Tower Tower
1170
1171 add Branch HCal/eflowTracks EFlowTrack Track
1172 add Branch PhotonEnergySmearing/eflowPhotons EFlowPhoton Tower
1173 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1174
1175 add Branch GenJetFinder/jets GenJet Jet
1176 add Branch GenMissingET/momentum GenMissingET MissingET
1177
1178 add Branch UniqueObjectFinder/jets Jet Jet
1179 add Branch UniqueObjectFinder/electrons Electron Electron
1180 add Branch UniqueObjectFinder/photons Photon Photon
1181 add Branch UniqueObjectFinder/muons Muon Muon
1182
1183 add Branch FatJetFinder/jets FatJet Jet
1184
1185 add Branch MissingET/momentum MissingET MissingET
1186 add Branch ScalarHT/energy ScalarHT ScalarHT
1187}
Note: See TracBrowser for help on using the repository browser.