Fork me on GitHub

source: git/cards/FCC/FCChh.tcl@ e5ea42e

Last change on this file since e5ea42e was 0c86fa5, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

align card with FCCSW version
  • Property mode set to 100644
File size: 25.1 KB
Line 
1#
2# Official Delphes card prepared by FCC-hh collaboration
3#
4# Main authors: Heather Gray (CERN)
5# Filip Moortgat (CERN)
6# Michele Selvaggi (CERN)
7#
8# Released on: Nov 14th, 2016
9#
10# Configuration: FCC-hh baseline detector
11#
12#######################################
13# Order of execution of various modules
14#######################################
15
16set ExecutionPath {
17
18 ParticlePropagator
19
20 ChargedHadronTrackingEfficiency
21 ElectronTrackingEfficiency
22 MuonTrackingEfficiency
23
24 ChargedHadronMomentumSmearing
25 ElectronMomentumSmearing
26 MuonMomentumSmearing
27
28 TrackMerger
29
30 ECal
31 HCal
32
33 Calorimeter
34 EFlowMerger
35 EFlowFilter
36
37 PhotonEfficiency
38 PhotonIsolation
39
40 ElectronFilter
41 ElectronIsolation
42
43 ChargedHadronFilter
44
45 MuonIsolation
46
47 NeutrinoFilter
48
49 MissingET
50 GenMissingET
51
52 GenJetFinder
53 FastJetFinder
54
55 JetEnergyScale
56
57 JetFlavorAssociation
58
59 BTagging
60 CTagging
61 TauTagging
62
63 ScalarHT
64
65 UniqueObjectFinder
66
67 TreeWriter
68}
69
70#################################
71# Propagate particles in cylinder
72#################################
73
74module ParticlePropagator ParticlePropagator {
75 set InputArray Delphes/stableParticles
76
77 set OutputArray stableParticles
78 set ChargedHadronOutputArray chargedHadrons
79 set ElectronOutputArray electrons
80 set MuonOutputArray muons
81
82 # radius of the magnetic field coverage, in m
83 set Radius 1.5
84 # half-length of the magnetic field coverage, in m
85 set HalfLength 5
86
87 # magnetic field
88 set Bz 4.0
89}
90
91####################################
92# Charged hadron tracking efficiency
93####################################
94
95module Efficiency ChargedHadronTrackingEfficiency {
96 set InputArray ParticlePropagator/chargedHadrons
97 set OutputArray chargedHadrons
98
99 # TBC (which eta_max ? which pT min?)
100
101 # tracking efficiency formula for charged hadrons
102
103 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
104(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
105(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
106(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
107(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
108(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
109(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
110(abs(eta) > 6.0) * (0.00)}
111
112}
113
114##############################
115# Electron tracking efficiency
116##############################
117
118module Efficiency ElectronTrackingEfficiency {
119 set InputArray ParticlePropagator/electrons
120 set OutputArray electrons
121
122# TBC (which eta_max ?)
123# putting same as charged hadrons for now...
124
125 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
126 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
127 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
128 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
129 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
130 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
131 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
132 (abs(eta) > 6.0) * (0.00)}
133
134}
135##########################
136# Muon tracking efficiency
137##########################
138
139module Efficiency MuonTrackingEfficiency {
140 set InputArray ParticlePropagator/muons
141 set OutputArray muons
142
143# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
144# what about high pT ?
145 # tracking efficiency formula for muons
146 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
147 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
148 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
149 (abs(eta) > 6.0) * (0.00)}
150
151}
152
153########################################
154# Momentum resolution for charged tracks
155########################################
156
157module MomentumSmearing ChargedHadronMomentumSmearing {
158 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
159 set OutputArray chargedHadrons
160
161 source momentumResolutionVsP.tcl
162}
163
164
165###################################
166# Momentum resolution for electrons
167###################################
168
169module MomentumSmearing ElectronMomentumSmearing {
170 set InputArray ElectronTrackingEfficiency/electrons
171 set OutputArray electrons
172
173 source momentumResolutionVsP.tcl
174}
175
176
177###############################
178# Momentum resolution for muons
179###############################
180
181module MomentumSmearing MuonMomentumSmearing {
182 set InputArray MuonTrackingEfficiency/muons
183 set OutputArray muons
184
185 # TBC for just putting tracker resolution/ need to add improvement at high pT
186
187 source muonMomentumResolutionVsP.tcl
188}
189
190##############
191# Track merger
192##############
193
194module Merger TrackMerger {
195# add InputArray InputArray
196 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
197 add InputArray ElectronMomentumSmearing/electrons
198 add InputArray MuonMomentumSmearing/muons
199 set OutputArray tracks
200}
201
202#############
203# ECAL
204#############
205
206# TBC : calos seems ok, check eta max value though.
207
208module SimpleCalorimeter ECal {
209 set ParticleInputArray ParticlePropagator/stableParticles
210 set TrackInputArray TrackMerger/tracks
211
212 set TowerOutputArray ecalTowers
213 set EFlowTrackOutputArray eflowTracks
214 set EFlowTowerOutputArray eflowPhotons
215
216 set IsEcal true
217
218 set EnergyMin 0.5
219 set EnergySignificanceMin 2.0
220
221 set SmearTowerCenter true
222
223 set pi [expr {acos(-1)}]
224
225 # lists of the edges of each tower in eta and phi
226 # each list starts with the lower edge of the first tower
227 # the list ends with the higher edged of the last tower
228
229 # 0.5 degree towers (5x5 mm^2)
230 set PhiBins {}
231 for {set i -256} {$i <= 256} {incr i} {
232 add PhiBins [expr {$i * $pi/256.0}]
233 }
234
235 # TBC
236 foreach eta {-2.4875 -2.475 -2.4625 -2.45 -2.4375 -2.425 -2.4125 -2.4 -2.3875 -2.375 -2.3625 -2.35 -2.3375 -2.325 -2.3125 -2.3 -2.2875 -2.275 -2.2625 -2.25 -2.2375 -2.225 -2.2125 -2.2 -2.1875 -2.175 -2.1625 -2.15 -2.1375 -2.125 -2.1125 -2.1 -2.0875 -2.075 -2.0625 -2.05 -2.0375 -2.025 -2.0125 -2 -1.9875 -1.975 -1.9625 -1.95 -1.9375 -1.925 -1.9125 -1.9 -1.8875 -1.875 -1.8625 -1.85 -1.8375 -1.825 -1.8125 -1.8 -1.7875 -1.775 -1.7625 -1.75 -1.7375 -1.725 -1.7125 -1.7 -1.6875 -1.675 -1.6625 -1.65 -1.6375 -1.625 -1.6125 -1.6 -1.5875 -1.575 -1.5625 -1.55 -1.5375 -1.525 -1.5125 -1.5 -1.4875 -1.475 -1.4625 -1.45 -1.4375 -1.425 -1.4125 -1.4 -1.3875 -1.375 -1.3625 -1.35 -1.3375 -1.325 -1.3125 -1.3 -1.2875 -1.275 -1.2625 -1.25 -1.2375 -1.225 -1.2125 -1.2 -1.1875 -1.175 -1.1625 -1.15 -1.1375 -1.125 -1.1125 -1.1 -1.0875 -1.075 -1.0625 -1.05 -1.0375 -1.025 -1.0125 -1.0 -0.9875 -0.975 -0.9625 -0.95 -0.9375 -0.925 -0.9125 -0.9 -0.8875 -0.875 -0.8625 -0.85 -0.8375 -0.825 -0.8125 -0.8 -0.7875 -0.775 -0.7625 -0.75 -0.7375 -0.725 -0.7125 -0.7 -0.6875 -0.675 -0.6625 -0.65 -0.6375 -0.625 -0.6125 -0.6 -0.5875 -0.575 -0.5625 -0.55 -0.5375 -0.525 -0.5125 -0.5 -0.4875 -0.475 -0.4625 -0.45 -0.4375 -0.425 -0.4125 -0.4 -0.3875 -0.375 -0.3625 -0.35 -0.3375 -0.325 -0.3125 -0.3 -0.2875 -0.275 -0.2625 -0.25 -0.2375 -0.225 -0.2125 -0.2 -0.1875 -0.175 -0.1625 -0.15 -0.1375 -0.125 -0.1125 -0.1 -0.0875 -0.075 -0.0625 -0.05 -0.0375 -0.025 -0.0125 0 0.0125 0.025 0.0375 0.05 0.0625 0.075 0.0875 0.1 0.1125 0.125 0.1375 0.15 0.1625 0.175 0.1875 0.2 0.2125 0.225 0.2375 0.25 0.2625 0.275 0.2875 0.3 0.3125 0.325 0.3375 0.35 0.3625 0.375 0.3875 0.4 0.4125 0.425 0.4375 0.45 0.4625 0.475 0.4875 0.5 0.5125 0.525 0.5375 0.55 0.5625 0.575 0.5875 0.6 0.6125 0.625 0.6375 0.65 0.6625 0.675 0.6875 0.7 0.7125 0.725 0.7375 0.75 0.7625 0.775 0.7875 0.8 0.8125 0.825 0.8375 0.85 0.8625 0.875 0.8875 0.9 0.9125 0.925 0.9375 0.95 0.9625 0.975 0.9875 1.0 1.0125 1.025 1.0375 1.05 1.0625 1.075 1.0875 1.1 1.1125 1.125 1.1375 1.15 1.1625 1.175 1.1875 1.2 1.2125 1.225 1.2375 1.25 1.2625 1.275 1.2875 1.3 1.3125 1.325 1.3375 1.35 1.3625 1.375 1.3875 1.4 1.4125 1.425 1.4375 1.45 1.4625 1.475 1.4875 1.5 1.5125 1.525 1.5375 1.55 1.5625 1.575 1.5875 1.6 1.6125 1.625 1.6375 1.65 1.6625 1.675 1.6875 1.7 1.7125 1.725 1.7375 1.75 1.7625 1.775 1.7875 1.8 1.8125 1.825 1.8375 1.85 1.8625 1.875 1.8875 1.9 1.9125 1.925 1.9375 1.95 1.9625 1.975 1.9875 2 2.0125 2.025 2.0375 2.05 2.0625 2.075 2.0875 2.1 2.1125 2.125 2.1375 2.15 2.1625 2.175 2.1875 2.2 2.2125 2.225 2.2375 2.25 2.2625 2.275 2.2875 2.3 2.3125 2.325 2.3375 2.35 2.3625 2.375 2.3875 2.4 2.4125 2.425 2.4375 2.45 2.4625 2.475 2.4875 2.5} {
237 add EtaPhiBins $eta $PhiBins
238 }
239 # 0.025 eta x 0.025 for eta between 2.5 and 4.0
240 set PhiBins {}
241 for {set i -128} {$i <= 128} {incr i} {
242 add PhiBins [expr {$i * $pi/128.0}]
243 }
244 foreach eta {-3.975 -3.95 -3.925 -3.9 -3.875 -3.85 -3.825 -3.8 -3.775 -3.75 -3.725 -3.7 -3.675 -3.65 -3.625 -3.6 -3.575 -3.55 -3.525 -3.5 -3.475 -3.45 -3.425 -3.4 -3.375 -3.35 -3.325 -3.3 -3.275 -3.25 -3.225 -3.2 -3.175 -3.15 -3.125 -3.1 -3.075 -3.05 -3.025 -3.0 -2.975 -2.95 -2.925 -2.9 -2.875 -2.85 -2.825 -2.8 -2.775 -2.75 -2.725 -2.7 -2.675 -2.65 -2.625 -2.6 -2.575 -2.55 -2.525 -2.5 2.525 2.55 2.575 2.6 2.625 2.65 2.675 2.7 2.725 2.75 2.775 2.8 2.825 2.85 2.875 2.9 2.925 2.95 2.975 3.0 3.025 3.05 3.075 3.1 3.125 3.15 3.175 3.2 3.225 3.25 3.275 3.3 3.325 3.35 3.375 3.4 3.425 3.45 3.475 3.5 3.525 3.55 3.575 3.6 3.625 3.65 3.675 3.7 3.725 3.75 3.775 3.8 3.825 3.85 3.875 3.9 3.925 3.95 3.975 4.0} {
245 add EtaPhiBins $eta $PhiBins
246 }
247
248 # 0.05 x 0.05 for eta between 4.0 and 6.0
249 set PhiBins {}
250 for {set i -64} {$i <= 64} {incr i} {
251 add PhiBins [expr {$i * $pi/64.0}]
252 }
253 foreach eta {-6.0 -5.95 -5.9 -5.85 -5.8 -5.75 -5.7 -5.65 -5.6 -5.55 -5.5 -5.45 -5.4 -5.35 -5.3 -5.25 -5.2 -5.15 -5.1 -5.05 -5.0 -4.95 -4.9 -4.85 -4.8 -4.75 -4.7 -4.65 -4.6 -4.55 -4.5 -4.45 -4.4 -4.35 -4.3 -4.25 -4.2 -4.15 -4.1 -4.05 -4 4.05 4.1 4.15 4.2 4.25 4.3 4.35 4.4 4.45 4.5 4.55 4.6 4.65 4.7 4.75 4.8 4.85 4.9 4.95 5.0 5.05 5.1 5.15 5.2 5.25 5.3 5.35 5.4 5.45 5.5 5.55 5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6.0} {
254 add EtaPhiBins $eta $PhiBins
255 }
256 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
257
258 add EnergyFraction {0} {0.0}
259 # energy fractions for e, gamma and pi0
260 add EnergyFraction {11} {1.0}
261 add EnergyFraction {22} {1.0}
262 add EnergyFraction {111} {1.0}
263 # energy fractions for muon, neutrinos and neutralinos
264 add EnergyFraction {12} {0.0}
265 add EnergyFraction {13} {0.0}
266 add EnergyFraction {14} {0.0}
267 add EnergyFraction {16} {0.0}
268 add EnergyFraction {1000022} {0.0}
269 add EnergyFraction {1000023} {0.0}
270 add EnergyFraction {1000025} {0.0}
271 add EnergyFraction {1000035} {0.0}
272 add EnergyFraction {1000045} {0.0}
273 # energy fractions for K0short and Lambda
274 # add EnergyFraction {310} {0.3}
275 # add EnergyFraction {3122} {0.3}
276
277 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
278 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
279 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
280
281
282}
283
284#############
285# HCAL
286#############
287
288module SimpleCalorimeter HCal {
289 set ParticleInputArray ParticlePropagator/stableParticles
290 set TrackInputArray ECal/eflowTracks
291
292 set TowerOutputArray hcalTowers
293 set EFlowTrackOutputArray eflowTracks
294 set EFlowTowerOutputArray eflowNeutralHadrons
295
296 set IsEcal false
297
298 set EnergyMin 1.0
299 set EnergySignificanceMin 2.0
300
301 set SmearTowerCenter true
302
303 set pi [expr {acos(-1)}]
304
305 # lists of the edges of each tower in eta and phi
306 # each list starts with the lower edge of the first tower
307 # the list ends with the higher edged of the last tower
308
309 # 6 degree towers
310 set PhiBins {}
311 for {set i -64} {$i <= 64} {incr i} {
312 add PhiBins [expr {$i * $pi/64.0}]
313 }
314
315 # TBC
316 foreach eta {-2.45 -2.4 -2.35 -2.3 -2.25 -2.2 -2.15 -2.1 -2.05 -2.0 -1.95 -1.9 -1.85 -1.8 -1.75 -1.7 -1.65 -1.6 -1.55 -1.5 -1.45 -1.4 -1.35 -1.3 -1.25 -1.2 -1.15 -1.1 -1.05 -1.0 -0.95 -0.9 -0.85 -0.8 -0.75 -0.7 -0.65 -0.6 -0.55 -0.5 -0.45 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0.0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1.0 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2.0 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5} {
317 add EtaPhiBins $eta $PhiBins
318 }
319
320 # 0.1 x 0.1 between eta 2.5 and 4.0
321 set PhiBins {}
322 for {set i -32} {$i <= 32} {incr i} {
323 add PhiBins [expr {$i * $pi/32.0}]
324 }
325 foreach eta {-3.9 -3.8 -3.7 -3.6 -3.5 -3.4 -3.3 -3.2 -3.1 -3.0 -2.9 -2.8 -2.7 -2.6 -2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0} {
326 add EtaPhiBins $eta $PhiBins
327 }
328
329 # 0.2 x 0.2 between eta 4.0 and 6.0
330 set PhiBins {}
331 for {set i -16} {$i <= 16} {incr i} {
332 add PhiBins [expr {$i * $pi/16.0}]
333 }
334 foreach eta {-6.0 -5.8 -5.6 -5.4 -5.2 -5.0 -4.8 -4.6 -4.4 -4.2 -4.0 4.2 4.4 4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0} {
335 add EtaPhiBins $eta $PhiBins
336 }
337
338 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
339 add EnergyFraction {0} {1.0}
340 # energy fractions for e, gamma and pi0
341 add EnergyFraction {11} {0.0}
342 add EnergyFraction {22} {0.0}
343 add EnergyFraction {111} {0.0}
344 # energy fractions for muon, neutrinos and neutralinos
345 add EnergyFraction {12} {0.0}
346 add EnergyFraction {13} {0.0}
347 add EnergyFraction {14} {0.0}
348 add EnergyFraction {16} {0.0}
349 add EnergyFraction {1000022} {0.0}
350 add EnergyFraction {1000023} {0.0}
351 add EnergyFraction {1000025} {0.0}
352 add EnergyFraction {1000035} {0.0}
353 add EnergyFraction {1000045} {0.0}
354 # energy fractions for K0short and Lambda
355 # add EnergyFraction {310} {0.7}
356 # add EnergyFraction {3122} {0.7}
357
358 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
359 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
360 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
361}
362
363#################
364# Electron filter
365#################
366
367module PdgCodeFilter ElectronFilter {
368 set InputArray HCal/eflowTracks
369 set OutputArray electrons
370 set Invert true
371 add PdgCode {11}
372 add PdgCode {-11}
373}
374
375
376######################
377# ChargedHadronFilter
378######################
379
380module PdgCodeFilter ChargedHadronFilter {
381 set InputArray HCal/eflowTracks
382 set OutputArray chargedHadrons
383
384 add PdgCode {11}
385 add PdgCode {-11}
386 add PdgCode {13}
387 add PdgCode {-13}
388}
389
390
391###################################################
392# Tower Merger (in case not using e-flow algorithm)
393###################################################
394
395module Merger Calorimeter {
396# add InputArray InputArray
397 add InputArray ECal/ecalTowers
398 add InputArray HCal/hcalTowers
399 set OutputArray towers
400}
401
402####################
403# Energy flow merger
404####################
405
406module Merger EFlowMerger {
407# add InputArray InputArray
408 add InputArray HCal/eflowTracks
409 add InputArray ECal/eflowPhotons
410 add InputArray HCal/eflowNeutralHadrons
411 set OutputArray eflow
412}
413
414######################
415# EFlowFilter
416######################
417
418module PdgCodeFilter EFlowFilter {
419 set InputArray EFlowMerger/eflow
420 set OutputArray eflow
421
422 add PdgCode {11}
423 add PdgCode {-11}
424 add PdgCode {13}
425 add PdgCode {-13}
426}
427
428
429###################
430# Missing ET merger
431###################
432
433module Merger MissingET {
434# add InputArray InputArray
435 add InputArray EFlowMerger/eflow
436 set MomentumOutputArray momentum
437}
438
439
440
441##################
442# Scalar HT merger
443##################
444
445module Merger ScalarHT {
446# add InputArray InputArray
447 add InputArray EFlowMerger/eflow
448 set EnergyOutputArray energy
449}
450
451#################
452# Neutrino Filter
453#################
454
455module PdgCodeFilter NeutrinoFilter {
456
457 set InputArray Delphes/stableParticles
458 set OutputArray filteredParticles
459
460 set PTMin 0.0
461
462 add PdgCode {12}
463 add PdgCode {14}
464 add PdgCode {16}
465 add PdgCode {-12}
466 add PdgCode {-14}
467 add PdgCode {-16}
468
469}
470
471
472#####################
473# MC truth jet finder
474#####################
475
476# TBC: is jet radius fine?
477
478module FastJetFinder GenJetFinder {
479# set InputArray NeutrinoFilter/filteredParticles
480 set InputArray Delphes/stableParticles
481
482 set OutputArray jets
483
484 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
485 set JetAlgorithm 6
486 set ParameterR 0.4
487
488 set JetPTMin 5.0
489}
490
491#########################
492# Gen Missing ET merger
493########################
494
495module Merger GenMissingET {
496
497# add InputArray InputArray
498 add InputArray NeutrinoFilter/filteredParticles
499 set MomentumOutputArray momentum
500}
501
502
503
504############
505# Jet finder
506############
507
508# TBC need to include jet substructure variables
509# TBC is jet radius fine?
510
511module FastJetFinder FastJetFinder {
512# set InputArray Calorimeter/towers
513 set InputArray EFlowMerger/eflow
514
515 set OutputArray jets
516
517 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
518 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness
519
520 set JetAlgorithm 6
521 set ParameterR 0.4
522
523 set JetPTMin 30.0
524}
525
526##################
527# Jet Energy Scale
528##################
529
530module EnergyScale JetEnergyScale {
531 set InputArray FastJetFinder/jets
532 set OutputArray jets
533
534 # scale formula for jets
535 set ScaleFormula {1.00}
536}
537
538
539########################
540# Jet Flavor Association
541########################
542
543module JetFlavorAssociation JetFlavorAssociation {
544
545 set PartonInputArray Delphes/partons
546 set ParticleInputArray Delphes/allParticles
547 set ParticleLHEFInputArray Delphes/allParticlesLHEF
548 set JetInputArray JetEnergyScale/jets
549
550 set DeltaR 0.5
551 set PartonPTMin 5.0
552 set PartonEtaMax 6.0
553
554}
555
556###################
557# Photon efficiency
558###################
559
560module Efficiency PhotonEfficiency {
561 set InputArray ECal/eflowPhotons
562 set OutputArray photons
563
564 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
565
566 set EfficiencyFormula {
567 (pt <= 1.0) * (0.00) + \
568 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
569 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
570 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
571
572 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
573 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
574 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
575
576 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
577 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
578 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
579 (abs(eta) > 6.0) * (0.00)}
580
581}
582
583##################
584# Photon isolation
585##################
586
587# TBC: check values for iso cuts
588
589module Isolation PhotonIsolation {
590 set CandidateInputArray PhotonEfficiency/photons
591 set IsolationInputArray EFlowFilter/eflow
592
593 set OutputArray photons
594
595 set DeltaRMax 0.3
596
597 set PTMin 0.5
598
599 set PTRatioMax 0.1
600}
601
602
603####################
604# Electron isolation
605####################
606
607# TBC: check values for iso cuts
608
609module Isolation ElectronIsolation {
610 set CandidateInputArray ElectronFilter/electrons
611 set IsolationInputArray EFlowFilter/eflow
612
613 set OutputArray electrons
614
615 set DeltaRMax 0.3
616
617 set PTMin 0.5
618
619 set PTRatioMax 0.1
620}
621
622
623################
624# Muon isolation
625################
626
627# TBC: check values for iso cuts
628
629module Isolation MuonIsolation {
630 set CandidateInputArray MuonMomentumSmearing/muons
631 set IsolationInputArray EFlowFilter/eflow
632
633 set OutputArray muons
634
635 set DeltaRMax 0.3
636
637 set PTMin 0.5
638
639 set PTRatioMax 0.2
640}
641
642
643###########
644# b-tagging
645###########
646
647module BTagging BTagging {
648 set JetInputArray JetEnergyScale/jets
649
650 set BitNumber 0
651
652 add EfficiencyFormula {0} {
653
654 (pt <= 10.0) * (0.00) +
655 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.001) + \
656 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
657 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
658 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \
659 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
660 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
661 (abs(eta) > 4.0) * (0.00)}
662
663 add EfficiencyFormula {4} {
664
665 (pt <= 10.0) * (0.00) +
666 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.04) + \
667 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
668 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
669 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
670 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
671 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
672 (abs(eta) > 4.0) * (0.00)}
673
674 add EfficiencyFormula {5} {
675
676 (pt <= 10.0) * (0.00) +
677 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
678 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
679 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
680 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
681 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
682 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
683 (abs(eta) >= 4.0) * (0.00)}
684
685}
686
687###########
688# c-tagging
689###########
690
691module BTagging CTagging {
692 set JetInputArray JetEnergyScale/jets
693
694 set BitNumber 1
695
696 add EfficiencyFormula {0} {
697
698 (pt <= 10.0) * (0.00) +
699 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
700 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
701
702 add EfficiencyFormula {4} {
703
704 (pt <= 10.0) * (0.00) +
705 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \
706 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
707
708 add EfficiencyFormula {5} {
709
710 (pt <= 10.0) * (0.00) +
711 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
712 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
713
714}
715
716
717#############
718# tau-tagging
719#############
720
721
722module TauTagging TauTagging {
723 set ParticleInputArray Delphes/allParticles
724 set PartonInputArray Delphes/partons
725 set JetInputArray JetEnergyScale/jets
726
727 set DeltaR 0.5
728
729 set TauPTMin 1.0
730
731 set TauEtaMax 4.0
732
733 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
734 add EfficiencyFormula {0} {
735
736 (pt <= 10.0) * (0.00) +
737 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
738 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
739 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
740 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
741 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
742 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
743 (abs(eta) > 4.0) * (0.00)}
744
745 add EfficiencyFormula {11} {
746
747 (pt <= 10.0) * (0.00) +
748 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
749 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
750 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
751 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
752 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
753 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
754 (abs(eta) > 4.0) * (0.00)}
755
756 add EfficiencyFormula {15} {
757
758 (pt <= 10.0) * (0.00) +
759 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.6) + \
760 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.6) *(8./9. - pt/30000.) + \
761 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
762 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
763 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
764 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
765 (abs(eta) > 4.0) * (0.00)}
766
767}
768
769#####################################################
770# Find uniquely identified photons/electrons/tau/jets
771#####################################################
772
773module UniqueObjectFinder UniqueObjectFinder {
774# earlier arrays take precedence over later ones
775# add InputArray InputArray OutputArray
776 add InputArray PhotonIsolation/photons photons
777 add InputArray ElectronIsolation/electrons electrons
778 add InputArray MuonIsolation/muons muons
779 add InputArray JetEnergyScale/jets jets
780
781}
782
783
784##################
785# ROOT tree writer
786##################
787
788module TreeWriter TreeWriter {
789# add Branch InputArray BranchName BranchClass
790 add Branch Delphes/allParticles Particle GenParticle
791
792 add Branch GenJetFinder/jets GenJet Jet
793 add Branch GenMissingET/momentum GenMissingET MissingET
794
795 add Branch TrackMerger/tracks Track Track
796 add Branch Calorimeter/towers Tower Tower
797
798 add Branch HCal/eflowTracks EFlowTrack Track
799 add Branch ECal/eflowPhotons EFlowPhoton Tower
800 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
801
802 add Branch UniqueObjectFinder/photons Photon Photon
803 add Branch UniqueObjectFinder/electrons Electron Electron
804 add Branch UniqueObjectFinder/muons Muon Muon
805 add Branch UniqueObjectFinder/jets Jet Jet
806
807 add Branch MissingET/momentum MissingET MissingET
808 add Branch ScalarHT/energy ScalarHT ScalarHT
809}
810
Note: See TracBrowser for help on using the repository browser.