Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ fdfac34

ImprovedOutputFile Timing dual_readout llp
Last change on this file since fdfac34 was a356d6d, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

corrected light mist rate in btag formulae

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