Fork me on GitHub

source: git/examples/delphes_card_CMS_PileUp.tcl@ 2e6a81b

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 2e6a81b was c526439, checked in by pavel <pavel@…>, 11 years ago

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