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 Michele Selvaggi, 10 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.