Fork me on GitHub

source: git/examples/delphes_card_ATLAS_PileUp.tcl@ 273e668

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 273e668 was f14bd6a, checked in by Michele <michele.selvaggi@…>, 10 years ago

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