Fork me on GitHub

source: git/cards/delphes_card_ILD.tcl@ 1f34dac

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

fixed GenMissingET

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