Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl@ 4de073b

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

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