Fork me on GitHub

source: git/examples/delphes_card_CMS_PileUp.tcl@ 83eff6d

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 83eff6d was 54b6dfc, checked in by mselvaggi <mselvaggi@…>, 11 years ago

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