Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ 6e91f77

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 6e91f77 was 59d6164, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

added GenMissingET to all cards
  • Property mode set to 100644
File size: 21.7 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 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
314 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
315
316 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
317 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
318 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
319}
320
321#################
322# Electron filter
323#################
324
325module PdgCodeFilter ElectronFilter {
326 set InputArray Calorimeter/eflowTracks
327 set OutputArray electrons
328 set Invert true
329 add PdgCode {11}
330 add PdgCode {-11}
331}
332
333##########################
334# Track pile-up subtractor
335##########################
336
337module TrackPileUpSubtractor TrackPileUpSubtractor {
338# add InputArray InputArray OutputArray
339 add InputArray Calorimeter/eflowTracks eflowTracks
340 add InputArray ElectronFilter/electrons electrons
341 add InputArray MuonMomentumSmearing/muons muons
342
343 set VertexInputArray PileUpMerger/vertices
344 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
345 # Z vertex resolution in m
346 set ZVertexResolution 0.0001
347}
348
349####################
350# Neutral tower merger
351####################
352
353module Merger NeutralTowerMerger {
354# add InputArray InputArray
355 add InputArray Calorimeter/eflowPhotons
356 add InputArray Calorimeter/eflowNeutralHadrons
357 set OutputArray eflowTowers
358}
359
360##################################
361# Energy flow merger (all tracks)
362##################################
363
364module Merger EFlowMergerAllTracks {
365# add InputArray InputArray
366 add InputArray TrackMerger/tracks
367 add InputArray Calorimeter/eflowPhotons
368 add InputArray Calorimeter/eflowNeutralHadrons
369 set OutputArray eflow
370}
371
372
373####################
374# Energy flow merger
375####################
376
377module Merger EFlowMerger {
378# add InputArray InputArray
379 add InputArray TrackPileUpSubtractor/eflowTracks
380 add InputArray Calorimeter/eflowPhotons
381 add InputArray Calorimeter/eflowNeutralHadrons
382 set OutputArray eflow
383}
384
385
386
387#############
388# Rho pile-up
389#############
390
391module FastJetGridMedianEstimator Rho {
392
393 set InputArray EFlowMerger/eflow
394 set RhoOutputArray rho
395
396 # add GridRange rapmin rapmax drap dphi
397 # rapmin - the minimum rapidity extent of the grid
398 # rapmax - the maximum rapidity extent of the grid
399 # drap - the grid spacing in rapidity
400 # dphi - the grid spacing in azimuth
401
402 add GridRange -5.0 -2.5 1.0 1.0
403 add GridRange -2.5 2.5 1.0 1.0
404 add GridRange 2.5 5.0 1.0 1.0
405
406}
407
408#####################
409# Neutrino Filter
410#####################
411
412module PdgCodeFilter NeutrinoFilter {
413
414 set InputArray Delphes/stableParticles
415 set OutputArray filteredParticles
416
417 set PTMin 0.0
418
419 add PdgCode {12}
420 add PdgCode {14}
421 add PdgCode {16}
422 add PdgCode {-12}
423 add PdgCode {-14}
424 add PdgCode {-16}
425
426}
427
428
429
430#####################
431# MC truth jet finder
432#####################
433
434module FastJetFinder GenJetFinder {
435 set InputArray NeutrinoFilter/filteredParticles
436
437 set OutputArray jets
438
439 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
440 set JetAlgorithm 6
441 set ParameterR 0.5
442
443 set JetPTMin 20.0
444}
445
446#########################
447# Gen Missing ET merger
448########################
449
450module Merger GenMissingET {
451# add InputArray InputArray
452 add InputArray NeutrinoFilter/filteredParticles
453 set MomentumOutputArray momentum
454}
455
456
457
458############
459# Jet finder
460############
461
462module FastJetFinder FastJetFinder {
463# set InputArray Calorimeter/towers
464 set InputArray EFlowMerger/eflow
465
466 set OutputArray jets
467
468 # 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
469 set AreaAlgorithm 5
470
471 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
472 set JetAlgorithm 6
473 set ParameterR 0.5
474
475 set JetPTMin 20.0
476}
477
478###########################
479# Jet Pile-Up ID
480###########################
481
482module PileUpJetID PileUpJetID {
483 set JetInputArray FastJetFinder/jets
484 set TrackInputArray Calorimeter/eflowTracks
485 set NeutralInputArray NeutralTowerMerger/eflowTowers
486
487 set VertexInputArray PileUpMerger/vertices
488 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
489 # Z vertex resolution in m
490 set ZVertexResolution 0.0001
491
492 set OutputArray jets
493
494 set UseConstituents 0
495 set ParameterR 0.5
496
497 set JetPTMin 20.0
498}
499
500###########################
501# Jet Pile-Up Subtraction
502###########################
503
504module JetPileUpSubtractor JetPileUpSubtractor {
505 set JetInputArray PileUpJetID/jets
506 set RhoInputArray Rho/rho
507
508 set OutputArray jets
509
510 set JetPTMin 20.0
511}
512
513##################
514# Jet Energy Scale
515##################
516
517module EnergyScale JetEnergyScale {
518 set InputArray JetPileUpSubtractor/jets
519 set OutputArray jets
520
521 # scale formula for jets
522 set ScaleFormula {1.0}
523}
524
525###################
526# Photon efficiency
527###################
528
529module Efficiency PhotonEfficiency {
530 set InputArray Calorimeter/eflowPhotons
531 set OutputArray photons
532
533 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
534
535 # efficiency formula for photons
536 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
537 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
538 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
539 (abs(eta) > 2.5) * (0.00)}
540}
541
542
543##################
544# Photon isolation
545##################
546
547module Isolation PhotonIsolation {
548 set CandidateInputArray PhotonEfficiency/photons
549 set IsolationInputArray EFlowMergerAllTracks/eflow
550 set RhoInputArray Rho/rho
551
552 set OutputArray photons
553
554 set DeltaRMax 0.5
555
556 set PTMin 0.5
557
558 set PTRatioMax 0.12
559}
560
561#####################
562# Electron efficiency
563#####################
564
565module Efficiency ElectronEfficiency {
566 set InputArray TrackPileUpSubtractor/electrons
567 set OutputArray electrons
568
569 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
570
571 # efficiency formula for electrons
572 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
573 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
574 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
575 (abs(eta) > 2.5) * (0.00)}
576}
577
578####################
579# Electron isolation
580####################
581
582module Isolation ElectronIsolation {
583 set CandidateInputArray ElectronEfficiency/electrons
584 set IsolationInputArray EFlowMergerAllTracks/eflow
585 set RhoInputArray Rho/rho
586
587 set OutputArray electrons
588
589 set DeltaRMax 0.5
590
591 set PTMin 0.5
592
593 set PTRatioMax 0.12
594}
595
596#################
597# Muon efficiency
598#################
599
600module Efficiency MuonEfficiency {
601 set InputArray TrackPileUpSubtractor/muons
602 set OutputArray muons
603
604 # set EfficiencyFormula {efficiency as a function of eta and pt}
605
606 # efficiency formula for muons
607 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
608 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
609 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
610 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
611 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
612 (abs(eta) > 2.4) * (0.00)}
613}
614
615################
616# Muon isolation
617################
618
619module Isolation MuonIsolation {
620 set CandidateInputArray MuonEfficiency/muons
621 set IsolationInputArray EFlowMergerAllTracks/eflow
622 set RhoInputArray Rho/rho
623
624 set OutputArray muons
625
626 set DeltaRMax 0.5
627
628 set PTMin 0.5
629
630 set PTRatioMax 0.25
631}
632
633###################
634# Missing ET merger
635###################
636
637module Merger MissingET {
638# add InputArray InputArray
639 add InputArray EFlowMergerAllTracks/eflow
640 set MomentumOutputArray momentum
641}
642
643
644
645##################
646# Scalar HT merger
647##################
648
649module Merger ScalarHT {
650# add InputArray InputArray
651 add InputArray UniqueObjectFinder/jets
652 add InputArray UniqueObjectFinder/electrons
653 add InputArray UniqueObjectFinder/photons
654 add InputArray UniqueObjectFinder/muons
655 set EnergyOutputArray energy
656}
657
658########################
659# Jet Flavor Association
660########################
661
662module JetFlavorAssociation JetFlavorAssociation {
663
664 set PartonInputArray Delphes/partons
665 set ParticleInputArray Delphes/allParticles
666 set ParticleLHEFInputArray Delphes/allParticlesLHEF
667 set JetInputArray JetEnergyScale/jets
668
669 set DeltaR 0.5
670 set PartonPTMin 1.0
671 set PartonEtaMax 2.5
672
673}
674
675###########
676# b-tagging
677###########
678
679module BTagging BTagging {
680 set JetInputArray JetEnergyScale/jets
681
682 set BitNumber 0
683
684 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
685 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
686 # gluon's PDG code has the lowest priority
687
688 # based on arXiv:1211.4462
689
690 # default efficiency formula (misidentification rate)
691 add EfficiencyFormula {0} {0.01+0.000038*pt}
692
693 # efficiency formula for c-jets (misidentification rate)
694 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
695
696 # efficiency formula for b-jets
697 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
698}
699
700#############
701# tau-tagging
702#############
703
704module TauTagging TauTagging {
705 set ParticleInputArray Delphes/allParticles
706 set PartonInputArray Delphes/partons
707 set JetInputArray JetEnergyScale/jets
708
709 set DeltaR 0.5
710
711 set TauPTMin 1.0
712
713 set TauEtaMax 2.5
714
715 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
716
717 # default efficiency formula (misidentification rate)
718 add EfficiencyFormula {0} {0.01}
719 # efficiency formula for tau-jets
720 add EfficiencyFormula {15} {0.6}
721}
722
723#####################################################
724# Find uniquely identified photons/electrons/tau/jets
725#####################################################
726
727module UniqueObjectFinder UniqueObjectFinder {
728# earlier arrays take precedence over later ones
729# add InputArray InputArray OutputArray
730 add InputArray PhotonIsolation/photons photons
731 add InputArray ElectronIsolation/electrons electrons
732 add InputArray MuonIsolation/muons muons
733 add InputArray JetEnergyScale/jets jets
734}
735
736##################
737# ROOT tree writer
738##################
739
740# tracks, towers and eflow objects are not stored by default in the output.
741# if needed (for jet constituent or other studies), uncomment the relevant
742# "add Branch ..." lines.
743
744module TreeWriter TreeWriter {
745# add Branch InputArray BranchName BranchClass
746 add Branch Delphes/allParticles Particle GenParticle
747
748# add Branch TrackMerger/tracks Track Track
749# add Branch Calorimeter/towers Tower Tower
750
751# add Branch Calorimeter/eflowTracks EFlowTrack Track
752# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
753# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
754
755 add Branch GenJetFinder/jets GenJet Jet
756 add Branch GenMissingET/momentum GenMissingET MissingET
757
758 add Branch UniqueObjectFinder/jets Jet Jet
759 add Branch UniqueObjectFinder/electrons Electron Electron
760 add Branch UniqueObjectFinder/photons Photon Photon
761 add Branch UniqueObjectFinder/muons Muon Muon
762 add Branch MissingET/momentum MissingET MissingET
763 add Branch ScalarHT/energy ScalarHT ScalarHT
764 add Branch Rho/rho Rho Rho
765 add Branch PileUpMerger/vertices Vertex Vertex
766}
Note: See TracBrowser for help on using the repository browser.