Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_noPU_HFnose_v0.tcl@ 2a2f310

Last change on this file since 2a2f310 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.1 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| < 4.175 (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 <= 154} {incr i} {
368 set eta [expr { -4.176 + $i * 0.0174}]
369 add EtaPhiBins $eta $PhiBins
370 }
371
372 for {set i 1} {$i <= 154} {incr i} {
373 set eta [expr { 1.4964 + $i * 0.0174}]
374 add EtaPhiBins $eta $PhiBins
375 }
376
377 # remaining HF part
378 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 4.2 < |eta| < 5.0
379
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.176 4.35 4.525 4.7 5} {
386 add EtaPhiBins $eta $PhiBins
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) <= 4.176) * sqrt(energy^2*0.01^2 + energy*0.67^2) + \
418 (abs(eta) >= 4.176 && 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| < 4.2 (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 <= 39} {incr i} {
466 set eta [expr { -4.2108 + $i * 0.0696}]
467 add EtaPhiBins $eta $PhiBins
468 }
469
470 for {set i 1} {$i <= 39} {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.2108 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) <= 4.2108) * sqrt(energy^2*0.05^2 + energy*1.96^2) + \
512 (abs(eta) > 4.2108 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.40^2)}
513
514}
515
516#################################
517# Energy resolution for electrons
518#################################
519
520module EnergySmearing PhotonEnergySmearing {
521 set InputArray ECal/eflowPhotons
522 set OutputArray eflowPhotons
523
524 # adding 1% extra photon smearing
525 set ResolutionFormula {energy*0.01}
526
527}
528
529
530
531#################
532# Electron filter
533#################
534
535module PdgCodeFilter ElectronFilter {
536 set InputArray HCal/eflowTracks
537 set OutputArray electrons
538 set Invert true
539 add PdgCode {11}
540 add PdgCode {-11}
541}
542
543
544######################
545# ChargedHadronFilter
546######################
547
548module PdgCodeFilter ChargedHadronFilter {
549 set InputArray HCal/eflowTracks
550 set OutputArray chargedHadrons
551
552 add PdgCode {11}
553 add PdgCode {-11}
554 add PdgCode {13}
555 add PdgCode {-13}
556}
557
558
559###################################################
560# Tower Merger (in case not using e-flow algorithm)
561###################################################
562
563module Merger Calorimeter {
564# add InputArray InputArray
565 add InputArray ECal/ecalTowers
566 add InputArray HCal/hcalTowers
567 set OutputArray towers
568}
569
570
571
572####################
573# Energy flow merger
574####################
575
576module Merger EFlowMerger {
577# add InputArray InputArray
578 add InputArray HCal/eflowTracks
579 add InputArray PhotonEnergySmearing/eflowPhotons
580 add InputArray HCal/eflowNeutralHadrons
581 set OutputArray eflow
582}
583
584######################
585# EFlowFilter
586######################
587
588module PdgCodeFilter EFlowFilter {
589 set InputArray EFlowMerger/eflow
590 set OutputArray eflow
591
592 add PdgCode {11}
593 add PdgCode {-11}
594 add PdgCode {13}
595 add PdgCode {-13}
596}
597
598
599###################
600# Photon efficiency
601###################
602
603module Efficiency PhotonEfficiency {
604 set InputArray PhotonEnergySmearing/eflowPhotons
605 set OutputArray photons
606
607 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
608
609 # efficiency formula for photons
610 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
611 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
612 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
613 (abs(eta) > 4.0) * (0.00)}
614}
615
616
617##################
618# Photon isolation
619##################
620
621module Isolation PhotonIsolation {
622 set CandidateInputArray PhotonEfficiency/photons
623 set IsolationInputArray EFlowFilter/eflow
624
625 set OutputArray photons
626
627 set DeltaRMax 0.5
628
629 set PTMin 0.5
630
631 set PTRatioMax 0.12
632}
633
634
635#####################
636# Electron efficiency
637#####################
638
639module Efficiency ElectronEfficiency {
640 set InputArray ElectronFilter/electrons
641 set OutputArray electrons
642
643 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
644 set EfficiencyFormula {
645 (pt <= 4.0) * (0.00) +
646 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.018) +
647 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.252) +
648 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.480) +
649 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.681) +
650 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.792) +
651 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.862) +
652 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.859) +
653 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.016) +
654 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.198) +
655 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.446) +
656 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.598) +
657 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.759) +
658 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.847) +
659 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.872) +
660 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 4.0 && pt <= 6.0) * (0.005) +
661 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 6.0 && pt <= 8.0) * (0.029) +
662 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 8.0 && pt <= 10.0) * (0.108) +
663 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 10.0 && pt <= 20.0) * (0.289) +
664 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 20.0 && pt <= 35.0) * (0.570) +
665 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 35.0 && pt <= 50.0) * (0.743) +
666 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 50.0 && pt <= 14000.0) * (0.828) +
667 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 6.0) * (0.026) +
668 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 6.0 && pt <= 8.0) * (0.045) +
669 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 8.0 && pt <= 10.0) * (0.133) +
670 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 20.0) * (0.411) +
671 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 20.0 && pt <= 35.0) * (0.629) +
672 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 35.0 && pt <= 50.0) * (0.761) +
673 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 50.0 && pt <= 14000.0) * (0.752) +
674 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.061) +
675 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.191) +
676 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.337) +
677 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.475) +
678 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.605) +
679 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.713) +
680 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.794) +
681 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.100) +
682 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.223) +
683 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.427) +
684 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.590) +
685 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.720) +
686 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.800) +
687 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.840) +
688 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 4.0 && pt <= 6.0) * (0.049) +
689 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 6.0 && pt <= 8.0) * (0.152) +
690 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 8.0 && pt <= 10.0) * (0.436) +
691 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.679) +
692 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 35.0) * (0.778) +
693 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 35.0 && pt <= 50.0) * (0.830) +
694 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 14000.0) * (0.919) +
695
696 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 4.0 && pt <= 6.0) * (0.049) +
697 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 6.0 && pt <= 8.0) * (0.152) +
698 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 8.0 && pt <= 10.0) * (0.436) +
699 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 20.0) * (0.679) +
700 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 20.0 && pt <= 35.0) * (0.778) +
701 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 35.0 && pt <= 50.0) * (0.830) +
702 (abs(eta) > 3.0 && abs(eta) <= 4.0) * (pt > 50.0 && pt <= 14000.0) * (0.919)}
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# tau-tagging
1127#############
1128
1129module TauTagging TauTagging {
1130 set ParticleInputArray Delphes/allParticles
1131 set PartonInputArray Delphes/partons
1132 set JetInputArray JetEnergyScale/jets
1133
1134 set DeltaR 0.5
1135
1136 set TauPTMin 20.0
1137
1138 set TauEtaMax 2.3
1139
1140 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1141
1142 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)) + \
1143 (abs(eta) > 2.3) * (0.000)
1144 }
1145 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) ) + \
1146 (abs(eta) > 2.3) * (0.000)
1147 }
1148}
1149
1150#####################################################
1151# Find uniquely identified photons/electrons/tau/jets
1152#####################################################
1153
1154module UniqueObjectFinder UniqueObjectFinder {
1155# earlier arrays take precedence over later ones
1156# add InputArray InputArray OutputArray
1157 add InputArray PhotonIsolation/photons photons
1158 add InputArray ElectronIsolation/electrons electrons
1159 add InputArray MuonIsolation/muons muons
1160 add InputArray JetEnergyScale/jets jets
1161}
1162
1163###############################################################################################################
1164# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1165###############################################################################################################
1166
1167module StatusPidFilter GenParticleFilter {
1168
1169 set InputArray Delphes/allParticles
1170 set OutputArray filteredParticles
1171 set PTMin 0.0
1172
1173}
1174
1175
1176##################
1177# ROOT tree writer
1178##################
1179
1180# tracks, towers and eflow objects are not stored by default in the output.
1181# if needed (for jet constituent or other studies), uncomment the relevant
1182# "add Branch ..." lines.
1183
1184module TreeWriter TreeWriter {
1185# add Branch InputArray BranchName BranchClass
1186 # add Branch Delphes/allParticles Particle GenParticle
1187 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1188
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.