Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl@ 5198eed

ImprovedOutputFile Timing dual_readout llp 3.3.3pre16
Last change on this file since 5198eed was 989ee8f, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

fixed typo in ECAL formula.
  • Property mode set to 100644
File size: 18.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 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 1.29
66 # half-length of the magnetic field coverage, in m
67 set HalfLength 3.00
68
69 # magnetic field
70 set Bz 3.8
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.1) * (0.00) +
85 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
86 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
87 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
88 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
89 (abs(eta) > 2.5) * (0.00)}
90}
91
92##############################
93# Electron tracking efficiency
94##############################
95
96module Efficiency ElectronTrackingEfficiency {
97 set InputArray ParticlePropagator/electrons
98 set OutputArray electrons
99
100 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
101
102 # tracking efficiency formula for electrons
103 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
104 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
105 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
106 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
107 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
108 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
109 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
110 (abs(eta) > 2.5) * (0.00)}
111}
112
113##########################
114# Muon tracking efficiency
115##########################
116
117module Efficiency MuonTrackingEfficiency {
118 set InputArray ParticlePropagator/muons
119 set OutputArray muons
120
121 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
122
123 # tracking efficiency formula for muons
124 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
125 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
126 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
127 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
128 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
129 (abs(eta) > 2.5) * (0.00)}
130}
131
132########################################
133# Momentum resolution for charged tracks
134########################################
135
136module MomentumSmearing ChargedHadronMomentumSmearing {
137 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
138 set OutputArray chargedHadrons
139
140 # set ResolutionFormula {resolution formula as a function of eta and pt}
141
142 # resolution formula for charged hadrons
143 # based on arXiv:1405.6569
144 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
145 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
147}
148
149###################################
150# Momentum resolution for electrons
151###################################
152
153module MomentumSmearing ElectronMomentumSmearing {
154 set InputArray ElectronTrackingEfficiency/electrons
155 set OutputArray electrons
156
157 # set ResolutionFormula {resolution formula as a function of eta and energy}
158
159 # resolution formula for electrons
160 # based on arXiv:1405.6569
161 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
162 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
163 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
164}
165
166###############################
167# Momentum resolution for muons
168###############################
169
170module MomentumSmearing MuonMomentumSmearing {
171 set InputArray MuonTrackingEfficiency/muons
172 set OutputArray muons
173
174 # set ResolutionFormula {resolution formula as a function of eta and pt}
175
176 # resolution formula for muons
177 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
178 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
179 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
180}
181
182##############
183# Track merger
184##############
185
186module Merger TrackMerger {
187# add InputArray InputArray
188 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
189 add InputArray ElectronMomentumSmearing/electrons
190 add InputArray MuonMomentumSmearing/muons
191 set OutputArray tracks
192}
193
194#############
195# Calorimeter
196#############
197
198module Calorimeter Calorimeter {
199 set ParticleInputArray ParticlePropagator/stableParticles
200 set TrackInputArray TrackMerger/tracks
201
202 set TowerOutputArray towers
203 set PhotonOutputArray photons
204
205 set EFlowTrackOutputArray eflowTracks
206 set EFlowPhotonOutputArray eflowPhotons
207 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
208
209 set ECalEnergyMin 0.5
210 set HCalEnergyMin 1.0
211
212 set ECalEnergySignificanceMin 1.0
213 set HCalEnergySignificanceMin 1.0
214
215 set SmearTowerCenter true
216
217 set pi [expr {acos(-1)}]
218
219 # lists of the edges of each tower in eta and phi
220 # each list starts with the lower edge of the first tower
221 # the list ends with the higher edged of the last tower
222
223 # 5 degrees towers
224 set PhiBins {}
225 for {set i -36} {$i <= 36} {incr i} {
226 add PhiBins [expr {$i * $pi/36.0}]
227 }
228 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} {
229 add EtaPhiBins $eta $PhiBins
230 }
231
232 # 10 degrees towers
233 set PhiBins {}
234 for {set i -18} {$i <= 18} {incr i} {
235 add PhiBins [expr {$i * $pi/18.0}]
236 }
237 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} {
238 add EtaPhiBins $eta $PhiBins
239 }
240
241 # 20 degrees towers
242 set PhiBins {}
243 for {set i -9} {$i <= 9} {incr i} {
244 add PhiBins [expr {$i * $pi/9.0}]
245 }
246 foreach eta {-5 -4.7 -4.525 4.7 5} {
247 add EtaPhiBins $eta $PhiBins
248 }
249
250 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
251 add EnergyFraction {0} {0.0 1.0}
252 # energy fractions for e, gamma and pi0
253 add EnergyFraction {11} {1.0 0.0}
254 add EnergyFraction {22} {1.0 0.0}
255 add EnergyFraction {111} {1.0 0.0}
256 # energy fractions for muon, neutrinos and neutralinos
257 add EnergyFraction {12} {0.0 0.0}
258 add EnergyFraction {13} {0.0 0.0}
259 add EnergyFraction {14} {0.0 0.0}
260 add EnergyFraction {16} {0.0 0.0}
261 add EnergyFraction {1000022} {0.0 0.0}
262 add EnergyFraction {1000023} {0.0 0.0}
263 add EnergyFraction {1000025} {0.0 0.0}
264 add EnergyFraction {1000035} {0.0 0.0}
265 add EnergyFraction {1000045} {0.0 0.0}
266 # energy fractions for K0short and Lambda
267 add EnergyFraction {310} {0.3 0.7}
268 add EnergyFraction {3122} {0.3 0.7}
269
270 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
271 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
272 set ECalResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
273 (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) +
274 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
275
276 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
277 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
278 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
279}
280
281####################
282# Energy flow merger
283####################
284
285module Merger EFlowMerger {
286# add InputArray InputArray
287 add InputArray Calorimeter/eflowTracks
288 add InputArray Calorimeter/eflowPhotons
289 add InputArray Calorimeter/eflowNeutralHadrons
290 set OutputArray eflow
291}
292
293###################
294# Photon efficiency
295###################
296
297module Efficiency PhotonEfficiency {
298 set InputArray Calorimeter/eflowPhotons
299 set OutputArray photons
300
301 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
302
303 # efficiency formula for photons
304 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
305 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
306 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
307 (abs(eta) > 2.5) * (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.5
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 <= 10.0) * (0.00) +
351 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
352 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
353 (abs(eta) > 2.5) * (0.00)}
354}
355
356####################
357# Electron isolation
358####################
359
360module Isolation ElectronIsolation {
361 set CandidateInputArray ElectronEfficiency/electrons
362 set IsolationInputArray EFlowMerger/eflow
363
364 set OutputArray electrons
365
366 set DeltaRMax 0.5
367
368 set PTMin 0.5
369
370 set PTRatioMax 0.12
371}
372
373#################
374# Muon efficiency
375#################
376
377module Efficiency MuonEfficiency {
378 set InputArray MuonMomentumSmearing/muons
379 set OutputArray muons
380
381 # set EfficiencyFormula {efficiency as a function of eta and pt}
382
383 # efficiency formula for muons
384 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
385 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
386 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
387 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
388 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
389 (abs(eta) > 2.4) * (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.5
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 20.0
468}
469
470#########################
471# Gen Missing ET merger
472########################
473
474module Merger GenMissingET {
475# add InputArray InputArray
476 add InputArray NeutrinoFilter/filteredParticles
477 set MomentumOutputArray momentum
478}
479
480
481
482############
483# Jet finder
484############
485
486module FastJetFinder FastJetFinder {
487# set InputArray Calorimeter/towers
488 set InputArray EFlowMerger/eflow
489
490 set OutputArray jets
491
492 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
493 set JetAlgorithm 6
494 set ParameterR 0.5
495
496 set JetPTMin 20.0
497}
498
499##################
500# Jet Energy Scale
501##################
502
503module EnergyScale JetEnergyScale {
504 set InputArray FastJetFinder/jets
505 set OutputArray jets
506
507 # scale formula for jets
508 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
509}
510
511########################
512# Jet Flavor Association
513########################
514
515module JetFlavorAssociation JetFlavorAssociation {
516
517 set PartonInputArray Delphes/partons
518 set ParticleInputArray Delphes/allParticles
519 set ParticleLHEFInputArray Delphes/allParticlesLHEF
520 set JetInputArray JetEnergyScale/jets
521
522 set DeltaR 0.5
523 set PartonPTMin 1.0
524 set PartonEtaMax 2.5
525
526}
527
528###########
529# b-tagging
530###########
531
532module BTagging BTagging {
533 set JetInputArray JetEnergyScale/jets
534
535 set BitNumber 0
536
537 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
538 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
539 # gluon's PDG code has the lowest priority
540
541 # based on arXiv:1211.4462
542
543 # default efficiency formula (misidentification rate)
544 add EfficiencyFormula {0} {0.01+0.000038*pt}
545
546 # efficiency formula for c-jets (misidentification rate)
547 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
548
549 # efficiency formula for b-jets
550 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
551}
552
553#############
554# tau-tagging
555#############
556
557module TauTagging TauTagging {
558 set ParticleInputArray Delphes/allParticles
559 set PartonInputArray Delphes/partons
560 set JetInputArray JetEnergyScale/jets
561
562 set DeltaR 0.5
563
564 set TauPTMin 1.0
565
566 set TauEtaMax 2.5
567
568 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
569
570 # default efficiency formula (misidentification rate)
571 add EfficiencyFormula {0} {0.01}
572 # efficiency formula for tau-jets
573 add EfficiencyFormula {15} {0.6}
574}
575
576#####################################################
577# Find uniquely identified photons/electrons/tau/jets
578#####################################################
579
580module UniqueObjectFinder UniqueObjectFinder {
581# earlier arrays take precedence over later ones
582# add InputArray InputArray OutputArray
583 add InputArray PhotonIsolation/photons photons
584 add InputArray ElectronIsolation/electrons electrons
585 add InputArray MuonIsolation/muons muons
586 add InputArray JetEnergyScale/jets jets
587}
588
589##################
590# ROOT tree writer
591##################
592
593# tracks, towers and eflow objects are not stored by default in the output.
594# if needed (for jet constituent or other studies), uncomment the relevant
595# "add Branch ..." lines.
596
597module TreeWriter TreeWriter {
598# add Branch InputArray BranchName BranchClass
599 add Branch Delphes/allParticles Particle GenParticle
600
601 add Branch TrackMerger/tracks Track Track
602 add Branch Calorimeter/towers Tower Tower
603
604 add Branch Calorimeter/eflowTracks EFlowTrack Track
605 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
606 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
607
608 add Branch GenJetFinder/jets GenJet Jet
609 add Branch GenMissingET/momentum GenMissingET MissingET
610
611 add Branch UniqueObjectFinder/jets Jet Jet
612 add Branch UniqueObjectFinder/electrons Electron Electron
613 add Branch UniqueObjectFinder/photons Photon Photon
614 add Branch UniqueObjectFinder/muons Muon Muon
615 add Branch MissingET/momentum MissingET MissingET
616 add Branch ScalarHT/energy ScalarHT ScalarHT
617}
Note: See TracBrowser for help on using the repository browser.