Fork me on GitHub

source: git/cards/delphes_card_ILD.tcl@ 7e227ae

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 7e227ae was 7442f9b, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

put back eflowtrack

  • Property mode set to 100644
File size: 16.6 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 ElectronMomentumSmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19
20 ECal
21 HCal
22
23 TowerMerger
24 EFlowMerger
25
26 PhotonEfficiency
27 PhotonIsolation
28
29 ElectronFilter
30 ElectronEfficiency
31 ElectronIsolation
32
33 MuonEfficiency
34 MuonIsolation
35
36 MissingET
37
38 NeutrinoFilter
39 GenJetFinder
40 FastJetFinder
41
42 JetEnergyScale
43
44 JetFlavorAssociation
45
46 BTagging
47
48 TauTagging
49
50 ScalarHT
51
52 UniqueObjectFinder
53
54 TreeWriter
55}
56
57#################################
58# Propagate particles in cylinder
59#################################
60
61module ParticlePropagator ParticlePropagator {
62 set InputArray Delphes/stableParticles
63
64 set OutputArray stableParticles
65 set ChargedHadronOutputArray chargedHadrons
66 set ElectronOutputArray electrons
67 set MuonOutputArray muons
68
69 # radius of the magnetic field coverage, in m
70 set Radius 1.8
71 # half-length of the magnetic field coverage, in m
72 set HalfLength 2.4
73
74 # magnetic field
75 set Bz 3.5
76}
77
78####################################
79# Charged hadron tracking efficiency
80####################################
81
82module Efficiency ChargedHadronTrackingEfficiency {
83 set InputArray ParticlePropagator/chargedHadrons
84 set OutputArray chargedHadrons
85
86 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
87
88 # tracking efficiency formula for charged hadrons
89 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
90 (abs(eta) <= 2.4) * (pt > 0.1) * (0.99) +
91 (abs(eta) > 2.4) * (0.00)}
92}
93
94##############################
95# Electron tracking efficiency
96##############################
97
98module Efficiency ElectronTrackingEfficiency {
99 set InputArray ParticlePropagator/electrons
100 set OutputArray electrons
101
102 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
103
104 # tracking efficiency formula for electrons
105 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
106 (abs(eta) <= 2.4) * (pt > 0.1) * (0.99) +
107 (abs(eta) > 2.4) * (0.00)}
108}
109
110##########################
111# Muon tracking efficiency
112##########################
113
114module Efficiency MuonTrackingEfficiency {
115 set InputArray ParticlePropagator/muons
116 set OutputArray muons
117
118 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
119
120 # tracking efficiency formula for muons
121 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
122 (abs(eta) <= 2.4) * (pt > 0.1) * (0.99) +
123 (abs(eta) > 2.4) * (0.00)}
124}
125
126########################################
127# Momentum resolution for charged tracks
128########################################
129
130module MomentumSmearing ChargedHadronMomentumSmearing {
131 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
132 set OutputArray chargedHadrons
133
134 # set ResolutionFormula {resolution formula as a function of eta and pt}
135
136 # resolution formula for charged hadrons
137 set ResolutionFormula { (abs(eta) <= 1.0) * sqrt(0.001^2 + pt^2*1.e-5^2) +
138 (abs(eta) > 1.0 && abs(eta) <= 2.4) * sqrt(0.01^2 + pt^2*1.e-4^2)}
139
140
141}
142
143###################################
144# Momentum resolution for electrons
145###################################
146
147module MomentumSmearing ElectronMomentumSmearing {
148 set InputArray ElectronTrackingEfficiency/electrons
149 set OutputArray electrons
150
151 # set ResolutionFormula {resolution formula as a function of eta and energy}
152
153 # resolution formula for charged hadrons
154 set ResolutionFormula { (abs(eta) <= 1.0) * sqrt(0.001^2 + pt^2*1.e-5^2) +
155 (abs(eta) > 1.0 && abs(eta) <= 2.4) * sqrt(0.01^2 + pt^2*1.e-4^2)}
156}
157
158###############################
159# Momentum resolution for muons
160###############################
161
162module MomentumSmearing MuonMomentumSmearing {
163 set InputArray MuonTrackingEfficiency/muons
164 set OutputArray muons
165
166 # set ResolutionFormula {resolution formula as a function of eta and pt}
167
168 # resolution formula for charged hadrons
169 set ResolutionFormula { (abs(eta) <= 1.0) * sqrt(0.001^2 + pt^2*1.e-5^2) +
170 (abs(eta) > 1.0 && abs(eta) <= 2.4) * sqrt(0.01^2 + pt^2*1.e-4^2)}
171
172}
173
174##############
175# Track merger
176##############
177
178module Merger TrackMerger {
179# add InputArray InputArray
180 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
181 add InputArray ElectronMomentumSmearing/electrons
182 add InputArray MuonMomentumSmearing/muons
183 set OutputArray tracks
184}
185
186#############
187# ECAL
188#############
189
190module SimpleCalorimeter ECal {
191 set ParticleInputArray ParticlePropagator/stableParticles
192 set TrackInputArray TrackMerger/tracks
193
194 set TowerOutputArray ecalTowers
195 set EFlowTrackOutputArray eflowTracks
196 set EFlowTowerOutputArray eflowPhotons
197
198 set IsEcal true
199
200 set EnergyMin 0.5
201 set EnergySignificanceMin 1.0
202
203 set SmearTowerCenter true
204
205 set pi [expr {acos(-1)}]
206
207 # lists of the edges of each tower in eta and phi
208 # each list starts with the lower edge of the first tower
209 # the list ends with the higher edged of the last tower
210
211 # 0.5 degree towers (5x5 mm^2)
212 set PhiBins {}
213 for {set i -360} {$i <= 360} {incr i} {
214 add PhiBins [expr {$i * $pi/360.0}]
215 }
216
217 # 0.01 unit in eta up to eta = 2.5
218 for {set i -500} {$i <= 500} {incr i} {
219 set eta [expr {$i * 0.005}]
220 add EtaPhiBins $eta $PhiBins
221 }
222
223 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
224
225 add EnergyFraction {0} {0.0}
226 # energy fractions for e, gamma and pi0
227 add EnergyFraction {11} {1.0}
228 add EnergyFraction {22} {1.0}
229 add EnergyFraction {111} {1.0}
230 # energy fractions for muon, neutrinos and neutralinos
231 add EnergyFraction {12} {0.0}
232 add EnergyFraction {13} {0.0}
233 add EnergyFraction {14} {0.0}
234 add EnergyFraction {16} {0.0}
235 add EnergyFraction {1000022} {0.0}
236 add EnergyFraction {1000023} {0.0}
237 add EnergyFraction {1000025} {0.0}
238 add EnergyFraction {1000035} {0.0}
239 add EnergyFraction {1000045} {0.0}
240 # energy fractions for K0short and Lambda
241 add EnergyFraction {310} {0.3}
242 add EnergyFraction {3122} {0.3}
243
244 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
245
246 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.01^2 + energy*0.15^2) }
247
248}
249
250#############
251# HCAL
252#############
253
254module SimpleCalorimeter HCal {
255 set ParticleInputArray ParticlePropagator/stableParticles
256 set TrackInputArray ECal/eflowTracks
257
258 set TowerOutputArray hcalTowers
259 set EFlowTrackOutputArray eflowTracks
260 set EFlowTowerOutputArray eflowNeutralHadrons
261
262 set IsEcal false
263
264 set EnergyMin 1.0
265 set EnergySignificanceMin 1.0
266
267 set SmearTowerCenter true
268
269 set pi [expr {acos(-1)}]
270
271 # lists of the edges of each tower in eta and phi
272 # each list starts with the lower edge of the first tower
273 # the list ends with the higher edged of the last tower
274
275
276 # 6 degree towers
277 set PhiBins {}
278 for {set i -60} {$i <= 60} {incr i} {
279 add PhiBins [expr {$i * $pi/60.0}]
280 }
281
282 # 0.5 unit in eta up to eta = 3
283 for {set i -60} {$i <= 60} {incr i} {
284 set eta [expr {$i * 0.05}]
285 add EtaPhiBins $eta $PhiBins
286 }
287
288
289 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
290 add EnergyFraction {0} {1.0}
291 # energy fractions for e, gamma and pi0
292 add EnergyFraction {11} {0.0}
293 add EnergyFraction {22} {0.0}
294 add EnergyFraction {111} {0.0}
295 # energy fractions for muon, neutrinos and neutralinos
296 add EnergyFraction {12} {0.0}
297 add EnergyFraction {13} {0.0}
298 add EnergyFraction {14} {0.0}
299 add EnergyFraction {16} {0.0}
300 add EnergyFraction {1000022} {0.0}
301 add EnergyFraction {1000023} {0.0}
302 add EnergyFraction {1000025} {0.0}
303 add EnergyFraction {1000035} {0.0}
304 add EnergyFraction {1000045} {0.0}
305 # energy fractions for K0short and Lambda
306 add EnergyFraction {310} {0.7}
307 add EnergyFraction {3122} {0.7}
308
309 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
310
311 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.015^2 + energy*0.50^2)}
312
313}
314
315#################
316# Electron filter
317#################
318
319module PdgCodeFilter ElectronFilter {
320 set InputArray HCal/eflowTracks
321 set OutputArray electrons
322 set Invert true
323 add PdgCode {11}
324 add PdgCode {-11}
325}
326
327###################################################
328# Tower Merger (in case not using e-flow algorithm)
329###################################################
330
331module Merger TowerMerger {
332# add InputArray InputArray
333 add InputArray ECal/ecalTowers
334 add InputArray HCal/hcalTowers
335 set OutputArray towers
336}
337
338####################
339# Energy flow merger
340####################
341
342module Merger EFlowMerger {
343# add InputArray InputArray
344 add InputArray HCal/eflowTracks
345 add InputArray ECal/eflowPhotons
346 add InputArray HCal/eflowNeutralHadrons
347 set OutputArray eflow
348}
349
350
351###################
352# Missing ET merger
353###################
354
355module Merger MissingET {
356# add InputArray InputArray
357 add InputArray EFlowMerger/eflow
358 set MomentumOutputArray momentum
359}
360
361
362##################
363# Scalar HT merger
364##################
365
366module Merger ScalarHT {
367# add InputArray InputArray
368 add InputArray EFlowMerger/eflow
369 set EnergyOutputArray energy
370}
371
372#################
373# Neutrino Filter
374#################
375
376module PdgCodeFilter NeutrinoFilter {
377
378 set InputArray Delphes/stableParticles
379 set OutputArray filteredParticles
380
381 set PTMin 0.0
382
383 add PdgCode {12}
384 add PdgCode {14}
385 add PdgCode {16}
386 add PdgCode {-12}
387 add PdgCode {-14}
388 add PdgCode {-16}
389
390}
391
392
393#####################
394# MC truth jet finder
395#####################
396
397module FastJetFinder GenJetFinder {
398 set InputArray NeutrinoFilter/filteredParticles
399
400 set OutputArray jets
401
402 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
403 set JetAlgorithm 6
404 set ParameterR 0.5
405
406 set JetPTMin 20.0
407}
408
409############
410# Jet finder
411############
412
413module FastJetFinder FastJetFinder {
414# set InputArray TowerMerger/towers
415 set InputArray EFlowMerger/eflow
416
417 set OutputArray jets
418
419 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
420 set JetAlgorithm 6
421 set ParameterR 0.5
422
423 set JetPTMin 20.0
424}
425
426##################
427# Jet Energy Scale
428##################
429
430module EnergyScale JetEnergyScale {
431 set InputArray FastJetFinder/jets
432 set OutputArray jets
433
434 # scale formula for jets
435 set ScaleFormula {1.00}
436}
437
438
439########################
440# Jet Flavor Association
441########################
442
443module JetFlavorAssociation JetFlavorAssociation {
444
445 set PartonInputArray Delphes/partons
446 set ParticleInputArray Delphes/allParticles
447 set ParticleLHEFInputArray Delphes/allParticlesLHEF
448 set JetInputArray JetEnergyScale/jets
449
450 set DeltaR 0.5
451 set PartonPTMin 1.0
452 set PartonEtaMax 2.5
453
454}
455
456###################
457# Photon efficiency
458###################
459
460module Efficiency PhotonEfficiency {
461 set InputArray ECal/eflowPhotons
462 set OutputArray photons
463
464 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
465
466 # efficiency formula for photons
467 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
468 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
469 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
470 (abs(eta) > 2.5) * (0.00)}
471}
472
473##################
474# Photon isolation
475##################
476
477module Isolation PhotonIsolation {
478 set CandidateInputArray PhotonEfficiency/photons
479 set IsolationInputArray EFlowMerger/eflow
480
481 set OutputArray photons
482
483 set DeltaRMax 0.5
484
485 set PTMin 0.5
486
487 set PTRatioMax 0.12
488}
489
490#####################
491# Electron efficiency
492#####################
493
494module Efficiency ElectronEfficiency {
495 set InputArray ElectronFilter/electrons
496 set OutputArray electrons
497
498 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
499
500 # efficiency formula for electrons
501 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
502 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
503 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
504 (abs(eta) > 2.5) * (0.00)}
505}
506
507####################
508# Electron isolation
509####################
510
511module Isolation ElectronIsolation {
512 set CandidateInputArray ElectronEfficiency/electrons
513 set IsolationInputArray EFlowMerger/eflow
514
515 set OutputArray electrons
516
517 set DeltaRMax 0.5
518
519 set PTMin 0.5
520
521 set PTRatioMax 0.12
522}
523
524#################
525# Muon efficiency
526#################
527
528module Efficiency MuonEfficiency {
529 set InputArray MuonMomentumSmearing/muons
530 set OutputArray muons
531
532 # set EfficiencyFormula {efficiency as a function of eta and pt}
533
534 # efficiency formula for muons
535 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
536 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
537 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
538 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
539 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
540 (abs(eta) > 2.4) * (0.00)}
541}
542
543################
544# Muon isolation
545################
546
547module Isolation MuonIsolation {
548 set CandidateInputArray MuonEfficiency/muons
549 set IsolationInputArray EFlowMerger/eflow
550
551 set OutputArray muons
552
553 set DeltaRMax 0.5
554
555 set PTMin 0.5
556
557 set PTRatioMax 0.25
558}
559
560
561###########
562# b-tagging
563###########
564
565module BTagging BTagging {
566 set JetInputArray JetEnergyScale/jets
567
568 set BitNumber 0
569
570 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
571 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
572 # gluon's PDG code has the lowest priority
573
574 # based on arXiv:1211.4462
575
576 # default efficiency formula (misidentification rate)
577 add EfficiencyFormula {0} {0.01+0.000038*pt}
578
579 # efficiency formula for c-jets (misidentification rate)
580 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
581
582 # efficiency formula for b-jets
583 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
584}
585
586#############
587# tau-tagging
588#############
589
590
591module TauTagging TauTagging {
592 set ParticleInputArray Delphes/allParticles
593 set PartonInputArray Delphes/partons
594 set JetInputArray JetEnergyScale/jets
595
596 set DeltaR 0.5
597
598 set TauPTMin 1.0
599
600 set TauEtaMax 4.0
601
602 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
603
604 # default efficiency formula (misidentification rate)
605 add EfficiencyFormula {0} {0.001}
606 # efficiency formula for tau-jets
607 add EfficiencyFormula {15} {0.4}
608}
609
610#####################################################
611# Find uniquely identified photons/electrons/tau/jets
612#####################################################
613
614module UniqueObjectFinder UniqueObjectFinder {
615# earlier arrays take precedence over later ones
616# add InputArray InputArray OutputArray
617 add InputArray PhotonIsolation/photons photons
618 add InputArray ElectronIsolation/electrons electrons
619 add InputArray MuonIsolation/muons muons
620 add InputArray JetEnergyScale/jets jets
621}
622
623
624##################
625# ROOT tree writer
626##################
627
628module TreeWriter TreeWriter {
629# add Branch InputArray BranchName BranchClass
630 add Branch Delphes/allParticles Particle GenParticle
631 add Branch GenJetFinder/jets GenJet Jet
632
633 add Branch TrackMerger/tracks Track Track
634 add Branch TowerMerger/towers Tower Tower
635
636 add Branch HCal/eflowTracks EFlowTrack Track
637 add Branch ECal/eflowPhotons EFlowPhoton Tower
638 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
639
640 add Branch UniqueObjectFinder/photons Photon Photon
641 add Branch UniqueObjectFinder/electrons Electron Electron
642 add Branch UniqueObjectFinder/muons Muon Muon
643 add Branch UniqueObjectFinder/jets Jet Jet
644
645 add Branch MissingET/momentum MissingET MissingET
646 add Branch ScalarHT/energy ScalarHT ScalarHT
647}
648
Note: See TracBrowser for help on using the repository browser.