Fork me on GitHub

source: svn/trunk/examples/delphes_card_CMS_PileUp.tcl@ 1396

Last change on this file since 1396 was 1378, checked in by Pavel Demin, 11 years ago

fix formatting
File size: 22.4 KB
Line 
1
2#######################################
3# Order of execution of various modules
4#######################################
5
6set ExecutionPath {
7
8 PileUpMerger
9 ParticlePropagator
10
11 ChargedHadronTrackingEfficiency
12 ElectronTrackingEfficiency
13 MuonTrackingEfficiency
14
15 ChargedHadronMomentumSmearing
16 ElectronEnergySmearing
17 MuonMomentumSmearing
18
19 TrackMerger
20 Calorimeter
21 TrackPileUpSubtractor
22 NeutralTowerMerger
23 EFlowMerger
24
25 GenJetFinder
26
27 Rho
28 FastJetFinder
29 PileUpJetID
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
62 set ParticleOutputArray stableParticles
63 set VertexOutputArray vertices
64
65 # pre-generated minbias input file
66 set PileUpFile MinBias.pileup
67
68 # average expected pile up
69 set MeanPileUp 10
70
71 # maximum spread in the beam direction in m
72 set ZVertexSpread 0.10
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.29
103 # half-length of the magnetic field coverage, in m
104 set HalfLength 3.00
105
106 # magnetic field
107 set Bz 3.8
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
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)}
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 # resolution formula for electrons
201 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) + \
202 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
203 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
204 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
205
206}
207
208###############################
209# Momentum resolution for muons
210###############################
211
212module MomentumSmearing MuonMomentumSmearing {
213 set InputArray MuonTrackingEfficiency/muons
214 set OutputArray muons
215
216 # set ResolutionFormula {resolution formula as a function of eta and pt}
217
218 # resolution formula for muons
219 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1 && pt <= 5.0) * (0.02) + \
220 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) + \
221 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
222 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
223 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) + \
224 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) + \
225 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
226 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) + \
227 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) + \
228 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) + \
229 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
230 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
231}
232
233##############
234# Track merger
235##############
236
237module Merger TrackMerger {
238# add InputArray InputArray
239 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
240 add InputArray ElectronEnergySmearing/electrons
241 add InputArray MuonMomentumSmearing/muons
242 set OutputArray tracks
243}
244
245#############
246# Calorimeter
247#############
248
249module Calorimeter Calorimeter {
250 set ParticleInputArray ParticlePropagator/stableParticles
251 set TrackInputArray TrackMerger/tracks
252
253 set TowerOutputArray towers
254 set PhotonOutputArray photons
255
256 set EFlowTrackOutputArray eflowTracks
257 set EFlowPhotonOutputArray eflowPhotons
258 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
259
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 # 5 degrees towers
267 set PhiBins {}
268 for {set i -36} {$i <= 36} {incr i} {
269 add PhiBins [expr {$i * $pi/36.0}]
270 }
271 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
272 add EtaPhiBins $eta $PhiBins
273 }
274
275 # 10 degrees towers
276 set PhiBins {}
277 for {set i -18} {$i <= 18} {incr i} {
278 add PhiBins [expr {$i * $pi/18.0}]
279 }
280 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
281 add EtaPhiBins $eta $PhiBins
282 }
283
284 # 20 degrees towers
285 set PhiBins {}
286 for {set i -9} {$i <= 9} {incr i} {
287 add PhiBins [expr {$i * $pi/9.0}]
288 }
289 foreach eta {-5 -4.7 -4.525 4.7 5} {
290 add EtaPhiBins $eta $PhiBins
291 }
292
293 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
294 add EnergyFraction {0} {0.0 1.0}
295 # energy fractions for e, gamma and pi0
296 add EnergyFraction {11} {1.0 0.0}
297 add EnergyFraction {22} {1.0 0.0}
298 add EnergyFraction {111} {1.0 0.0}
299 # energy fractions for muon, neutrinos and neutralinos
300 add EnergyFraction {12} {0.0 0.0}
301 add EnergyFraction {13} {0.0 0.0}
302 add EnergyFraction {14} {0.0 0.0}
303 add EnergyFraction {16} {0.0 0.0}
304 add EnergyFraction {1000022} {0.0 0.0}
305 add EnergyFraction {1000023} {0.0 0.0}
306 add EnergyFraction {1000025} {0.0 0.0}
307 add EnergyFraction {1000035} {0.0 0.0}
308 add EnergyFraction {1000045} {0.0 0.0}
309 # energy fractions for K0short and Lambda
310 add EnergyFraction {310} {0.3 0.7}
311 add EnergyFraction {3122} {0.3 0.7}
312
313 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
314 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
315 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
316
317 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
318 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
319 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
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 set VertexInputArray PileUpMerger/vertices
333 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
334 # Z vertex resolution in m
335 set ZVertexResolution 0.0001
336}
337
338####################
339# Neutral tower merger
340####################
341
342module Merger NeutralTowerMerger {
343# add InputArray InputArray
344 add InputArray Calorimeter/eflowPhotons
345 add InputArray Calorimeter/eflowNeutralHadrons
346 set OutputArray eflowTowers
347}
348
349
350####################
351# Energy flow merger
352####################
353
354module Merger EFlowMerger {
355# add InputArray InputArray
356 add InputArray TrackPileUpSubtractor/eflowTracks
357 add InputArray Calorimeter/eflowPhotons
358 add InputArray Calorimeter/eflowNeutralHadrons
359 set OutputArray eflow
360}
361
362
363#############
364# Rho pile-up
365#############
366
367module FastJetFinder Rho {
368# set InputArray Calorimeter/towers
369 set InputArray EFlowMerger/eflow
370
371 set ComputeRho true
372 set RhoOutputArray rho
373
374 # 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
375 set AreaAlgorithm 5
376
377 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
378 set JetAlgorithm 4
379 set ParameterR 0.6
380 set GhostEtaMax 5.0
381
382 add RhoEtaRange 0.0 2.5
383 add RhoEtaRange 2.5 5.0
384
385 set JetPTMin 0.0
386}
387
388#####################
389# MC truth jet finder
390#####################
391
392module FastJetFinder GenJetFinder {
393 set InputArray Delphes/stableParticles
394
395 set OutputArray jets
396
397 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
398 set JetAlgorithm 6
399 set ParameterR 0.5
400
401 set JetPTMin 20.0
402}
403
404############
405# Jet finder
406############
407
408module FastJetFinder FastJetFinder {
409# set InputArray Calorimeter/towers
410 set InputArray EFlowMerger/eflow
411
412 set OutputArray jets
413
414 # 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
415 set AreaAlgorithm 5
416
417 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
418 set JetAlgorithm 6
419 set ParameterR 0.5
420
421 set JetPTMin 20.0
422}
423
424###########################
425# Jet Pile-Up ID
426###########################
427
428module PileUpJetID PileUpJetID {
429 set JetInputArray FastJetFinder/jets
430 set TrackInputArray Calorimeter/eflowTracks
431 set NeutralInputArray NeutralTowerMerger/eflowTowers
432
433 set VertexInputArray PileUpMerger/vertices
434 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
435 # Z vertex resolution in m
436 set ZVertexResolution 0.0001
437
438 set OutputArray jets
439
440 set UseConstituents 0
441 set ParameterR 0.5
442
443 set JetPTMin 20.0
444}
445
446###########################
447# Jet Pile-Up Subtraction
448###########################
449
450module JetPileUpSubtractor JetPileUpSubtractor {
451 set JetInputArray PileUpJetID/jets
452 set RhoInputArray Rho/rho
453
454 set OutputArray jets
455
456 set JetPTMin 20.0
457}
458
459##################
460# Jet Energy Scale
461##################
462
463module EnergyScale JetEnergyScale {
464 set InputArray JetPileUpSubtractor/jets
465 set OutputArray jets
466
467 # scale formula for jets
468 set ScaleFormula {1.0}
469}
470
471###################
472# Photon efficiency
473###################
474
475module Efficiency PhotonEfficiency {
476 set InputArray Calorimeter/eflowPhotons
477 set OutputArray photons
478
479 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
480
481 # efficiency formula for photons
482 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
483 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
484 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
485 (abs(eta) > 2.5) * (0.00)}
486}
487
488
489##################
490# Photon isolation
491##################
492
493module Isolation PhotonIsolation {
494 set CandidateInputArray PhotonEfficiency/photons
495 set IsolationInputArray EFlowMerger/eflow
496 set RhoInputArray Rho/rho
497
498 set OutputArray photons
499
500 set DeltaRMax 0.5
501
502 set PTMin 0.5
503
504 set PTRatioMax 0.1
505}
506
507#####################
508# Electron efficiency
509#####################
510
511module Efficiency ElectronEfficiency {
512 set InputArray TrackPileUpSubtractor/electrons
513 set OutputArray electrons
514
515 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
516
517 # efficiency formula for electrons
518 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
519 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
520 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
521 (abs(eta) > 2.5) * (0.00)}
522}
523
524####################
525# Electron isolation
526####################
527
528module Isolation ElectronIsolation {
529 set CandidateInputArray ElectronEfficiency/electrons
530 set IsolationInputArray EFlowMerger/eflow
531 set RhoInputArray Rho/rho
532
533 set OutputArray electrons
534
535 set DeltaRMax 0.5
536
537 set PTMin 0.5
538
539 set PTRatioMax 0.1
540}
541
542#################
543# Muon efficiency
544#################
545
546module Efficiency MuonEfficiency {
547 set InputArray TrackPileUpSubtractor/muons
548 set OutputArray muons
549
550 # set EfficiencyFormula {efficiency as a function of eta and pt}
551
552 # efficiency formula for muons
553 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
554 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
555 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
556 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
557 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
558 (abs(eta) > 2.4) * (0.00)}
559}
560
561################
562# Muon isolation
563################
564
565module Isolation MuonIsolation {
566 set CandidateInputArray MuonEfficiency/muons
567 set IsolationInputArray EFlowMerger/eflow
568 set RhoInputArray Rho/rho
569
570 set OutputArray muons
571
572 set DeltaRMax 0.5
573
574 set PTMin 0.5
575
576 set PTRatioMax 0.1
577}
578
579###################
580# Missing ET merger
581###################
582
583module Merger MissingET {
584# add InputArray InputArray
585 add InputArray EFlowMerger/eflow
586 set MomentumOutputArray momentum
587}
588
589
590
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 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
625 # default efficiency formula (misidentification rate)
626 add EfficiencyFormula {0} {0.001}
627
628 # efficiency formula for c-jets (misidentification rate)
629 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
630 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
631 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
632 (abs(eta) > 2.5) * (0.000)}
633
634 # efficiency formula for b-jets
635 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
636 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
637 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
638 (abs(eta) > 2.5) * (0.000)}
639}
640
641module TauTagging TauTagging {
642 set ParticleInputArray Delphes/allParticles
643 set PartonInputArray Delphes/partons
644 set JetInputArray JetEnergyScale/jets
645
646 set DeltaR 0.5
647
648 set TauPTMin 1.0
649
650 set TauEtaMax 2.5
651
652 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
653
654 # default efficiency formula (misidentification rate)
655 add EfficiencyFormula {0} {0.001}
656 # efficiency formula for tau-jets
657 add EfficiencyFormula {15} {0.4}
658}
659
660#####################################################
661# Find uniquely identified photons/electrons/tau/jets
662#####################################################
663
664module UniqueObjectFinder UniqueObjectFinder {
665# earlier arrays take precedence over later ones
666# add InputArray InputArray OutputArray
667 add InputArray PhotonIsolation/photons photons
668 add InputArray ElectronIsolation/electrons electrons
669 add InputArray MuonIsolation/muons muons
670 add InputArray JetEnergyScale/jets jets
671}
672
673##################
674# ROOT tree writer
675##################
676
677# tracks, towers and eflow objects are not stored by default in the output.
678# if needed (for jet constituent or other studies), uncomment the relevant
679# "add Branch ..." lines.
680
681module TreeWriter TreeWriter {
682# add Branch InputArray BranchName BranchClass
683 add Branch Delphes/allParticles Particle GenParticle
684
685# add Branch TrackMerger/tracks Track Track
686# add Branch Calorimeter/towers Tower Tower
687
688# add Branch Calorimeter/eflowTracks EFlowTrack Track
689# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
690# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
691
692 add Branch GenJetFinder/jets GenJet Jet
693 add Branch UniqueObjectFinder/jets Jet Jet
694 add Branch UniqueObjectFinder/electrons Electron Electron
695 add Branch UniqueObjectFinder/photons Photon Photon
696 add Branch UniqueObjectFinder/muons Muon Muon
697 add Branch MissingET/momentum MissingET MissingET
698 add Branch ScalarHT/energy ScalarHT ScalarHT
699 add Branch Rho/rho Rho Rho
700 add Branch PileUpMerger/vertices Vertex Vertex
701}
702
Note: See TracBrowser for help on using the repository browser.