Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl@ d2cf547

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

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