Fork me on GitHub

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

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

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