Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl@ 21eab4f

Last change on this file since 21eab4f was 1270747, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

remove maxEvent flag
  • Property mode set to 100644
File size: 39.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
17 PileUpMerger
18 ParticlePropagator
19
20 ChargedHadronTrackingEfficiency
21 ElectronTrackingEfficiency
22 MuonTrackingEfficiency
23
24 ChargedHadronMomentumSmearing
25 ElectronEnergySmearing
26 MuonMomentumSmearing
27
28 TrackMerger
29
30 ECal
31 HCal
32
33 PhotonEnergySmearing
34 ElectronFilter
35
36 TrackPileUpSubtractor
37 RecoPuFilter
38
39 TowerMerger
40 NeutralEFlowMerger
41
42 EFlowMerger
43 EFlowMergerCHS
44 Rho
45
46 LeptonFilterNoLep
47 LeptonFilterLep
48 RunPUPPIBase
49 RunPUPPI
50
51 PhotonFilter
52
53 PhotonIsolation
54 PhotonIsolationCHS
55 PhotonEfficiency
56 PhotonEfficiencyCHS
57
58 ElectronIsolation
59 ElectronIsolationCHS
60
61 ElectronEfficiency
62 ElectronEfficiencyCHS
63
64 MuonIsolation
65 MuonIsolationCHS
66
67 MuonLooseIdEfficiency
68 MuonTightIdEfficiency
69
70 MuonLooseIdEfficiencyCHS
71 MuonTightIdEfficiencyCHS
72
73 NeutrinoFilter
74
75 MissingET
76 PuppiMissingET
77 GenMissingET
78 GenPileUpMissingET
79
80 GenJetFinder
81 GenJetFinderAK8
82 FastJetFinder
83 FastJetFinderAK8
84 JetPileUpSubtractor
85 JetPileUpSubtractorAK8
86 FastJetFinderPUPPI
87 FastJetFinderPUPPIAK8
88
89 ScalarHT
90
91 JetEnergyScale
92 JetEnergyScaleAK8
93 JetEnergyScalePUPPI
94 JetEnergyScalePUPPIAK8
95
96 JetFlavorAssociation
97 JetFlavorAssociationAK8
98 JetFlavorAssociationPUPPI
99 JetFlavorAssociationPUPPIAK8
100
101 BTaggingLoose
102 BTaggingMedium
103 BTaggingTight
104 BTaggingLooseAK8
105 BTaggingMediumAK8
106 BTaggingTightAK8
107 BTaggingLoosePUPPI
108 BTaggingMediumPUPPI
109 BTaggingTightPUPPI
110 BTaggingLoosePUPPIAK8
111 BTaggingMediumPUPPIAK8
112 BTaggingTightPUPPIAK8
113
114 TauTagging
115
116 GenParticleFilter
117
118 TreeWriter
119}
120
121###############
122# PileUp Merger
123###############
124
125module PileUpMerger PileUpMerger {
126 set InputArray Delphes/stableParticles
127
128 set ParticleOutputArray stableParticles
129 set VertexOutputArray vertices
130
131 # pre-generated minbias input file
132 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
133 #set PileUpFile MinBias.pileup
134
135 # average expected pile up
136 set MeanPileUp 0
137
138 # maximum spread in the beam direction in m
139 set ZVertexSpread 0.25
140
141 # maximum spread in time in s
142 set TVertexSpread 800E-12
143
144 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) - {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
145 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
146
147}
148
149
150#################################
151# Propagate particles in cylinder
152#################################
153
154module ParticlePropagator ParticlePropagator {
155 #set InputArray PileUpMerger/stableParticles
156 set InputArray Delphes/stableParticles
157
158 set OutputArray stableParticles
159 set ChargedHadronOutputArray chargedHadrons
160 set ElectronOutputArray electrons
161 set MuonOutputArray muons
162
163 # radius of the magnetic field coverage, in m
164 set Radius 1.29
165 # half-length of the magnetic field coverage, in m
166 set HalfLength 3.0
167
168 # magnetic field
169 set Bz 3.8
170}
171
172
173####################################
174# Charged hadron tracking efficiency
175####################################
176
177module Efficiency ChargedHadronTrackingEfficiency {
178 ## particles after propagation
179 set InputArray ParticlePropagator/chargedHadrons
180 set OutputArray chargedHadrons
181 # tracking efficiency formula for charged hadrons
182 set EfficiencyFormula {
183 (pt <= 0.2) * (0.00) + \
184 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
185 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
186 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
187 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
188 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
189 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
190 (abs(eta) > 4.0) * (0.00)
191 }
192}
193
194
195#####################################
196# Electron tracking efficiency - ID
197####################################
198
199module Efficiency ElectronTrackingEfficiency {
200 set InputArray ParticlePropagator/electrons
201 set OutputArray electrons
202 # tracking efficiency formula for electrons
203 set EfficiencyFormula {
204 (pt <= 0.2) * (0.00) + \
205 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
206 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
207 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
208 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
209 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
210 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
211 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
212 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
213 (abs(eta) > 4.0) * (0.00)
214
215 }
216}
217
218##########################
219# Muon tracking efficiency
220##########################
221
222module Efficiency MuonTrackingEfficiency {
223 set InputArray ParticlePropagator/muons
224 set OutputArray muons
225 # tracking efficiency formula for muons
226 set EfficiencyFormula {
227 (pt <= 0.2) * (0.00) + \
228 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
229 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
230 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
231 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
232 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
233 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
234 (abs(eta) > 4.0) * (0.00)
235
236 }
237}
238
239
240########################################
241# Momentum resolution for charged tracks
242########################################
243
244module MomentumSmearing ChargedHadronMomentumSmearing {
245 ## hadrons after having applied the tracking efficiency
246 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
247 set OutputArray chargedHadrons
248 # resolution formula for charged hadrons ,
249
250 source trackMomentumResolution_vs_p_PIX4022.tcl
251}
252
253#################################
254# Energy resolution for electrons
255#################################
256
257module EnergySmearing ElectronEnergySmearing {
258 set InputArray ElectronTrackingEfficiency/electrons
259 set OutputArray electrons
260
261 # set ResolutionFormula {resolution formula as a function of eta and energy}
262
263 # resolution formula for electrons
264
265 # taking something flat in energy for now, ECAL will take over at high energy anyway.
266 # inferred from hep-ex/1306.2016 and 1502.02701
267 set ResolutionFormula {
268
269 (abs(eta) <= 1.5) * (energy*0.028) +
270 (abs(eta) > 1.5 && abs(eta) <= 1.75) * (energy*0.037) +
271 (abs(eta) > 1.75 && abs(eta) <= 2.15) * (energy*0.038) +
272 (abs(eta) > 2.15 && abs(eta) <= 3.00) * (energy*0.044) +
273 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)}
274
275}
276
277###############################
278# Momentum resolution for muons
279###############################
280
281module MomentumSmearing MuonMomentumSmearing {
282 set InputArray MuonTrackingEfficiency/muons
283 set OutputArray muons
284 # resolution formula for muons
285
286 # up to |eta| < 2.8 take measurement from tracking + muon chambers
287 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
288 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
289 source muonMomentumResolution.tcl
290}
291
292##############
293# Track merger
294##############
295
296module Merger TrackMerger {
297# add InputArray InputArray
298 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
299 add InputArray ElectronEnergySmearing/electrons
300 add InputArray MuonMomentumSmearing/muons
301 set OutputArray tracks
302}
303
304#############
305# ECAL
306#############
307
308module SimpleCalorimeter ECal {
309 set ParticleInputArray ParticlePropagator/stableParticles
310 set TrackInputArray TrackMerger/tracks
311
312 set TowerOutputArray ecalTowers
313 set EFlowTrackOutputArray eflowTracks
314 set EFlowTowerOutputArray eflowPhotons
315
316 set IsEcal true
317
318 set EnergyMin 0.5
319 set EnergySignificanceMin 1.0
320
321 set SmearTowerCenter true
322
323 set pi [expr {acos(-1)}]
324
325 # lists of the edges of each tower in eta and phi
326 # each list starts with the lower edge of the first tower
327 # the list ends with the higher edged of the last tower
328
329 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
330
331 set PhiBins {}
332 for {set i -180} {$i <= 180} {incr i} {
333 add PhiBins [expr {$i * $pi/180.0}]
334 }
335
336 # 0.02 unit in eta up to eta = 1.5 (barrel)
337 for {set i -85} {$i <= 86} {incr i} {
338 set eta [expr {$i * 0.0174}]
339 add EtaPhiBins $eta $PhiBins
340 }
341
342 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
343
344 set PhiBins {}
345 for {set i -180} {$i <= 180} {incr i} {
346 add PhiBins [expr {$i * $pi/180.0}]
347 }
348
349 # 0.02 unit in eta up to eta = 3
350 for {set i 1} {$i <= 84} {incr i} {
351 set eta [expr { -2.958 + $i * 0.0174}]
352 add EtaPhiBins $eta $PhiBins
353 }
354
355 for {set i 1} {$i <= 84} {incr i} {
356 set eta [expr { 1.4964 + $i * 0.0174}]
357 add EtaPhiBins $eta $PhiBins
358 }
359
360 # take present CMS granularity for HF
361
362 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
363 set PhiBins {}
364 for {set i -18} {$i <= 18} {incr i} {
365 add PhiBins [expr {$i * $pi/18.0}]
366 }
367
368 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} {
369 add EtaPhiBins $eta $PhiBins
370 }
371
372
373 add EnergyFraction {0} {0.0}
374 # energy fractions for e, gamma and pi0
375 add EnergyFraction {11} {1.0}
376 add EnergyFraction {22} {1.0}
377 add EnergyFraction {111} {1.0}
378 # energy fractions for muon, neutrinos and neutralinos
379 add EnergyFraction {12} {0.0}
380 add EnergyFraction {13} {0.0}
381 add EnergyFraction {14} {0.0}
382 add EnergyFraction {16} {0.0}
383 add EnergyFraction {1000022} {0.0}
384 add EnergyFraction {1000023} {0.0}
385 add EnergyFraction {1000025} {0.0}
386 add EnergyFraction {1000035} {0.0}
387 add EnergyFraction {1000045} {0.0}
388 # energy fractions for K0short and Lambda
389 add EnergyFraction {310} {0.3}
390 add EnergyFraction {3122} {0.3}
391
392 # set ResolutionFormula {resolution formula as a function of eta and energy}
393
394 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
395 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
396
397 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
398 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
399 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
400 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
401 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
402
403}
404
405#############
406# HCAL
407#############
408
409module SimpleCalorimeter HCal {
410 set ParticleInputArray ParticlePropagator/stableParticles
411 set TrackInputArray ECal/eflowTracks
412
413 set TowerOutputArray hcalTowers
414 set EFlowTrackOutputArray eflowTracks
415 set EFlowTowerOutputArray eflowNeutralHadrons
416
417 set IsEcal false
418
419 set EnergyMin 1.0
420 set EnergySignificanceMin 1.0
421
422 set SmearTowerCenter true
423
424 set pi [expr {acos(-1)}]
425
426 # lists of the edges of each tower in eta and phi
427 # each list starts with the lower edge of the first tower
428 # the list ends with the higher edged of the last tower
429
430 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
431
432 set PhiBins {}
433 for {set i -36} {$i <= 36} {incr i} {
434 add PhiBins [expr {$i * $pi/36.0}]
435 }
436 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} {
437 add EtaPhiBins $eta $PhiBins
438 }
439
440 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
441
442 set PhiBins {}
443 for {set i -180} {$i <= 180} {incr i} {
444 add PhiBins [expr {$i * $pi/180.0}]
445 }
446
447 # 0.02 unit in eta up to eta = 3
448 for {set i 1} {$i <= 84} {incr i} {
449 set eta [expr { -2.958 + $i * 0.0174}]
450 add EtaPhiBins $eta $PhiBins
451 }
452
453 for {set i 1} {$i <= 84} {incr i} {
454 set eta [expr { 1.4964 + $i * 0.0174}]
455 add EtaPhiBins $eta $PhiBins
456 }
457
458 # take present CMS granularity for HF
459
460 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
461 set PhiBins {}
462 for {set i -18} {$i <= 18} {incr i} {
463 add PhiBins [expr {$i * $pi/18.0}]
464 }
465
466 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} {
467 add EtaPhiBins $eta $PhiBins
468 }
469
470
471 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
472 add EnergyFraction {0} {1.0}
473 # energy fractions for e, gamma and pi0
474 add EnergyFraction {11} {0.0}
475 add EnergyFraction {22} {0.0}
476 add EnergyFraction {111} {0.0}
477 # energy fractions for muon, neutrinos and neutralinos
478 add EnergyFraction {12} {0.0}
479 add EnergyFraction {13} {0.0}
480 add EnergyFraction {14} {0.0}
481 add EnergyFraction {16} {0.0}
482 add EnergyFraction {1000022} {0.0}
483 add EnergyFraction {1000023} {0.0}
484 add EnergyFraction {1000025} {0.0}
485 add EnergyFraction {1000035} {0.0}
486 add EnergyFraction {1000045} {0.0}
487 # energy fractions for K0short and Lambda
488 add EnergyFraction {310} {0.7}
489 add EnergyFraction {3122} {0.7}
490
491# set ResolutionFormula {resolution formula as a function of eta and energy}
492 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
493 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
494 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
495
496}
497
498#################################
499# Energy resolution for electrons
500#################################
501
502module EnergySmearing PhotonEnergySmearing {
503 set InputArray ECal/eflowPhotons
504 set OutputArray eflowPhotons
505
506 # adding 1% extra photon smearing
507 set ResolutionFormula {energy*0.01}
508
509}
510
511
512
513#################
514# Electron filter
515#################
516
517module PdgCodeFilter ElectronFilter {
518 set InputArray HCal/eflowTracks
519 set OutputArray electrons
520 set Invert true
521 add PdgCode {11}
522 add PdgCode {-11}
523}
524
525
526##########################
527# Track pile-up subtractor
528##########################
529
530module TrackPileUpSubtractor TrackPileUpSubtractor {
531# add InputArray InputArray OutputArray
532 add InputArray HCal/eflowTracks eflowTracks
533 add InputArray ElectronFilter/electrons electrons
534 add InputArray MuonMomentumSmearing/muons muons
535
536 set VertexInputArray PileUpMerger/vertices
537 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
538 # Z vertex resolution in m
539 set ZVertexResolution 0.0001
540}
541
542########################
543# Reco PU filter
544########################
545
546module RecoPuFilter RecoPuFilter {
547 set InputArray HCal/eflowTracks
548 set OutputArray eflowTracks
549}
550
551###################################################
552# Tower Merger (in case not using e-flow algorithm)
553###################################################
554
555module Merger TowerMerger {
556# add InputArray InputArray
557 add InputArray ECal/ecalTowers
558 add InputArray HCal/hcalTowers
559 set OutputArray towers
560}
561
562####################
563# Neutral eflow erger
564####################
565
566module Merger NeutralEFlowMerger {
567# add InputArray InputArray
568 add InputArray PhotonEnergySmearing/eflowPhotons
569 add InputArray HCal/eflowNeutralHadrons
570 set OutputArray eflowTowers
571}
572
573#####################
574# Energy flow merger
575#####################
576
577module Merger EFlowMerger {
578# add InputArray InputArray
579 add InputArray HCal/eflowTracks
580 add InputArray PhotonEnergySmearing/eflowPhotons
581 add InputArray HCal/eflowNeutralHadrons
582 set OutputArray eflow
583}
584
585############################
586# Energy flow merger no PU
587############################
588
589module Merger EFlowMergerCHS {
590# add InputArray InputArray
591 add InputArray RecoPuFilter/eflowTracks
592 add InputArray PhotonEnergySmearing/eflowPhotons
593 add InputArray HCal/eflowNeutralHadrons
594 set OutputArray eflow
595}
596
597#########################################
598### Run the puppi code (to be tuned) ###
599#########################################
600
601module PdgCodeFilter LeptonFilterNoLep {
602 set InputArray HCal/eflowTracks
603 set OutputArray eflowTracksNoLeptons
604 set Invert false
605 add PdgCode {13}
606 add PdgCode {-13}
607 add PdgCode {11}
608 add PdgCode {-11}
609}
610
611module PdgCodeFilter LeptonFilterLep {
612 set InputArray HCal/eflowTracks
613 set OutputArray eflowTracksLeptons
614 set Invert true
615 add PdgCode {11}
616 add PdgCode {-11}
617 add PdgCode {13}
618 add PdgCode {-13}
619}
620
621module RunPUPPI RunPUPPIBase {
622 ## input information
623 set TrackInputArray LeptonFilterNoLep/eflowTracksNoLeptons
624 set NeutralInputArray NeutralEFlowMerger/eflowTowers
625 set PVInputArray PileUpMerger/vertices
626 set MinPuppiWeight 0.05
627 set UseExp false
628 set UseNoLep false
629
630 ## define puppi algorithm parameters (more than one for the same eta region is possible)
631 add EtaMinBin 0.0 1.5 4.0
632 add EtaMaxBin 1.5 4.0 10.0
633 add PtMinBin 0.0 0.0 0.0
634 add ConeSizeBin 0.2 0.2 0.2
635 add RMSPtMinBin 0.1 0.5 0.5
636 add RMSScaleFactorBin 1.0 1.0 1.0
637 add NeutralMinEBin 0.2 0.2 0.5
638 add NeutralPtSlope 0.006 0.013 0.067
639 add ApplyCHS true true true
640 add UseCharged true true false
641 add ApplyLowPUCorr true true true
642 add MetricId 5 5 5
643 add CombId 0 0 0
644
645 ## output name
646 set OutputArray PuppiParticles
647 set OutputArrayTracks puppiTracks
648 set OutputArrayNeutrals puppiNeutrals
649}
650
651module Merger RunPUPPI {
652 add InputArray RunPUPPIBase/PuppiParticles
653 add InputArray LeptonFilterLep/eflowTracksLeptons
654 set OutputArray PuppiParticles
655}
656
657###################
658# Missing ET merger
659###################
660
661module Merger MissingET {
662# add InputArray InputArray
663# add InputArray RunPUPPI/PuppiParticles
664 add InputArray EFlowMerger/eflow
665 set MomentumOutputArray momentum
666}
667
668module Merger PuppiMissingET {
669 #add InputArray InputArray
670 add InputArray RunPUPPI/PuppiParticles
671 #add InputArray EFlowMerger/eflow
672 set MomentumOutputArray momentum
673}
674
675###################
676# Ger PileUp Missing ET
677###################
678
679module Merger GenPileUpMissingET {
680# add InputArray InputArray
681# add InputArray RunPUPPI/PuppiParticles
682 add InputArray ParticlePropagator/stableParticles
683 set MomentumOutputArray momentum
684}
685
686##################
687# Scalar HT merger
688##################
689
690module Merger ScalarHT {
691# add InputArray InputArray
692 add InputArray RunPUPPI/PuppiParticles
693 set EnergyOutputArray energy
694}
695
696#################
697# Neutrino Filter
698#################
699
700module PdgCodeFilter NeutrinoFilter {
701
702 set InputArray Delphes/stableParticles
703 set OutputArray filteredParticles
704
705 set PTMin 0.0
706
707 add PdgCode {12}
708 add PdgCode {14}
709 add PdgCode {16}
710 add PdgCode {-12}
711 add PdgCode {-14}
712 add PdgCode {-16}
713
714}
715
716#####################
717# MC truth jet finder
718#####################
719
720module FastJetFinder GenJetFinder {
721 set InputArray NeutrinoFilter/filteredParticles
722
723 set OutputArray jets
724
725 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
726 set JetAlgorithm 6
727 set ParameterR 0.4
728
729 set JetPTMin 15.0
730}
731
732module FastJetFinder GenJetFinderAK8 {
733 set InputArray NeutrinoFilter/filteredParticles
734
735 set OutputArray jetsAK8
736
737 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
738 set JetAlgorithm 6
739 set ParameterR 0.8
740
741 set JetPTMin 200.0
742}
743
744#########################
745# Gen Missing ET merger
746########################
747
748module Merger GenMissingET {
749
750# add InputArray InputArray
751 add InputArray NeutrinoFilter/filteredParticles
752 set MomentumOutputArray momentum
753}
754
755
756#############
757# Rho pile-up
758#############
759
760module FastJetFinder Rho {
761# set InputArray Calorimeter/towers
762 set InputArray EFlowMergerCHS/eflow
763
764 set ComputeRho true
765 set RhoOutputArray rho
766
767 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
768 set AreaAlgorithm 5
769
770 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
771 set JetAlgorithm 4
772 set ParameterR 0.4
773 set GhostEtaMax 5.0
774
775 add RhoEtaRange -5.0 -4.0
776 add RhoEtaRange -4.0 -1.5
777 add RhoEtaRange -1.5 1.5
778 add RhoEtaRange 1.5 4.0
779 add RhoEtaRange 4.0 5.0
780
781 set JetPTMin 0.0
782}
783
784
785##############
786# Jet finder
787##############
788
789module FastJetFinder FastJetFinder {
790# set InputArray TowerMerger/towers
791 set InputArray EFlowMergerCHS/eflow
792
793 set OutputArray jets
794
795 set AreaAlgorithm 5
796
797 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
798 set JetAlgorithm 6
799 set ParameterR 0.4
800
801 set JetPTMin 15.0
802}
803
804#module Class Name
805module FastJetFinder FastJetFinderAK8 {
806# set InputArray TowerMerger/towers
807 set InputArray EFlowMergerCHS/eflow
808
809 set OutputArray jets
810
811 set AreaAlgorithm 5
812
813 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
814 set JetAlgorithm 6
815 set ParameterR 0.8
816
817 set ComputeNsubjettiness 1
818 set Beta 1.0
819 set AxisMode 4
820
821 set ComputeTrimming 1
822 set RTrim 0.2
823 set PtFracTrim 0.05
824
825 set ComputePruning 1
826 set ZcutPrun 0.1
827 set RcutPrun 0.5
828 set RPrun 0.8
829
830 set ComputeSoftDrop 1
831 set BetaSoftDrop 0.0
832 set SymmetryCutSoftDrop 0.1
833 set R0SoftDrop 0.8
834
835 set JetPTMin 200.0
836}
837
838###########################
839# Jet Pile-Up Subtraction
840###########################
841
842module JetPileUpSubtractor JetPileUpSubtractor {
843 set JetInputArray FastJetFinder/jets
844 set RhoInputArray Rho/rho
845
846 set OutputArray jets
847
848 set JetPTMin 15.0
849}
850
851##############################
852# Jet Pile-Up Subtraction AK8
853##############################
854
855module JetPileUpSubtractor JetPileUpSubtractorAK8 {
856 set JetInputArray FastJetFinderAK8/jets
857 set RhoInputArray Rho/rho
858
859 set OutputArray jets
860
861 set JetPTMin 15.0
862}
863
864module FastJetFinder FastJetFinderPUPPI {
865# set InputArray TowerMerger/towers
866 set InputArray RunPUPPI/PuppiParticles
867
868 set OutputArray jets
869
870 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
871 set JetAlgorithm 6
872 set ParameterR 0.4
873
874 set JetPTMin 15.0
875}
876
877
878module FastJetFinder FastJetFinderPUPPIAK8 {
879# set InputArray TowerMerger/towers
880 set InputArray RunPUPPI/PuppiParticles
881
882 set OutputArray jets
883
884 set JetAlgorithm 6
885 set ParameterR 0.8
886
887 set ComputeNsubjettiness 1
888 set Beta 1.0
889 set AxisMode 4
890
891 set ComputeTrimming 1
892 set RTrim 0.2
893 set PtFracTrim 0.05
894
895 set ComputePruning 1
896 set ZcutPrun 0.1
897 set RcutPrun 0.5
898 set RPrun 0.8
899
900 set ComputeSoftDrop 1
901 set BetaSoftDrop 0.0
902 set SymmetryCutSoftDrop 0.1
903 set R0SoftDrop 0.8
904
905 set JetPTMin 200.0
906}
907
908##################
909# Jet Energy Scale
910##################
911
912module EnergyScale JetEnergyScale {
913 set InputArray JetPileUpSubtractor/jets
914 set OutputArray jets
915
916 # scale formula for jets
917 set ScaleFormula {1.00}
918}
919
920module EnergyScale JetEnergyScaleAK8 {
921 set InputArray JetPileUpSubtractorAK8/jets
922 set OutputArray jets
923
924 # scale formula for jets
925 set ScaleFormula {1.00}
926}
927
928module EnergyScale JetEnergyScalePUPPI {
929 set InputArray FastJetFinderPUPPI/jets
930 set OutputArray jets
931
932 # scale formula for jets
933 set ScaleFormula {1.00}
934}
935
936module EnergyScale JetEnergyScalePUPPIAK8 {
937 set InputArray FastJetFinderPUPPIAK8/jets
938 set OutputArray jets
939
940 # scale formula for jets
941 set ScaleFormula {1.00}
942}
943
944#################
945# Photon filter
946#################
947
948module PdgCodeFilter PhotonFilter {
949 set InputArray PhotonEnergySmearing/eflowPhotons
950 set OutputArray photons
951 set Invert true
952 set PTMin 5.0
953 add PdgCode {22}
954}
955
956
957####################
958# Photon isolation #
959####################
960
961module Isolation PhotonIsolation {
962
963 # particle for which calculate the isolation
964 set CandidateInputArray PhotonFilter/photons
965
966 # isolation collection
967 set IsolationInputArray RunPUPPI/PuppiParticles
968
969 # output array
970 set OutputArray photons
971
972 # veto isolation cand. based on proximity to input cand.
973 set DeltaRMin 0.01
974 set UseMiniCone true
975
976 # isolation cone
977 set DeltaRMax 0.3
978
979 # minimum pT
980 set PTMin 1.0
981
982 # iso ratio to cut
983 set PTRatioMax 9999.
984
985}
986
987
988########################
989# Photon isolation CHS #
990########################
991
992module Isolation PhotonIsolationCHS {
993
994 # particle for which calculate the isolation
995 set CandidateInputArray PhotonFilter/photons
996
997 # isolation collection
998 set IsolationInputArray EFlowMerger/eflow
999
1000 # output array
1001 set OutputArray photons
1002
1003 # isolation cone
1004 set DeltaRMax 0.3
1005
1006 # minimum pT
1007 set PTMin 1.0
1008
1009 # iso ratio to cut
1010 set PTRatioMax 9999.
1011
1012}
1013
1014
1015#####################
1016# Photon efficiency #
1017#####################
1018
1019module Efficiency PhotonEfficiency {
1020
1021 ## input particles
1022 set InputArray PhotonIsolation/photons
1023 ## output particles
1024 set OutputArray photons
1025 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1026 # efficiency formula for photons
1027 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1028 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1029 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1030 (abs(eta) > 4.0) * (0.00)}
1031
1032}
1033
1034
1035#####################
1036# Photon efficiency #
1037#####################
1038
1039module Efficiency PhotonEfficiencyCHS {
1040
1041 ## input particles
1042 set InputArray PhotonIsolationCHS/photons
1043 ## output particles
1044 set OutputArray photons
1045 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1046 # efficiency formula for photons
1047 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1048 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1049 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1050 (abs(eta) > 4.0) * (0.00)}
1051
1052}
1053
1054######################
1055# Electron isolation #
1056######################
1057
1058module Isolation ElectronIsolation {
1059
1060 set CandidateInputArray ElectronFilter/electrons
1061
1062 # isolation collection
1063 set IsolationInputArray RunPUPPI/PuppiParticles
1064 #set IsolationInputArray EFlowMerger/eflow
1065
1066 set OutputArray electrons
1067
1068 set DeltaRMax 0.3
1069 set PTMin 1.0
1070 set PTRatioMax 9999.
1071
1072}
1073
1074
1075##########################
1076# Electron isolation CHS #
1077##########################
1078
1079module Isolation ElectronIsolationCHS {
1080
1081 set CandidateInputArray ElectronFilter/electrons
1082
1083 # isolation collection
1084 set IsolationInputArray EFlowMerger/eflow
1085
1086 set OutputArray electrons
1087
1088 set DeltaRMax 0.3
1089 set PTMin 1.0
1090 set PTRatioMax 9999.
1091
1092}
1093
1094
1095#######################
1096# Electron efficiency #
1097#######################
1098
1099module Efficiency ElectronEfficiency {
1100
1101 set InputArray ElectronIsolation/electrons
1102 set OutputArray electrons
1103
1104 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1105 # efficiency formula for electrons
1106 set EfficiencyFormula {
1107 (pt <= 4.0) * (0.00) + \
1108 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1109 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1110 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1111 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1112 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1113 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1114 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1115 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1116 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1117 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1118 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1119 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1120 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1121 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1122 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1123 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1124 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1125 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1126 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1127 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1128 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1129 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1130 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1131 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1132 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1133 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1134 (abs(eta) > 4.0) * (0.00)
1135
1136 }
1137}
1138
1139###########################
1140# Electron efficiency CHS #
1141###########################
1142
1143module Efficiency ElectronEfficiencyCHS {
1144
1145 set InputArray ElectronIsolationCHS/electrons
1146 set OutputArray electrons
1147
1148 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1149 # efficiency formula for electrons
1150 set EfficiencyFormula {
1151 (pt <= 4.0) * (0.00) + \
1152 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1153 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1154 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1155 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1156 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1157 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1158 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1159 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1160 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1161 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1162 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1163 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1164 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1165 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1166 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1167 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1168 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1169 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1170 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1171 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1172 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1173 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1174 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1175 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1176 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1177 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1178 (abs(eta) > 4.0) * (0.00)
1179
1180 }
1181}
1182
1183
1184##################
1185# Muon isolation #
1186##################
1187
1188module Isolation MuonIsolation {
1189 set CandidateInputArray MuonMomentumSmearing/muons
1190
1191 # isolation collection
1192 set IsolationInputArray RunPUPPI/PuppiParticles
1193
1194 set OutputArray muons
1195
1196 set DeltaRMax 0.3
1197 set PTMin 1.0
1198 set PTRatioMax 9999.
1199
1200}
1201
1202######################
1203# Muon isolation CHS #
1204######################
1205
1206module Isolation MuonIsolationCHS {
1207 set CandidateInputArray MuonMomentumSmearing/muons
1208
1209 # isolation collection
1210 set IsolationInputArray EFlowMerger/eflow
1211
1212 set OutputArray muons
1213
1214 set DeltaRMax 0.3
1215 set PTMin 1.0
1216 set PTRatioMax 9999.
1217
1218}
1219
1220
1221#####################
1222# Muon Loose Id #
1223#####################
1224
1225module Efficiency MuonLooseIdEfficiency {
1226 set InputArray MuonIsolation/muons
1227 set OutputArray muons
1228 # tracking + LooseID efficiency formula for muons
1229 source muonLooseId.tcl
1230
1231}
1232
1233
1234##################
1235# Muon Tight Id #
1236##################
1237
1238module Efficiency MuonTightIdEfficiency {
1239 set InputArray MuonIsolation/muons
1240 set OutputArray muons
1241 # tracking + TightID efficiency formula for muons
1242 source muonTightId.tcl
1243}
1244
1245
1246#####################
1247# Muon Loose Id CHS #
1248#####################
1249
1250module Efficiency MuonLooseIdEfficiencyCHS {
1251 set InputArray MuonIsolationCHS/muons
1252 set OutputArray muons
1253 # tracking + LooseID efficiency formula for muons
1254 source muonLooseId.tcl
1255
1256}
1257
1258
1259######################
1260# Muon Tight Id CHS #
1261######################
1262
1263module Efficiency MuonTightIdEfficiencyCHS {
1264 set InputArray MuonIsolationCHS/muons
1265 set OutputArray muons
1266 # tracking + TightID efficiency formula for muons
1267 source muonTightId.tcl
1268}
1269
1270
1271########################
1272# Jet Flavor Association
1273########################
1274
1275module JetFlavorAssociation JetFlavorAssociation {
1276
1277 set PartonInputArray Delphes/partons
1278 set ParticleInputArray Delphes/allParticles
1279 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1280 set JetInputArray JetEnergyScale/jets
1281
1282 set DeltaR 0.5
1283 set PartonPTMin 10.0
1284 set PartonEtaMax 4.0
1285
1286}
1287
1288module JetFlavorAssociation JetFlavorAssociationAK8 {
1289
1290 set PartonInputArray Delphes/partons
1291 set ParticleInputArray Delphes/allParticles
1292 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1293 set JetInputArray JetEnergyScaleAK8/jets
1294
1295 set DeltaR 0.8
1296 set PartonPTMin 100.0
1297 set PartonEtaMax 4.0
1298
1299}
1300
1301module JetFlavorAssociation JetFlavorAssociationPUPPI {
1302
1303 set PartonInputArray Delphes/partons
1304 set ParticleInputArray Delphes/allParticles
1305 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1306 set JetInputArray JetEnergyScalePUPPI/jets
1307
1308 set DeltaR 0.5
1309 set PartonPTMin 10.0
1310 set PartonEtaMax 4.0
1311
1312}
1313
1314module JetFlavorAssociation JetFlavorAssociationPUPPIAK8 {
1315
1316 set PartonInputArray Delphes/partons
1317 set ParticleInputArray Delphes/allParticles
1318 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1319 set JetInputArray JetEnergyScalePUPPIAK8/jets
1320
1321 set DeltaR 0.8
1322 set PartonPTMin 100.0
1323 set PartonEtaMax 4.0
1324
1325}
1326
1327
1328#############
1329# b-tagging #
1330#############
1331module BTagging BTaggingLoose {
1332
1333 set JetInputArray JetEnergyScale/jets
1334
1335 set BitNumber 0
1336
1337 source btagLoose.tcl
1338}
1339
1340module BTagging BTaggingLooseAK8 {
1341
1342 set JetInputArray JetEnergyScaleAK8/jets
1343
1344 set BitNumber 0
1345
1346 source btagLoose.tcl
1347}
1348
1349module BTagging BTaggingLoosePUPPI {
1350
1351 set JetInputArray JetEnergyScalePUPPI/jets
1352
1353 set BitNumber 0
1354
1355 source btagLoose.tcl
1356}
1357
1358module BTagging BTaggingLoosePUPPIAK8 {
1359
1360 set JetInputArray JetEnergyScalePUPPIAK8/jets
1361
1362 set BitNumber 0
1363
1364 source btagLoose.tcl
1365}
1366
1367
1368#############
1369# b-tagging #
1370#############
1371module BTagging BTaggingMedium {
1372
1373 set JetInputArray JetEnergyScale/jets
1374
1375 set BitNumber 1
1376
1377 source btagMedium.tcl
1378}
1379
1380module BTagging BTaggingMediumAK8 {
1381
1382 set JetInputArray JetEnergyScaleAK8/jets
1383
1384 set BitNumber 1
1385
1386 source btagMedium.tcl
1387}
1388
1389module BTagging BTaggingMediumPUPPI {
1390
1391 set JetInputArray JetEnergyScalePUPPI/jets
1392
1393 set BitNumber 1
1394
1395 source btagMedium.tcl
1396}
1397
1398module BTagging BTaggingMediumPUPPIAK8 {
1399
1400 set JetInputArray JetEnergyScalePUPPIAK8/jets
1401
1402 set BitNumber 1
1403
1404 source btagMedium.tcl
1405}
1406
1407
1408
1409#############
1410# b-tagging #
1411#############
1412module BTagging BTaggingTight {
1413
1414 set JetInputArray JetEnergyScale/jets
1415
1416 set BitNumber 2
1417
1418 source btagTight.tcl
1419}
1420
1421module BTagging BTaggingTightAK8 {
1422
1423 set JetInputArray JetEnergyScaleAK8/jets
1424
1425 set BitNumber 2
1426
1427 source btagTight.tcl
1428}
1429
1430module BTagging BTaggingTightPUPPI {
1431
1432 set JetInputArray JetEnergyScalePUPPI/jets
1433
1434 set BitNumber 2
1435
1436 source btagTight.tcl
1437}
1438
1439module BTagging BTaggingTightPUPPIAK8 {
1440
1441 set JetInputArray JetEnergyScalePUPPIAK8/jets
1442
1443 set BitNumber 2
1444
1445 source btagTight.tcl
1446}
1447
1448#############
1449# tau-tagging
1450#############
1451
1452
1453module TauTagging TauTagging {
1454 set ParticleInputArray Delphes/allParticles
1455 set PartonInputArray Delphes/partons
1456 set JetInputArray JetEnergyScale/jets
1457
1458 set DeltaR 0.5
1459
1460 set TauPTMin 20.0
1461
1462 set TauEtaMax 2.3
1463
1464 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1465
1466 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)) + \
1467 (abs(eta) > 2.3) * (0.000)
1468 }
1469 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) ) + \
1470 (abs(eta) > 2.3) * (0.000)
1471 }
1472}
1473
1474
1475###############################################################################################################
1476# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1477###############################################################################################################
1478
1479module StatusPidFilter GenParticleFilter {
1480
1481 set InputArray Delphes/allParticles
1482 set OutputArray filteredParticles
1483 set PTMin 5.0
1484
1485}
1486
1487
1488##################
1489# ROOT tree writer
1490##################
1491
1492module TreeWriter TreeWriter {
1493# add Branch InputArray BranchName BranchClass
1494 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1495 add Branch PileUpMerger/vertices Vertex Vertex
1496
1497 add Branch GenJetFinder/jets GenJet Jet
1498 add Branch GenJetFinderAK8/jetsAK8 GenJetAK8 Jet
1499 add Branch GenMissingET/momentum GenMissingET MissingET
1500
1501# add Branch HCal/eflowTracks EFlowTrack Track
1502# add Branch ECal/eflowPhotons EFlowPhoton Tower
1503# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1504
1505 add Branch PhotonEfficiency/photons Photon Photon
1506 add Branch ElectronEfficiency/electrons Electron Electron
1507 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
1508 add Branch MuonTightIdEfficiency/muons MuonTight Muon
1509
1510 add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
1511 add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
1512 add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
1513 add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon
1514
1515 add Branch JetEnergyScale/jets Jet Jet
1516# add Branch FatJetFinder/jets FatJet Jet
1517 add Branch JetEnergyScalePUPPI/jets JetPUPPI Jet
1518 add Branch JetEnergyScaleAK8/jets JetAK8 Jet
1519 add Branch JetEnergyScalePUPPIAK8/jets JetPUPPIAK8 Jet
1520
1521 add Branch MissingET/momentum MissingET MissingET
1522 add Branch PuppiMissingET/momentum PuppiMissingET MissingET
1523 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1524 add Branch ScalarHT/energy ScalarHT ScalarHT
1525
1526}
Note: See TracBrowser for help on using the repository browser.