Fork me on GitHub

source: git/cards/delphes_card_ILD.tcl@ 914fb04

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

fixed eflow at high pt

  • Property mode set to 100644
File size: 13.0 KB
Line 
1# based on arXiv:1306.6329
2
3#######################################
4# Order of execution of various modules
5#######################################
6
7set ExecutionPath {
8 ParticlePropagator
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronEnergySmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19 AngularSmearing
20 ImpactParameterSmearing
21
22 ECal
23 HCal
24
25 TowerMerger
26 EFlowMerger
27
28 MissingET
29
30 NeutrinoFilter
31 GenJetFinder
32 FastJetFinder
33
34 JetEnergyScale
35
36 TrackCountingBTagging
37 TauTagging
38
39 ScalarHT
40
41 TreeWriter
42}
43
44#################################
45# Propagate particles in cylinder
46#################################
47
48module ParticlePropagator ParticlePropagator {
49 set InputArray Delphes/stableParticles
50
51 set OutputArray stableParticles
52 set ChargedHadronOutputArray chargedHadrons
53 set ElectronOutputArray electrons
54 set MuonOutputArray muons
55
56 # radius of the magnetic field coverage, in m
57 set Radius 1.8
58 # half-length of the magnetic field coverage, in m
59 set HalfLength 2.4
60
61 # magnetic field
62 set Bz 3.5
63}
64
65####################################
66# Charged hadron tracking efficiency
67####################################
68
69module Efficiency ChargedHadronTrackingEfficiency {
70 set InputArray ParticlePropagator/chargedHadrons
71 set OutputArray chargedHadrons
72
73 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
74
75 # tracking efficiency formula for charged hadrons
76 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
77 (abs(eta) <= 2.4) * (pt > 0.1) * (0.99) +
78 (abs(eta) > 2.4) * (0.00)}
79}
80
81##############################
82# Electron tracking efficiency
83##############################
84
85module Efficiency ElectronTrackingEfficiency {
86 set InputArray ParticlePropagator/electrons
87 set OutputArray electrons
88
89 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
90
91 # tracking efficiency formula for electrons
92 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
93 (abs(eta) <= 2.4) * (pt > 0.1) * (0.99) +
94 (abs(eta) > 2.4) * (0.00)}
95}
96
97##########################
98# Muon tracking efficiency
99##########################
100
101module Efficiency MuonTrackingEfficiency {
102 set InputArray ParticlePropagator/muons
103 set OutputArray muons
104
105 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
106
107 # tracking efficiency formula for muons
108 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
109 (abs(eta) <= 2.4) * (pt > 0.1) * (0.99) +
110 (abs(eta) > 2.4) * (0.00)}
111}
112
113########################################
114# Momentum resolution for charged tracks
115########################################
116
117module MomentumSmearing ChargedHadronMomentumSmearing {
118 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
119 set OutputArray chargedHadrons
120
121 # set ResolutionFormula {resolution formula as a function of eta and pt}
122
123 # resolution formula for charged hadrons
124 set ResolutionFormula { (abs(eta) <= 1.0) * (0.001 + pt*1.e-5) +
125 (abs(eta) > 1.0 && abs(eta) <= 2.4) * (0.01 + pt*1.e-4)}
126
127
128}
129
130#################################
131# Energy resolution for electrons
132#################################
133
134module EnergySmearing ElectronEnergySmearing {
135 set InputArray ElectronTrackingEfficiency/electrons
136 set OutputArray electrons
137
138 # set ResolutionFormula {resolution formula as a function of eta and energy}
139
140 # resolution formula for electrons (we keep nominal tracking reso since ECAL can't beat this, even at high E)
141 set ResolutionFormula { (abs(eta) <= 1.0) * (0.001 + pt*1.e-5) +
142 (abs(eta) > 1.0 && abs(eta) <= 2.4) * (0.01 + pt*1.e-4)}
143}
144
145###############################
146# Momentum resolution for muons
147###############################
148
149module MomentumSmearing MuonMomentumSmearing {
150 set InputArray MuonTrackingEfficiency/muons
151 set OutputArray muons
152
153 # set ResolutionFormula {resolution formula as a function of eta and pt}
154
155 # resolution formula for muons
156 set ResolutionFormula { (abs(eta) <= 1.0) * (0.001 + pt*1.e-5) +
157 (abs(eta) > 1.0 && abs(eta) <= 2.4) * (0.01 + pt*1.e-4)}
158
159
160}
161
162##############
163# Track merger
164##############
165
166module Merger TrackMerger {
167# add InputArray InputArray
168 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
169 add InputArray ElectronEnergySmearing/electrons
170 add InputArray MuonMomentumSmearing/muons
171 set OutputArray tracks
172}
173
174
175################################
176# Track angular smearing
177################################
178
179module AngularSmearing AngularSmearing {
180 set InputArray TrackMerger/tracks
181 set OutputArray tracks
182
183
184 # angular smearing in eta formula as a function of pt and eta
185 set EtaResolutionFormula { 0.001 }
186
187 # angular smearing in phi formula as a function of pt and eta
188 set PhiResolutionFormula { 0.001 }
189
190}
191
192################################
193# Track impact parameter smearing
194################################
195
196module ImpactParameterSmearing ImpactParameterSmearing {
197 set InputArray AngularSmearing/tracks
198 set OutputArray tracks
199
200
201 # absolute impact parameter smearing formula (in mm) as a function of pt and eta
202 set ResolutionFormula {(pt > 0.1 && pt <= 5.0) * (0.010) +
203 (pt > 5.0) * (0.005)}
204
205}
206
207#############
208# ECAL
209#############
210
211module SimpleCalorimeter ECal {
212 set ParticleInputArray ParticlePropagator/stableParticles
213 set TrackInputArray ImpactParameterSmearing/tracks
214
215 set TowerOutputArray ecalTowers
216 set EFlowTrackOutputArray eflowTracks
217 set EFlowTowerOutputArray eflowPhotons
218
219 set IsEcal true
220
221 set EnergyMin 0.5
222 set EnergySignificanceMin 1.0
223
224 set SmearTowerCenter true
225
226 set pi [expr {acos(-1)}]
227
228 # lists of the edges of each tower in eta and phi
229 # each list starts with the lower edge of the first tower
230 # the list ends with the higher edged of the last tower
231
232 # 0.5 degree towers (5x5 mm^2)
233 set PhiBins {}
234 for {set i -360} {$i <= 360} {incr i} {
235 add PhiBins [expr {$i * $pi/360.0}]
236 }
237
238 # 0.01 unit in eta up to eta = 2.5
239 for {set i -500} {$i <= 500} {incr i} {
240 set eta [expr {$i * 0.005}]
241 add EtaPhiBins $eta $PhiBins
242 }
243
244 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
245
246 add EnergyFraction {0} {0.0}
247 # energy fractions for e, gamma and pi0
248 add EnergyFraction {11} {1.0}
249 add EnergyFraction {22} {1.0}
250 add EnergyFraction {111} {1.0}
251 # energy fractions for muon, neutrinos and neutralinos
252 add EnergyFraction {12} {0.0}
253 add EnergyFraction {13} {0.0}
254 add EnergyFraction {14} {0.0}
255 add EnergyFraction {16} {0.0}
256 add EnergyFraction {1000022} {0.0}
257 add EnergyFraction {1000023} {0.0}
258 add EnergyFraction {1000025} {0.0}
259 add EnergyFraction {1000035} {0.0}
260 add EnergyFraction {1000045} {0.0}
261 # energy fractions for K0short and Lambda
262 add EnergyFraction {310} {0.3}
263 add EnergyFraction {3122} {0.3}
264
265 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
266
267 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.15^2) }
268
269}
270
271#############
272# HCAL
273#############
274
275module SimpleCalorimeter HCal {
276 set ParticleInputArray ParticlePropagator/stableParticles
277 set TrackInputArray ECal/eflowTracks
278
279 set TowerOutputArray hcalTowers
280 set EFlowTrackOutputArray eflowTracks
281 set EFlowTowerOutputArray eflowNeutralHadrons
282
283 set IsEcal false
284
285 set EnergyMin 1.0
286 set EnergySignificanceMin 1.0
287
288 set SmearTowerCenter true
289
290 set pi [expr {acos(-1)}]
291
292 # lists of the edges of each tower in eta and phi
293 # each list starts with the lower edge of the first tower
294 # the list ends with the higher edged of the last tower
295
296
297 # 6 degree towers
298 set PhiBins {}
299 for {set i -60} {$i <= 60} {incr i} {
300 add PhiBins [expr {$i * $pi/60.0}]
301 }
302
303 # 0.5 unit in eta up to eta = 3
304 for {set i -60} {$i <= 60} {incr i} {
305 set eta [expr {$i * 0.05}]
306 add EtaPhiBins $eta $PhiBins
307 }
308
309
310 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
311 add EnergyFraction {0} {1.0}
312 # energy fractions for e, gamma and pi0
313 add EnergyFraction {11} {0.0}
314 add EnergyFraction {22} {0.0}
315 add EnergyFraction {111} {0.0}
316 # energy fractions for muon, neutrinos and neutralinos
317 add EnergyFraction {12} {0.0}
318 add EnergyFraction {13} {0.0}
319 add EnergyFraction {14} {0.0}
320 add EnergyFraction {16} {0.0}
321 add EnergyFraction {1000022} {0.0}
322 add EnergyFraction {1000023} {0.0}
323 add EnergyFraction {1000025} {0.0}
324 add EnergyFraction {1000035} {0.0}
325 add EnergyFraction {1000045} {0.0}
326 # energy fractions for K0short and Lambda
327 add EnergyFraction {310} {0.7}
328 add EnergyFraction {3122} {0.7}
329
330 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
331
332 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.015^2 + energy*0.50^2)}
333
334}
335
336####################
337# Tower Merger (in case not using e-flow algorithm)
338####################
339
340module Merger TowerMerger {
341# add InputArray InputArray
342 add InputArray ECal/ecalTowers
343 add InputArray HCal/hcalTowers
344 set OutputArray towers
345}
346
347####################
348# Energy flow merger
349####################
350
351module Merger EFlowMerger {
352# add InputArray InputArray
353 add InputArray HCal/eflowTracks
354 add InputArray ECal/eflowPhotons
355 add InputArray HCal/eflowNeutralHadrons
356 set OutputArray eflow
357}
358
359
360###################
361# Missing ET merger
362###################
363
364module Merger MissingET {
365# add InputArray InputArray
366 add InputArray EFlowMerger/eflow
367 set MomentumOutputArray momentum
368}
369
370
371##################
372# Scalar HT merger
373##################
374
375module Merger ScalarHT {
376# add InputArray InputArray
377 add InputArray EFlowMerger/eflow
378 set EnergyOutputArray energy
379}
380
381#####################
382# Neutrino Filter
383#####################
384
385module PdgCodeFilter NeutrinoFilter {
386
387 set InputArray Delphes/stableParticles
388 set OutputArray filteredParticles
389
390 set PTMin 0.0
391
392 add PdgCode {12}
393 add PdgCode {14}
394 add PdgCode {16}
395 add PdgCode {-12}
396 add PdgCode {-14}
397 add PdgCode {-16}
398
399}
400
401
402#####################
403# MC truth jet finder
404#####################
405
406module FastJetFinder GenJetFinder {
407 set InputArray NeutrinoFilter/filteredParticles
408
409 set OutputArray jets
410
411 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
412 set JetAlgorithm 6
413 set ParameterR 0.5
414
415 set JetPTMin 20.0
416}
417
418############
419# Jet finder
420############
421
422module FastJetFinder FastJetFinder {
423# set InputArray TowerMerger/towers
424 set InputArray EFlowMerger/eflow
425
426 set OutputArray jets
427
428 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
429 set JetAlgorithm 6
430 set ParameterR 0.5
431
432 set JetPTMin 20.0
433}
434
435##################
436# Jet Energy Scale
437##################
438
439module EnergyScale JetEnergyScale {
440 set InputArray FastJetFinder/jets
441 set OutputArray jets
442
443 # scale formula for jets
444 set ScaleFormula {1.00}
445}
446
447##########################
448# Track Counting b-tagging
449##########################
450
451module TrackCountingBTagging TrackCountingBTagging {
452 set TrackInputArray ImpactParameterSmearing/tracks
453 set JetInputArray JetEnergyScale/jets
454
455 set BitNumber 0
456
457 # maximum distance between jet and track
458 set DeltaR 0.3
459
460 # minimum pt of tracks
461 set TrackPTMin 1.0
462
463 # minimum transverse impact parameter (in mm)
464 set TrackIPMax 2.0
465
466 # minimum ip significance for the track to be counted
467 set SigMin 6.5
468
469 # minimum number of tracks (high efficiency n=2, high purity n=3)
470 set Ntracks 3
471}
472
473
474##########################
475# tau-tagging
476##########################
477
478
479module TauTagging TauTagging {
480 set ParticleInputArray Delphes/allParticles
481 set PartonInputArray Delphes/partons
482 set JetInputArray JetEnergyScale/jets
483
484 set DeltaR 0.5
485
486 set TauPTMin 1.0
487
488 set TauEtaMax 4.0
489
490 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
491
492 # default efficiency formula (misidentification rate)
493 add EfficiencyFormula {0} {0.001}
494 # efficiency formula for tau-jets
495 add EfficiencyFormula {15} {0.4}
496}
497
498##################
499# ROOT tree writer
500##################
501
502module TreeWriter TreeWriter {
503# add Branch InputArray BranchName BranchClass
504 add Branch Delphes/allParticles Particle GenParticle
505 add Branch GenJetFinder/jets GenJet Jet
506
507 add Branch ChargedHadronMomentumSmearing/chargedHadrons ChargedHadron Track
508 add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
509 add Branch ECal/eflowPhotons Photon Photon
510
511 add Branch ElectronEnergySmearing/electrons Electron Electron
512 add Branch MuonMomentumSmearing/muons Muon Muon
513 add Branch JetEnergyScale/jets Jet Jet
514 add Branch MissingET/momentum MissingET MissingET
515 add Branch ScalarHT/energy ScalarHT ScalarHT
516}
517
Note: See TracBrowser for help on using the repository browser.