Fork me on GitHub

source: git/cards/delphes_card_FCC.tcl@ 16217da

Last change on this file since 16217da was 679a87e, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

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