Fork me on GitHub

source: git/cards/delphes_card_ATLAS_PileUp.tcl@ 8707eeb

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 8707eeb was 7b50bf8, checked in by Michele Selvaggi <michele.selvaggi@…>, 10 years ago

Changed grid in the fwd region for pile-up subtraction.
  • Property mode set to 100644
File size: 21.7 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
[a2983ec]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) +
[22dc7fd]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
[a2983ec]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) +
[89c5be4]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
[a2983ec]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) +
[89c5be4]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
[a2983ec]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) +
[89c5be4]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
[a2983ec]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) +
[ce1cbf1]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
[a2983ec]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) +
[89c5be4]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
[a2983ec]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) +
[89c5be4]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
[a2983ec]251 set EFlowTrackOutputArray eflowTracks
252 set EFlowPhotonOutputArray eflowPhotons
253 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
254
255 set ECalEnergyMin 0.5
256 set HCalEnergyMin 1.0
[1c8d9db]257
[38bf1ae]258 set ECalEnergySignificanceMin 1.0
259 set HCalEnergySignificanceMin 1.0
[1c8d9db]260
[a2983ec]261 set SmearTowerCenter true
[282f591]262
[89c5be4]263 set pi [expr {acos(-1)}]
264
265 # lists of the edges of each tower in eta and phi
266 # each list starts with the lower edge of the first tower
267 # the list ends with the higher edged of the last tower
268
269 # 10 degrees towers
270 set PhiBins {}
271 for {set i -18} {$i <= 18} {incr i} {
272 add PhiBins [expr {$i * $pi/18.0}]
273 }
274 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} {
275 add EtaPhiBins $eta $PhiBins
276 }
277
278 # 20 degrees towers
279 set PhiBins {}
280 for {set i -9} {$i <= 9} {incr i} {
281 add PhiBins [expr {$i * $pi/9.0}]
282 }
283 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} {
284 add EtaPhiBins $eta $PhiBins
285 }
286
287 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
288 add EnergyFraction {0} {0.0 1.0}
289 # energy fractions for e, gamma and pi0
290 add EnergyFraction {11} {1.0 0.0}
291 add EnergyFraction {22} {1.0 0.0}
292 add EnergyFraction {111} {1.0 0.0}
293 # energy fractions for muon, neutrinos and neutralinos
294 add EnergyFraction {12} {0.0 0.0}
295 add EnergyFraction {13} {0.0 0.0}
296 add EnergyFraction {14} {0.0 0.0}
297 add EnergyFraction {16} {0.0 0.0}
298 add EnergyFraction {1000022} {0.0 0.0}
299 add EnergyFraction {1000023} {0.0 0.0}
300 add EnergyFraction {1000025} {0.0 0.0}
301 add EnergyFraction {1000035} {0.0 0.0}
302 add EnergyFraction {1000045} {0.0 0.0}
303 # energy fractions for K0short and Lambda
304 add EnergyFraction {310} {0.3 0.7}
305 add EnergyFraction {3122} {0.3 0.7}
306
307 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
308 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
309 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
310 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
[a2983ec]311 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
[89c5be4]312 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
313
314 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
315 # http://arxiv.org/pdf/hep-ex/0004009v1
316 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
[a2983ec]317 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
318 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
[f42fdd9]319 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
[89c5be4]320}
321
322##########################
323# Track pile-up subtractor
324##########################
325
326module TrackPileUpSubtractor TrackPileUpSubtractor {
327# add InputArray InputArray OutputArray
328 add InputArray Calorimeter/eflowTracks eflowTracks
329 add InputArray ElectronEnergySmearing/electrons electrons
330 add InputArray MuonMomentumSmearing/muons muons
331
332 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
333 # Z vertex resolution in m
334 set ZVertexResolution 0.0001
335}
336
337####################
[27bf162]338# Neutral tower merger
[89c5be4]339####################
340
[27bf162]341module Merger NeutralTowerMerger {
342# add InputArray InputArray
343 add InputArray Calorimeter/eflowPhotons
344 add InputArray Calorimeter/eflowNeutralHadrons
345 set OutputArray eflowTowers
[282f591]346}
[27bf162]347
[6fb1a5d]348##################################
349# Energy flow merger (all tracks)
350##################################
351
352module Merger EFlowMergerAllTracks {
353# add InputArray InputArray
354 add InputArray TrackMerger/tracks
355 add InputArray Calorimeter/eflowPhotons
356 add InputArray Calorimeter/eflowNeutralHadrons
357 set OutputArray eflow
358}
359
[27bf162]360
361####################
362# Energy flow merger
363####################
[282f591]364
[89c5be4]365module Merger EFlowMerger {
366# add InputArray InputArray
[27bf162]367 add InputArray Calorimeter/eflowTracks
368 add InputArray Calorimeter/eflowPhotons
369 add InputArray Calorimeter/eflowNeutralHadrons
[89c5be4]370 set OutputArray eflow
[282f591]371}
[27bf162]372
[89c5be4]373#############
374# Rho pile-up
375#############
376
[6fb1a5d]377module FastJetGridMedianEstimator Rho {
[d4b9697]378
[89c5be4]379 set InputArray Calorimeter/towers
380 set RhoOutputArray rho
381
[d4b9697]382 # add GridRange rapmin rapmax drap dphi
383 # rapmin - the minimum rapidity extent of the grid
384 # rapmax - the maximum rapidity extent of the grid
385 # drap - the grid spacing in rapidity
386 # dphi - the grid spacing in azimuth
387
[7b50bf8]388 add GridRange -5.0 -2.5 1.0 1.0
[914a1f5]389 add GridRange -2.5 2.5 0.5 0.5
[7b50bf8]390 add GridRange 2.5 5.0 1.0 1.0
[984cd31]391
[89c5be4]392}
393
[be2222c]394
395#####################
396# Neutrino Filter
397#####################
398
399module PdgCodeFilter NeutrinoFilter {
[d4b9697]400
[be2222c]401 set InputArray Delphes/stableParticles
402 set OutputArray filteredParticles
403
404 set PTMin 0.0
[d4b9697]405
[be2222c]406 add PdgCode {12}
407 add PdgCode {14}
408 add PdgCode {16}
409 add PdgCode {-12}
410 add PdgCode {-14}
411 add PdgCode {-16}
412
413}
414
[89c5be4]415#####################
416# MC truth jet finder
417#####################
418
419module FastJetFinder GenJetFinder {
[be2222c]420 set InputArray NeutrinoFilter/filteredParticles
[89c5be4]421
422 set OutputArray jets
423
424 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
425 set JetAlgorithm 6
426 set ParameterR 0.6
427
428 set JetPTMin 20.0
429}
430
431############
432# Jet finder
433############
434
435module FastJetFinder FastJetFinder {
436 set InputArray Calorimeter/towers
437
438 set OutputArray jets
439
440 # 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
441 set AreaAlgorithm 5
442
443 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
444 set JetAlgorithm 6
445 set ParameterR 0.6
446
447 set JetPTMin 20.0
448}
449
450###########################
451# Jet Pile-Up Subtraction
452###########################
453
454module JetPileUpSubtractor JetPileUpSubtractor {
455 set JetInputArray FastJetFinder/jets
456 set RhoInputArray Rho/rho
457
458 set OutputArray jets
459
460 set JetPTMin 20.0
461}
462
463##################
464# Jet Energy Scale
465##################
466
467module EnergyScale JetEnergyScale {
468 set InputArray JetPileUpSubtractor/jets
469 set OutputArray jets
470
471 # scale formula for jets
472 set ScaleFormula {1.0}
473}
474
475###################
476# Photon efficiency
477###################
478
479module Efficiency PhotonEfficiency {
480 set InputArray Calorimeter/photons
481 set OutputArray photons
482
483 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
484
485 # efficiency formula for photons
[a2983ec]486 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
487 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
488 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
[89c5be4]489 (abs(eta) > 2.5) * (0.00)}
490}
491
492##################
493# Photon isolation
494##################
495
496module Isolation PhotonIsolation {
497 set CandidateInputArray PhotonEfficiency/photons
498 set IsolationInputArray EFlowMerger/eflow
499 set RhoInputArray Rho/rho
500
501 set OutputArray photons
502
503 set DeltaRMax 0.5
504
505 set PTMin 0.5
506
507 set PTRatioMax 0.1
508}
509
510#####################
511# Electron efficiency
512#####################
513
514module Efficiency ElectronEfficiency {
515 set InputArray TrackPileUpSubtractor/electrons
516 set OutputArray electrons
517
518 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
519
520 # efficiency formula for electrons
[a2983ec]521 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
522 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
523 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
[89c5be4]524 (abs(eta) > 2.5) * (0.00)}
525}
526
527####################
528# Electron isolation
529####################
530
531module Isolation ElectronIsolation {
532 set CandidateInputArray ElectronEfficiency/electrons
533 set IsolationInputArray EFlowMerger/eflow
534 set RhoInputArray Rho/rho
535
536 set OutputArray electrons
537
538 set DeltaRMax 0.5
539
540 set PTMin 0.5
541
542 set PTRatioMax 0.1
543}
544
545#################
546# Muon efficiency
547#################
548
549module Efficiency MuonEfficiency {
550 set InputArray TrackPileUpSubtractor/muons
551 set OutputArray muons
552
553 # set EfficiencyFormula {efficiency as a function of eta and pt}
554
555 # efficiency formula for muons
[a2983ec]556 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
557 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
558 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) +
[89c5be4]559 (abs(eta) > 2.7) * (0.00)}
560}
561
562################
563# Muon isolation
564################
565
566module Isolation MuonIsolation {
567 set CandidateInputArray MuonEfficiency/muons
568 set IsolationInputArray EFlowMerger/eflow
569 set RhoInputArray Rho/rho
570
571 set OutputArray muons
572
573 set DeltaRMax 0.5
574
575 set PTMin 0.5
576
577 set PTRatioMax 0.1
578}
579
580###################
581# Missing ET merger
582###################
583
584module Merger MissingET {
585# add InputArray InputArray
[6fb1a5d]586 add InputArray EFlowMergerAllTracks/eflow
[89c5be4]587 set MomentumOutputArray momentum
588}
589
[27bf162]590
[89c5be4]591##################
592# Scalar HT merger
593##################
594
595module Merger ScalarHT {
596# add InputArray InputArray
597 add InputArray UniqueObjectFinder/jets
598 add InputArray UniqueObjectFinder/electrons
599 add InputArray UniqueObjectFinder/photons
600 add InputArray UniqueObjectFinder/muons
601 set EnergyOutputArray energy
602}
603
604###########
605# b-tagging
606###########
607
608module BTagging BTagging {
609 set PartonInputArray Delphes/partons
610 set JetInputArray JetEnergyScale/jets
611
612 set BitNumber 0
613
614 set DeltaR 0.5
615
616 set PartonPTMin 1.0
617
618 set PartonEtaMax 2.5
619
620 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
621 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
622 # gluon's PDG code has the lowest priority
623
624 # default efficiency formula (misidentification rate)
625 add EfficiencyFormula {0} {0.001}
626
627 # efficiency formula for c-jets (misidentification rate)
[a2983ec]628 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) +
629 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
630 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
[89c5be4]631 (abs(eta) > 2.5) * (0.000)}
632
633 # efficiency formula for b-jets
[a2983ec]634 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) +
635 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
636 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
[89c5be4]637 (abs(eta) > 2.5) * (0.000)}
638}
639
640module TauTagging TauTagging {
641 set ParticleInputArray Delphes/allParticles
642 set PartonInputArray Delphes/partons
643 set JetInputArray JetEnergyScale/jets
644
645 set DeltaR 0.5
646
647 set TauPTMin 1.0
648
649 set TauEtaMax 2.5
650
651 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
652
653 # default efficiency formula (misidentification rate)
654 add EfficiencyFormula {0} {0.001}
655 # efficiency formula for tau-jets
656 add EfficiencyFormula {15} {0.4}
657}
658
659#####################################################
660# Find uniquely identified photons/electrons/tau/jets
661#####################################################
662
663module UniqueObjectFinder UniqueObjectFinder {
664# earlier arrays take precedence over later ones
665# add InputArray InputArray OutputArray
666 add InputArray PhotonIsolation/photons photons
667 add InputArray ElectronIsolation/electrons electrons
668 add InputArray MuonIsolation/muons muons
669 add InputArray JetEnergyScale/jets jets
670}
671
672##################
673# ROOT tree writer
674##################
675
[27bf162]676# tracks, towers and eflow objects are not stored by default in the output.
[282f591]677# if needed (for jet constituent or other studies), uncomment the relevant
678# "add Branch ..." lines.
[27bf162]679
[89c5be4]680module TreeWriter TreeWriter {
681# add Branch InputArray BranchName BranchClass
682 add Branch Delphes/allParticles Particle GenParticle
[27bf162]683
684# add Branch TrackMerger/tracks Track Track
[1c8d9db]685 add Branch Calorimeter/towers Tower Tower
[282f591]686
[89c5be4]687# add Branch Calorimeter/eflowTracks EFlowTrack Track
[27bf162]688# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
689# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
690
[89c5be4]691 add Branch GenJetFinder/jets GenJet Jet
692 add Branch UniqueObjectFinder/jets Jet Jet
693 add Branch UniqueObjectFinder/electrons Electron Electron
694 add Branch UniqueObjectFinder/photons Photon Photon
695 add Branch UniqueObjectFinder/muons Muon Muon
696 add Branch MissingET/momentum MissingET MissingET
697 add Branch ScalarHT/energy ScalarHT ScalarHT
698 add Branch Rho/rho Rho Rho
[fb21bc8]699 add Branch PileUpMerger/vertices Vertex Vertex
[27bf162]700
[89c5be4]701}
Note: See TracBrowser for help on using the repository browser.