Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ d67b86b

ImprovedOutputFile Timing dual_readout llp
Last change on this file since d67b86b was 934d037, checked in by Pavel Demin <pavel.demin@…>, 9 years ago

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