Fork me on GitHub

source: git/examples/delphes_card_ATLAS_PileUp.tcl@ 6fb1a5d

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

implemented FastJetGridMedianEstimator module, added 1k events MinBias file, put back pileup tracks in MET calculation, updated CMS and ATLAS pile-up cards.
  • Property mode set to 100644
File size: 21.6 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
[be2222c]24
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
[1c8d9db]251 set EcalTowerMinEnergy 0.50
252 set HcalTowerMinEnergy 1.00
253
254 set EcalTowerMinSignificance 1.0
255 set HcalTowerMinSignificance 1.0
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 {
376
[89c5be4]377 set InputArray Calorimeter/towers
378 set RhoOutputArray rho
379
[6fb1a5d]380 # etamin etamax gridsize_eta gridsize_phi
381
382 add GridRange 0.0 2.5 0.5 0.5
383 add GridRange 2.5 5.0 0.5 0.5
[984cd31]384
[89c5be4]385}
386
[be2222c]387
388#####################
389# Neutrino Filter
390#####################
391
392module PdgCodeFilter NeutrinoFilter {
393
394 set InputArray Delphes/stableParticles
395 set OutputArray filteredParticles
396
397 set PTMin 0.0
398
399 add PdgCode {12}
400 add PdgCode {14}
401 add PdgCode {16}
402 add PdgCode {-12}
403 add PdgCode {-14}
404 add PdgCode {-16}
405
406}
407
[89c5be4]408#####################
409# MC truth jet finder
410#####################
411
412module FastJetFinder GenJetFinder {
[be2222c]413 set InputArray NeutrinoFilter/filteredParticles
[89c5be4]414
415 set OutputArray jets
416
417 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
418 set JetAlgorithm 6
419 set ParameterR 0.6
420
421 set JetPTMin 20.0
422}
423
424############
425# Jet finder
426############
427
428module FastJetFinder FastJetFinder {
429 set InputArray Calorimeter/towers
430
431 set OutputArray jets
432
433 # 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
434 set AreaAlgorithm 5
435
436 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
437 set JetAlgorithm 6
438 set ParameterR 0.6
439
440 set JetPTMin 20.0
441}
442
443###########################
444# Jet Pile-Up Subtraction
445###########################
446
447module JetPileUpSubtractor JetPileUpSubtractor {
448 set JetInputArray FastJetFinder/jets
449 set RhoInputArray Rho/rho
450
451 set OutputArray jets
452
453 set JetPTMin 20.0
454}
455
456##################
457# Jet Energy Scale
458##################
459
460module EnergyScale JetEnergyScale {
461 set InputArray JetPileUpSubtractor/jets
462 set OutputArray jets
463
464 # scale formula for jets
465 set ScaleFormula {1.0}
466}
467
468###################
469# Photon efficiency
470###################
471
472module Efficiency PhotonEfficiency {
473 set InputArray Calorimeter/photons
474 set OutputArray photons
475
476 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
477
478 # efficiency formula for photons
479 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
480 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
481 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
482 (abs(eta) > 2.5) * (0.00)}
483}
484
485##################
486# Photon isolation
487##################
488
489module Isolation PhotonIsolation {
490 set CandidateInputArray PhotonEfficiency/photons
491 set IsolationInputArray EFlowMerger/eflow
492 set RhoInputArray Rho/rho
493
494 set OutputArray photons
495
496 set DeltaRMax 0.5
497
498 set PTMin 0.5
499
500 set PTRatioMax 0.1
501}
502
503#####################
504# Electron efficiency
505#####################
506
507module Efficiency ElectronEfficiency {
508 set InputArray TrackPileUpSubtractor/electrons
509 set OutputArray electrons
510
511 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
512
513 # efficiency formula for electrons
514 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
515 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
516 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
517 (abs(eta) > 2.5) * (0.00)}
518}
519
520####################
521# Electron isolation
522####################
523
524module Isolation ElectronIsolation {
525 set CandidateInputArray ElectronEfficiency/electrons
526 set IsolationInputArray EFlowMerger/eflow
527 set RhoInputArray Rho/rho
528
529 set OutputArray electrons
530
531 set DeltaRMax 0.5
532
533 set PTMin 0.5
534
535 set PTRatioMax 0.1
536}
537
538#################
539# Muon efficiency
540#################
541
542module Efficiency MuonEfficiency {
543 set InputArray TrackPileUpSubtractor/muons
544 set OutputArray muons
545
546 # set EfficiencyFormula {efficiency as a function of eta and pt}
547
548 # efficiency formula for muons
549 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
550 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
551 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
552 (abs(eta) > 2.7) * (0.00)}
553}
554
555################
556# Muon isolation
557################
558
559module Isolation MuonIsolation {
560 set CandidateInputArray MuonEfficiency/muons
561 set IsolationInputArray EFlowMerger/eflow
562 set RhoInputArray Rho/rho
563
564 set OutputArray muons
565
566 set DeltaRMax 0.5
567
568 set PTMin 0.5
569
570 set PTRatioMax 0.1
571}
572
573###################
574# Missing ET merger
575###################
576
577module Merger MissingET {
578# add InputArray InputArray
[6fb1a5d]579 add InputArray EFlowMergerAllTracks/eflow
[89c5be4]580 set MomentumOutputArray momentum
581}
582
[27bf162]583
[89c5be4]584##################
585# Scalar HT merger
586##################
587
588module Merger ScalarHT {
589# add InputArray InputArray
590 add InputArray UniqueObjectFinder/jets
591 add InputArray UniqueObjectFinder/electrons
592 add InputArray UniqueObjectFinder/photons
593 add InputArray UniqueObjectFinder/muons
594 set EnergyOutputArray energy
595}
596
597###########
598# b-tagging
599###########
600
601module BTagging BTagging {
602 set PartonInputArray Delphes/partons
603 set JetInputArray JetEnergyScale/jets
604
605 set BitNumber 0
606
607 set DeltaR 0.5
608
609 set PartonPTMin 1.0
610
611 set PartonEtaMax 2.5
612
613 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
614 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
615 # gluon's PDG code has the lowest priority
616
617 # default efficiency formula (misidentification rate)
618 add EfficiencyFormula {0} {0.001}
619
620 # efficiency formula for c-jets (misidentification rate)
621 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
622 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
623 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
624 (abs(eta) > 2.5) * (0.000)}
625
626 # efficiency formula for b-jets
627 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
628 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
629 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
630 (abs(eta) > 2.5) * (0.000)}
631}
632
633module TauTagging TauTagging {
634 set ParticleInputArray Delphes/allParticles
635 set PartonInputArray Delphes/partons
636 set JetInputArray JetEnergyScale/jets
637
638 set DeltaR 0.5
639
640 set TauPTMin 1.0
641
642 set TauEtaMax 2.5
643
644 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
645
646 # default efficiency formula (misidentification rate)
647 add EfficiencyFormula {0} {0.001}
648 # efficiency formula for tau-jets
649 add EfficiencyFormula {15} {0.4}
650}
651
652#####################################################
653# Find uniquely identified photons/electrons/tau/jets
654#####################################################
655
656module UniqueObjectFinder UniqueObjectFinder {
657# earlier arrays take precedence over later ones
658# add InputArray InputArray OutputArray
659 add InputArray PhotonIsolation/photons photons
660 add InputArray ElectronIsolation/electrons electrons
661 add InputArray MuonIsolation/muons muons
662 add InputArray JetEnergyScale/jets jets
663}
664
665##################
666# ROOT tree writer
667##################
668
[27bf162]669# tracks, towers and eflow objects are not stored by default in the output.
[282f591]670# if needed (for jet constituent or other studies), uncomment the relevant
671# "add Branch ..." lines.
[27bf162]672
[89c5be4]673module TreeWriter TreeWriter {
674# add Branch InputArray BranchName BranchClass
675 add Branch Delphes/allParticles Particle GenParticle
[27bf162]676
677# add Branch TrackMerger/tracks Track Track
[1c8d9db]678 add Branch Calorimeter/towers Tower Tower
[282f591]679
[89c5be4]680# add Branch Calorimeter/eflowTracks EFlowTrack Track
[27bf162]681# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
682# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
683
[89c5be4]684 add Branch GenJetFinder/jets GenJet Jet
685 add Branch UniqueObjectFinder/jets Jet Jet
686 add Branch UniqueObjectFinder/electrons Electron Electron
687 add Branch UniqueObjectFinder/photons Photon Photon
688 add Branch UniqueObjectFinder/muons Muon Muon
689 add Branch MissingET/momentum MissingET MissingET
690 add Branch ScalarHT/energy ScalarHT ScalarHT
691 add Branch Rho/rho Rho Rho
[fb21bc8]692 add Branch PileUpMerger/vertices Vertex Vertex
[27bf162]693
[89c5be4]694}
695
[27bf162]696
697
Note: See TracBrowser for help on using the repository browser.