Fork me on GitHub

source: svn/trunk/examples/delphes_card_CMS_PileUp.tcl@ 1349

Last change on this file since 1349 was 1349, checked in by Michele Selvaggi, 11 years ago

JetId v2
File size: 21.9 KB
RevLine 
[1009]1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
[1031]6
[1009]7 PileUpMerger
8 ParticlePropagator
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronEnergySmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19 Calorimeter
[1072]20 TrackPileUpSubtractor
[1009]21 EFlowMerger
22
[1074]23 GenJetFinder
24
[1031]25 Rho
[1009]26 FastJetFinder
[1349]27 PileUpJetID
[1009]28 JetPileUpSubtractor
[1031]29
[1256]30 JetEnergyScale
[1074]31
[1009]32 PhotonEfficiency
33 PhotonIsolation
[1031]34
[1009]35 ElectronEfficiency
36 ElectronIsolation
[1031]37
[1009]38 MuonEfficiency
39 MuonIsolation
[1031]40
[1009]41 MissingET
[1031]42
[1009]43 BTagging
44 TauTagging
[1031]45
[1009]46 UniqueObjectFinder
[1031]47
[1009]48 ScalarHT
[1031]49
[1009]50 TreeWriter
51}
52
[1016]53###############
[1009]54# PileUp Merger
[1016]55###############
[1009]56
57module PileUpMerger PileUpMerger {
58 set InputArray Delphes/stableParticles
[1016]59
[1326]60 set ParticleOutputArray stableParticles
61 set VertexOutputArray vertices
[1031]62
[1016]63 # pre-generated minbias input file
[1349]64 set PileUpFile ../../Delphes/MinBias.pileup
[1009]65
[1031]66 # average expected pile up
[1345]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
[1009]85}
[1026]86
[1009]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 {
[1072]113 set InputArray ParticlePropagator/chargedHadrons
[1009]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 {
[1072]132 set InputArray ParticlePropagator/electrons
[1009]133 set OutputArray electrons
134
135 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
[1031]136
[1009]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
[1031]148##########################
[1009]149# Muon tracking efficiency
[1031]150##########################
[1009]151
152module Efficiency MuonTrackingEfficiency {
[1072]153 set InputArray ParticlePropagator/muons
[1009]154 set OutputArray muons
155
156 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
[1031]157
[1009]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
[1278]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)}
[1009]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
[1256]198 # resolution formula for electrons
[1301]199 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \
[1256]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) + \
[1009]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
[1256]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) + \
[1320]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) + \
[1256]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) + \
[1320]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) + \
[1256]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) + \
[1320]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)}
[1009]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
[1301]239 add InputArray MuonMomentumSmearing/muons
[1009]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)}]
[1031]258
[1009]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 {}
[1031]265 for {set i -36} {$i <= 36} {incr i} {
[1009]266 add PhiBins [expr {$i * $pi/36.0}]
267 }
[1031]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} {
[1009]269 add EtaPhiBins $eta $PhiBins
270 }
271
272 # 10 degrees towers
273 set PhiBins {}
[1031]274 for {set i -18} {$i <= 18} {incr i} {
[1009]275 add PhiBins [expr {$i * $pi/18.0}]
276 }
[1031]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} {
[1009]278 add EtaPhiBins $eta $PhiBins
279 }
[1031]280
[1009]281 # 20 degrees towers
282 set PhiBins {}
[1031]283 for {set i -9} {$i <= 9} {incr i} {
[1009]284 add PhiBins [expr {$i * $pi/9.0}]
[1031]285 }
286 foreach eta {-5 -4.7 -4.525 4.7 5} {
[1009]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}
[1076]296 # energy fractions for muon, neutrinos and neutralinos
[1009]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}
[1076]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}
[1009]306 # energy fractions for K0short and Lambda
307 add EnergyFraction {310} {0.3 0.7}
308 add EnergyFraction {3122} {0.3 0.7}
[1031]309
[1009]310 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
[1256]311 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
[1009]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}
[1278]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)}
[1009]317}
318
[1072]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
[1349]328
329 set VertexInputArray PileUpMerger/vertices
[1072]330 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
331 # Z vertex resolution in m
332 set ZVertexResolution 0.0001
333}
334
[1009]335####################
336# Energy flow merger
337####################
338
339module Merger EFlowMerger {
340# add InputArray InputArray
[1072]341 add InputArray TrackPileUpSubtractor/eflowTracks
[1009]342 add InputArray Calorimeter/eflowTowers
343 set OutputArray eflow
344}
345
[1031]346#############
[1009]347# Rho pile-up
[1031]348#############
[1009]349
350module FastJetFinder Rho {
351# set InputArray Calorimeter/towers
352 set InputArray EFlowMerger/eflow
[1031]353
[1009]354 set ComputeRho true
[1031]355 set RhoOutputArray rho
[1009]356
[1023]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
[1009]358 set AreaAlgorithm 5
[1023]359
[1009]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
[1316]364
365 add RhoEtaRange 0.0 2.5
366 add RhoEtaRange 2.5 5.0
[1009]367
368 set JetPTMin 0.0
369}
370
[1074]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
[1009]387############
388# Jet finder
389############
390
391module FastJetFinder FastJetFinder {
392# set InputArray Calorimeter/towers
393 set InputArray EFlowMerger/eflow
394
395 set OutputArray jets
[1071]396
[1036]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
[1071]399
[1009]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
407###########################
[1349]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
429###########################
[1009]430# Jet Pile-Up Subtraction
431###########################
432
433module JetPileUpSubtractor JetPileUpSubtractor {
[1349]434 set JetInputArray PileUpJetID/jets
[1009]435 set RhoInputArray Rho/rho
[1031]436
[1009]437 set OutputArray jets
[1031]438
[1009]439 set JetPTMin 20.0
440}
441
[1256]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
[1009]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
[1023]478 set RhoInputArray Rho/rho
[1031]479
[1009]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 {
[1072]494 set InputArray TrackPileUpSubtractor/electrons
[1009]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
[1031]514
[1009]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 {
[1072]529 set InputArray TrackPileUpSubtractor/muons
[1009]530 set OutputArray muons
531
532 # set EfficiencyFormula {efficiency as a function of eta and pt}
533
534 # efficiency formula for muons
[1320]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)}
[1009]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
[1031]551
[1009]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
[1301]567 add InputArray Calorimeter/eflowTracks
568 add InputArray Calorimeter/eflowTowers
[1009]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
[1256]581 add InputArray UniqueObjectFinder/muons
[1009]582 set EnergyOutputArray energy
583}
584
585###########
586# b-tagging
587###########
588
589module BTagging BTagging {
590 set PartonInputArray Delphes/partons
[1256]591 set JetInputArray JetEnergyScale/jets
[1009]592
[1099]593 set BitNumber 0
594
[1009]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)) + \
[1031]613 (abs(eta) > 2.5) * (0.000)}
614
[1009]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
[1256]625 set JetInputArray JetEnergyScale/jets
[1009]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
[1256]650 add InputArray MuonIsolation/muons muons
651 add InputArray JetEnergyScale/jets jets
[1009]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
[1257]663 add Branch Calorimeter/eflowTracks EFlowTrack Track
664 add Branch Calorimeter/eflowTowers EFlowTower Tower
[1074]665 add Branch GenJetFinder/jets GenJet Jet
[1009]666 add Branch UniqueObjectFinder/jets Jet Jet
667 add Branch UniqueObjectFinder/electrons Electron Electron
668 add Branch UniqueObjectFinder/photons Photon Photon
[1256]669 add Branch UniqueObjectFinder/muons Muon Muon
[1009]670 add Branch MissingET/momentum MissingET MissingET
671 add Branch ScalarHT/energy ScalarHT ScalarHT
[1114]672 add Branch Rho/rho Rho Rho
[1326]673 add Branch PileUpMerger/vertices Vertex Vertex
[1009]674}
675
Note: See TracBrowser for help on using the repository browser.