Fork me on GitHub

source: git/cards/delphes_card_FCC.tcl@ 21eab4f

Last change on this file since 21eab4f was 59d6164, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

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