Fork me on GitHub

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

Last change on this file since b03c0d0 was 94f8f5f, checked in by Michele Selvaggi <michele.selvaggi@…>, 5 years ago

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