Fork me on GitHub

source: git/cards/delphes_card_ILD.tcl@ 0413f44

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

input collection for ElectronFilter

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