Fork me on GitHub

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

ImprovedOutputFile Timing dual_readout llp
Last change on this file since ed333b3 was ed333b3, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

added gen particle skimmer to FCChh card
  • Property mode set to 100644
File size: 25.5 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: Oct 02th, 2016
9#
10# Configuration: FCC-hh generic detector without dipole in the forward region, no pile-up
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 GenParticleFilter
68
69 TreeWriter
70}
71
72#################################
73# Propagate particles in cylinder
74#################################
75
76module ParticlePropagator ParticlePropagator {
77 set InputArray Delphes/stableParticles
78
79 set OutputArray stableParticles
80 set ChargedHadronOutputArray chargedHadrons
81 set ElectronOutputArray electrons
82 set MuonOutputArray muons
83
84 # radius of the magnetic field coverage, in m
85 set Radius 1.5
86 # half-length of the magnetic field coverage, in m
87 set HalfLength 5
88
89 # magnetic field
90 set Bz 4.0
91}
92
93####################################
94# Charged hadron tracking efficiency
95####################################
96
97module Efficiency ChargedHadronTrackingEfficiency {
98 set InputArray ParticlePropagator/chargedHadrons
99 set OutputArray chargedHadrons
100
101 # TBC (which eta_max ? which pT min?)
102
103 # tracking efficiency formula for charged hadrons
104
105 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
106(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
107(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
108(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
109(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
110(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
111(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
112(abs(eta) > 6.0) * (0.00)}
113
114}
115
116##############################
117# Electron tracking efficiency
118##############################
119
120module Efficiency ElectronTrackingEfficiency {
121 set InputArray ParticlePropagator/electrons
122 set OutputArray electrons
123
124# TBC (which eta_max ?)
125# putting same as charged hadrons for now...
126
127 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
128 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
129 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
130 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
131 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
132 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
133 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
134 (abs(eta) > 6.0) * (0.00)}
135
136}
137##########################
138# Muon tracking efficiency
139##########################
140
141module Efficiency MuonTrackingEfficiency {
142 set InputArray ParticlePropagator/muons
143 set OutputArray muons
144
145# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
146# what about high pT ?
147 # tracking efficiency formula for muons
148 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
149 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
150 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
151 (abs(eta) > 6.0) * (0.00)}
152
153}
154
155########################################
156# Momentum resolution for charged tracks
157########################################
158
159module MomentumSmearing ChargedHadronMomentumSmearing {
160 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
161 set OutputArray chargedHadrons
162
163 source momentumResolutionVsP.tcl
164}
165
166
167###################################
168# Momentum resolution for electrons
169###################################
170
171module MomentumSmearing ElectronMomentumSmearing {
172 set InputArray ElectronTrackingEfficiency/electrons
173 set OutputArray electrons
174
175 source momentumResolutionVsP.tcl
176}
177
178
179###############################
180# Momentum resolution for muons
181###############################
182
183module MomentumSmearing MuonMomentumSmearing {
184 set InputArray MuonTrackingEfficiency/muons
185 set OutputArray muons
186
187 # TBC for just putting tracker resolution/ need to add improvement at high pT
188
189 source muonMomentumResolutionVsP.tcl
190}
191
192##############
193# Track merger
194##############
195
196module Merger TrackMerger {
197# add InputArray InputArray
198 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
199 add InputArray ElectronMomentumSmearing/electrons
200 add InputArray MuonMomentumSmearing/muons
201 set OutputArray tracks
202}
203
204#############
205# ECAL
206#############
207
208# TBC : calos seems ok, check eta max value though.
209
210module SimpleCalorimeter ECal {
211 set ParticleInputArray ParticlePropagator/stableParticles
212 set TrackInputArray TrackMerger/tracks
213
214 set TowerOutputArray ecalTowers
215 set EFlowTrackOutputArray eflowTracks
216 set EFlowTowerOutputArray eflowPhotons
217
218 set IsEcal true
219
220 set EnergyMin 0.5
221 set EnergySignificanceMin 2.0
222
223 set SmearTowerCenter true
224
225 set pi [expr {acos(-1)}]
226
227 # lists of the edges of each tower in eta and phi
228 # each list starts with the lower edge of the first tower
229 # the list ends with the higher edged of the last tower
230
231 # 0.5 degree towers (5x5 mm^2)
232 set PhiBins {}
233 for {set i -256} {$i <= 256} {incr i} {
234 add PhiBins [expr {$i * $pi/256.0}]
235 }
236
237 # TBC
238 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} {
239 add EtaPhiBins $eta $PhiBins
240 }
241 # 0.025 eta x 0.025 for eta between 2.5 and 4.0
242 set PhiBins {}
243 for {set i -128} {$i <= 128} {incr i} {
244 add PhiBins [expr {$i * $pi/128.0}]
245 }
246 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} {
247 add EtaPhiBins $eta $PhiBins
248 }
249
250 # 0.05 x 0.05 for eta between 4.0 and 6.0
251 set PhiBins {}
252 for {set i -64} {$i <= 64} {incr i} {
253 add PhiBins [expr {$i * $pi/64.0}]
254 }
255 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} {
256 add EtaPhiBins $eta $PhiBins
257 }
258 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
259
260 add EnergyFraction {0} {0.0}
261 # energy fractions for e, gamma and pi0
262 add EnergyFraction {11} {1.0}
263 add EnergyFraction {22} {1.0}
264 add EnergyFraction {111} {1.0}
265 # energy fractions for muon, neutrinos and neutralinos
266 add EnergyFraction {12} {0.0}
267 add EnergyFraction {13} {0.0}
268 add EnergyFraction {14} {0.0}
269 add EnergyFraction {16} {0.0}
270 add EnergyFraction {1000022} {0.0}
271 add EnergyFraction {1000023} {0.0}
272 add EnergyFraction {1000025} {0.0}
273 add EnergyFraction {1000035} {0.0}
274 add EnergyFraction {1000045} {0.0}
275 # energy fractions for K0short and Lambda
276 # add EnergyFraction {310} {0.3}
277 # add EnergyFraction {3122} {0.3}
278
279 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
280 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.01^2 + energy*0.10^2) + \
281 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.01^2 + energy*0.10^2)}
282
283
284}
285
286#############
287# HCAL
288#############
289
290module SimpleCalorimeter HCal {
291 set ParticleInputArray ParticlePropagator/stableParticles
292 set TrackInputArray ECal/eflowTracks
293
294 set TowerOutputArray hcalTowers
295 set EFlowTrackOutputArray eflowTracks
296 set EFlowTowerOutputArray eflowNeutralHadrons
297
298 set IsEcal false
299
300 set EnergyMin 1.0
301 set EnergySignificanceMin 2.0
302
303 set SmearTowerCenter true
304
305 set pi [expr {acos(-1)}]
306
307 # lists of the edges of each tower in eta and phi
308 # each list starts with the lower edge of the first tower
309 # the list ends with the higher edged of the last tower
310
311 # 6 degree towers
312 set PhiBins {}
313 for {set i -64} {$i <= 64} {incr i} {
314 add PhiBins [expr {$i * $pi/64.0}]
315 }
316
317 # TBC
318 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} {
319 add EtaPhiBins $eta $PhiBins
320 }
321
322 # 0.1 x 0.1 between eta 2.5 and 4.0
323 set PhiBins {}
324 for {set i -32} {$i <= 32} {incr i} {
325 add PhiBins [expr {$i * $pi/32.0}]
326 }
327 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} {
328 add EtaPhiBins $eta $PhiBins
329 }
330
331 # 0.2 x 0.2 between eta 4.0 and 6.0
332 set PhiBins {}
333 for {set i -16} {$i <= 16} {incr i} {
334 add PhiBins [expr {$i * $pi/16.0}]
335 }
336 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} {
337 add EtaPhiBins $eta $PhiBins
338 }
339
340 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
341 add EnergyFraction {0} {1.0}
342 # energy fractions for e, gamma and pi0
343 add EnergyFraction {11} {0.0}
344 add EnergyFraction {22} {0.0}
345 add EnergyFraction {111} {0.0}
346 # energy fractions for muon, neutrinos and neutralinos
347 add EnergyFraction {12} {0.0}
348 add EnergyFraction {13} {0.0}
349 add EnergyFraction {14} {0.0}
350 add EnergyFraction {16} {0.0}
351 add EnergyFraction {1000022} {0.0}
352 add EnergyFraction {1000023} {0.0}
353 add EnergyFraction {1000025} {0.0}
354 add EnergyFraction {1000035} {0.0}
355 add EnergyFraction {1000045} {0.0}
356 # energy fractions for K0short and Lambda
357 # add EnergyFraction {310} {0.7}
358 # add EnergyFraction {3122} {0.7}
359
360 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
361 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
362 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
363}
364
365#################
366# Electron filter
367#################
368
369module PdgCodeFilter ElectronFilter {
370 set InputArray HCal/eflowTracks
371 set OutputArray electrons
372 set Invert true
373 add PdgCode {11}
374 add PdgCode {-11}
375}
376
377
378######################
379# ChargedHadronFilter
380######################
381
382module PdgCodeFilter ChargedHadronFilter {
383 set InputArray HCal/eflowTracks
384 set OutputArray chargedHadrons
385
386 add PdgCode {11}
387 add PdgCode {-11}
388 add PdgCode {13}
389 add PdgCode {-13}
390}
391
392
393###################################################
394# Tower Merger (in case not using e-flow algorithm)
395###################################################
396
397module Merger Calorimeter {
398# add InputArray InputArray
399 add InputArray ECal/ecalTowers
400 add InputArray HCal/hcalTowers
401 set OutputArray towers
402}
403
404####################
405# Energy flow merger
406####################
407
408module Merger EFlowMerger {
409# add InputArray InputArray
410 add InputArray HCal/eflowTracks
411 add InputArray ECal/eflowPhotons
412 add InputArray HCal/eflowNeutralHadrons
413 set OutputArray eflow
414}
415
416######################
417# EFlowFilter
418######################
419
420module PdgCodeFilter EFlowFilter {
421 set InputArray EFlowMerger/eflow
422 set OutputArray eflow
423
424 add PdgCode {11}
425 add PdgCode {-11}
426 add PdgCode {13}
427 add PdgCode {-13}
428}
429
430
431###################
432# Missing ET merger
433###################
434
435module Merger MissingET {
436# add InputArray InputArray
437 add InputArray EFlowMerger/eflow
438 set MomentumOutputArray momentum
439}
440
441
442
443##################
444# Scalar HT merger
445##################
446
447module Merger ScalarHT {
448# add InputArray InputArray
449 add InputArray EFlowMerger/eflow
450 set EnergyOutputArray energy
451}
452
453#################
454# Neutrino Filter
455#################
456
457module PdgCodeFilter NeutrinoFilter {
458
459 set InputArray Delphes/stableParticles
460 set OutputArray filteredParticles
461
462 set PTMin 0.0
463
464 add PdgCode {12}
465 add PdgCode {14}
466 add PdgCode {16}
467 add PdgCode {-12}
468 add PdgCode {-14}
469 add PdgCode {-16}
470
471}
472
473
474#####################
475# MC truth jet finder
476#####################
477
478# TBC: is jet radius fine?
479
480module FastJetFinder GenJetFinder {
481# set InputArray NeutrinoFilter/filteredParticles
482 set InputArray Delphes/stableParticles
483
484 set OutputArray jets
485
486 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
487 set JetAlgorithm 6
488 set ParameterR 0.4
489
490 set JetPTMin 5.0
491}
492
493#########################
494# Gen Missing ET merger
495########################
496
497module Merger GenMissingET {
498
499# add InputArray InputArray
500 add InputArray NeutrinoFilter/filteredParticles
501 set MomentumOutputArray momentum
502}
503
504
505
506############
507# Jet finder
508############
509
510# TBC need to include jet substructure variables
511# TBC is jet radius fine?
512
513module FastJetFinder FastJetFinder {
514# set InputArray Calorimeter/towers
515 set InputArray EFlowMerger/eflow
516
517 set OutputArray jets
518
519 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
520 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness
521
522 set JetAlgorithm 6
523 set ParameterR 0.4
524
525 set JetPTMin 30.0
526}
527
528##################
529# Jet Energy Scale
530##################
531
532module EnergyScale JetEnergyScale {
533 set InputArray FastJetFinder/jets
534 set OutputArray jets
535
536 # scale formula for jets
537 set ScaleFormula {1.00}
538}
539
540
541########################
542# Jet Flavor Association
543########################
544
545module JetFlavorAssociation JetFlavorAssociation {
546
547 set PartonInputArray Delphes/partons
548 set ParticleInputArray Delphes/allParticles
549 set ParticleLHEFInputArray Delphes/allParticlesLHEF
550 set JetInputArray JetEnergyScale/jets
551
552 set DeltaR 0.5
553 set PartonPTMin 5.0
554 set PartonEtaMax 6.0
555
556}
557
558###################
559# Photon efficiency
560###################
561
562module Efficiency PhotonEfficiency {
563 set InputArray ECal/eflowPhotons
564 set OutputArray photons
565
566 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
567
568 set EfficiencyFormula {
569 (pt <= 1.0) * (0.00) + \
570 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
571 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
572 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
573
574 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
575 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
576 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
577
578 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
579 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
580 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
581 (abs(eta) > 6.0) * (0.00)}
582
583}
584
585##################
586# Photon isolation
587##################
588
589# TBC: check values for iso cuts
590
591module Isolation PhotonIsolation {
592 set CandidateInputArray PhotonEfficiency/photons
593 set IsolationInputArray EFlowFilter/eflow
594
595 set OutputArray photons
596
597 set DeltaRMax 0.3
598
599 set PTMin 0.5
600
601 set PTRatioMax 0.1
602}
603
604
605####################
606# Electron isolation
607####################
608
609# TBC: check values for iso cuts
610
611module Isolation ElectronIsolation {
612 set CandidateInputArray ElectronFilter/electrons
613 set IsolationInputArray EFlowFilter/eflow
614
615 set OutputArray electrons
616
617 set DeltaRMax 0.3
618
619 set PTMin 0.5
620
621 set PTRatioMax 0.1
622}
623
624
625################
626# Muon isolation
627################
628
629# TBC: check values for iso cuts
630
631module Isolation MuonIsolation {
632 set CandidateInputArray MuonMomentumSmearing/muons
633 set IsolationInputArray EFlowFilter/eflow
634
635 set OutputArray muons
636
637 set DeltaRMax 0.3
638
639 set PTMin 0.5
640
641 set PTRatioMax 0.2
642}
643
644
645###########
646# b-tagging
647###########
648
649module BTagging BTagging {
650 set JetInputArray JetEnergyScale/jets
651
652 set BitNumber 0
653
654 add EfficiencyFormula {0} {
655
656 (pt <= 10.0) * (0.00) +
657 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.001) + \
658 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
659 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
660 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \
661 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
662 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
663 (abs(eta) > 4.0) * (0.00)}
664
665 add EfficiencyFormula {4} {
666
667 (pt <= 10.0) * (0.00) +
668 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.04) + \
669 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
670 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
671 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
672 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
673 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
674 (abs(eta) > 4.0) * (0.00)}
675
676 add EfficiencyFormula {5} {
677
678 (pt <= 10.0) * (0.00) +
679 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
680 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
681 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
682 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
683 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
684 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
685 (abs(eta) >= 4.0) * (0.00)}
686
687}
688
689###########
690# c-tagging
691###########
692
693module BTagging CTagging {
694 set JetInputArray JetEnergyScale/jets
695
696 set BitNumber 1
697
698 add EfficiencyFormula {0} {
699
700 (pt <= 10.0) * (0.00) +
701 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
702 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
703
704 add EfficiencyFormula {4} {
705
706 (pt <= 10.0) * (0.00) +
707 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \
708 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
709
710 add EfficiencyFormula {5} {
711
712 (pt <= 10.0) * (0.00) +
713 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
714 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
715
716}
717
718
719#############
720# tau-tagging
721#############
722
723
724module TauTagging TauTagging {
725 set ParticleInputArray Delphes/allParticles
726 set PartonInputArray Delphes/partons
727 set JetInputArray JetEnergyScale/jets
728
729 set DeltaR 0.5
730
731 set TauPTMin 1.0
732
733 set TauEtaMax 4.0
734
735 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
736 add EfficiencyFormula {0} {
737
738 (pt <= 10.0) * (0.00) +
739 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
740 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
741 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
742 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
743 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
744 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
745 (abs(eta) > 4.0) * (0.00)}
746
747 add EfficiencyFormula {11} {
748
749 (pt <= 10.0) * (0.00) +
750 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
751 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
752 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
753 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
754 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
755 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
756 (abs(eta) > 4.0) * (0.00)}
757
758 add EfficiencyFormula {15} {
759
760 (pt <= 10.0) * (0.00) +
761 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.6) + \
762 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.6) *(8./9. - pt/30000.) + \
763 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
764 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
765 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
766 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
767 (abs(eta) > 4.0) * (0.00)}
768
769}
770
771#####################################################
772# Find uniquely identified photons/electrons/tau/jets
773#####################################################
774
775module UniqueObjectFinder UniqueObjectFinder {
776# earlier arrays take precedence over later ones
777# add InputArray InputArray OutputArray
778 add InputArray PhotonIsolation/photons photons
779 add InputArray ElectronIsolation/electrons electrons
780 add InputArray MuonIsolation/muons muons
781 add InputArray JetEnergyScale/jets jets
782
783}
784
785######################
786# GenParticleFilter
787######################
788
789# store only interesting particles
790module StatusPidFilter GenParticleFilter {
791 set InputArray Delphes/allParticles
792 set OutputArray filteredParticles
793}
794
795##################
796# ROOT tree writer
797##################
798
799module TreeWriter TreeWriter {
800# add Branch InputArray BranchName BranchClass
801 add Branch Delphes/allParticles Particle GenParticle
802 add Branch GenParticleFilter/filteredParticles FilteredParticle GenParticle
803
804 add Branch GenJetFinder/jets GenJet Jet
805 add Branch GenMissingET/momentum GenMissingET MissingET
806
807 add Branch TrackMerger/tracks Track Track
808 add Branch Calorimeter/towers Tower Tower
809
810 add Branch HCal/eflowTracks EFlowTrack Track
811 add Branch ECal/eflowPhotons EFlowPhoton Tower
812 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
813
814 add Branch UniqueObjectFinder/photons Photon Photon
815 add Branch UniqueObjectFinder/electrons Electron Electron
816 add Branch UniqueObjectFinder/muons Muon Muon
817 add Branch UniqueObjectFinder/jets Jet Jet
818
819 add Branch MissingET/momentum MissingET MissingET
820 add Branch ScalarHT/energy ScalarHT ScalarHT
821}
822
Note: See TracBrowser for help on using the repository browser.