Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_0PU.tcl@ b0ba2b4

ImprovedOutputFile Timing dual_readout llp 3.3.3pre07
Last change on this file since b0ba2b4 was b0ba2b4, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

removed eta dependence from e/gamma resolution in central region
  • Property mode set to 100644
File size: 22.5 KB
RevLine 
[f7f33af]1#
2# Phase II - No Pile-Up
3#
4# Main authors: Michele Selvaggi (UCL)
5#
6# Released on:
7#
8# Version: v01
9#
10#
11#######################################
12# Order of execution of various modules
13#######################################
14
15set ExecutionPath {
16 ParticlePropagator
17
18 ChargedHadronTrackingEfficiency
19 ElectronTrackingEfficiency
20 MuonTrackingEfficiency
21
22 ChargedHadronMomentumSmearing
[df7986a]23 ElectronEnergySmearing
[f7f33af]24 MuonMomentumSmearing
25
26 TrackMerger
27
28 ECal
29 HCal
30
31 TowerMerger
32 EFlowMerger
33
34 PhotonIsolation
35 PhotonEfficiency
36
37 ElectronFilter
38 ElectronIsolation
39 ElectronEfficiency
40
41 MuonIsolation
42 MuonLooseIdEfficiency
43 MuonTightIdEfficiency
44
45 NeutrinoFilter
46
47 MissingET
48 GenMissingET
49
50 GenJetFinder
51 FastJetFinder
52
53 ScalarHT
54
55 JetEnergyScale
56
57 JetFlavorAssociation
58
59 BTaggingLoose
60 BTaggingMedium
61 BTaggingTight
62
63 TauTagging
64
65 TreeWriter
66}
67
68#################################
69# Propagate particles in cylinder
70#################################
71
72module ParticlePropagator ParticlePropagator {
73 set InputArray Delphes/stableParticles
74
75 set OutputArray stableParticles
76 set ChargedHadronOutputArray chargedHadrons
77 set ElectronOutputArray electrons
78 set MuonOutputArray muons
79
80 # radius of the magnetic field coverage, in m
81 set Radius 1.29
82 # half-length of the magnetic field coverage, in m
83 set HalfLength 3.0
84
85 # magnetic field
86 set Bz 3.8
87}
88
89
90####################################
91# Charged hadron tracking efficiency
92####################################
93
94module Efficiency ChargedHadronTrackingEfficiency {
95 ## particles after propagation
96 set InputArray ParticlePropagator/chargedHadrons
97 set OutputArray chargedHadrons
98 # tracking efficiency formula for charged hadrons
99 set EfficiencyFormula {
100 (pt <= 0.2) * (0.00) + \
101 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
102 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
103 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
104 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
105 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
106 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
107 (abs(eta) > 4.0) * (0.00)
108 }
109}
110
111
112#####################################
113# Electron tracking efficiency - ID
114####################################
115
116module Efficiency ElectronTrackingEfficiency {
117 set InputArray ParticlePropagator/electrons
118 set OutputArray electrons
119 # tracking efficiency formula for electrons
120 set EfficiencyFormula {
121 (pt <= 0.2) * (0.00) + \
122 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
123 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
124 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
125 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
126 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
127 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
128 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
129 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
130 (abs(eta) > 4.0) * (0.00)
131
132 }
133}
134
135##########################
136# Muon tracking efficiency
137##########################
138
139module Efficiency MuonTrackingEfficiency {
140 set InputArray ParticlePropagator/muons
141 set OutputArray muons
142 # tracking efficiency formula for muons
143 set EfficiencyFormula {
144 (pt <= 0.2) * (0.00) + \
145 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
146 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
147 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
148 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
149 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
150 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
151 (abs(eta) > 4.0) * (0.00)
152
153 }
154}
155
156########################################
157# Momentum resolution for charged tracks
158########################################
159
160module MomentumSmearing ChargedHadronMomentumSmearing {
161 ## hadrons after having applied the tracking efficiency
162 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
163 set OutputArray chargedHadrons
164 # resolution formula for charged hadrons ,
165
166 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
167 source trackMomentumResolution.tcl
168}
169
[df7986a]170#################################
171# Energy resolution for electrons
172#################################
[f7f33af]173
[df7986a]174module EnergySmearing ElectronEnergySmearing {
175 set InputArray ElectronTrackingEfficiency/electrons
[f7f33af]176 set OutputArray electrons
[df7986a]177
178 # set ResolutionFormula {resolution formula as a function of eta and energy}
179
[f7f33af]180 # resolution formula for electrons
[df7986a]181
182 # taking something flat in energy for now, ECAL will take over at high energy anyway.
183 # inferred from hep-ex/1306.2016 and 1502.02701
184 set ResolutionFormula {
185
[b0ba2b4]186 (abs(eta) <= 1.5) * (energy*0.028) +
[df7986a]187 (abs(eta) > 1.5 && abs(eta) <= 1.75) * (energy*0.037) +
188 (abs(eta) > 1.75 && abs(eta) <= 2.15) * (energy*0.038) +
189 (abs(eta) > 2.15 && abs(eta) <= 3.00) * (energy*0.044) +
190 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)}
[f7f33af]191
192}
193
194###############################
195# Momentum resolution for muons
196###############################
197
198module MomentumSmearing MuonMomentumSmearing {
199 set InputArray MuonTrackingEfficiency/muons
200 set OutputArray muons
201 # resolution formula for muons
202
203 # up to |eta| < 2.8 take measurement from tracking + muon chambers
204 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
205 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
206 source muonMomentumResolution.tcl
207}
208
209
210##############
211# Track merger
212##############
213
214module Merger TrackMerger {
215# add InputArray InputArray
216 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
[df7986a]217 add InputArray ElectronEnergySmearing/electrons
[f7f33af]218 add InputArray MuonMomentumSmearing/muons
219 set OutputArray tracks
220}
221
222#############
223# ECAL
224#############
225
226module SimpleCalorimeter ECal {
227 set ParticleInputArray ParticlePropagator/stableParticles
228 set TrackInputArray TrackMerger/tracks
229
230 set TowerOutputArray ecalTowers
231 set EFlowTrackOutputArray eflowTracks
232 set EFlowTowerOutputArray eflowPhotons
233
234 set IsEcal true
235
236 set EnergyMin 0.5
237 set EnergySignificanceMin 1.0
238
239 set SmearTowerCenter true
240
241 set pi [expr {acos(-1)}]
242
243 # lists of the edges of each tower in eta and phi
244 # each list starts with the lower edge of the first tower
245 # the list ends with the higher edged of the last tower
246
247 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
248
249 set PhiBins {}
250 for {set i -180} {$i <= 180} {incr i} {
251 add PhiBins [expr {$i * $pi/180.0}]
252 }
253
254 # 0.02 unit in eta up to eta = 1.5 (barrel)
255 for {set i -85} {$i <= 86} {incr i} {
256 set eta [expr {$i * 0.0174}]
257 add EtaPhiBins $eta $PhiBins
258 }
259
260 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
261
262 set PhiBins {}
263 for {set i -180} {$i <= 180} {incr i} {
264 add PhiBins [expr {$i * $pi/180.0}]
265 }
266
267 # 0.02 unit in eta up to eta = 3
268 for {set i 1} {$i <= 84} {incr i} {
269 set eta [expr { -2.958 + $i * 0.0174}]
270 add EtaPhiBins $eta $PhiBins
271 }
272
273 for {set i 1} {$i <= 84} {incr i} {
274 set eta [expr { 1.4964 + $i * 0.0174}]
275 add EtaPhiBins $eta $PhiBins
276 }
277
278 # take present CMS granularity for HF
279
280 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
281 set PhiBins {}
282 for {set i -18} {$i <= 18} {incr i} {
283 add PhiBins [expr {$i * $pi/18.0}]
284 }
285
286 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
287 add EtaPhiBins $eta $PhiBins
288 }
289
290
291 add EnergyFraction {0} {0.0}
292 # energy fractions for e, gamma and pi0
293 add EnergyFraction {11} {1.0}
294 add EnergyFraction {22} {1.0}
295 add EnergyFraction {111} {1.0}
296 # energy fractions for muon, neutrinos and neutralinos
297 add EnergyFraction {12} {0.0}
298 add EnergyFraction {13} {0.0}
299 add EnergyFraction {14} {0.0}
300 add EnergyFraction {16} {0.0}
301 add EnergyFraction {1000022} {0.0}
302 add EnergyFraction {1000023} {0.0}
303 add EnergyFraction {1000025} {0.0}
304 add EnergyFraction {1000035} {0.0}
305 add EnergyFraction {1000045} {0.0}
306 # energy fractions for K0short and Lambda
307 add EnergyFraction {310} {0.3}
308 add EnergyFraction {3122} {0.3}
309
310 # set ResolutionFormula {resolution formula as a function of eta and energy}
311
[df7986a]312 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
[f7f33af]313 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
314
[b0ba2b4]315 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
[db7f5e4]316 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
317 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
318 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
319 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
[f7f33af]320}
321
322#############
323# HCAL
324#############
325
326module SimpleCalorimeter HCal {
327 set ParticleInputArray ParticlePropagator/stableParticles
328 set TrackInputArray ECal/eflowTracks
329
330 set TowerOutputArray hcalTowers
331 set EFlowTrackOutputArray eflowTracks
332 set EFlowTowerOutputArray eflowNeutralHadrons
333
334 set IsEcal false
335
336 set EnergyMin 1.0
337 set EnergySignificanceMin 1.0
338
339 set SmearTowerCenter true
340
341 set pi [expr {acos(-1)}]
342
343 # lists of the edges of each tower in eta and phi
344 # each list starts with the lower edge of the first tower
345 # the list ends with the higher edged of the last tower
346
347 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
348
349 set PhiBins {}
350 for {set i -36} {$i <= 36} {incr i} {
351 add PhiBins [expr {$i * $pi/36.0}]
352 }
353 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.65} {
354 add EtaPhiBins $eta $PhiBins
355 }
356
357 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
358
359 set PhiBins {}
360 for {set i -180} {$i <= 180} {incr i} {
361 add PhiBins [expr {$i * $pi/180.0}]
362 }
363
364 # 0.02 unit in eta up to eta = 3
365 for {set i 1} {$i <= 84} {incr i} {
366 set eta [expr { -2.958 + $i * 0.0174}]
367 add EtaPhiBins $eta $PhiBins
368 }
369
370 for {set i 1} {$i <= 84} {incr i} {
371 set eta [expr { 1.4964 + $i * 0.0174}]
372 add EtaPhiBins $eta $PhiBins
373 }
374
375 # take present CMS granularity for HF
376
377 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
378 set PhiBins {}
379 for {set i -18} {$i <= 18} {incr i} {
380 add PhiBins [expr {$i * $pi/18.0}]
381 }
382
383 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
384 add EtaPhiBins $eta $PhiBins
385 }
386
387 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
388 add EnergyFraction {0} {1.0}
389 # energy fractions for e, gamma and pi0
390 add EnergyFraction {11} {0.0}
391 add EnergyFraction {22} {0.0}
392 add EnergyFraction {111} {0.0}
393 # energy fractions for muon, neutrinos and neutralinos
394 add EnergyFraction {12} {0.0}
395 add EnergyFraction {13} {0.0}
396 add EnergyFraction {14} {0.0}
397 add EnergyFraction {16} {0.0}
398 add EnergyFraction {1000022} {0.0}
399 add EnergyFraction {1000023} {0.0}
400 add EnergyFraction {1000025} {0.0}
401 add EnergyFraction {1000035} {0.0}
402 add EnergyFraction {1000045} {0.0}
403 # energy fractions for K0short and Lambda
404 add EnergyFraction {310} {0.7}
405 add EnergyFraction {3122} {0.7}
406
407# set ResolutionFormula {resolution formula as a function of eta and energy}
408 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
409 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
410 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
411
412}
413
414#################
415# Electron filter
416#################
417
418module PdgCodeFilter ElectronFilter {
419 set InputArray HCal/eflowTracks
420 set OutputArray electrons
421 set Invert true
422 add PdgCode {11}
423 add PdgCode {-11}
424}
425
426###################################################
427# Tower Merger (in case not using e-flow algorithm)
428###################################################
429
430module Merger TowerMerger {
431# add InputArray InputArray
432 add InputArray ECal/ecalTowers
433 add InputArray HCal/hcalTowers
434 set OutputArray towers
435}
436
437####################
438# Energy flow merger
439####################
440
441module Merger EFlowMerger {
442# add InputArray InputArray
443 add InputArray HCal/eflowTracks
444 add InputArray ECal/eflowPhotons
445 add InputArray HCal/eflowNeutralHadrons
446 set OutputArray eflow
447}
448
449
450###################
451# Missing ET merger
452###################
453
454module Merger MissingET {
455# add InputArray InputArray
456 add InputArray EFlowMerger/eflow
457 set MomentumOutputArray momentum
458}
459
460
461
462##################
463# Scalar HT merger
464##################
465
466module Merger ScalarHT {
467# add InputArray InputArray
468 add InputArray EFlowMerger/eflow
469 set EnergyOutputArray energy
470}
471
472#################
473# Neutrino Filter
474#################
475
476module PdgCodeFilter NeutrinoFilter {
477
478 set InputArray Delphes/stableParticles
479 set OutputArray filteredParticles
480
481 set PTMin 0.0
482
483 add PdgCode {12}
484 add PdgCode {14}
485 add PdgCode {16}
486 add PdgCode {-12}
487 add PdgCode {-14}
488 add PdgCode {-16}
489
490}
491
492
493#####################
494# MC truth jet finder
495#####################
496
497module FastJetFinder GenJetFinder {
498 set InputArray NeutrinoFilter/filteredParticles
499
500 set OutputArray jets
501
502 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
503 set JetAlgorithm 6
504 set ParameterR 0.4
505
506 set JetPTMin 15.0
507}
508
509#########################
510# Gen Missing ET merger
511########################
512
513module Merger GenMissingET {
514
515# add InputArray InputArray
516 add InputArray NeutrinoFilter/filteredParticles
517 set MomentumOutputArray momentum
518}
519
520
521
522############
523# Jet finder
524############
525
526module FastJetFinder FastJetFinder {
527# set InputArray TowerMerger/towers
528 set InputArray EFlowMerger/eflow
529
530 set OutputArray jets
531
532 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
533 set JetAlgorithm 6
534 set ParameterR 0.4
535
536 set JetPTMin 15.0
537}
538
539##################
540# Jet Energy Scale
541##################
542
543module EnergyScale JetEnergyScale {
544 set InputArray FastJetFinder/jets
545 set OutputArray jets
546
547 # scale formula for jets
548 set ScaleFormula {1.00}
549}
550
551####################
552# Photon isolation #
553####################
554
555module Isolation PhotonIsolation {
556
557 # particle for which calculate the isolation
558 set CandidateInputArray ECal/eflowPhotons
559
560 # isolation collection
561 set IsolationInputArray EFlowMerger/eflow
562
563 # output array
564 set OutputArray photons
565
566 # isolation cone
567 set DeltaRMax 0.3
568
569 # minimum pT
570 set PTMin 1.0
571
572 # iso ratio to cut
573 set PTRatioMax 9999.
574}
575
576
577
578#####################
579# Photon efficiency #
580#####################
581
582module Efficiency PhotonEfficiency {
583
584 ## input particles
585 set InputArray PhotonIsolation/photons
586 ## output particles
587 set OutputArray photons
588 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
589 # efficiency formula for photons
590 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
591 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
592 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
593 (abs(eta) > 4.0) * (0.00)}
594
595}
596
597
598######################
599# Electron isolation #
600######################
601
602module Isolation ElectronIsolation {
603
604 set CandidateInputArray ElectronFilter/electrons
605
606 # isolation collection
607 set IsolationInputArray EFlowMerger/eflow
608
609 set OutputArray electrons
610
611 set DeltaRMax 0.3
612 set PTMin 1.0
613 set PTRatioMax 9999.
614
615}
616
617
618
619#######################
620# Electron efficiency #
621#######################
622
623module Efficiency ElectronEfficiency {
624
625 set InputArray ElectronIsolation/electrons
626 set OutputArray electrons
627
628 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
629 # efficiency formula for electrons
630 set EfficiencyFormula {
631 (pt <= 4.0) * (0.00) + \
632 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
633 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
634 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
635 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
636 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
637 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
638 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
639 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
640 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
641 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
642 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
643 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
644 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
645 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
646 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
647 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
648 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
649 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
650 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
651 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
652 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
653 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
654 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
655 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
656 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
657 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
658 (abs(eta) > 4.0) * (0.00)
659
660 }
661}
662
663##################
664# Muon isolation #
665##################
666
667module Isolation MuonIsolation {
668 set CandidateInputArray MuonMomentumSmearing/muons
669
670 # isolation collection
671 set IsolationInputArray EFlowMerger/eflow
672
673 set OutputArray muons
674
675 set DeltaRMax 0.3
676 set PTMin 1.0
677 set PTRatioMax 9999.
678
679}
680
681##################
682# Muon Loose Id #
683##################
684
685module Efficiency MuonLooseIdEfficiency {
686 set InputArray MuonIsolation/muons
687 set OutputArray muons
688 # tracking + LooseID efficiency formula for muons
689 source muonLooseId.tcl
690
691}
692
693
694##################
695# Muon Tight Id #
696##################
697
698module Efficiency MuonTightIdEfficiency {
699 set InputArray MuonIsolation/muons
700 set OutputArray muons
701 # tracking + TightID efficiency formula for muons
702 source muonTightId.tcl
703}
704
705########################
706# Jet Flavor Association
707########################
708
709module JetFlavorAssociation JetFlavorAssociation {
710
711 set PartonInputArray Delphes/partons
712 set ParticleInputArray Delphes/allParticles
713 set ParticleLHEFInputArray Delphes/allParticlesLHEF
714 set JetInputArray JetEnergyScale/jets
715
716 set DeltaR 0.5
717 set PartonPTMin 10.0
718 set PartonEtaMax 4.0
719
720}
721
722
723#############
724# b-tagging #
725#############
726module BTagging BTaggingLoose {
727
728 set JetInputArray JetEnergyScale/jets
729
730 set BitNumber 0
731
732 source btagLoose.tcl
733}
734
735
736#############
737# b-tagging #
738#############
739module BTagging BTaggingMedium {
740
741 set JetInputArray JetEnergyScale/jets
742
743 set BitNumber 1
744
745 source btagMedium.tcl
746}
747
748
749
750#############
751# b-tagging #
752#############
753module BTagging BTaggingTight {
754
755 set JetInputArray JetEnergyScale/jets
756
757 set BitNumber 2
758
759 source btagTight.tcl
760}
761
762
763
764#############
765# tau-tagging
766#############
767
768
769module TauTagging TauTagging {
770 set ParticleInputArray Delphes/allParticles
771 set PartonInputArray Delphes/partons
772 set JetInputArray JetEnergyScale/jets
773
774 set DeltaR 0.5
775
776 set TauPTMin 20.0
777
778 set TauEtaMax 2.3
779
780 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
781
782 add EfficiencyFormula {0} { (abs(eta) < 2.3) * ((( -0.00621816+0.00130097*pt-2.19642e-5*pt^2+1.49393e-7*pt^3-4.58972e-10*pt^4+5.27983e-13*pt^5 )) * (pt<250) + 0.0032*(pt>250)) + \
783 (abs(eta) > 2.3) * (0.000)
784 }
785 add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.97*0.77*( (0.32 + 0.01*pt - 0.000054*pt*pt )*(pt<100)+0.78*(pt>100) ) + \
786 (abs(eta) > 2.3) * (0.000)
787 }
788}
789
790
791##################
792# ROOT tree writer
793##################
794
795module TreeWriter TreeWriter {
796# add Branch InputArray BranchName BranchClass
797 add Branch Delphes/allParticles Particle GenParticle
798
799 add Branch GenJetFinder/jets GenJet Jet
800 add Branch GenMissingET/momentum GenMissingET MissingET
801
802 add Branch HCal/eflowTracks EFlowTrack Track
803 add Branch ECal/eflowPhotons EFlowPhoton Tower
804 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
805
806 add Branch PhotonEfficiency/photons Photon Photon
807 add Branch ElectronEfficiency/electrons Electron Electron
808 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
809 add Branch MuonTightIdEfficiency/muons MuonTight Muon
810
811 add Branch JetEnergyScale/jets Jet Jet
812
813 add Branch MissingET/momentum MissingET MissingET
814 add Branch ScalarHT/energy ScalarHT ScalarHT
815}
816
Note: See TracBrowser for help on using the repository browser.