Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ 8f0b34c

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 8f0b34c was eab6d46, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

apply muon eff drop in tracking
  • Property mode set to 100644
File size: 21.6 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.25
78
79 # maximum spread in time in s
80 set TVertexSpread 800E-12
81
82 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
83 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
84
85
86}
87
88#################################
89# Propagate particles in cylinder
90#################################
91
92module ParticlePropagator ParticlePropagator {
93 set InputArray PileUpMerger/stableParticles
94
95 set OutputArray stableParticles
96 set ChargedHadronOutputArray chargedHadrons
97 set ElectronOutputArray electrons
98 set MuonOutputArray muons
99
100 # radius of the magnetic field coverage, in m
101 set Radius 1.29
102 # half-length of the magnetic field coverage, in m
103 set HalfLength 3.00
104
105 # magnetic field
106 set Bz 3.8
107}
108
109####################################
110# Charged hadron tracking efficiency
111####################################
112
113module Efficiency ChargedHadronTrackingEfficiency {
114 set InputArray ParticlePropagator/chargedHadrons
115 set OutputArray chargedHadrons
116
117 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
118
119 # tracking efficiency formula for charged hadrons
120 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
121 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
122 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
124 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
125 (abs(eta) > 2.5) * (0.00)}
126}
127
128##############################
129# Electron tracking efficiency
130##############################
131
132module Efficiency ElectronTrackingEfficiency {
133 set InputArray ParticlePropagator/electrons
134 set OutputArray electrons
135
136 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
137
138 # tracking efficiency formula for electrons
139 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
140 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
141 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
142 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
146 (abs(eta) > 2.5) * (0.00)}
147}
148
149##########################
150# Muon tracking efficiency
151##########################
152
153module Efficiency MuonTrackingEfficiency {
154 set InputArray ParticlePropagator/muons
155 set OutputArray muons
156
157 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
158
159 # tracking efficiency formula for muons
160 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
161 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
162 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
163 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
164
165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
166 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
167 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
168 (abs(eta) > 2.5) * (0.00)}
169}
170
171########################################
172# Momentum resolution for charged tracks
173########################################
174
175module MomentumSmearing ChargedHadronMomentumSmearing {
176 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
177 set OutputArray chargedHadrons
178
179 # set ResolutionFormula {resolution formula as a function of eta and pt}
180
181 # resolution formula for charged hadrons
182 # based on arXiv:1405.6569
183 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
184 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
186}
187
188###################################
189# Momentum resolution for electrons
190###################################
191
192module MomentumSmearing ElectronMomentumSmearing {
193 set InputArray ElectronTrackingEfficiency/electrons
194 set OutputArray electrons
195
196 # set ResolutionFormula {resolution formula as a function of eta and energy}
197
198 # resolution formula for electrons
199 # based on arXiv:1405.6569
200 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
201 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
202 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
203}
204
205###############################
206# Momentum resolution for muons
207###############################
208
209module MomentumSmearing MuonMomentumSmearing {
210 set InputArray MuonTrackingEfficiency/muons
211 set OutputArray muons
212
213 # set ResolutionFormula {resolution formula as a function of eta and pt}
214
215 # resolution formula for muons
216 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
217 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
218 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
219}
220
221##############
222# Track merger
223##############
224
225module Merger TrackMerger {
226# add InputArray InputArray
227 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
228 add InputArray ElectronMomentumSmearing/electrons
229 add InputArray MuonMomentumSmearing/muons
230 set OutputArray tracks
231}
232
233#############
234# Calorimeter
235#############
236
237module Calorimeter Calorimeter {
238 set ParticleInputArray ParticlePropagator/stableParticles
239 set TrackInputArray TrackMerger/tracks
240
241 set TowerOutputArray towers
242 set PhotonOutputArray photons
243
244 set EFlowTrackOutputArray eflowTracks
245 set EFlowPhotonOutputArray eflowPhotons
246 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
247
248 set ECalEnergyMin 0.5
249 set HCalEnergyMin 1.0
250
251 set ECalEnergySignificanceMin 1.0
252 set HCalEnergySignificanceMin 1.0
253
254 set SmearTowerCenter true
255
256 set pi [expr {acos(-1)}]
257
258 # lists of the edges of each tower in eta and phi
259 # each list starts with the lower edge of the first tower
260 # the list ends with the higher edged of the last tower
261
262 # 5 degrees towers
263 set PhiBins {}
264 for {set i -36} {$i <= 36} {incr i} {
265 add PhiBins [expr {$i * $pi/36.0}]
266 }
267 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} {
268 add EtaPhiBins $eta $PhiBins
269 }
270
271 # 10 degrees towers
272 set PhiBins {}
273 for {set i -18} {$i <= 18} {incr i} {
274 add PhiBins [expr {$i * $pi/18.0}]
275 }
276 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} {
277 add EtaPhiBins $eta $PhiBins
278 }
279
280 # 20 degrees towers
281 set PhiBins {}
282 for {set i -9} {$i <= 9} {incr i} {
283 add PhiBins [expr {$i * $pi/9.0}]
284 }
285 foreach eta {-5 -4.7 -4.525 4.7 5} {
286 add EtaPhiBins $eta $PhiBins
287 }
288
289 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
290 add EnergyFraction {0} {0.0 1.0}
291 # energy fractions for e, gamma and pi0
292 add EnergyFraction {11} {1.0 0.0}
293 add EnergyFraction {22} {1.0 0.0}
294 add EnergyFraction {111} {1.0 0.0}
295 # energy fractions for muon, neutrinos and neutralinos
296 add EnergyFraction {12} {0.0 0.0}
297 add EnergyFraction {13} {0.0 0.0}
298 add EnergyFraction {14} {0.0 0.0}
299 add EnergyFraction {16} {0.0 0.0}
300 add EnergyFraction {1000022} {0.0 0.0}
301 add EnergyFraction {1000023} {0.0 0.0}
302 add EnergyFraction {1000025} {0.0 0.0}
303 add EnergyFraction {1000035} {0.0 0.0}
304 add EnergyFraction {1000045} {0.0 0.0}
305 # energy fractions for K0short and Lambda
306 add EnergyFraction {310} {0.3 0.7}
307 add EnergyFraction {3122} {0.3 0.7}
308
309 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
310 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
311 set ECalResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
312 (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) +
313 (abs(eta) > 2.5 && 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# Gen Missing ET merger
447########################
448
449module Merger GenMissingET {
450# add InputArray InputArray
451 add InputArray NeutrinoFilter/filteredParticles
452 set MomentumOutputArray momentum
453}
454
455
456
457############
458# Jet finder
459############
460
461module FastJetFinder FastJetFinder {
462# set InputArray Calorimeter/towers
463 set InputArray EFlowMerger/eflow
464
465 set OutputArray jets
466
467 # 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
468 set AreaAlgorithm 5
469
470 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
471 set JetAlgorithm 6
472 set ParameterR 0.5
473
474 set JetPTMin 20.0
475}
476
477###########################
478# Jet Pile-Up ID
479###########################
480
481module PileUpJetID PileUpJetID {
482 set JetInputArray FastJetFinder/jets
483 set TrackInputArray Calorimeter/eflowTracks
484 set NeutralInputArray NeutralTowerMerger/eflowTowers
485
486 set VertexInputArray PileUpMerger/vertices
487 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
488 # Z vertex resolution in m
489 set ZVertexResolution 0.0001
490
491 set OutputArray jets
492
493 set UseConstituents 0
494 set ParameterR 0.5
495
496 set JetPTMin 20.0
497}
498
499###########################
500# Jet Pile-Up Subtraction
501###########################
502
503module JetPileUpSubtractor JetPileUpSubtractor {
504 set JetInputArray PileUpJetID/jets
505 set RhoInputArray Rho/rho
506
507 set OutputArray jets
508
509 set JetPTMin 20.0
510}
511
512##################
513# Jet Energy Scale
514##################
515
516module EnergyScale JetEnergyScale {
517 set InputArray JetPileUpSubtractor/jets
518 set OutputArray jets
519
520 # scale formula for jets
521 set ScaleFormula {1.0}
522}
523
524###################
525# Photon efficiency
526###################
527
528module Efficiency PhotonEfficiency {
529 set InputArray Calorimeter/eflowPhotons
530 set OutputArray photons
531
532 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
533
534 # efficiency formula for photons
535 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
536 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
537 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
538 (abs(eta) > 2.5) * (0.00)}
539}
540
541
542##################
543# Photon isolation
544##################
545
546module Isolation PhotonIsolation {
547 set CandidateInputArray PhotonEfficiency/photons
548 set IsolationInputArray EFlowMergerAllTracks/eflow
549 set RhoInputArray Rho/rho
550
551 set OutputArray photons
552
553 set DeltaRMax 0.5
554
555 set PTMin 0.5
556
557 set PTRatioMax 0.12
558}
559
560#####################
561# Electron efficiency
562#####################
563
564module Efficiency ElectronEfficiency {
565 set InputArray TrackPileUpSubtractor/electrons
566 set OutputArray electrons
567
568 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
569
570 # efficiency formula for electrons
571 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
572 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
573 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
574 (abs(eta) > 2.5) * (0.00)}
575}
576
577####################
578# Electron isolation
579####################
580
581module Isolation ElectronIsolation {
582 set CandidateInputArray ElectronEfficiency/electrons
583 set IsolationInputArray EFlowMergerAllTracks/eflow
584 set RhoInputArray Rho/rho
585
586 set OutputArray electrons
587
588 set DeltaRMax 0.5
589
590 set PTMin 0.5
591
592 set PTRatioMax 0.12
593}
594
595#################
596# Muon efficiency
597#################
598
599module Efficiency MuonEfficiency {
600 set InputArray TrackPileUpSubtractor/muons
601 set OutputArray muons
602
603 # set EfficiencyFormula {efficiency as a function of eta and pt}
604
605 # efficiency formula for muons
606 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
607 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
608 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
609 (abs(eta) > 2.4) * (0.00)}
610
611}
612
613################
614# Muon isolation
615################
616
617module Isolation MuonIsolation {
618 set CandidateInputArray MuonEfficiency/muons
619 set IsolationInputArray EFlowMergerAllTracks/eflow
620 set RhoInputArray Rho/rho
621
622 set OutputArray muons
623
624 set DeltaRMax 0.5
625
626 set PTMin 0.5
627
628 set PTRatioMax 0.25
629}
630
631###################
632# Missing ET merger
633###################
634
635module Merger MissingET {
636# add InputArray InputArray
637 add InputArray EFlowMergerAllTracks/eflow
638 set MomentumOutputArray momentum
639}
640
641
642
643##################
644# Scalar HT merger
645##################
646
647module Merger ScalarHT {
648# add InputArray InputArray
649 add InputArray UniqueObjectFinder/jets
650 add InputArray UniqueObjectFinder/electrons
651 add InputArray UniqueObjectFinder/photons
652 add InputArray UniqueObjectFinder/muons
653 set EnergyOutputArray energy
654}
655
656########################
657# Jet Flavor Association
658########################
659
660module JetFlavorAssociation JetFlavorAssociation {
661
662 set PartonInputArray Delphes/partons
663 set ParticleInputArray Delphes/allParticles
664 set ParticleLHEFInputArray Delphes/allParticlesLHEF
665 set JetInputArray JetEnergyScale/jets
666
667 set DeltaR 0.5
668 set PartonPTMin 1.0
669 set PartonEtaMax 2.5
670
671}
672
673###########
674# b-tagging
675###########
676
677module BTagging BTagging {
678 set JetInputArray JetEnergyScale/jets
679
680 set BitNumber 0
681
682 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
683 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
684 # gluon's PDG code has the lowest priority
685
686 # based on arXiv:1211.4462
687
688 # default efficiency formula (misidentification rate)
689 add EfficiencyFormula {0} {0.01+0.000038*pt}
690
691 # efficiency formula for c-jets (misidentification rate)
692 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
693
694 # efficiency formula for b-jets
695 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
696}
697
698#############
699# tau-tagging
700#############
701
702module TauTagging TauTagging {
703 set ParticleInputArray Delphes/allParticles
704 set PartonInputArray Delphes/partons
705 set JetInputArray JetEnergyScale/jets
706
707 set DeltaR 0.5
708
709 set TauPTMin 1.0
710
711 set TauEtaMax 2.5
712
713 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
714
715 # default efficiency formula (misidentification rate)
716 add EfficiencyFormula {0} {0.01}
717 # efficiency formula for tau-jets
718 add EfficiencyFormula {15} {0.6}
719}
720
721#####################################################
722# Find uniquely identified photons/electrons/tau/jets
723#####################################################
724
725module UniqueObjectFinder UniqueObjectFinder {
726# earlier arrays take precedence over later ones
727# add InputArray InputArray OutputArray
728 add InputArray PhotonIsolation/photons photons
729 add InputArray ElectronIsolation/electrons electrons
730 add InputArray MuonIsolation/muons muons
731 add InputArray JetEnergyScale/jets jets
732}
733
734##################
735# ROOT tree writer
736##################
737
738# tracks, towers and eflow objects are not stored by default in the output.
739# if needed (for jet constituent or other studies), uncomment the relevant
740# "add Branch ..." lines.
741
742module TreeWriter TreeWriter {
743# add Branch InputArray BranchName BranchClass
744 add Branch Delphes/allParticles Particle GenParticle
745
746# add Branch TrackMerger/tracks Track Track
747# add Branch Calorimeter/towers Tower Tower
748
749# add Branch Calorimeter/eflowTracks EFlowTrack Track
750# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
751# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
752
753 add Branch GenJetFinder/jets GenJet Jet
754 add Branch GenMissingET/momentum GenMissingET MissingET
755
756 add Branch UniqueObjectFinder/jets Jet Jet
757 add Branch UniqueObjectFinder/electrons Electron Electron
758 add Branch UniqueObjectFinder/photons Photon Photon
759 add Branch UniqueObjectFinder/muons Muon Muon
760 add Branch MissingET/momentum MissingET MissingET
761 add Branch ScalarHT/energy ScalarHT ScalarHT
762 add Branch Rho/rho Rho Rho
763 add Branch PileUpMerger/vertices Vertex Vertex
764}
Note: See TracBrowser for help on using the repository browser.