Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl@ e55f5b0

ImprovedOutputFile Timing dual_readout llp
Last change on this file since e55f5b0 was fe3908f, checked in by Alexandre Mertens <alexandre.mertens@…>, 8 years ago

working with calos
  • Property mode set to 100644
File size: 18.6 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 FastCaloJetFinder
38
39 JetEnergyScale
40
41 JetFlavorAssociation
42
43 BTagging
44 TauTagging
45
46 UniqueObjectFinder
47
48 ScalarHT
49
50 TreeWriter
51}
52
53#################################
54# Propagate particles in cylinder
55#################################
56
57module ParticlePropagator ParticlePropagator {
58 set InputArray Delphes/stableParticles
59
60 set OutputArray stableParticles
61 set ChargedHadronOutputArray chargedHadrons
62 set ElectronOutputArray electrons
63 set MuonOutputArray muons
64
65 # radius of the magnetic field coverage, in m
66 set Radius 1.29
67 # half-length of the magnetic field coverage, in m
68 set HalfLength 3.00
69
70 # magnetic field
71 set Bz 3.8
72}
73
74####################################
75# Charged hadron tracking efficiency
76####################################
77
78module Efficiency ChargedHadronTrackingEfficiency {
79 set InputArray ParticlePropagator/chargedHadrons
80 set OutputArray chargedHadrons
81
82 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
83
84 # tracking efficiency formula for charged hadrons
85 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
86 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
87 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
88 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
89 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
90 (abs(eta) > 2.5) * (0.00)}
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.1) * (0.00) +
105 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
106 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
107 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
108 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
109 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
110 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
111 (abs(eta) > 2.5) * (0.00)}
112}
113
114##########################
115# Muon tracking efficiency
116##########################
117
118module Efficiency MuonTrackingEfficiency {
119 set InputArray ParticlePropagator/muons
120 set OutputArray muons
121
122 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
123
124 # tracking efficiency formula for muons
125 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
126 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
127 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
128 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
129 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
130 (abs(eta) > 2.5) * (0.00)}
131}
132
133########################################
134# Momentum resolution for charged tracks
135########################################
136
137module MomentumSmearing ChargedHadronMomentumSmearing {
138 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
139 set OutputArray chargedHadrons
140
141 # set ResolutionFormula {resolution formula as a function of eta and pt}
142
143 # resolution formula for charged hadrons
144 # based on arXiv:1405.6569
145 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
146 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
148}
149
150###################################
151# Momentum resolution for electrons
152###################################
153
154module MomentumSmearing ElectronMomentumSmearing {
155 set InputArray ElectronTrackingEfficiency/electrons
156 set OutputArray electrons
157
158 # set ResolutionFormula {resolution formula as a function of eta and energy}
159
160 # resolution formula for electrons
161 # based on arXiv:1405.6569
162 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
163 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
164 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
165}
166
167###############################
168# Momentum resolution for muons
169###############################
170
171module MomentumSmearing MuonMomentumSmearing {
172 set InputArray MuonTrackingEfficiency/muons
173 set OutputArray muons
174
175 # set ResolutionFormula {resolution formula as a function of eta and pt}
176
177 # resolution formula for muons
178 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
179 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
180 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
181}
182
183##############
184# Track merger
185##############
186
187module Merger TrackMerger {
188# add InputArray InputArray
189 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
190 add InputArray ElectronMomentumSmearing/electrons
191 add InputArray MuonMomentumSmearing/muons
192 set OutputArray tracks
193}
194
195#############
196# Calorimeter
197#############
198
199module Calorimeter Calorimeter {
200 set ParticleInputArray ParticlePropagator/stableParticles
201 set TrackInputArray TrackMerger/tracks
202
203 set TowerOutputArray towers
204 set PhotonOutputArray photons
205
206 set EFlowTrackOutputArray eflowTracks
207 set EFlowPhotonOutputArray eflowPhotons
208 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
209
210 set ECalEnergyMin 0.5
211 set HCalEnergyMin 1.0
212
213 set ECalEnergySignificanceMin 2.0
214 set HCalEnergySignificanceMin 2.0
215
216 set SmearTowerCenter true
217
218 set pi [expr {acos(-1)}]
219
220 # lists of the edges of each tower in eta and phi
221 # each list starts with the lower edge of the first tower
222 # the list ends with the higher edged of the last tower
223
224 # 5 degrees towers
225 set PhiBins {}
226 for {set i -36} {$i <= 36} {incr i} {
227 add PhiBins [expr {$i * $pi/36.0}]
228 }
229 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
230 add EtaPhiBins $eta $PhiBins
231 }
232
233 # 10 degrees towers
234 set PhiBins {}
235 for {set i -18} {$i <= 18} {incr i} {
236 add PhiBins [expr {$i * $pi/18.0}]
237 }
238 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
239 add EtaPhiBins $eta $PhiBins
240 }
241
242 # 20 degrees towers
243 set PhiBins {}
244 for {set i -9} {$i <= 9} {incr i} {
245 add PhiBins [expr {$i * $pi/9.0}]
246 }
247 foreach eta {-5 -4.7 -4.525 4.7 5} {
248 add EtaPhiBins $eta $PhiBins
249 }
250
251 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
252 add EnergyFraction {0} {0.0 1.0}
253 # energy fractions for e, gamma and pi0
254 add EnergyFraction {11} {1.0 0.0}
255 add EnergyFraction {22} {1.0 0.0}
256 add EnergyFraction {111} {1.0 0.0}
257 # energy fractions for muon, neutrinos and neutralinos
258 add EnergyFraction {12} {0.0 0.0}
259 add EnergyFraction {13} {0.0 0.0}
260 add EnergyFraction {14} {0.0 0.0}
261 add EnergyFraction {16} {0.0 0.0}
262 add EnergyFraction {1000022} {0.0 0.0}
263 add EnergyFraction {1000023} {0.0 0.0}
264 add EnergyFraction {1000025} {0.0 0.0}
265 add EnergyFraction {1000035} {0.0 0.0}
266 add EnergyFraction {1000045} {0.0 0.0}
267 # energy fractions for K0short and Lambda
268 add EnergyFraction {310} {0.3 0.7}
269 add EnergyFraction {3122} {0.3 0.7}
270
271 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
272 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
273 set ECalResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
274 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
275 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
276
277 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
278 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
279 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
280}
281
282####################
283# Energy flow merger
284####################
285
286module Merger EFlowMerger {
287# add InputArray InputArray
288 add InputArray Calorimeter/eflowTracks
289 add InputArray Calorimeter/eflowPhotons
290 add InputArray Calorimeter/eflowNeutralHadrons
291 set OutputArray eflow
292}
293
294###################
295# Photon efficiency
296###################
297
298module Efficiency PhotonEfficiency {
299 set InputArray Calorimeter/eflowPhotons
300 set OutputArray photons
301
302 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
303
304 # efficiency formula for photons
305 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
306 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
307 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
308 (abs(eta) > 2.5) * (0.00)}
309}
310
311##################
312# Photon isolation
313##################
314
315module Isolation PhotonIsolation {
316 set CandidateInputArray PhotonEfficiency/photons
317 set IsolationInputArray EFlowMerger/eflow
318
319 set OutputArray photons
320
321 set DeltaRMax 0.5
322
323 set PTMin 0.5
324
325 set PTRatioMax 0.12
326}
327
328#################
329# Electron filter
330#################
331
332module PdgCodeFilter ElectronFilter {
333 set InputArray Calorimeter/eflowTracks
334 set OutputArray electrons
335 set Invert true
336 add PdgCode {11}
337 add PdgCode {-11}
338}
339
340#####################
341# Electron efficiency
342#####################
343
344module Efficiency ElectronEfficiency {
345 set InputArray ElectronFilter/electrons
346 set OutputArray electrons
347
348 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
349
350 # efficiency formula for electrons
351 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
352 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
353 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
354 (abs(eta) > 2.5) * (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.5
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 <= 10.0) * (0.00) +
386 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
387 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
388 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
389 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
390 (abs(eta) > 2.4) * (0.00)}
391}
392
393################
394# Muon isolation
395################
396
397module Isolation MuonIsolation {
398 set CandidateInputArray MuonEfficiency/muons
399 set IsolationInputArray EFlowMerger/eflow
400
401 set OutputArray muons
402
403 set DeltaRMax 0.5
404
405 set PTMin 0.5
406
407 set PTRatioMax 0.25
408}
409
410###################
411# Missing ET merger
412###################
413
414module Merger MissingET {
415# add InputArray InputArray
416 add InputArray EFlowMerger/eflow
417 set MomentumOutputArray momentum
418}
419
420##################
421# Scalar HT merger
422##################
423
424module Merger ScalarHT {
425# add InputArray InputArray
426 add InputArray UniqueObjectFinder/jets
427 add InputArray UniqueObjectFinder/electrons
428 add InputArray UniqueObjectFinder/photons
429 add InputArray UniqueObjectFinder/muons
430 set EnergyOutputArray energy
431}
432
433
434#####################
435# Neutrino Filter
436#####################
437
438module PdgCodeFilter NeutrinoFilter {
439
440 set InputArray Delphes/stableParticles
441 set OutputArray filteredParticles
442
443 set PTMin 0.0
444
445 add PdgCode {12}
446 add PdgCode {14}
447 add PdgCode {16}
448 add PdgCode {-12}
449 add PdgCode {-14}
450 add PdgCode {-16}
451
452}
453
454
455#####################
456# MC truth jet finder
457#####################
458
459module FastJetFinder GenJetFinder {
460 set InputArray NeutrinoFilter/filteredParticles
461
462 set OutputArray jets
463
464 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
465 set JetAlgorithm 6
466 set ParameterR 0.5
467
468 set JetPTMin 20.0
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 20.0
498}
499
500################
501# caloJet finder
502################
503
504module FastJetFinder FastCaloJetFinder {
505 set InputArray Calorimeter/towers
506
507 set OutputArray jets
508
509 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
510 set JetAlgorithm 6
511 set ParameterR 0.5
512
513 set JetPTMin 20.0
514}
515
516##################
517# Jet Energy Scale
518##################
519
520module EnergyScale JetEnergyScale {
521 set InputArray FastJetFinder/jets
522 set OutputArray jets
523
524 # scale formula for jets
525 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
526}
527
528########################
529# Jet Flavor Association
530########################
531
532module JetFlavorAssociation JetFlavorAssociation {
533
534 set PartonInputArray Delphes/partons
535 set ParticleInputArray Delphes/allParticles
536 set ParticleLHEFInputArray Delphes/allParticlesLHEF
537 set JetInputArray JetEnergyScale/jets
538
539 set DeltaR 0.5
540 set PartonPTMin 1.0
541 set PartonEtaMax 2.5
542
543}
544
545###########
546# b-tagging
547###########
548
549module BTagging BTagging {
550 set JetInputArray JetEnergyScale/jets
551
552 set BitNumber 0
553
554 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
555 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
556 # gluon's PDG code has the lowest priority
557
558 # based on arXiv:1211.4462
559
560 # default efficiency formula (misidentification rate)
561 add EfficiencyFormula {0} {0.01+0.000038*pt}
562
563 # efficiency formula for c-jets (misidentification rate)
564 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
565
566 # efficiency formula for b-jets
567 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
568}
569
570#############
571# tau-tagging
572#############
573
574module TauTagging TauTagging {
575 set ParticleInputArray Delphes/allParticles
576 set PartonInputArray Delphes/partons
577 set JetInputArray JetEnergyScale/jets
578
579 set DeltaR 0.5
580
581 set TauPTMin 1.0
582
583 set TauEtaMax 2.5
584
585 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
586
587 # default efficiency formula (misidentification rate)
588 add EfficiencyFormula {0} {0.01}
589 # efficiency formula for tau-jets
590 add EfficiencyFormula {15} {0.6}
591}
592
593#####################################################
594# Find uniquely identified photons/electrons/tau/jets
595#####################################################
596
597module UniqueObjectFinder UniqueObjectFinder {
598# earlier arrays take precedence over later ones
599# add InputArray InputArray OutputArray
600 add InputArray PhotonIsolation/photons photons
601 add InputArray ElectronIsolation/electrons electrons
602 add InputArray MuonIsolation/muons muons
603 add InputArray JetEnergyScale/jets jets
604}
605
606##################
607# ROOT tree writer
608##################
609
610# tracks, towers and eflow objects are not stored by default in the output.
611# if needed (for jet constituent or other studies), uncomment the relevant
612# "add Branch ..." lines.
613
614module TreeWriter TreeWriter {
615# add Branch InputArray BranchName BranchClass
616 add Branch Delphes/allParticles Particle GenParticle
617
618 add Branch TrackMerger/tracks Track Track
619 add Branch Calorimeter/towers Tower Tower
620
621 add Branch Calorimeter/eflowTracks EFlowTrack Track
622 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
623 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
624
625 add Branch GenJetFinder/jets GenJet Jet
626 add Branch GenMissingET/momentum GenMissingET MissingET
627
628 add Branch UniqueObjectFinder/jets Jet Jet
629 add Branch FastCaloJetFinder/jets CaloJet Jet
630 add Branch UniqueObjectFinder/electrons Electron Electron
631 add Branch UniqueObjectFinder/photons Photon Photon
632 add Branch UniqueObjectFinder/muons Muon Muon
633 add Branch MissingET/momentum MissingET MissingET
634 add Branch ScalarHT/energy ScalarHT ScalarHT
635}
Note: See TracBrowser for help on using the repository browser.