Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ a3a556c

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

put correct bs in pu cards
  • 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) * (0.99) +
163 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
164 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
165 (abs(eta) > 2.5) * (0.00)}
166}
167
168########################################
169# Momentum resolution for charged tracks
170########################################
171
172module MomentumSmearing ChargedHadronMomentumSmearing {
173 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
174 set OutputArray chargedHadrons
175
176 # set ResolutionFormula {resolution formula as a function of eta and pt}
177
178 # resolution formula for charged hadrons
179 # based on arXiv:1405.6569
180 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
181 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
183}
184
185###################################
186# Momentum resolution for electrons
187###################################
188
189module MomentumSmearing ElectronMomentumSmearing {
190 set InputArray ElectronTrackingEfficiency/electrons
191 set OutputArray electrons
192
193 # set ResolutionFormula {resolution formula as a function of eta and energy}
194
195 # resolution formula for electrons
196 # based on arXiv:1405.6569
197 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
198 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
199 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
200}
201
202###############################
203# Momentum resolution for muons
204###############################
205
206module MomentumSmearing MuonMomentumSmearing {
207 set InputArray MuonTrackingEfficiency/muons
208 set OutputArray muons
209
210 # set ResolutionFormula {resolution formula as a function of eta and pt}
211
212 # resolution formula for muons
213 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
214 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
215 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
216}
217
218##############
219# Track merger
220##############
221
222module Merger TrackMerger {
223# add InputArray InputArray
224 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
225 add InputArray ElectronMomentumSmearing/electrons
226 add InputArray MuonMomentumSmearing/muons
227 set OutputArray tracks
228}
229
230#############
231# Calorimeter
232#############
233
234module Calorimeter Calorimeter {
235 set ParticleInputArray ParticlePropagator/stableParticles
236 set TrackInputArray TrackMerger/tracks
237
238 set TowerOutputArray towers
239 set PhotonOutputArray photons
240
241 set EFlowTrackOutputArray eflowTracks
242 set EFlowPhotonOutputArray eflowPhotons
243 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
244
245 set ECalEnergyMin 0.5
246 set HCalEnergyMin 1.0
247
248 set ECalEnergySignificanceMin 1.0
249 set HCalEnergySignificanceMin 1.0
250
251 set SmearTowerCenter true
252
253 set pi [expr {acos(-1)}]
254
255 # lists of the edges of each tower in eta and phi
256 # each list starts with the lower edge of the first tower
257 # the list ends with the higher edged of the last tower
258
259 # 5 degrees towers
260 set PhiBins {}
261 for {set i -36} {$i <= 36} {incr i} {
262 add PhiBins [expr {$i * $pi/36.0}]
263 }
264 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} {
265 add EtaPhiBins $eta $PhiBins
266 }
267
268 # 10 degrees towers
269 set PhiBins {}
270 for {set i -18} {$i <= 18} {incr i} {
271 add PhiBins [expr {$i * $pi/18.0}]
272 }
273 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} {
274 add EtaPhiBins $eta $PhiBins
275 }
276
277 # 20 degrees towers
278 set PhiBins {}
279 for {set i -9} {$i <= 9} {incr i} {
280 add PhiBins [expr {$i * $pi/9.0}]
281 }
282 foreach eta {-5 -4.7 -4.525 4.7 5} {
283 add EtaPhiBins $eta $PhiBins
284 }
285
286 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
287 add EnergyFraction {0} {0.0 1.0}
288 # energy fractions for e, gamma and pi0
289 add EnergyFraction {11} {1.0 0.0}
290 add EnergyFraction {22} {1.0 0.0}
291 add EnergyFraction {111} {1.0 0.0}
292 # energy fractions for muon, neutrinos and neutralinos
293 add EnergyFraction {12} {0.0 0.0}
294 add EnergyFraction {13} {0.0 0.0}
295 add EnergyFraction {14} {0.0 0.0}
296 add EnergyFraction {16} {0.0 0.0}
297 add EnergyFraction {1000022} {0.0 0.0}
298 add EnergyFraction {1000023} {0.0 0.0}
299 add EnergyFraction {1000025} {0.0 0.0}
300 add EnergyFraction {1000035} {0.0 0.0}
301 add EnergyFraction {1000045} {0.0 0.0}
302 # energy fractions for K0short and Lambda
303 add EnergyFraction {310} {0.3 0.7}
304 add EnergyFraction {3122} {0.3 0.7}
305
306 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
307 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
308 set ECalResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
309 (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) +
310 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
311
312 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
313 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
314 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
315}
316
317#################
318# Electron filter
319#################
320
321module PdgCodeFilter ElectronFilter {
322 set InputArray Calorimeter/eflowTracks
323 set OutputArray electrons
324 set Invert true
325 add PdgCode {11}
326 add PdgCode {-11}
327}
328
329##########################
330# Track pile-up subtractor
331##########################
332
333module TrackPileUpSubtractor TrackPileUpSubtractor {
334# add InputArray InputArray OutputArray
335 add InputArray Calorimeter/eflowTracks eflowTracks
336 add InputArray ElectronFilter/electrons electrons
337 add InputArray MuonMomentumSmearing/muons muons
338
339 set VertexInputArray PileUpMerger/vertices
340 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
341 # Z vertex resolution in m
342 set ZVertexResolution 0.0001
343}
344
345####################
346# Neutral tower merger
347####################
348
349module Merger NeutralTowerMerger {
350# add InputArray InputArray
351 add InputArray Calorimeter/eflowPhotons
352 add InputArray Calorimeter/eflowNeutralHadrons
353 set OutputArray eflowTowers
354}
355
356##################################
357# Energy flow merger (all tracks)
358##################################
359
360module Merger EFlowMergerAllTracks {
361# add InputArray InputArray
362 add InputArray TrackMerger/tracks
363 add InputArray Calorimeter/eflowPhotons
364 add InputArray Calorimeter/eflowNeutralHadrons
365 set OutputArray eflow
366}
367
368
369####################
370# Energy flow merger
371####################
372
373module Merger EFlowMerger {
374# add InputArray InputArray
375 add InputArray TrackPileUpSubtractor/eflowTracks
376 add InputArray Calorimeter/eflowPhotons
377 add InputArray Calorimeter/eflowNeutralHadrons
378 set OutputArray eflow
379}
380
381
382
383#############
384# Rho pile-up
385#############
386
387module FastJetGridMedianEstimator Rho {
388
389 set InputArray EFlowMerger/eflow
390 set RhoOutputArray rho
391
392 # add GridRange rapmin rapmax drap dphi
393 # rapmin - the minimum rapidity extent of the grid
394 # rapmax - the maximum rapidity extent of the grid
395 # drap - the grid spacing in rapidity
396 # dphi - the grid spacing in azimuth
397
398 add GridRange -5.0 -2.5 1.0 1.0
399 add GridRange -2.5 2.5 1.0 1.0
400 add GridRange 2.5 5.0 1.0 1.0
401
402}
403
404#####################
405# Neutrino Filter
406#####################
407
408module PdgCodeFilter NeutrinoFilter {
409
410 set InputArray Delphes/stableParticles
411 set OutputArray filteredParticles
412
413 set PTMin 0.0
414
415 add PdgCode {12}
416 add PdgCode {14}
417 add PdgCode {16}
418 add PdgCode {-12}
419 add PdgCode {-14}
420 add PdgCode {-16}
421
422}
423
424
425
426#####################
427# MC truth jet finder
428#####################
429
430module FastJetFinder GenJetFinder {
431 set InputArray NeutrinoFilter/filteredParticles
432
433 set OutputArray jets
434
435 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
436 set JetAlgorithm 6
437 set ParameterR 0.5
438
439 set JetPTMin 20.0
440}
441
442#########################
443# Gen Missing ET merger
444########################
445
446module Merger GenMissingET {
447# add InputArray InputArray
448 add InputArray NeutrinoFilter/filteredParticles
449 set MomentumOutputArray momentum
450}
451
452
453
454############
455# Jet finder
456############
457
458module FastJetFinder FastJetFinder {
459# set InputArray Calorimeter/towers
460 set InputArray EFlowMerger/eflow
461
462 set OutputArray jets
463
464 # 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
465 set AreaAlgorithm 5
466
467 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
468 set JetAlgorithm 6
469 set ParameterR 0.5
470
471 set JetPTMin 20.0
472}
473
474###########################
475# Jet Pile-Up ID
476###########################
477
478module PileUpJetID PileUpJetID {
479 set JetInputArray FastJetFinder/jets
480 set TrackInputArray Calorimeter/eflowTracks
481 set NeutralInputArray NeutralTowerMerger/eflowTowers
482
483 set VertexInputArray PileUpMerger/vertices
484 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
485 # Z vertex resolution in m
486 set ZVertexResolution 0.0001
487
488 set OutputArray jets
489
490 set UseConstituents 0
491 set ParameterR 0.5
492
493 set JetPTMin 20.0
494}
495
496###########################
497# Jet Pile-Up Subtraction
498###########################
499
500module JetPileUpSubtractor JetPileUpSubtractor {
501 set JetInputArray PileUpJetID/jets
502 set RhoInputArray Rho/rho
503
504 set OutputArray jets
505
506 set JetPTMin 20.0
507}
508
509##################
510# Jet Energy Scale
511##################
512
513module EnergyScale JetEnergyScale {
514 set InputArray JetPileUpSubtractor/jets
515 set OutputArray jets
516
517 # scale formula for jets
518 set ScaleFormula {1.0}
519}
520
521###################
522# Photon efficiency
523###################
524
525module Efficiency PhotonEfficiency {
526 set InputArray Calorimeter/eflowPhotons
527 set OutputArray photons
528
529 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
530
531 # efficiency formula for photons
532 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
533 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
534 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
535 (abs(eta) > 2.5) * (0.00)}
536}
537
538
539##################
540# Photon isolation
541##################
542
543module Isolation PhotonIsolation {
544 set CandidateInputArray PhotonEfficiency/photons
545 set IsolationInputArray EFlowMergerAllTracks/eflow
546 set RhoInputArray Rho/rho
547
548 set OutputArray photons
549
550 set DeltaRMax 0.5
551
552 set PTMin 0.5
553
554 set PTRatioMax 0.12
555}
556
557#####################
558# Electron efficiency
559#####################
560
561module Efficiency ElectronEfficiency {
562 set InputArray TrackPileUpSubtractor/electrons
563 set OutputArray electrons
564
565 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
566
567 # efficiency formula for electrons
568 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
569 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
570 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
571 (abs(eta) > 2.5) * (0.00)}
572}
573
574####################
575# Electron isolation
576####################
577
578module Isolation ElectronIsolation {
579 set CandidateInputArray ElectronEfficiency/electrons
580 set IsolationInputArray EFlowMergerAllTracks/eflow
581 set RhoInputArray Rho/rho
582
583 set OutputArray electrons
584
585 set DeltaRMax 0.5
586
587 set PTMin 0.5
588
589 set PTRatioMax 0.12
590}
591
592#################
593# Muon efficiency
594#################
595
596module Efficiency MuonEfficiency {
597 set InputArray TrackPileUpSubtractor/muons
598 set OutputArray muons
599
600 # set EfficiencyFormula {efficiency as a function of eta and pt}
601
602 # efficiency formula for muons
603 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
604 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
605 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
606 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
607 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
608 (abs(eta) > 2.4) * (0.00)}
609}
610
611################
612# Muon isolation
613################
614
615module Isolation MuonIsolation {
616 set CandidateInputArray MuonEfficiency/muons
617 set IsolationInputArray EFlowMergerAllTracks/eflow
618 set RhoInputArray Rho/rho
619
620 set OutputArray muons
621
622 set DeltaRMax 0.5
623
624 set PTMin 0.5
625
626 set PTRatioMax 0.25
627}
628
629###################
630# Missing ET merger
631###################
632
633module Merger MissingET {
634# add InputArray InputArray
635 add InputArray EFlowMergerAllTracks/eflow
636 set MomentumOutputArray momentum
637}
638
639
640
641##################
642# Scalar HT merger
643##################
644
645module Merger ScalarHT {
646# add InputArray InputArray
647 add InputArray UniqueObjectFinder/jets
648 add InputArray UniqueObjectFinder/electrons
649 add InputArray UniqueObjectFinder/photons
650 add InputArray UniqueObjectFinder/muons
651 set EnergyOutputArray energy
652}
653
654########################
655# Jet Flavor Association
656########################
657
658module JetFlavorAssociation JetFlavorAssociation {
659
660 set PartonInputArray Delphes/partons
661 set ParticleInputArray Delphes/allParticles
662 set ParticleLHEFInputArray Delphes/allParticlesLHEF
663 set JetInputArray JetEnergyScale/jets
664
665 set DeltaR 0.5
666 set PartonPTMin 1.0
667 set PartonEtaMax 2.5
668
669}
670
671###########
672# b-tagging
673###########
674
675module BTagging BTagging {
676 set JetInputArray JetEnergyScale/jets
677
678 set BitNumber 0
679
680 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
681 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
682 # gluon's PDG code has the lowest priority
683
684 # based on arXiv:1211.4462
685
686 # default efficiency formula (misidentification rate)
687 add EfficiencyFormula {0} {0.01+0.000038*pt}
688
689 # efficiency formula for c-jets (misidentification rate)
690 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
691
692 # efficiency formula for b-jets
693 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
694}
695
696#############
697# tau-tagging
698#############
699
700module TauTagging TauTagging {
701 set ParticleInputArray Delphes/allParticles
702 set PartonInputArray Delphes/partons
703 set JetInputArray JetEnergyScale/jets
704
705 set DeltaR 0.5
706
707 set TauPTMin 1.0
708
709 set TauEtaMax 2.5
710
711 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
712
713 # default efficiency formula (misidentification rate)
714 add EfficiencyFormula {0} {0.01}
715 # efficiency formula for tau-jets
716 add EfficiencyFormula {15} {0.6}
717}
718
719#####################################################
720# Find uniquely identified photons/electrons/tau/jets
721#####################################################
722
723module UniqueObjectFinder UniqueObjectFinder {
724# earlier arrays take precedence over later ones
725# add InputArray InputArray OutputArray
726 add InputArray PhotonIsolation/photons photons
727 add InputArray ElectronIsolation/electrons electrons
728 add InputArray MuonIsolation/muons muons
729 add InputArray JetEnergyScale/jets jets
730}
731
732##################
733# ROOT tree writer
734##################
735
736# tracks, towers and eflow objects are not stored by default in the output.
737# if needed (for jet constituent or other studies), uncomment the relevant
738# "add Branch ..." lines.
739
740module TreeWriter TreeWriter {
741# add Branch InputArray BranchName BranchClass
742 add Branch Delphes/allParticles Particle GenParticle
743
744# add Branch TrackMerger/tracks Track Track
745# add Branch Calorimeter/towers Tower Tower
746
747# add Branch Calorimeter/eflowTracks EFlowTrack Track
748# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
749# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
750
751 add Branch GenJetFinder/jets GenJet Jet
752 add Branch GenMissingET/momentum GenMissingET MissingET
753
754 add Branch UniqueObjectFinder/jets Jet Jet
755 add Branch UniqueObjectFinder/electrons Electron Electron
756 add Branch UniqueObjectFinder/photons Photon Photon
757 add Branch UniqueObjectFinder/muons Muon Muon
758 add Branch MissingET/momentum MissingET MissingET
759 add Branch ScalarHT/energy ScalarHT ScalarHT
760 add Branch Rho/rho Rho Rho
761 add Branch PileUpMerger/vertices Vertex Vertex
762}
Note: See TracBrowser for help on using the repository browser.