Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_0PU.tcl@ 6e91f77

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 6e91f77 was f7f33af, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

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