Fork me on GitHub

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

Last change on this file since 1397 was 1397, checked in by mselvaggi, 6 years ago

added Set RandomSeed? 0 to all data cards

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