Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_200PU.tcl@ b1fd38d

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

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