Fork me on GitHub

source: git/examples/delphes_card_CMS_PileUp.tcl@ 24d005f

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 24d005f was 22dc7fd, checked in by mselvaggi <mselvaggi@…>, 11 years ago

timing implemented
  • Property mode set to 100644
File size: 21.3 KB
RevLine 
[d7d2da3]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 ElectronEnergySmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19 Calorimeter
20 TrackPileUpSubtractor
21 EFlowMerger
22
23 GenJetFinder
24
25 Rho
26 FastJetFinder
27 JetPileUpSubtractor
28
[3c9e6d2]29 JetEnergyScale
[d7d2da3]30
31 PhotonEfficiency
32 PhotonIsolation
33
34 ElectronEfficiency
35 ElectronIsolation
36
37 MuonEfficiency
38 MuonIsolation
39
40 MissingET
41
42 BTagging
43 TauTagging
44
45 UniqueObjectFinder
46
47 ScalarHT
48
49 TreeWriter
50}
51
52###############
53# PileUp Merger
54###############
55
56module PileUpMerger PileUpMerger {
57 set InputArray Delphes/stableParticles
58
[fb21bc8]59 set ParticleOutputArray stableParticles
60 set VertexOutputArray vertices
[d7d2da3]61
62 # pre-generated minbias input file
63 set PileUpFile MinBias.pileup
64
65 # average expected pile up
[22dc7fd]66 set MeanPileUp 10
67
68 # maximum spread in the beam direction in m
69 set ZVertexSpread 0.10
70
71 # maximum spread in time in s
72 set TVertexSpread 1.5E-09
73
74 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
75
76 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
77
78 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
79 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
80 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + \
81 # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
82
83
[d7d2da3]84}
85
86#################################
87# Propagate particles in cylinder
88#################################
89
90module ParticlePropagator ParticlePropagator {
91 set InputArray PileUpMerger/stableParticles
92
93 set OutputArray stableParticles
94 set ChargedHadronOutputArray chargedHadrons
95 set ElectronOutputArray electrons
96 set MuonOutputArray muons
97
98 # radius of the magnetic field coverage, in m
99 set Radius 1.29
100 # half-length of the magnetic field coverage, in m
101 set HalfLength 3.00
102
103 # magnetic field
104 set Bz 3.8
105}
106
107####################################
108# Charged hadron tracking efficiency
109####################################
110
111module Efficiency ChargedHadronTrackingEfficiency {
112 set InputArray ParticlePropagator/chargedHadrons
113 set OutputArray chargedHadrons
114
115 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
116
117 # tracking efficiency formula for charged hadrons
118 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
119 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
120 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
121 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
123 (abs(eta) > 2.5) * (0.00)}
124}
125
126##############################
127# Electron tracking efficiency
128##############################
129
130module Efficiency ElectronTrackingEfficiency {
131 set InputArray ParticlePropagator/electrons
132 set OutputArray electrons
133
134 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
135
136 # tracking efficiency formula for electrons
137 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
138 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
139 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
140 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
141 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
144 (abs(eta) > 2.5) * (0.00)}
145}
146
147##########################
148# Muon tracking efficiency
149##########################
150
151module Efficiency MuonTrackingEfficiency {
152 set InputArray ParticlePropagator/muons
153 set OutputArray muons
154
155 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
156
157 # tracking efficiency formula for muons
158 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
159 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
160 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
161 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
162 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
163 (abs(eta) > 2.5) * (0.00)}
164}
165
166########################################
167# Momentum resolution for charged tracks
168########################################
169
170module MomentumSmearing ChargedHadronMomentumSmearing {
171 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
172 set OutputArray chargedHadrons
173
174 # set ResolutionFormula {resolution formula as a function of eta and pt}
175
176 # resolution formula for charged hadrons
[a1ac20e]177 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
178 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
179 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
180 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
183 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
184 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
[d7d2da3]185}
186
187#################################
188# Energy resolution for electrons
189#################################
190
191module EnergySmearing ElectronEnergySmearing {
192 set InputArray ElectronTrackingEfficiency/electrons
193 set OutputArray electrons
194
195 # set ResolutionFormula {resolution formula as a function of eta and energy}
196
[3c9e6d2]197 # resolution formula for electrons
[a0ff736]198 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \
[3c9e6d2]199 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
200 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
[d7d2da3]201 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
202
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
[3c9e6d2]216 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
217 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
[1697699]218 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
219 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
[3c9e6d2]220 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
221 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
[1697699]222 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
223 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
[3c9e6d2]224 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
225 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
[1697699]226 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
227 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
[d7d2da3]228}
229
230##############
231# Track merger
232##############
233
234module Merger TrackMerger {
235# add InputArray InputArray
236 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
237 add InputArray ElectronEnergySmearing/electrons
[a0ff736]238 add InputArray MuonMomentumSmearing/muons
[d7d2da3]239 set OutputArray tracks
240}
241
242#############
243# Calorimeter
244#############
245
246module Calorimeter Calorimeter {
247 set ParticleInputArray ParticlePropagator/stableParticles
248 set TrackInputArray TrackMerger/tracks
249
250 set TowerOutputArray towers
251 set PhotonOutputArray photons
252
253 set EFlowTrackOutputArray eflowTracks
254 set EFlowTowerOutputArray eflowTowers
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}
[3c9e6d2]310 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
[d7d2da3]311 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
312
313 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
[a1ac20e]314 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
315 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
[d7d2da3]316}
317
318##########################
319# Track pile-up subtractor
320##########################
321
322module TrackPileUpSubtractor TrackPileUpSubtractor {
323# add InputArray InputArray OutputArray
324 add InputArray Calorimeter/eflowTracks eflowTracks
325 add InputArray ElectronEnergySmearing/electrons electrons
326 add InputArray MuonMomentumSmearing/muons muons
327
328 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
329 # Z vertex resolution in m
330 set ZVertexResolution 0.0001
331}
332
333####################
334# Energy flow merger
335####################
336
337module Merger EFlowMerger {
338# add InputArray InputArray
339 add InputArray TrackPileUpSubtractor/eflowTracks
340 add InputArray Calorimeter/eflowTowers
341 set OutputArray eflow
342}
343
344#############
345# Rho pile-up
346#############
347
348module FastJetFinder Rho {
349# set InputArray Calorimeter/towers
350 set InputArray EFlowMerger/eflow
351
352 set ComputeRho true
353 set RhoOutputArray rho
354
355 # 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
356 set AreaAlgorithm 5
357
358 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
359 set JetAlgorithm 4
360 set ParameterR 0.6
361 set GhostEtaMax 5.0
[e9065e7]362
363 add RhoEtaRange 0.0 2.5
364 add RhoEtaRange 2.5 5.0
[d7d2da3]365
366 set JetPTMin 0.0
367}
368
369#####################
370# MC truth jet finder
371#####################
372
373module FastJetFinder GenJetFinder {
374 set InputArray Delphes/stableParticles
375
376 set OutputArray jets
377
378 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
379 set JetAlgorithm 6
380 set ParameterR 0.5
381
382 set JetPTMin 20.0
383}
384
385############
386# Jet finder
387############
388
389module FastJetFinder FastJetFinder {
390# set InputArray Calorimeter/towers
391 set InputArray EFlowMerger/eflow
392
393 set OutputArray jets
394
395 # 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
396 set AreaAlgorithm 5
397
398 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
399 set JetAlgorithm 6
400 set ParameterR 0.5
401
402 set JetPTMin 20.0
403}
404
405###########################
406# Jet Pile-Up Subtraction
407###########################
408
409module JetPileUpSubtractor JetPileUpSubtractor {
410 set JetInputArray FastJetFinder/jets
411 set RhoInputArray Rho/rho
412
413 set OutputArray jets
414
415 set JetPTMin 20.0
416}
417
[3c9e6d2]418##################
419# Jet Energy Scale
420##################
421
422module EnergyScale JetEnergyScale {
423 set InputArray JetPileUpSubtractor/jets
424 set OutputArray jets
425
426 # scale formula for jets
427 set ScaleFormula {1.0}
428}
429
[d7d2da3]430###################
431# Photon efficiency
432###################
433
434module Efficiency PhotonEfficiency {
435 set InputArray Calorimeter/photons
436 set OutputArray photons
437
438 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
439
440 # efficiency formula for photons
441 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
442 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
443 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
444 (abs(eta) > 2.5) * (0.00)}
445}
446
447##################
448# Photon isolation
449##################
450
451module Isolation PhotonIsolation {
452 set CandidateInputArray PhotonEfficiency/photons
453 set IsolationInputArray EFlowMerger/eflow
454 set RhoInputArray Rho/rho
455
456 set OutputArray photons
457
458 set DeltaRMax 0.5
459
460 set PTMin 0.5
461
462 set PTRatioMax 0.1
463}
464
465#####################
466# Electron efficiency
467#####################
468
469module Efficiency ElectronEfficiency {
470 set InputArray TrackPileUpSubtractor/electrons
471 set OutputArray electrons
472
473 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
474
475 # efficiency formula for electrons
476 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
477 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
478 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
479 (abs(eta) > 2.5) * (0.00)}
480}
481
482####################
483# Electron isolation
484####################
485
486module Isolation ElectronIsolation {
487 set CandidateInputArray ElectronEfficiency/electrons
488 set IsolationInputArray EFlowMerger/eflow
489 set RhoInputArray Rho/rho
490
491 set OutputArray electrons
492
493 set DeltaRMax 0.5
494
495 set PTMin 0.5
496
497 set PTRatioMax 0.1
498}
499
500#################
501# Muon efficiency
502#################
503
504module Efficiency MuonEfficiency {
505 set InputArray TrackPileUpSubtractor/muons
506 set OutputArray muons
507
508 # set EfficiencyFormula {efficiency as a function of eta and pt}
509
510 # efficiency formula for muons
[1697699]511 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
512 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
513 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
514 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
515 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
516 (abs(eta) > 2.4) * (0.00)}
[d7d2da3]517}
518
519################
520# Muon isolation
521################
522
523module Isolation MuonIsolation {
524 set CandidateInputArray MuonEfficiency/muons
525 set IsolationInputArray EFlowMerger/eflow
526 set RhoInputArray Rho/rho
527
528 set OutputArray muons
529
530 set DeltaRMax 0.5
531
532 set PTMin 0.5
533
534 set PTRatioMax 0.1
535}
536
537###################
538# Missing ET merger
539###################
540
541module Merger MissingET {
542# add InputArray InputArray
[a0ff736]543 add InputArray Calorimeter/eflowTracks
544 add InputArray Calorimeter/eflowTowers
[d7d2da3]545 set MomentumOutputArray momentum
546}
547
548##################
549# Scalar HT merger
550##################
551
552module Merger ScalarHT {
553# add InputArray InputArray
554 add InputArray UniqueObjectFinder/jets
555 add InputArray UniqueObjectFinder/electrons
556 add InputArray UniqueObjectFinder/photons
[3c9e6d2]557 add InputArray UniqueObjectFinder/muons
[d7d2da3]558 set EnergyOutputArray energy
559}
560
561###########
562# b-tagging
563###########
564
565module BTagging BTagging {
566 set PartonInputArray Delphes/partons
[3c9e6d2]567 set JetInputArray JetEnergyScale/jets
[d7d2da3]568
[264bf40]569 set BitNumber 0
570
[d7d2da3]571 set DeltaR 0.5
572
573 set PartonPTMin 1.0
574
575 set PartonEtaMax 2.5
576
577 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
578 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
579 # gluon's PDG code has the lowest priority
580
581 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
582 # default efficiency formula (misidentification rate)
583 add EfficiencyFormula {0} {0.001}
584
585 # efficiency formula for c-jets (misidentification rate)
586 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
587 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
588 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
589 (abs(eta) > 2.5) * (0.000)}
590
591 # efficiency formula for b-jets
592 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
593 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
594 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
595 (abs(eta) > 2.5) * (0.000)}
596}
597
598module TauTagging TauTagging {
599 set ParticleInputArray Delphes/allParticles
600 set PartonInputArray Delphes/partons
[3c9e6d2]601 set JetInputArray JetEnergyScale/jets
[d7d2da3]602
603 set DeltaR 0.5
604
605 set TauPTMin 1.0
606
607 set TauEtaMax 2.5
608
609 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
610
611 # default efficiency formula (misidentification rate)
612 add EfficiencyFormula {0} {0.001}
613 # efficiency formula for tau-jets
614 add EfficiencyFormula {15} {0.4}
615}
616
617#####################################################
618# Find uniquely identified photons/electrons/tau/jets
619#####################################################
620
621module UniqueObjectFinder UniqueObjectFinder {
622# earlier arrays take precedence over later ones
623# add InputArray InputArray OutputArray
624 add InputArray PhotonIsolation/photons photons
625 add InputArray ElectronIsolation/electrons electrons
[3c9e6d2]626 add InputArray MuonIsolation/muons muons
627 add InputArray JetEnergyScale/jets jets
[d7d2da3]628}
629
630##################
631# ROOT tree writer
632##################
633
634module TreeWriter TreeWriter {
635# add Branch InputArray BranchName BranchClass
636 add Branch Delphes/allParticles Particle GenParticle
637 add Branch TrackMerger/tracks Track Track
638 add Branch Calorimeter/towers Tower Tower
[3bb5979]639 add Branch Calorimeter/eflowTracks EFlowTrack Track
640 add Branch Calorimeter/eflowTowers EFlowTower Tower
[d7d2da3]641 add Branch GenJetFinder/jets GenJet Jet
642 add Branch UniqueObjectFinder/jets Jet Jet
643 add Branch UniqueObjectFinder/electrons Electron Electron
644 add Branch UniqueObjectFinder/photons Photon Photon
[3c9e6d2]645 add Branch UniqueObjectFinder/muons Muon Muon
[d7d2da3]646 add Branch MissingET/momentum MissingET MissingET
647 add Branch ScalarHT/energy ScalarHT ScalarHT
[71648c2]648 add Branch Rho/rho Rho Rho
[fb21bc8]649 add Branch PileUpMerger/vertices Vertex Vertex
[d7d2da3]650}
651
Note: See TracBrowser for help on using the repository browser.