Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_0PU_noHGCAL.tcl@ 83e77ee

Last change on this file since 83e77ee was a393b57, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

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