Fork me on GitHub

source: git/cards/delphes_card_ATLAS_PileUp.tcl@ d4b9697

ImprovedOutputFile Timing dual_readout llp
Last change on this file since d4b9697 was d4b9697, checked in by Pavel Demin <pavel.demin@…>, 10 years ago

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