Fork me on GitHub

Ticket #419: ITK-like_Delphes_ATLAS_card_withPileup.tcl

File ITK-like_Delphes_ATLAS_card_withPileup.tcl, 23.6 KB (added by Jessica Metcalfe, 10 years ago)
Line 
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
21 NeutralTowerMerger
22 EFlowMerger
23
24 GenJetFinder
25
26 Rho
27 FastJetFinder
28 JetPileUpSubtractor
29
30 JetEnergyScale
31
32 PhotonEfficiency
33 PhotonIsolation
34
35 ElectronEfficiency
36 ElectronIsolation
37
38 MuonEfficiency
39 MuonIsolation
40
41 MissingET
42
43 BTagging
44 TauTagging
45
46 UniqueObjectFinder
47
48 ScalarHT
49
50 TreeWriter
51}
52
53###############
54# PileUp Merger
55###############
56
57module PileUpMerger PileUpMerger {
58 set InputArray Delphes/stableParticles
59
60 set ParticleOutputArray stableParticles
61 set VertexOutputArray vertices
62
63 # pre-generated minbias input file
64 set PileUpFile MinBias100K_14TeV.pileup
65
66 # average expected pile up
67 set MeanPileUp 80
68 #80
69
70 # maximum spread in the beam direction in m
71 set ZVertexSpread 0.10
72
73
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)
78
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
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# These efficiencies are combined from Delphes-3.2.1 ATLAS_PileUP Card and LoI for Phase-II upgrade.
114# LoI: http://cds.cern.ch/record/1502664/files/LHCC-I-023.pdf from p.66, Figure 6.9.
115# For both eta regions and pT<0.1 GeV, the efficiencies are taken from Delphes-3.2.1 ATLAS_PileUP Card.
116# For both eta regions and pT>0.1 GeV the efficiencies are taken from LoI and represent Pion tracking efficiencies.
117# -> efficiencies for 5GeV Pions are taken for 1.0 < pT < 100 range and efficiencies for 100GeV Pions are taken for pT>100GeV range
118
119module Efficiency ChargedHadronTrackingEfficiency {
120 set InputArray ParticlePropagator/chargedHadrons
121 set OutputArray chargedHadrons
122
123 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
124
125 # tracking efficiency formula for charged hadrons
126 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
127 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
128 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.92) + \
129 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.94) + \
130 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
131 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 1.0e2) * (0.85) + \
132 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0e2) * (0.89) + \
133 (abs(eta) > 4.0) * (0.00)}
134}
135
136##############################
137# Electron tracking efficiency
138##############################
139# These efficiencies are combined from Delphes-3.2.1 ATLAS_PileUP Card and LoI for Phase-II upgrade.
140# LoI: http://cds.cern.ch/record/1502664/files/LHCC-I-023.pdf from p.66, Figure 6.9.
141# For both eta regions and pT<0.1 GeV, the efficiencies are taken from Delphes-3.2.1 ATLAS_PileUP Card.
142# For both eta regions and pT>0.1 GeV the efficiencies are taken from LoI and represent Pion tracking efficiencies.
143# -> efficiencies for 5GeV Electrons are taken for 1.0 < pT < 100 range and efficiencies for 100GeV Electrons are taken for pT>100GeV range
144
145
146module Efficiency ElectronTrackingEfficiency {
147 set InputArray ParticlePropagator/electrons
148 set OutputArray electrons
149
150 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
151
152 # tracking efficiency formula for electrons
153 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
154 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
155 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.90) + \
156 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.96) + \
157 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
158 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 1.0e2) * (0.80) + \
159 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0e2) * (0.94) + \
160 (abs(eta) > 4.0) * (0.00)}
161}
162
163##########################
164# Muon tracking efficiency
165##########################
166# These efficiencies are combined from Delphes-3.2.1 ATLAS_PileUP Card and LoI for Phase-II upgrade.
167# LoI: http://cds.cern.ch/record/1502664/files/LHCC-I-023.pdf from p.66, Figure 6.9.
168# For both eta regions and pT<0.1 GeV, the efficiencies are taken from Delphes-3.2.1 ATLAS_PileUP Card.
169# For both eta regions and pT>0.1 GeV the efficiencies are taken from LoI and represent Pion tracking efficiencies.
170# -> efficiencies for 5GeV Muons are taken for 1.0 < pT < 100 range and efficiencies for 100GeV Muons are taken for pT>100GeV range
171
172module Efficiency MuonTrackingEfficiency {
173 set InputArray ParticlePropagator/muons
174 set OutputArray muons
175
176 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
177
178 # tracking efficiency formula for muons
179 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
180 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
181 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.99) + \
182 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
183 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
184 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 1.0e2) * (0.99) + \
185 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0e2) * (0.99) + \
186 (abs(eta) > 4.0) * (0.00)}
187}
188
189########################################
190# Momentum resolution for charged tracks
191########################################
192# Taken from LoI for Phase-II upgrade:
193# LoI: http://cds.cern.ch/record/1502664/files/LHCC-I-023.pdf from p.63, Figure 6.5.
194# Assuming that the momentum resolution for charged tracks is the same as for the muons.
195
196module MomentumSmearing ChargedHadronMomentumSmearing {
197 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
198 set OutputArray chargedHadrons
199
200 # set ResolutionFormula {resolution formula as a function of eta and pt}
201
202 # resolution formula for charged hadrons
203 set ResolutionFormula {
204 (abs(eta) <= 1.5) * (pt <= 1.0e2) * (0.015) + \
205 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.02) + \
206 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt <= 1.0e2) * (0.05) + \
207 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0e2) * (0.07)}
208}
209
210#################################
211# Energy resolution for electrons
212#################################
213# Taken from ES note: http://cds.cern.ch/record/1527529/files/ATL-PHYS-PUB-2013-004.pdf
214# p7, section 3.7
215
216module EnergySmearing ElectronEnergySmearing {
217 set InputArray ElectronTrackingEfficiency/electrons
218 set OutputArray electrons
219
220 # set ResolutionFormula {resolution formula as a function of eta }
221
222 set ResolutionFormula { (abs(eta) <= 1.4) * sqrt( 0.3^2 + energy*0.10^2 + energy^2*0.010^2) + \
223 (abs(eta) > 1.4 && abs(eta) <= 4.0) * sqrt( 0.3^2 + energy*0.15^2 + energy^2*0.015^2)}
224
225}
226
227###############################
228# Momentum resolution for muons
229###############################
230# Taken from LoI for Phase-II upgrade:
231# LoI: http://cds.cern.ch/record/1502664/files/LHCC-I-023.pdf from p.63, Figure 6.5.
232
233
234module MomentumSmearing MuonMomentumSmearing {
235 set InputArray MuonTrackingEfficiency/muons
236 set OutputArray muons
237
238 # set ResolutionFormula {resolution formula as a function of eta and pt}
239
240 # resolution formula for muons
241 set ResolutionFormula { (abs(eta) <= 1.5) * (pt <= 1.0e2) * (0.015) + \
242 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.02) + \
243 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt <= 1.0e2) * (0.05) + \
244 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 1.0e2) * (0.07)}
245}
246
247##############
248# Track merger
249##############
250
251module Merger TrackMerger {
252# add InputArray InputArray
253 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
254 add InputArray ElectronEnergySmearing/electrons
255 add InputArray MuonMomentumSmearing/muons
256 set OutputArray tracks
257}
258
259#############
260# Calorimeter
261#############
262# No changes here
263
264module Calorimeter Calorimeter {
265 set ParticleInputArray ParticlePropagator/stableParticles
266 set TrackInputArray TrackMerger/tracks
267
268 set TowerOutputArray towers
269 set PhotonOutputArray photons
270
271 set EFlowTrackOutputArray eflowTracks
272 set EFlowPhotonOutputArray eflowPhotons
273 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
274
275 set pi [expr {acos(-1)}]
276
277 # lists of the edges of each tower in eta and phi
278 # each list starts with the lower edge of the first tower
279 # the list ends with the higher edged of the last tower
280
281 # 10 degrees towers
282 set PhiBins {}
283 for {set i -18} {$i <= 18} {incr i} {
284 add PhiBins [expr {$i * $pi/18.0}]
285 }
286 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} {
287 add EtaPhiBins $eta $PhiBins
288 }
289
290 # 20 degrees towers
291 set PhiBins {}
292 for {set i -9} {$i <= 9} {incr i} {
293 add PhiBins [expr {$i * $pi/9.0}]
294 }
295 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} {
296 add EtaPhiBins $eta $PhiBins
297 }
298
299 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
300 add EnergyFraction {0} {0.0 1.0}
301 # energy fractions for e, gamma and pi0
302 add EnergyFraction {11} {1.0 0.0}
303 add EnergyFraction {22} {1.0 0.0}
304 add EnergyFraction {111} {1.0 0.0}
305 # energy fractions for muon, neutrinos and neutralinos
306 add EnergyFraction {12} {0.0 0.0}
307 add EnergyFraction {13} {0.0 0.0}
308 add EnergyFraction {14} {0.0 0.0}
309 add EnergyFraction {16} {0.0 0.0}
310 add EnergyFraction {1000022} {0.0 0.0}
311 add EnergyFraction {1000023} {0.0 0.0}
312 add EnergyFraction {1000025} {0.0 0.0}
313 add EnergyFraction {1000035} {0.0 0.0}
314 add EnergyFraction {1000045} {0.0 0.0}
315 # energy fractions for K0short and Lambda
316 add EnergyFraction {310} {0.3 0.7}
317 add EnergyFraction {3122} {0.3 0.7}
318
319 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
320 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
321 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
322 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
323 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
324 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
325
326 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
327 # http://arxiv.org/pdf/hep-ex/0004009v1
328 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
329 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
330 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
331 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
332}
333
334##########################
335# Track pile-up subtractor
336##########################
337
338module TrackPileUpSubtractor TrackPileUpSubtractor {
339# add InputArray InputArray OutputArray
340 add InputArray Calorimeter/eflowTracks eflowTracks
341 add InputArray ElectronEnergySmearing/electrons electrons
342 add InputArray MuonMomentumSmearing/muons muons
343
344 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
345 # Z vertex resolution in m
346 set ZVertexResolution 0.0001
347 #0.0001
348}
349
350####################
351# Neutral tower merger
352####################
353
354module Merger NeutralTowerMerger {
355# add InputArray InputArray
356 add InputArray Calorimeter/eflowPhotons
357 add InputArray Calorimeter/eflowNeutralHadrons
358 set OutputArray eflowTowers
359}
360
361
362####################
363# Energy flow merger
364####################
365
366module Merger EFlowMerger {
367# add InputArray InputArray
368 add InputArray Calorimeter/eflowTracks
369 add InputArray Calorimeter/eflowPhotons
370 add InputArray Calorimeter/eflowNeutralHadrons
371 set OutputArray eflow
372}
373
374
375#############
376# Rho pile-up
377#############
378
379module FastJetFinder Rho {
380 set InputArray Calorimeter/towers
381
382 set ComputeRho true
383 set RhoOutputArray rho
384
385 # 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
386 set AreaAlgorithm 5
387
388 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
389 set JetAlgorithm 4
390 set ParameterR 0.4
391 set GhostEtaMax 5.0
392 set RhoEtaMax 5.0
393
394 add RhoEtaRange 0.0 5.0
395
396 set JetPTMin 0.0
397}
398
399#####################
400# MC truth jet finder
401#####################
402
403module FastJetFinder GenJetFinder {
404 set InputArray Delphes/stableParticles
405
406 set OutputArray jets
407
408 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
409 set JetAlgorithm 6
410 set ParameterR 0.4
411
412 set JetPTMin 20.0
413}
414
415############
416# Jet finder
417############
418
419module FastJetFinder FastJetFinder {
420 set InputArray Calorimeter/towers
421
422 set OutputArray jets
423
424 # 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
425 set AreaAlgorithm 5
426
427 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
428 set JetAlgorithm 6
429 set ParameterR 0.4
430
431 set JetPTMin 20.0
432}
433
434###########################
435# Jet Pile-Up Subtraction
436###########################
437
438module JetPileUpSubtractor JetPileUpSubtractor {
439 set JetInputArray FastJetFinder/jets
440 set RhoInputArray Rho/rho
441
442 set OutputArray jets
443
444 set JetPTMin 20.0
445}
446
447##################
448# Jet Energy Scale
449##################
450
451module EnergyScale JetEnergyScale {
452 set InputArray JetPileUpSubtractor/jets
453 set OutputArray jets
454
455 # scale formula for jets
456 set ScaleFormula {1.0}
457}
458
459###################
460# Photon efficiency
461###################
462# Taken from ECFA: http://cds.cern.ch/record/1604420/files/ATL-PHYS-PUB-2013-009.pdf
463# p.2
464
465module Efficiency PhotonEfficiency {
466 set InputArray Calorimeter/photons
467 set OutputArray photons
468
469 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
470
471 # efficiency formula for photons
472 set EfficiencyFormula { (pt <= 15.0) * (0.00) + \
473 (abs(eta) <= 4.0) * (pt > 15.0) * ( 0.76 - 1.98 * exp( (-1.0*pt)/16.1 ) ) + \
474 (abs(eta) > 4.0) * (0.00)}
475}
476
477##################
478# Photon isolation
479##################
480# All isolation criteria were optimised to dR=0.3, pTmin=1.0GeV and pTratio=0.1
481
482module Isolation PhotonIsolation {
483 set CandidateInputArray PhotonEfficiency/photons
484 set IsolationInputArray EFlowMerger/eflow
485 set RhoInputArray Rho/rho
486
487 set OutputArray photons
488
489 set DeltaRMax 0.3
490
491 set PTMin 1.0
492
493 set PTRatioMax 0.1
494}
495
496#####################
497# Electron efficiency
498#####################
499# Taken from ES: http://cds.cern.ch/record/1604420/files/ATL-PHYS-PUB-2013-009.pdf
500# p.6, tight efficiency
501
502module Efficiency ElectronEfficiency {
503 set InputArray TrackPileUpSubtractor/electrons
504 set OutputArray electrons
505
506 # set EfficiencyFormula {efficiency formula as a function of eta}
507
508 # efficiency formula for electrons
509 set EfficiencyFormula {
510 (abs(eta) <= 4.0) * ( 0.85 - 0.191 * exp( 1.0 - pt/20.0)) + \
511 (abs(eta) > 4.0) * (0.00)}
512}
513
514
515####################
516# Electron isolation
517####################
518
519module Isolation ElectronIsolation {
520 set CandidateInputArray ElectronEfficiency/electrons
521 set IsolationInputArray EFlowMerger/eflow
522 set RhoInputArray Rho/rho
523
524 set OutputArray electrons
525
526 set DeltaRMax 0.3
527
528 set PTMin 1.0
529
530 set PTRatioMax 0.1
531}
532
533#################
534# Muon efficiency
535#################
536# Taken from ECFA: http://cds.cern.ch/record/1604420/files/ATL-PHYS-PUB-2013-009.pdf
537# p.3, table 1, for pT>10
538
539module Efficiency MuonEfficiency {
540 set InputArray TrackPileUpSubtractor/muons
541 set OutputArray muons
542
543 # set EfficiencyFormula {efficiency as a function of eta and pt}
544
545 # efficiency formula for muons
546 set EfficiencyFormula {
547 (abs(eta) <= 0.1) * (0.54) + \
548 (abs(eta) > 0.1 && abs(eta) <= 4.0) * (0.97) + \
549 (abs(eta) > 4.0) * (0.00)}
550}
551
552################
553# Muon isolation
554################
555
556module Isolation MuonIsolation {
557 set CandidateInputArray MuonEfficiency/muons
558 set IsolationInputArray EFlowMerger/eflow
559 set RhoInputArray Rho/rho
560
561 set OutputArray muons
562
563 set DeltaRMax 0.3
564
565 set PTMin 1.0
566
567 set PTRatioMax 0.1
568}
569
570###################
571# Missing ET merger
572###################
573
574module Merger MissingET {
575# add InputArray InputArray
576 add InputArray EFlowMerger/eflow
577 set MomentumOutputArray momentum
578}
579
580
581##################
582# Scalar HT merger
583##################
584
585module Merger ScalarHT {
586# add InputArray InputArray
587 add InputArray UniqueObjectFinder/jets
588 add InputArray UniqueObjectFinder/electrons
589 add InputArray UniqueObjectFinder/photons
590 add InputArray UniqueObjectFinder/muons
591 set EnergyOutputArray energy
592}
593
594###########
595# b-tagging
596###########
597
598module BTagging BTagging {
599 set PartonInputArray Delphes/partons
600 set JetInputArray JetEnergyScale/jets
601
602 set BitNumber 0
603
604 set DeltaR 0.5
605
606 set PartonPTMin 1.0
607
608 set PartonEtaMax 4.0
609
610 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
611 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
612 # gluon's PDG code has the lowest priority
613
614 # default efficiency formula (misidentification rate)
615 add EfficiencyFormula {0} {0.001}
616
617 # efficiency formula for c-jets (misidentification rate)
618 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
619 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
620 (abs(eta) > 1.2 && abs(eta) <= 4.0) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
621 (abs(eta) > 4.0) * (0.000)}
622
623 # efficiency formula for b-jets
624 # Taken from TruthToRecoFunctions, from ECFA: http://cds.cern.ch/record/1604420/files/ATL-PHYS-PUB-2013-009.pdf
625 # if (jet_pt<20) jet_pt = 20;
626 # if (jet_pt>300) jet_pt = 300;
627
628 add EfficiencyFormula {5} {
629 (abs(eta) <= 4.0) * (pt < 20.0) * ( 0.394 - 0.342*log(pt) - 0.298*abs(20.0) + 0.194*(log(pt))^2 - 0.060*(abs(20.0))^2 + 0.082*log(pt)*abs(20.0) + 0.035*(abs(20.0))^3 -0.022*(log(pt))^3 -0.003*(log(pt)*abs(20.0))^2 ) + \
630
631 (abs(eta) <= 4.0) * (pt >= 20.0 && pt <= 300.0) * ( 0.394 - 0.342*log(pt) - 0.298*abs(eta) + 0.194*(log(pt))^2 - 0.060*(abs(eta))^2 + 0.082*log(pt)*abs(eta) + 0.035*(abs(eta))^3 -0.022*(log(pt))^3 -0.003*(log(pt)*abs(eta))^2 ) + \
632
633 (abs(eta) <= 4.0) * (pt > 300.0) * ( 0.394 - 0.342*log(pt) - 0.298*abs(300.0) + 0.194*(log(pt))^2 - 0.060*(abs(300.0))^2 + 0.082*log(pt)*abs(300.0) + 0.035*(abs(300.0))^3 -0.022*(log(pt))^3 -0.003*(log(pt)*abs(300.0))^2 ) + \
634
635 (abs(eta) > 4.0) * (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
674# tracks, towers and eflow objects are not stored by default in the output.
675# if needed (for jet constituent or other studies), uncomment the relevant
676# "add Branch ..." lines.
677
678module TreeWriter TreeWriter {
679# add Branch InputArray BranchName BranchClass
680 add Branch Delphes/allParticles Particle GenParticle
681
682# add Branch TrackMerger/tracks Track Track
683# add Branch Calorimeter/towers Tower Tower
684
685# add Branch Calorimeter/eflowTracks EFlowTrack Track
686# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
687# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
688
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
697 add Branch PileUpMerger/vertices Vertex Vertex
698
699 # add Branch PileUpMerger/stableParticles PileUpMix GenParticle
700
701}
702
703
704