Fork me on GitHub

source: git/cards/delphes_card_ATLAS_PileUp.tcl@ c627b07

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

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