Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ ab3bdd9

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

adding + in the ECAL formula
  • Property mode set to 100644
File size: 21.9 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6
7 PileUpMerger
8 ParticlePropagator
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronMomentumSmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19 Calorimeter
20 ElectronFilter
21 TrackPileUpSubtractor
22 NeutralTowerMerger
23 EFlowMergerAllTracks
24 EFlowMerger
25
26 NeutrinoFilter
27 GenJetFinder
28 GenMissingET
29
30 Rho
31 FastJetFinder
32 PileUpJetID
33 JetPileUpSubtractor
34
35 JetEnergyScale
36
37 PhotonEfficiency
38 PhotonIsolation
39
40 ElectronEfficiency
41 ElectronIsolation
42
43 MuonEfficiency
44 MuonIsolation
45
46 MissingET
47
48 JetFlavorAssociation
49
50 BTagging
51 TauTagging
52
53 UniqueObjectFinder
54
55 ScalarHT
56
57 TreeWriter
58}
59
60###############
61# PileUp Merger
62###############
63
64module PileUpMerger PileUpMerger {
65 set InputArray Delphes/stableParticles
66
67 set ParticleOutputArray stableParticles
68 set VertexOutputArray vertices
69
70 # pre-generated minbias input file
71 set PileUpFile MinBias.pileup
72
73 # average expected pile up
74 set MeanPileUp 50
75
76 # maximum spread in the beam direction in m
77 set ZVertexSpread 0.10
78
79 # maximum spread in time in s
80 set TVertexSpread 1.5E-09
81
82 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
83
84 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
85
86 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
87 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) +
88 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) +
89 # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
90
91
92}
93
94#################################
95# Propagate particles in cylinder
96#################################
97
98module ParticlePropagator ParticlePropagator {
99 set InputArray PileUpMerger/stableParticles
100
101 set OutputArray stableParticles
102 set ChargedHadronOutputArray chargedHadrons
103 set ElectronOutputArray electrons
104 set MuonOutputArray muons
105
106 # radius of the magnetic field coverage, in m
107 set Radius 1.29
108 # half-length of the magnetic field coverage, in m
109 set HalfLength 3.00
110
111 # magnetic field
112 set Bz 3.8
113}
114
115####################################
116# Charged hadron tracking efficiency
117####################################
118
119module Efficiency ChargedHadronTrackingEfficiency {
120 set InputArray ParticlePropagator/chargedHadrons
121 set OutputArray chargedHadrons
122
123 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
124
125 # tracking efficiency formula for charged hadrons
126 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
127 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
128 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
129 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
130 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
131 (abs(eta) > 2.5) * (0.00)}
132}
133
134##############################
135# Electron tracking efficiency
136##############################
137
138module Efficiency ElectronTrackingEfficiency {
139 set InputArray ParticlePropagator/electrons
140 set OutputArray electrons
141
142 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
143
144 # tracking efficiency formula for electrons
145 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
146 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
147 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
148 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
150 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
151 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
152 (abs(eta) > 2.5) * (0.00)}
153}
154
155##########################
156# Muon tracking efficiency
157##########################
158
159module Efficiency MuonTrackingEfficiency {
160 set InputArray ParticlePropagator/muons
161 set OutputArray muons
162
163 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
164
165 # tracking efficiency formula for muons
166 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
167 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
168 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
169 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
170 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
171 (abs(eta) > 2.5) * (0.00)}
172}
173
174########################################
175# Momentum resolution for charged tracks
176########################################
177
178module MomentumSmearing ChargedHadronMomentumSmearing {
179 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
180 set OutputArray chargedHadrons
181
182 # set ResolutionFormula {resolution formula as a function of eta and pt}
183
184 # resolution formula for charged hadrons
185 # based on arXiv:1405.6569
186 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
187 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
188 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
189}
190
191###################################
192# Momentum resolution for electrons
193###################################
194
195module MomentumSmearing ElectronMomentumSmearing {
196 set InputArray ElectronTrackingEfficiency/electrons
197 set OutputArray electrons
198
199 # set ResolutionFormula {resolution formula as a function of eta and energy}
200
201 # resolution formula for electrons
202 # based on arXiv:1405.6569
203 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
204 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
205 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
206}
207
208###############################
209# Momentum resolution for muons
210###############################
211
212module MomentumSmearing MuonMomentumSmearing {
213 set InputArray MuonTrackingEfficiency/muons
214 set OutputArray muons
215
216 # set ResolutionFormula {resolution formula as a function of eta and pt}
217
218 # resolution formula for muons
219 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.0e-4^2) +
220 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.0e-4^2) +
221 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
222}
223
224##############
225# Track merger
226##############
227
228module Merger TrackMerger {
229# add InputArray InputArray
230 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
231 add InputArray ElectronMomentumSmearing/electrons
232 add InputArray MuonMomentumSmearing/muons
233 set OutputArray tracks
234}
235
236#############
237# Calorimeter
238#############
239
240module Calorimeter Calorimeter {
241 set ParticleInputArray ParticlePropagator/stableParticles
242 set TrackInputArray TrackMerger/tracks
243
244 set TowerOutputArray towers
245 set PhotonOutputArray photons
246
247 set EFlowTrackOutputArray eflowTracks
248 set EFlowPhotonOutputArray eflowPhotons
249 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
250
251 set ECalEnergyMin 0.5
252 set HCalEnergyMin 1.0
253
254 set ECalEnergySignificanceMin 1.0
255 set HCalEnergySignificanceMin 1.0
256
257 set SmearTowerCenter true
258
259 set pi [expr {acos(-1)}]
260
261 # lists of the edges of each tower in eta and phi
262 # each list starts with the lower edge of the first tower
263 # the list ends with the higher edged of the last tower
264
265 # 5 degrees towers
266 set PhiBins {}
267 for {set i -36} {$i <= 36} {incr i} {
268 add PhiBins [expr {$i * $pi/36.0}]
269 }
270 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} {
271 add EtaPhiBins $eta $PhiBins
272 }
273
274 # 10 degrees towers
275 set PhiBins {}
276 for {set i -18} {$i <= 18} {incr i} {
277 add PhiBins [expr {$i * $pi/18.0}]
278 }
279 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} {
280 add EtaPhiBins $eta $PhiBins
281 }
282
283 # 20 degrees towers
284 set PhiBins {}
285 for {set i -9} {$i <= 9} {incr i} {
286 add PhiBins [expr {$i * $pi/9.0}]
287 }
288 foreach eta {-5 -4.7 -4.525 4.7 5} {
289 add EtaPhiBins $eta $PhiBins
290 }
291
292 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
293 add EnergyFraction {0} {0.0 1.0}
294 # energy fractions for e, gamma and pi0
295 add EnergyFraction {11} {1.0 0.0}
296 add EnergyFraction {22} {1.0 0.0}
297 add EnergyFraction {111} {1.0 0.0}
298 # energy fractions for muon, neutrinos and neutralinos
299 add EnergyFraction {12} {0.0 0.0}
300 add EnergyFraction {13} {0.0 0.0}
301 add EnergyFraction {14} {0.0 0.0}
302 add EnergyFraction {16} {0.0 0.0}
303 add EnergyFraction {1000022} {0.0 0.0}
304 add EnergyFraction {1000023} {0.0 0.0}
305 add EnergyFraction {1000025} {0.0 0.0}
306 add EnergyFraction {1000035} {0.0 0.0}
307 add EnergyFraction {1000045} {0.0 0.0}
308 # energy fractions for K0short and Lambda
309 add EnergyFraction {310} {0.3 0.7}
310 add EnergyFraction {3122} {0.3 0.7}
311
312 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
313 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
314 set ECalResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
315 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(eta-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
316 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
317
318 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
319 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
320 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
321}
322
323#################
324# Electron filter
325#################
326
327module PdgCodeFilter ElectronFilter {
328 set InputArray Calorimeter/eflowTracks
329 set OutputArray electrons
330 set Invert true
331 add PdgCode {11}
332 add PdgCode {-11}
333}
334
335##########################
336# Track pile-up subtractor
337##########################
338
339module TrackPileUpSubtractor TrackPileUpSubtractor {
340# add InputArray InputArray OutputArray
341 add InputArray Calorimeter/eflowTracks eflowTracks
342 add InputArray ElectronFilter/electrons electrons
343 add InputArray MuonMomentumSmearing/muons muons
344
345 set VertexInputArray PileUpMerger/vertices
346 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
347 # Z vertex resolution in m
348 set ZVertexResolution 0.0001
349}
350
351####################
352# Neutral tower merger
353####################
354
355module Merger NeutralTowerMerger {
356# add InputArray InputArray
357 add InputArray Calorimeter/eflowPhotons
358 add InputArray Calorimeter/eflowNeutralHadrons
359 set OutputArray eflowTowers
360}
361
362##################################
363# Energy flow merger (all tracks)
364##################################
365
366module Merger EFlowMergerAllTracks {
367# add InputArray InputArray
368 add InputArray TrackMerger/tracks
369 add InputArray Calorimeter/eflowPhotons
370 add InputArray Calorimeter/eflowNeutralHadrons
371 set OutputArray eflow
372}
373
374
375####################
376# Energy flow merger
377####################
378
379module Merger EFlowMerger {
380# add InputArray InputArray
381 add InputArray TrackPileUpSubtractor/eflowTracks
382 add InputArray Calorimeter/eflowPhotons
383 add InputArray Calorimeter/eflowNeutralHadrons
384 set OutputArray eflow
385}
386
387
388
389#############
390# Rho pile-up
391#############
392
393module FastJetGridMedianEstimator Rho {
394
395 set InputArray EFlowMerger/eflow
396 set RhoOutputArray rho
397
398 # add GridRange rapmin rapmax drap dphi
399 # rapmin - the minimum rapidity extent of the grid
400 # rapmax - the maximum rapidity extent of the grid
401 # drap - the grid spacing in rapidity
402 # dphi - the grid spacing in azimuth
403
404 add GridRange -5.0 -2.5 1.0 1.0
405 add GridRange -2.5 2.5 1.0 1.0
406 add GridRange 2.5 5.0 1.0 1.0
407
408}
409
410#####################
411# Neutrino Filter
412#####################
413
414module PdgCodeFilter NeutrinoFilter {
415
416 set InputArray Delphes/stableParticles
417 set OutputArray filteredParticles
418
419 set PTMin 0.0
420
421 add PdgCode {12}
422 add PdgCode {14}
423 add PdgCode {16}
424 add PdgCode {-12}
425 add PdgCode {-14}
426 add PdgCode {-16}
427
428}
429
430
431
432#####################
433# MC truth jet finder
434#####################
435
436module FastJetFinder GenJetFinder {
437 set InputArray NeutrinoFilter/filteredParticles
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# Gen Missing ET merger
450########################
451
452module Merger GenMissingET {
453# add InputArray InputArray
454 add InputArray NeutrinoFilter/filteredParticles
455 set MomentumOutputArray momentum
456}
457
458
459
460############
461# Jet finder
462############
463
464module FastJetFinder FastJetFinder {
465# set InputArray Calorimeter/towers
466 set InputArray EFlowMerger/eflow
467
468 set OutputArray jets
469
470 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
471 set AreaAlgorithm 5
472
473 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
474 set JetAlgorithm 6
475 set ParameterR 0.5
476
477 set JetPTMin 20.0
478}
479
480###########################
481# Jet Pile-Up ID
482###########################
483
484module PileUpJetID PileUpJetID {
485 set JetInputArray FastJetFinder/jets
486 set TrackInputArray Calorimeter/eflowTracks
487 set NeutralInputArray NeutralTowerMerger/eflowTowers
488
489 set VertexInputArray PileUpMerger/vertices
490 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
491 # Z vertex resolution in m
492 set ZVertexResolution 0.0001
493
494 set OutputArray jets
495
496 set UseConstituents 0
497 set ParameterR 0.5
498
499 set JetPTMin 20.0
500}
501
502###########################
503# Jet Pile-Up Subtraction
504###########################
505
506module JetPileUpSubtractor JetPileUpSubtractor {
507 set JetInputArray PileUpJetID/jets
508 set RhoInputArray Rho/rho
509
510 set OutputArray jets
511
512 set JetPTMin 20.0
513}
514
515##################
516# Jet Energy Scale
517##################
518
519module EnergyScale JetEnergyScale {
520 set InputArray JetPileUpSubtractor/jets
521 set OutputArray jets
522
523 # scale formula for jets
524 set ScaleFormula {1.0}
525}
526
527###################
528# Photon efficiency
529###################
530
531module Efficiency PhotonEfficiency {
532 set InputArray Calorimeter/eflowPhotons
533 set OutputArray photons
534
535 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
536
537 # efficiency formula for photons
538 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
539 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
540 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
541 (abs(eta) > 2.5) * (0.00)}
542}
543
544
545##################
546# Photon isolation
547##################
548
549module Isolation PhotonIsolation {
550 set CandidateInputArray PhotonEfficiency/photons
551 set IsolationInputArray EFlowMergerAllTracks/eflow
552 set RhoInputArray Rho/rho
553
554 set OutputArray photons
555
556 set DeltaRMax 0.5
557
558 set PTMin 0.5
559
560 set PTRatioMax 0.12
561}
562
563#####################
564# Electron efficiency
565#####################
566
567module Efficiency ElectronEfficiency {
568 set InputArray TrackPileUpSubtractor/electrons
569 set OutputArray electrons
570
571 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
572
573 # efficiency formula for electrons
574 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
575 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
576 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
577 (abs(eta) > 2.5) * (0.00)}
578}
579
580####################
581# Electron isolation
582####################
583
584module Isolation ElectronIsolation {
585 set CandidateInputArray ElectronEfficiency/electrons
586 set IsolationInputArray EFlowMergerAllTracks/eflow
587 set RhoInputArray Rho/rho
588
589 set OutputArray electrons
590
591 set DeltaRMax 0.5
592
593 set PTMin 0.5
594
595 set PTRatioMax 0.12
596}
597
598#################
599# Muon efficiency
600#################
601
602module Efficiency MuonEfficiency {
603 set InputArray TrackPileUpSubtractor/muons
604 set OutputArray muons
605
606 # set EfficiencyFormula {efficiency as a function of eta and pt}
607
608 # efficiency formula for muons
609 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
610 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
611 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
612 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
613 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
614 (abs(eta) > 2.4) * (0.00)}
615}
616
617################
618# Muon isolation
619################
620
621module Isolation MuonIsolation {
622 set CandidateInputArray MuonEfficiency/muons
623 set IsolationInputArray EFlowMergerAllTracks/eflow
624 set RhoInputArray Rho/rho
625
626 set OutputArray muons
627
628 set DeltaRMax 0.5
629
630 set PTMin 0.5
631
632 set PTRatioMax 0.25
633}
634
635###################
636# Missing ET merger
637###################
638
639module Merger MissingET {
640# add InputArray InputArray
641 add InputArray EFlowMergerAllTracks/eflow
642 set MomentumOutputArray momentum
643}
644
645
646
647##################
648# Scalar HT merger
649##################
650
651module Merger ScalarHT {
652# add InputArray InputArray
653 add InputArray UniqueObjectFinder/jets
654 add InputArray UniqueObjectFinder/electrons
655 add InputArray UniqueObjectFinder/photons
656 add InputArray UniqueObjectFinder/muons
657 set EnergyOutputArray energy
658}
659
660########################
661# Jet Flavor Association
662########################
663
664module JetFlavorAssociation JetFlavorAssociation {
665
666 set PartonInputArray Delphes/partons
667 set ParticleInputArray Delphes/allParticles
668 set ParticleLHEFInputArray Delphes/allParticlesLHEF
669 set JetInputArray JetEnergyScale/jets
670
671 set DeltaR 0.5
672 set PartonPTMin 1.0
673 set PartonEtaMax 2.5
674
675}
676
677###########
678# b-tagging
679###########
680
681module BTagging BTagging {
682 set JetInputArray JetEnergyScale/jets
683
684 set BitNumber 0
685
686 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
687 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
688 # gluon's PDG code has the lowest priority
689
690 # based on arXiv:1211.4462
691
692 # default efficiency formula (misidentification rate)
693 add EfficiencyFormula {0} {0.01+0.000038*pt}
694
695 # efficiency formula for c-jets (misidentification rate)
696 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
697
698 # efficiency formula for b-jets
699 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
700}
701
702#############
703# tau-tagging
704#############
705
706module TauTagging TauTagging {
707 set ParticleInputArray Delphes/allParticles
708 set PartonInputArray Delphes/partons
709 set JetInputArray JetEnergyScale/jets
710
711 set DeltaR 0.5
712
713 set TauPTMin 1.0
714
715 set TauEtaMax 2.5
716
717 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
718
719 # default efficiency formula (misidentification rate)
720 add EfficiencyFormula {0} {0.01}
721 # efficiency formula for tau-jets
722 add EfficiencyFormula {15} {0.6}
723}
724
725#####################################################
726# Find uniquely identified photons/electrons/tau/jets
727#####################################################
728
729module UniqueObjectFinder UniqueObjectFinder {
730# earlier arrays take precedence over later ones
731# add InputArray InputArray OutputArray
732 add InputArray PhotonIsolation/photons photons
733 add InputArray ElectronIsolation/electrons electrons
734 add InputArray MuonIsolation/muons muons
735 add InputArray JetEnergyScale/jets jets
736}
737
738##################
739# ROOT tree writer
740##################
741
742# tracks, towers and eflow objects are not stored by default in the output.
743# if needed (for jet constituent or other studies), uncomment the relevant
744# "add Branch ..." lines.
745
746module TreeWriter TreeWriter {
747# add Branch InputArray BranchName BranchClass
748 add Branch Delphes/allParticles Particle GenParticle
749
750# add Branch TrackMerger/tracks Track Track
751# add Branch Calorimeter/towers Tower Tower
752
753# add Branch Calorimeter/eflowTracks EFlowTrack Track
754# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
755# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
756
757 add Branch GenJetFinder/jets GenJet Jet
758 add Branch GenMissingET/momentum GenMissingET MissingET
759
760 add Branch UniqueObjectFinder/jets Jet Jet
761 add Branch UniqueObjectFinder/electrons Electron Electron
762 add Branch UniqueObjectFinder/photons Photon Photon
763 add Branch UniqueObjectFinder/muons Muon Muon
764 add Branch MissingET/momentum MissingET MissingET
765 add Branch ScalarHT/energy ScalarHT ScalarHT
766 add Branch Rho/rho Rho Rho
767 add Branch PileUpMerger/vertices Vertex Vertex
768}
Note: See TracBrowser for help on using the repository browser.