Fork me on GitHub

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

Last change on this file since 1322 was 1320, checked in by Pavel Demin, 11 years ago

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