Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_0PU.tcl@ 0a122b5

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

degrade by hands photon reso in CMS UPG cards
  • Property mode set to 100644
File size: 23.5 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# ECAL
227#############
228
229module SimpleCalorimeter ECal {
230 set ParticleInputArray ParticlePropagator/stableParticles
231 set TrackInputArray TrackMerger/tracks
232
233 set TowerOutputArray ecalTowers
234 set EFlowTrackOutputArray eflowTracks
235 set EFlowTowerOutputArray eflowPhotons
236
237 set IsEcal true
238
239 set EnergyMin 0.5
240 set EnergySignificanceMin 1.0
241
242 set SmearTowerCenter true
243
244 set pi [expr {acos(-1)}]
245
246 # lists of the edges of each tower in eta and phi
247 # each list starts with the lower edge of the first tower
248 # the list ends with the higher edged of the last tower
249
250 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
251
252 set PhiBins {}
253 for {set i -180} {$i <= 180} {incr i} {
254 add PhiBins [expr {$i * $pi/180.0}]
255 }
256
257 # 0.02 unit in eta up to eta = 1.5 (barrel)
258 for {set i -85} {$i <= 86} {incr i} {
259 set eta [expr {$i * 0.0174}]
260 add EtaPhiBins $eta $PhiBins
261 }
262
263 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
264
265 set PhiBins {}
266 for {set i -180} {$i <= 180} {incr i} {
267 add PhiBins [expr {$i * $pi/180.0}]
268 }
269
270 # 0.02 unit in eta up to eta = 3
271 for {set i 1} {$i <= 84} {incr i} {
272 set eta [expr { -2.958 + $i * 0.0174}]
273 add EtaPhiBins $eta $PhiBins
274 }
275
276 for {set i 1} {$i <= 84} {incr i} {
277 set eta [expr { 1.4964 + $i * 0.0174}]
278 add EtaPhiBins $eta $PhiBins
279 }
280
281 # take present CMS granularity for HF
282
283 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
284 set PhiBins {}
285 for {set i -18} {$i <= 18} {incr i} {
286 add PhiBins [expr {$i * $pi/18.0}]
287 }
288
289 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} {
290 add EtaPhiBins $eta $PhiBins
291 }
292
293
294 add EnergyFraction {0} {0.0}
295 # energy fractions for e, gamma and pi0
296 add EnergyFraction {11} {1.0}
297 add EnergyFraction {22} {1.0}
298 add EnergyFraction {111} {1.0}
299 # energy fractions for muon, neutrinos and neutralinos
300 add EnergyFraction {12} {0.0}
301 add EnergyFraction {13} {0.0}
302 add EnergyFraction {14} {0.0}
303 add EnergyFraction {16} {0.0}
304 add EnergyFraction {1000022} {0.0}
305 add EnergyFraction {1000023} {0.0}
306 add EnergyFraction {1000025} {0.0}
307 add EnergyFraction {1000035} {0.0}
308 add EnergyFraction {1000045} {0.0}
309 # energy fractions for K0short and Lambda
310 add EnergyFraction {310} {0.3}
311 add EnergyFraction {3122} {0.3}
312
313 # set ResolutionFormula {resolution formula as a function of eta and energy}
314
315 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
316 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
317
318 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
319 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
320 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
321 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
322 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
323}
324
325#############
326# HCAL
327#############
328
329module SimpleCalorimeter HCal {
330 set ParticleInputArray ParticlePropagator/stableParticles
331 set TrackInputArray ECal/eflowTracks
332
333 set TowerOutputArray hcalTowers
334 set EFlowTrackOutputArray eflowTracks
335 set EFlowTowerOutputArray eflowNeutralHadrons
336
337 set IsEcal false
338
339 set EnergyMin 1.0
340 set EnergySignificanceMin 1.0
341
342 set SmearTowerCenter true
343
344 set pi [expr {acos(-1)}]
345
346 # lists of the edges of each tower in eta and phi
347 # each list starts with the lower edge of the first tower
348 # the list ends with the higher edged of the last tower
349
350 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
351
352 set PhiBins {}
353 for {set i -36} {$i <= 36} {incr i} {
354 add PhiBins [expr {$i * $pi/36.0}]
355 }
356 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} {
357 add EtaPhiBins $eta $PhiBins
358 }
359
360 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
361
362 set PhiBins {}
363 for {set i -180} {$i <= 180} {incr i} {
364 add PhiBins [expr {$i * $pi/180.0}]
365 }
366
367 # 0.02 unit in eta up to eta = 3
368 for {set i 1} {$i <= 84} {incr i} {
369 set eta [expr { -2.958 + $i * 0.0174}]
370 add EtaPhiBins $eta $PhiBins
371 }
372
373 for {set i 1} {$i <= 84} {incr i} {
374 set eta [expr { 1.4964 + $i * 0.0174}]
375 add EtaPhiBins $eta $PhiBins
376 }
377
378 # take present CMS granularity for HF
379
380 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
381 set PhiBins {}
382 for {set i -18} {$i <= 18} {incr i} {
383 add PhiBins [expr {$i * $pi/18.0}]
384 }
385
386 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} {
387 add EtaPhiBins $eta $PhiBins
388 }
389
390 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
391 add EnergyFraction {0} {1.0}
392 # energy fractions for e, gamma and pi0
393 add EnergyFraction {11} {0.0}
394 add EnergyFraction {22} {0.0}
395 add EnergyFraction {111} {0.0}
396 # energy fractions for muon, neutrinos and neutralinos
397 add EnergyFraction {12} {0.0}
398 add EnergyFraction {13} {0.0}
399 add EnergyFraction {14} {0.0}
400 add EnergyFraction {16} {0.0}
401 add EnergyFraction {1000022} {0.0}
402 add EnergyFraction {1000023} {0.0}
403 add EnergyFraction {1000025} {0.0}
404 add EnergyFraction {1000035} {0.0}
405 add EnergyFraction {1000045} {0.0}
406 # energy fractions for K0short and Lambda
407 add EnergyFraction {310} {0.7}
408 add EnergyFraction {3122} {0.7}
409
410# set ResolutionFormula {resolution formula as a function of eta and energy}
411 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
412 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
413 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
414
415}
416
417#################
418# Electron filter
419#################
420
421module PdgCodeFilter ElectronFilter {
422 set InputArray HCal/eflowTracks
423 set OutputArray electrons
424 set Invert true
425 add PdgCode {11}
426 add PdgCode {-11}
427}
428
429###################################################
430# Tower Merger (in case not using e-flow algorithm)
431###################################################
432
433module Merger TowerMerger {
434# add InputArray InputArray
435 add InputArray ECal/ecalTowers
436 add InputArray HCal/hcalTowers
437 set OutputArray towers
438}
439
440#################################
441# Energy resolution for electrons
442#################################
443
444module EnergySmearing PhotonEnergySmearing {
445 set InputArray ECal/eflowPhotons
446 set OutputArray eflowPhotons
447
448 # adding 1% extra photon smearing
449 set ResolutionFormula {energy*0.01}
450
451}
452
453
454####################
455# Energy flow merger
456####################
457
458module Merger EFlowMerger {
459# add InputArray InputArray
460 add InputArray HCal/eflowTracks
461 add InputArray PhotonEnergySmearing/eflowPhotons
462 add InputArray HCal/eflowNeutralHadrons
463 set OutputArray eflow
464}
465
466
467###################
468# Missing ET merger
469###################
470
471module Merger MissingET {
472# add InputArray InputArray
473 add InputArray EFlowMerger/eflow
474 set MomentumOutputArray momentum
475}
476
477
478
479##################
480# Scalar HT merger
481##################
482
483module Merger ScalarHT {
484# add InputArray InputArray
485 add InputArray EFlowMerger/eflow
486 set EnergyOutputArray energy
487}
488
489#################
490# Neutrino Filter
491#################
492
493module PdgCodeFilter NeutrinoFilter {
494
495 set InputArray Delphes/stableParticles
496 set OutputArray filteredParticles
497
498 set PTMin 0.0
499
500 add PdgCode {12}
501 add PdgCode {14}
502 add PdgCode {16}
503 add PdgCode {-12}
504 add PdgCode {-14}
505 add PdgCode {-16}
506
507}
508
509
510#####################
511# MC truth jet finder
512#####################
513
514module FastJetFinder GenJetFinder {
515 set InputArray NeutrinoFilter/filteredParticles
516
517 set OutputArray jets
518
519 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
520 set JetAlgorithm 6
521 set ParameterR 0.4
522
523 set JetPTMin 15.0
524}
525
526#########################
527# Gen Missing ET merger
528########################
529
530module Merger GenMissingET {
531
532# add InputArray InputArray
533 add InputArray NeutrinoFilter/filteredParticles
534 set MomentumOutputArray momentum
535}
536
537
538
539############
540# Jet finder
541############
542
543module FastJetFinder FastJetFinder {
544# set InputArray TowerMerger/towers
545 set InputArray EFlowMerger/eflow
546
547 set OutputArray jets
548
549 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
550 set JetAlgorithm 6
551 set ParameterR 0.4
552
553 set JetPTMin 15.0
554}
555
556
557
558
559############
560# Jet finder
561############
562
563module FastJetFinder FatJetFinder {
564# set InputArray TowerMerger/towers
565 set InputArray EFlowMerger/eflow
566
567 set OutputArray jets
568
569 set JetAlgorithm 5
570 set ParameterR 0.8
571
572 set ComputeNsubjettiness 1
573 set Beta 1.0
574 set AxisMode 4
575
576 set ComputeTrimming 1
577 set RTrim 0.2
578 set PtFracTrim 0.05
579
580 set ComputePruning 1
581 set ZcutPrun 0.1
582 set RcutPrun 0.5
583 set RPrun 0.8
584
585 set ComputeSoftDrop 1
586 set BetaSoftDrop 0.0
587 set SymmetryCutSoftDrop 0.1
588 set R0SoftDrop 0.8
589
590 set JetPTMin 200.0
591}
592
593
594
595##################
596# Jet Energy Scale
597##################
598
599module EnergyScale JetEnergyScale {
600 set InputArray FastJetFinder/jets
601 set OutputArray jets
602
603 # scale formula for jets
604 set ScaleFormula {1.00}
605}
606
607####################
608# Photon isolation #
609####################
610
611module Isolation PhotonIsolation {
612
613 # particle for which calculate the isolation
614 set CandidateInputArray PhotonEnergySmearing/eflowPhotons
615
616 # isolation collection
617 set IsolationInputArray EFlowMerger/eflow
618
619 # output array
620 set OutputArray photons
621
622 # isolation cone
623 set DeltaRMax 0.3
624
625 # minimum pT
626 set PTMin 1.0
627
628 # iso ratio to cut
629 set PTRatioMax 9999.
630}
631
632
633
634#####################
635# Photon efficiency #
636#####################
637
638module Efficiency PhotonEfficiency {
639
640 ## input particles
641 set InputArray PhotonIsolation/photons
642 ## output particles
643 set OutputArray photons
644 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
645 # efficiency formula for photons
646 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
647 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
648 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
649 (abs(eta) > 4.0) * (0.00)}
650
651}
652
653
654######################
655# Electron isolation #
656######################
657
658module Isolation ElectronIsolation {
659
660 set CandidateInputArray ElectronFilter/electrons
661
662 # isolation collection
663 set IsolationInputArray EFlowMerger/eflow
664
665 set OutputArray electrons
666
667 set DeltaRMax 0.3
668 set PTMin 1.0
669 set PTRatioMax 9999.
670
671}
672
673
674
675#######################
676# Electron efficiency #
677#######################
678
679module Efficiency ElectronEfficiency {
680
681 set InputArray ElectronIsolation/electrons
682 set OutputArray electrons
683
684 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
685 # efficiency formula for electrons
686 set EfficiencyFormula {
687 (pt <= 4.0) * (0.00) + \
688 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
689 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
690 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
691 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
692 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
693 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
694 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
695 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
696 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
697 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
698 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
699 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
700 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
701 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
702 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
703 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
704 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
705 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
706 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
707 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
708 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
709 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
710 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
711 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
712 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
713 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
714 (abs(eta) > 4.0) * (0.00)
715
716 }
717}
718
719##################
720# Muon isolation #
721##################
722
723module Isolation MuonIsolation {
724 set CandidateInputArray MuonMomentumSmearing/muons
725
726 # isolation collection
727 set IsolationInputArray EFlowMerger/eflow
728
729 set OutputArray muons
730
731 set DeltaRMax 0.3
732 set PTMin 1.0
733 set PTRatioMax 9999.
734
735}
736
737##################
738# Muon Loose Id #
739##################
740
741module Efficiency MuonLooseIdEfficiency {
742 set InputArray MuonIsolation/muons
743 set OutputArray muons
744 # tracking + LooseID efficiency formula for muons
745 source muonLooseId.tcl
746
747}
748
749
750##################
751# Muon Tight Id #
752##################
753
754module Efficiency MuonTightIdEfficiency {
755 set InputArray MuonIsolation/muons
756 set OutputArray muons
757 # tracking + TightID efficiency formula for muons
758 source muonTightId.tcl
759}
760
761########################
762# Jet Flavor Association
763########################
764
765module JetFlavorAssociation JetFlavorAssociation {
766
767 set PartonInputArray Delphes/partons
768 set ParticleInputArray Delphes/allParticles
769 set ParticleLHEFInputArray Delphes/allParticlesLHEF
770 set JetInputArray JetEnergyScale/jets
771
772 set DeltaR 0.5
773 set PartonPTMin 10.0
774 set PartonEtaMax 4.0
775
776}
777
778
779#############
780# b-tagging #
781#############
782module BTagging BTaggingLoose {
783
784 set JetInputArray JetEnergyScale/jets
785
786 set BitNumber 0
787
788 source btagLoose.tcl
789}
790
791
792#############
793# b-tagging #
794#############
795module BTagging BTaggingMedium {
796
797 set JetInputArray JetEnergyScale/jets
798
799 set BitNumber 1
800
801 source btagMedium.tcl
802}
803
804
805
806#############
807# b-tagging #
808#############
809module BTagging BTaggingTight {
810
811 set JetInputArray JetEnergyScale/jets
812
813 set BitNumber 2
814
815 source btagTight.tcl
816}
817
818
819
820#############
821# tau-tagging
822#############
823
824
825module TauTagging TauTagging {
826 set ParticleInputArray Delphes/allParticles
827 set PartonInputArray Delphes/partons
828 set JetInputArray JetEnergyScale/jets
829
830 set DeltaR 0.5
831
832 set TauPTMin 20.0
833
834 set TauEtaMax 2.3
835
836 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
837
838 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)) + \
839 (abs(eta) > 2.3) * (0.000)
840 }
841 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) ) + \
842 (abs(eta) > 2.3) * (0.000)
843 }
844}
845
846
847##################
848# ROOT tree writer
849##################
850
851module TreeWriter TreeWriter {
852# add Branch InputArray BranchName BranchClass
853 add Branch Delphes/allParticles Particle GenParticle
854
855 add Branch GenJetFinder/jets GenJet Jet
856 add Branch GenMissingET/momentum GenMissingET MissingET
857
858 add Branch HCal/eflowTracks EFlowTrack Track
859 add Branch PhotonEnergySmearing/eflowPhotons EFlowPhoton Tower
860 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
861
862 add Branch PhotonEfficiency/photons Photon Photon
863 add Branch ElectronEfficiency/electrons Electron Electron
864 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
865 add Branch MuonTightIdEfficiency/muons MuonTight Muon
866
867 add Branch JetEnergyScale/jets Jet Jet
868 add Branch FatJetFinder/jets FatJet Jet
869
870 add Branch MissingET/momentum MissingET MissingET
871 add Branch ScalarHT/energy ScalarHT ScalarHT
872}
873
Note: See TracBrowser for help on using the repository browser.