Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_200PU.tcl@ 8cbe720

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

cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl
  • Property mode set to 100644
File size: 39.3 KB
Line 
1set MaxEvents 100
2
3#
4# Phase II - No Pile-Up
5#
6# Main authors: Michele Selvaggi (UCL)
7#
8# Released on:
9#
10# Version: v01
11#
12#
13#######################################
14# Order of execution of various modules
15#######################################
16
17set ExecutionPath {
18
19 PileUpMerger
20 ParticlePropagator
21
22 ChargedHadronTrackingEfficiency
23 ElectronTrackingEfficiency
24 MuonTrackingEfficiency
25
26 ChargedHadronMomentumSmearing
27 ElectronEnergySmearing
28 MuonMomentumSmearing
29
30 TrackMerger
31
32 ECal
33 HCal
34
35 PhotonEnergySmearing
36 ElectronFilter
37
38 TrackPileUpSubtractor
39 RecoPuFilter
40
41 TowerMerger
42 NeutralEFlowMerger
43
44 EFlowMerger
45 EFlowMergerCHS
46 Rho
47
48 LeptonFilterNoLep
49 LeptonFilterLep
50 RunPUPPIBase
51 RunPUPPI
52
53 PhotonFilter
54
55 PhotonIsolation
56 PhotonIsolationCHS
57 PhotonEfficiency
58 PhotonEfficiencyCHS
59
60 ElectronIsolation
61 ElectronIsolationCHS
62
63 ElectronEfficiency
64 ElectronEfficiencyCHS
65
66 MuonIsolation
67 MuonIsolationCHS
68
69 MuonLooseIdEfficiency
70 MuonTightIdEfficiency
71
72 MuonLooseIdEfficiencyCHS
73 MuonTightIdEfficiencyCHS
74
75 NeutrinoFilter
76
77 MissingET
78 PuppiMissingET
79 GenMissingET
80 GenPileUpMissingET
81
82 GenJetFinder
83 GenJetFinderAK8
84 FastJetFinder
85 FastJetFinderAK8
86 JetPileUpSubtractor
87 JetPileUpSubtractorAK8
88 FastJetFinderPUPPI
89 FastJetFinderPUPPIAK8
90
91 ScalarHT
92
93 JetEnergyScale
94 JetEnergyScaleAK8
95 JetEnergyScalePUPPI
96 JetEnergyScalePUPPIAK8
97
98 JetFlavorAssociation
99 JetFlavorAssociationAK8
100 JetFlavorAssociationPUPPI
101 JetFlavorAssociationPUPPIAK8
102
103 BTaggingLoose
104 BTaggingMedium
105 BTaggingTight
106 BTaggingLooseAK8
107 BTaggingMediumAK8
108 BTaggingTightAK8
109 BTaggingLoosePUPPI
110 BTaggingMediumPUPPI
111 BTaggingTightPUPPI
112 BTaggingLoosePUPPIAK8
113 BTaggingMediumPUPPIAK8
114 BTaggingTightPUPPIAK8
115
116 TauTagging
117
118 GenParticleFilter
119
120 TreeWriter
121}
122
123###############
124# PileUp Merger
125###############
126
127module PileUpMerger PileUpMerger {
128 set InputArray Delphes/stableParticles
129
130 set ParticleOutputArray stableParticles
131 set VertexOutputArray vertices
132
133 # pre-generated minbias input file
134 #set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
135 set PileUpFile MinBias.pileup
136
137 # average expected pile up
138 set MeanPileUp 200
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# Photon isolation #
961####################
962
963module Isolation PhotonIsolation {
964
965 # particle for which calculate the isolation
966 set CandidateInputArray PhotonFilter/photons
967
968 # isolation collection
969 set IsolationInputArray RunPUPPI/PuppiParticles
970
971 # output array
972 set OutputArray photons
973
974 # veto isolation cand. based on proximity to input cand.
975 set DeltaRMin 0.01
976 set UseMiniCone true
977
978 # isolation cone
979 set DeltaRMax 0.3
980
981 # minimum pT
982 set PTMin 1.0
983
984 # iso ratio to cut
985 set PTRatioMax 9999.
986
987}
988
989
990########################
991# Photon isolation CHS #
992########################
993
994module Isolation PhotonIsolationCHS {
995
996 # particle for which calculate the isolation
997 set CandidateInputArray PhotonFilter/photons
998
999 # isolation collection
1000 set IsolationInputArray EFlowMerger/eflow
1001
1002 # output array
1003 set OutputArray photons
1004
1005 # isolation cone
1006 set DeltaRMax 0.3
1007
1008 # minimum pT
1009 set PTMin 1.0
1010
1011 # iso ratio to cut
1012 set PTRatioMax 9999.
1013
1014}
1015
1016
1017#####################
1018# Photon efficiency #
1019#####################
1020
1021module Efficiency PhotonEfficiency {
1022
1023 ## input particles
1024 set InputArray PhotonIsolation/photons
1025 ## output particles
1026 set OutputArray photons
1027 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1028 # efficiency formula for photons
1029 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1030 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1031 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1032 (abs(eta) > 4.0) * (0.00)}
1033
1034}
1035
1036
1037#####################
1038# Photon efficiency #
1039#####################
1040
1041module Efficiency PhotonEfficiencyCHS {
1042
1043 ## input particles
1044 set InputArray PhotonIsolationCHS/photons
1045 ## output particles
1046 set OutputArray photons
1047 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1048 # efficiency formula for photons
1049 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1050 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1051 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1052 (abs(eta) > 4.0) * (0.00)}
1053
1054}
1055
1056######################
1057# Electron isolation #
1058######################
1059
1060module Isolation ElectronIsolation {
1061
1062 set CandidateInputArray ElectronFilter/electrons
1063
1064 # isolation collection
1065 set IsolationInputArray RunPUPPI/PuppiParticles
1066 #set IsolationInputArray EFlowMerger/eflow
1067
1068 set OutputArray electrons
1069
1070 set DeltaRMax 0.3
1071 set PTMin 1.0
1072 set PTRatioMax 9999.
1073
1074}
1075
1076
1077##########################
1078# Electron isolation CHS #
1079##########################
1080
1081module Isolation ElectronIsolationCHS {
1082
1083 set CandidateInputArray ElectronFilter/electrons
1084
1085 # isolation collection
1086 set IsolationInputArray EFlowMerger/eflow
1087
1088 set OutputArray electrons
1089
1090 set DeltaRMax 0.3
1091 set PTMin 1.0
1092 set PTRatioMax 9999.
1093
1094}
1095
1096
1097#######################
1098# Electron efficiency #
1099#######################
1100
1101module Efficiency ElectronEfficiency {
1102
1103 set InputArray ElectronIsolation/electrons
1104 set OutputArray electrons
1105
1106 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1107 # efficiency formula for electrons
1108 set EfficiencyFormula {
1109 (pt <= 4.0) * (0.00) + \
1110 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1111 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1112 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1113 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1114 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1115 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1116 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1117 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1118 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1119 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1120 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1121 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1122 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1123 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1124 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1125 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1126 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1127 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1128 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1129 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1130 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1131 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1132 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1133 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1134 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1135 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1136 (abs(eta) > 4.0) * (0.00)
1137
1138 }
1139}
1140
1141###########################
1142# Electron efficiency CHS #
1143###########################
1144
1145module Efficiency ElectronEfficiencyCHS {
1146
1147 set InputArray ElectronIsolationCHS/electrons
1148 set OutputArray electrons
1149
1150 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1151 # efficiency formula for electrons
1152 set EfficiencyFormula {
1153 (pt <= 4.0) * (0.00) + \
1154 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1155 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1156 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1157 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1158 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1159 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1160 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1161 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1162 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1163 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1164 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1165 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1166 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1167 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1168 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1169 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1170 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1171 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1172 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1173 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1174 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1175 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1176 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1177 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1178 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1179 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1180 (abs(eta) > 4.0) * (0.00)
1181
1182 }
1183}
1184
1185
1186##################
1187# Muon isolation #
1188##################
1189
1190module Isolation MuonIsolation {
1191 set CandidateInputArray MuonMomentumSmearing/muons
1192
1193 # isolation collection
1194 set IsolationInputArray RunPUPPI/PuppiParticles
1195
1196 set OutputArray muons
1197
1198 set DeltaRMax 0.3
1199 set PTMin 1.0
1200 set PTRatioMax 9999.
1201
1202}
1203
1204######################
1205# Muon isolation CHS #
1206######################
1207
1208module Isolation MuonIsolationCHS {
1209 set CandidateInputArray MuonMomentumSmearing/muons
1210
1211 # isolation collection
1212 set IsolationInputArray EFlowMerger/eflow
1213
1214 set OutputArray muons
1215
1216 set DeltaRMax 0.3
1217 set PTMin 1.0
1218 set PTRatioMax 9999.
1219
1220}
1221
1222
1223#####################
1224# Muon Loose Id #
1225#####################
1226
1227module Efficiency MuonLooseIdEfficiency {
1228 set InputArray MuonIsolation/muons
1229 set OutputArray muons
1230 # tracking + LooseID efficiency formula for muons
1231 source muonLooseId.tcl
1232
1233}
1234
1235
1236##################
1237# Muon Tight Id #
1238##################
1239
1240module Efficiency MuonTightIdEfficiency {
1241 set InputArray MuonIsolation/muons
1242 set OutputArray muons
1243 # tracking + TightID efficiency formula for muons
1244 source muonTightId.tcl
1245}
1246
1247
1248#####################
1249# Muon Loose Id CHS #
1250#####################
1251
1252module Efficiency MuonLooseIdEfficiencyCHS {
1253 set InputArray MuonIsolationCHS/muons
1254 set OutputArray muons
1255 # tracking + LooseID efficiency formula for muons
1256 source muonLooseId.tcl
1257
1258}
1259
1260
1261######################
1262# Muon Tight Id CHS #
1263######################
1264
1265module Efficiency MuonTightIdEfficiencyCHS {
1266 set InputArray MuonIsolationCHS/muons
1267 set OutputArray muons
1268 # tracking + TightID efficiency formula for muons
1269 source muonTightId.tcl
1270}
1271
1272
1273########################
1274# Jet Flavor Association
1275########################
1276
1277module JetFlavorAssociation JetFlavorAssociation {
1278
1279 set PartonInputArray Delphes/partons
1280 set ParticleInputArray Delphes/allParticles
1281 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1282 set JetInputArray JetEnergyScale/jets
1283
1284 set DeltaR 0.5
1285 set PartonPTMin 10.0
1286 set PartonEtaMax 4.0
1287
1288}
1289
1290module JetFlavorAssociation JetFlavorAssociationAK8 {
1291
1292 set PartonInputArray Delphes/partons
1293 set ParticleInputArray Delphes/allParticles
1294 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1295 set JetInputArray JetEnergyScaleAK8/jets
1296
1297 set DeltaR 0.8
1298 set PartonPTMin 100.0
1299 set PartonEtaMax 4.0
1300
1301}
1302
1303module JetFlavorAssociation JetFlavorAssociationPUPPI {
1304
1305 set PartonInputArray Delphes/partons
1306 set ParticleInputArray Delphes/allParticles
1307 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1308 set JetInputArray JetEnergyScalePUPPI/jets
1309
1310 set DeltaR 0.5
1311 set PartonPTMin 10.0
1312 set PartonEtaMax 4.0
1313
1314}
1315
1316module JetFlavorAssociation JetFlavorAssociationPUPPIAK8 {
1317
1318 set PartonInputArray Delphes/partons
1319 set ParticleInputArray Delphes/allParticles
1320 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1321 set JetInputArray JetEnergyScalePUPPIAK8/jets
1322
1323 set DeltaR 0.8
1324 set PartonPTMin 100.0
1325 set PartonEtaMax 4.0
1326
1327}
1328
1329
1330#############
1331# b-tagging #
1332#############
1333module BTagging BTaggingLoose {
1334
1335 set JetInputArray JetEnergyScale/jets
1336
1337 set BitNumber 0
1338
1339 source btagLoose.tcl
1340}
1341
1342module BTagging BTaggingLooseAK8 {
1343
1344 set JetInputArray JetEnergyScaleAK8/jets
1345
1346 set BitNumber 0
1347
1348 source btagLoose.tcl
1349}
1350
1351module BTagging BTaggingLoosePUPPI {
1352
1353 set JetInputArray JetEnergyScalePUPPI/jets
1354
1355 set BitNumber 0
1356
1357 source btagLoose.tcl
1358}
1359
1360module BTagging BTaggingLoosePUPPIAK8 {
1361
1362 set JetInputArray JetEnergyScalePUPPIAK8/jets
1363
1364 set BitNumber 0
1365
1366 source btagLoose.tcl
1367}
1368
1369
1370#############
1371# b-tagging #
1372#############
1373module BTagging BTaggingMedium {
1374
1375 set JetInputArray JetEnergyScale/jets
1376
1377 set BitNumber 1
1378
1379 source btagMedium.tcl
1380}
1381
1382module BTagging BTaggingMediumAK8 {
1383
1384 set JetInputArray JetEnergyScaleAK8/jets
1385
1386 set BitNumber 1
1387
1388 source btagMedium.tcl
1389}
1390
1391module BTagging BTaggingMediumPUPPI {
1392
1393 set JetInputArray JetEnergyScalePUPPI/jets
1394
1395 set BitNumber 1
1396
1397 source btagMedium.tcl
1398}
1399
1400module BTagging BTaggingMediumPUPPIAK8 {
1401
1402 set JetInputArray JetEnergyScalePUPPIAK8/jets
1403
1404 set BitNumber 1
1405
1406 source btagMedium.tcl
1407}
1408
1409
1410
1411#############
1412# b-tagging #
1413#############
1414module BTagging BTaggingTight {
1415
1416 set JetInputArray JetEnergyScale/jets
1417
1418 set BitNumber 2
1419
1420 source btagTight.tcl
1421}
1422
1423module BTagging BTaggingTightAK8 {
1424
1425 set JetInputArray JetEnergyScaleAK8/jets
1426
1427 set BitNumber 2
1428
1429 source btagTight.tcl
1430}
1431
1432module BTagging BTaggingTightPUPPI {
1433
1434 set JetInputArray JetEnergyScalePUPPI/jets
1435
1436 set BitNumber 2
1437
1438 source btagTight.tcl
1439}
1440
1441module BTagging BTaggingTightPUPPIAK8 {
1442
1443 set JetInputArray JetEnergyScalePUPPIAK8/jets
1444
1445 set BitNumber 2
1446
1447 source btagTight.tcl
1448}
1449
1450#############
1451# tau-tagging
1452#############
1453
1454
1455module TauTagging TauTagging {
1456 set ParticleInputArray Delphes/allParticles
1457 set PartonInputArray Delphes/partons
1458 set JetInputArray JetEnergyScale/jets
1459
1460 set DeltaR 0.5
1461
1462 set TauPTMin 20.0
1463
1464 set TauEtaMax 2.3
1465
1466 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1467
1468 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)) + \
1469 (abs(eta) > 2.3) * (0.000)
1470 }
1471 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) ) + \
1472 (abs(eta) > 2.3) * (0.000)
1473 }
1474}
1475
1476
1477###############################################################################################################
1478# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1479###############################################################################################################
1480
1481module StatusPidFilter GenParticleFilter {
1482
1483 set InputArray Delphes/allParticles
1484 set OutputArray filteredParticles
1485 set PTMin 5.0
1486
1487}
1488
1489
1490##################
1491# ROOT tree writer
1492##################
1493
1494module TreeWriter TreeWriter {
1495# add Branch InputArray BranchName BranchClass
1496 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1497 add Branch PileUpMerger/vertices Vertex Vertex
1498
1499 add Branch GenJetFinder/jets GenJet Jet
1500 add Branch GenJetFinderAK8/jetsAK8 GenJetAK8 Jet
1501 add Branch GenMissingET/momentum GenMissingET MissingET
1502
1503# add Branch HCal/eflowTracks EFlowTrack Track
1504# add Branch ECal/eflowPhotons EFlowPhoton Tower
1505# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1506
1507 add Branch PhotonEfficiency/photons Photon Photon
1508 add Branch ElectronEfficiency/electrons Electron Electron
1509 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
1510 add Branch MuonTightIdEfficiency/muons MuonTight Muon
1511
1512 add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
1513 add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
1514 add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
1515 add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon
1516
1517 add Branch JetEnergyScale/jets Jet Jet
1518# add Branch FatJetFinder/jets FatJet Jet
1519 add Branch JetEnergyScalePUPPI/jets JetPUPPI Jet
1520 add Branch JetEnergyScaleAK8/jets JetAK8 Jet
1521 add Branch JetEnergyScalePUPPIAK8/jets JetPUPPIAK8 Jet
1522
1523 add Branch MissingET/momentum MissingET MissingET
1524 add Branch PuppiMissingET/momentum PuppiMissingET MissingET
1525 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1526 add Branch ScalarHT/energy ScalarHT ScalarHT
1527
1528}
Note: See TracBrowser for help on using the repository browser.