Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl@ 479f808

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

fixed photon isolation with puppi
  • Property mode set to 100644
File size: 27.2 KB
RevLine 
[740d430]1set MaxEvents 10
2
[f7f33af]3#
4# Phase II - Pile-Up
5#
6# Main authors: Michele Selvaggi (UCL)
[740d430]7#
8# Released on:
[f7f33af]9#
10# Version: v01
11#
12#
13#######################################
14# Order of execution of various modules
15#######################################
16
17set ExecutionPath {
[740d430]18
19 PileUpMerger
[f7f33af]20 ParticlePropagator
21
22 ChargedHadronTrackingEfficiency
23 ElectronTrackingEfficiency
24 MuonTrackingEfficiency
25
26 ChargedHadronMomentumSmearing
[df7986a]27 ElectronEnergySmearing
[f7f33af]28 MuonMomentumSmearing
29
30 TrackMerger
[740d430]31
[f7f33af]32 ECal
33 HCal
34
35 ElectronFilter
36 TrackPileUpSubtractor
[740d430]37
[f7f33af]38 TowerMerger
39 NeutralEFlowMerger
40 EFlowMergerAllTracks
41 EFlowMerger
42
[e921a28]43 LeptonFilterNoLep
44 LeptonFilterLep
45 RunPUPPIBase
[f7f33af]46 RunPUPPI
47
[740d430]48 PhotonFilter
[e921a28]49
[f7f33af]50 PhotonIsolation
51 PhotonEfficiency
[740d430]52
[f7f33af]53 ElectronIsolation
54 ElectronEfficiency
[740d430]55
[f7f33af]56 MuonIsolation
57 MuonLooseIdEfficiency
58 MuonTightIdEfficiency
[740d430]59
[f7f33af]60 NeutrinoFilter
61
62 MissingET
[e921a28]63 PuppiMissingET
[f7f33af]64 GenMissingET
65 GenPileUpMissingET
66
67 GenJetFinder
68 FastJetFinder
69
70 ScalarHT
[740d430]71
[f7f33af]72 JetEnergyScale
73
74 JetFlavorAssociation
75
76 BTaggingLoose
77 BTaggingMedium
78 BTaggingTight
79
80 TauTagging
[740d430]81
[f7f33af]82 GenParticleFilter
[740d430]83
[f7f33af]84 TreeWriter
85}
86
87
88
89###############
90# PileUp Merger
91###############
92
93module PileUpMerger PileUpMerger {
94 set InputArray Delphes/stableParticles
95
96 set ParticleOutputArray stableParticles
97 set VertexOutputArray vertices
98
99 # pre-generated minbias input file
[740d430]100 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
[f7f33af]101
102 # average expected pile up
103 set MeanPileUp 200
[740d430]104
[f7f33af]105 # maximum spread in the beam direction in m
[94c5375]106 set ZVertexSpread 0.25
[f7f33af]107
108 # maximum spread in time in s
[94c5375]109 set TVertexSpread 800E-12
[f7f33af]110
111 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
[94c5375]112 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
[740d430]113
[f7f33af]114}
115
116
117#################################
118# Propagate particles in cylinder
119#################################
120
121module ParticlePropagator ParticlePropagator {
122 set InputArray PileUpMerger/stableParticles
123
124 set OutputArray stableParticles
125 set ChargedHadronOutputArray chargedHadrons
126 set ElectronOutputArray electrons
127 set MuonOutputArray muons
128
129 # radius of the magnetic field coverage, in m
130 set Radius 1.29
131 # half-length of the magnetic field coverage, in m
132 set HalfLength 3.0
133
134 # magnetic field
135 set Bz 3.8
136}
137
138
139####################################
140# Charged hadron tracking efficiency
141####################################
142
143module Efficiency ChargedHadronTrackingEfficiency {
144 ## particles after propagation
145 set InputArray ParticlePropagator/chargedHadrons
146 set OutputArray chargedHadrons
147 # tracking efficiency formula for charged hadrons
148 set EfficiencyFormula {
149 (pt <= 0.2) * (0.00) + \
150 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
151 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
152 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
153 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
154 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
155 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
156 (abs(eta) > 4.0) * (0.00)
157 }
158}
159
160
161#####################################
162# Electron tracking efficiency - ID
163####################################
164
165module Efficiency ElectronTrackingEfficiency {
166 set InputArray ParticlePropagator/electrons
167 set OutputArray electrons
168 # tracking efficiency formula for electrons
[740d430]169 set EfficiencyFormula {
[f7f33af]170 (pt <= 0.2) * (0.00) + \
171 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
172 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
173 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
174 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
175 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
176 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
177 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
178 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
179 (abs(eta) > 4.0) * (0.00)
180
181 }
182}
183
184##########################
185# Muon tracking efficiency
186##########################
187
188module Efficiency MuonTrackingEfficiency {
189 set InputArray ParticlePropagator/muons
190 set OutputArray muons
191 # tracking efficiency formula for muons
192 set EfficiencyFormula {
193 (pt <= 0.2) * (0.00) + \
194 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
195 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
196 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
197 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
198 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
199 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
200 (abs(eta) > 4.0) * (0.00)
[740d430]201
[f7f33af]202 }
203}
204
205
206########################################
207# Momentum resolution for charged tracks
208########################################
209
210module MomentumSmearing ChargedHadronMomentumSmearing {
211 ## hadrons after having applied the tracking efficiency
212 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
213 set OutputArray chargedHadrons
[740d430]214 # resolution formula for charged hadrons ,
215
[f7f33af]216 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
217 source trackMomentumResolution.tcl
218}
219
[df7986a]220#################################
221# Energy resolution for electrons
222#################################
[f7f33af]223
[df7986a]224module EnergySmearing ElectronEnergySmearing {
225 set InputArray ElectronTrackingEfficiency/electrons
[f7f33af]226 set OutputArray electrons
[df7986a]227
228 # set ResolutionFormula {resolution formula as a function of eta and energy}
229
[f7f33af]230 # resolution formula for electrons
[740d430]231
[df7986a]232 # taking something flat in energy for now, ECAL will take over at high energy anyway.
233 # inferred from hep-ex/1306.2016 and 1502.02701
[740d430]234 set ResolutionFormula {
235
[df7986a]236 (abs(eta) <= 1.5) * (1+0.64*abs(eta)^2)*(energy*0.028) +
237 (abs(eta) > 1.5 && abs(eta) <= 1.75) * (energy*0.037) +
238 (abs(eta) > 1.75 && abs(eta) <= 2.15) * (energy*0.038) +
239 (abs(eta) > 2.15 && abs(eta) <= 3.00) * (energy*0.044) +
[740d430]240 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)}
[f7f33af]241
242}
243
244###############################
245# Momentum resolution for muons
246###############################
247
248module MomentumSmearing MuonMomentumSmearing {
249 set InputArray MuonTrackingEfficiency/muons
250 set OutputArray muons
251 # resolution formula for muons
[740d430]252
253 # up to |eta| < 2.8 take measurement from tracking + muon chambers
254 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
[f7f33af]255 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
[740d430]256 source muonMomentumResolution.tcl
[f7f33af]257}
258
259
260
261
262##############
263# Track merger
264##############
265
266module Merger TrackMerger {
267# add InputArray InputArray
268 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
[df7986a]269 add InputArray ElectronEnergySmearing/electrons
[f7f33af]270 add InputArray MuonMomentumSmearing/muons
271 set OutputArray tracks
272}
273
274#############
275# ECAL
276#############
277
278module SimpleCalorimeter ECal {
279 set ParticleInputArray ParticlePropagator/stableParticles
280 set TrackInputArray TrackMerger/tracks
281
282 set TowerOutputArray ecalTowers
283 set EFlowTrackOutputArray eflowTracks
284 set EFlowTowerOutputArray eflowPhotons
285
[740d430]286 set IsEcal true
287
[f7f33af]288 set EnergyMin 0.5
289 set EnergySignificanceMin 1.0
290
291 set SmearTowerCenter true
292
293 set pi [expr {acos(-1)}]
294
295 # lists of the edges of each tower in eta and phi
296 # each list starts with the lower edge of the first tower
297 # the list ends with the higher edged of the last tower
298
299 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
[740d430]300
[f7f33af]301 set PhiBins {}
302 for {set i -180} {$i <= 180} {incr i} {
303 add PhiBins [expr {$i * $pi/180.0}]
304 }
305
306 # 0.02 unit in eta up to eta = 1.5 (barrel)
307 for {set i -85} {$i <= 86} {incr i} {
308 set eta [expr {$i * 0.0174}]
309 add EtaPhiBins $eta $PhiBins
310 }
311
312 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
[740d430]313
[f7f33af]314 set PhiBins {}
315 for {set i -180} {$i <= 180} {incr i} {
316 add PhiBins [expr {$i * $pi/180.0}]
317 }
318
319 # 0.02 unit in eta up to eta = 3
320 for {set i 1} {$i <= 84} {incr i} {
321 set eta [expr { -2.958 + $i * 0.0174}]
322 add EtaPhiBins $eta $PhiBins
323 }
324
325 for {set i 1} {$i <= 84} {incr i} {
326 set eta [expr { 1.4964 + $i * 0.0174}]
327 add EtaPhiBins $eta $PhiBins
328 }
329
[740d430]330 # take present CMS granularity for HF
331
[f7f33af]332 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
333 set PhiBins {}
334 for {set i -18} {$i <= 18} {incr i} {
335 add PhiBins [expr {$i * $pi/18.0}]
336 }
[740d430]337
[f7f33af]338 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} {
339 add EtaPhiBins $eta $PhiBins
340 }
341
342
343 add EnergyFraction {0} {0.0}
344 # energy fractions for e, gamma and pi0
345 add EnergyFraction {11} {1.0}
346 add EnergyFraction {22} {1.0}
347 add EnergyFraction {111} {1.0}
348 # energy fractions for muon, neutrinos and neutralinos
349 add EnergyFraction {12} {0.0}
350 add EnergyFraction {13} {0.0}
351 add EnergyFraction {14} {0.0}
352 add EnergyFraction {16} {0.0}
353 add EnergyFraction {1000022} {0.0}
354 add EnergyFraction {1000023} {0.0}
355 add EnergyFraction {1000025} {0.0}
356 add EnergyFraction {1000035} {0.0}
357 add EnergyFraction {1000045} {0.0}
358 # energy fractions for K0short and Lambda
359 add EnergyFraction {310} {0.3}
360 add EnergyFraction {3122} {0.3}
361
362 # set ResolutionFormula {resolution formula as a function of eta and energy}
[740d430]363
[df7986a]364 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
[f7f33af]365 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
366
[740d430]367 set ResolutionFormula { (abs(eta) <= 1.50) * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
[db7f5e4]368 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
369 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
370 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
371 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
[f7f33af]372
373}
374
375#############
376# HCAL
377#############
378
379module SimpleCalorimeter HCal {
380 set ParticleInputArray ParticlePropagator/stableParticles
381 set TrackInputArray ECal/eflowTracks
382
383 set TowerOutputArray hcalTowers
384 set EFlowTrackOutputArray eflowTracks
385 set EFlowTowerOutputArray eflowNeutralHadrons
386
[740d430]387 set IsEcal false
388
[f7f33af]389 set EnergyMin 1.0
390 set EnergySignificanceMin 1.0
391
392 set SmearTowerCenter true
393
394 set pi [expr {acos(-1)}]
395
396 # lists of the edges of each tower in eta and phi
397 # each list starts with the lower edge of the first tower
398 # the list ends with the higher edged of the last tower
[740d430]399
[f7f33af]400 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
401
402 set PhiBins {}
403 for {set i -36} {$i <= 36} {incr i} {
404 add PhiBins [expr {$i * $pi/36.0}]
405 }
406 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} {
407 add EtaPhiBins $eta $PhiBins
408 }
409
410 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
[740d430]411
[f7f33af]412 set PhiBins {}
413 for {set i -180} {$i <= 180} {incr i} {
414 add PhiBins [expr {$i * $pi/180.0}]
415 }
416
417 # 0.02 unit in eta up to eta = 3
418 for {set i 1} {$i <= 84} {incr i} {
419 set eta [expr { -2.958 + $i * 0.0174}]
420 add EtaPhiBins $eta $PhiBins
421 }
422
423 for {set i 1} {$i <= 84} {incr i} {
424 set eta [expr { 1.4964 + $i * 0.0174}]
425 add EtaPhiBins $eta $PhiBins
426 }
427
[740d430]428 # take present CMS granularity for HF
429
[f7f33af]430 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
431 set PhiBins {}
432 for {set i -18} {$i <= 18} {incr i} {
433 add PhiBins [expr {$i * $pi/18.0}]
434 }
[740d430]435
[f7f33af]436 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} {
437 add EtaPhiBins $eta $PhiBins
438 }
439
440
441 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
442 add EnergyFraction {0} {1.0}
443 # energy fractions for e, gamma and pi0
444 add EnergyFraction {11} {0.0}
445 add EnergyFraction {22} {0.0}
446 add EnergyFraction {111} {0.0}
447 # energy fractions for muon, neutrinos and neutralinos
448 add EnergyFraction {12} {0.0}
449 add EnergyFraction {13} {0.0}
450 add EnergyFraction {14} {0.0}
451 add EnergyFraction {16} {0.0}
452 add EnergyFraction {1000022} {0.0}
453 add EnergyFraction {1000023} {0.0}
454 add EnergyFraction {1000025} {0.0}
455 add EnergyFraction {1000035} {0.0}
456 add EnergyFraction {1000045} {0.0}
457 # energy fractions for K0short and Lambda
458 add EnergyFraction {310} {0.7}
459 add EnergyFraction {3122} {0.7}
460
461# set ResolutionFormula {resolution formula as a function of eta and energy}
462 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
463 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
464 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
465
466}
467
468#################
469# Electron filter
470#################
471
472module PdgCodeFilter ElectronFilter {
473 set InputArray HCal/eflowTracks
474 set OutputArray electrons
475 set Invert true
476 add PdgCode {11}
477 add PdgCode {-11}
478}
479
[740d430]480
[f7f33af]481##########################
482# Track pile-up subtractor
483##########################
484
485module TrackPileUpSubtractor TrackPileUpSubtractor {
486# add InputArray InputArray OutputArray
487 add InputArray HCal/eflowTracks eflowTracks
488 add InputArray ElectronFilter/electrons electrons
489 add InputArray MuonMomentumSmearing/muons muons
490
491 set VertexInputArray PileUpMerger/vertices
492 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
493 # Z vertex resolution in m
494 set ZVertexResolution 0.0001
495}
496
497
498###################################################
499# Tower Merger (in case not using e-flow algorithm)
500###################################################
501
502module Merger TowerMerger {
503# add InputArray InputArray
504 add InputArray ECal/ecalTowers
505 add InputArray HCal/hcalTowers
506 set OutputArray towers
507}
508
509
510####################
511# Neutral eflow erger
512####################
513
514module Merger NeutralEFlowMerger {
515# add InputArray InputArray
516 add InputArray ECal/eflowPhotons
517 add InputArray HCal/eflowNeutralHadrons
518 set OutputArray eflowTowers
519}
520
521
522####################
523# Energy flow merger
524####################
525
526module Merger EFlowMerger {
527# add InputArray InputArray
528 add InputArray HCal/eflowTracks
529 add InputArray ECal/eflowPhotons
530 add InputArray HCal/eflowNeutralHadrons
531 set OutputArray eflow
532}
533
534##################################
535# Energy flow merger (all tracks)
536##################################
537
538module Merger EFlowMergerAllTracks {
539# add InputArray InputArray
540 add InputArray TrackMerger/tracks
541 add InputArray ECal/eflowPhotons
542 add InputArray HCal/eflowNeutralHadrons
543 set OutputArray eflow
544}
545
546#########################################
547### Run the puppi code (to be tuned) ###
548#########################################
549
[e921a28]550module PdgCodeFilter LeptonFilterNoLep {
551 set InputArray HCal/eflowTracks
552 set OutputArray eflowTracksNoLeptons
553 set Invert false
554 add PdgCode {13}
555 add PdgCode {-13}
556 add PdgCode {11}
557 add PdgCode {-11}
558}
559
560module PdgCodeFilter LeptonFilterLep {
561 set InputArray HCal/eflowTracks
562 set OutputArray eflowTracksLeptons
563 set Invert true
564 add PdgCode {11}
565 add PdgCode {-11}
566 add PdgCode {13}
567 add PdgCode {-13}
568}
569
570module RunPUPPI RunPUPPIBase {
[f7f33af]571 ## input information
[e921a28]572 set TrackInputArray LeptonFilterNoLep/eflowTracksNoLeptons
[f7f33af]573 set NeutralInputArray NeutralEFlowMerger/eflowTowers
574 set PVInputArray PileUpMerger/vertices
575 set MinPuppiWeight 0.05
576 set UseExp false
[e921a28]577 set UseNoLep false
[740d430]578
579 ## define puppi algorithm parameters (more than one for the same eta region is possible)
[e921a28]580 add EtaMinBin 0.0 1.5 4.0
581 add EtaMaxBin 1.5 4.0 10.0
582 add PtMinBin 0.0 0.0 0.0
583 add ConeSizeBin 0.2 0.2 0.2
584 add RMSPtMinBin 0.1 0.5 0.5
585 add RMSScaleFactorBin 1.0 1.0 1.0
586 add NeutralMinEBin 0.2 0.2 0.5
587 add NeutralPtSlope 0.006 0.013 0.067
588 add ApplyCHS true true true
589 add UseCharged true true false
[740d430]590 add ApplyLowPUCorr true true true
[e921a28]591 add MetricId 5 5 5
592 add CombId 0 0 0
[f7f33af]593
594 ## output name
595 set OutputArray PuppiParticles
596 set OutputArrayTracks puppiTracks
597 set OutputArrayNeutrals puppiNeutrals
[740d430]598}
[f7f33af]599
[e921a28]600module Merger RunPUPPI {
601 add InputArray RunPUPPIBase/PuppiParticles
602 add InputArray LeptonFilterLep/eflowTracksLeptons
603 set OutputArray PuppiParticles
604}
[f7f33af]605
606###################
607# Missing ET merger
608###################
609
610module Merger MissingET {
611# add InputArray InputArray
612# add InputArray RunPUPPI/PuppiParticles
613 add InputArray EFlowMerger/eflow
614 set MomentumOutputArray momentum
615}
616
[e921a28]617module Merger PuppiMissingET {
618 #add InputArray InputArray
619 add InputArray RunPUPPI/PuppiParticles
620 #add InputArray EFlowMerger/eflow
621 set MomentumOutputArray momentum
622}
623
[f7f33af]624###################
[740d430]625# Ger PileUp Missing ET
[f7f33af]626###################
627
628module Merger GenPileUpMissingET {
629# add InputArray InputArray
630# add InputArray RunPUPPI/PuppiParticles
631 add InputArray ParticlePropagator/stableParticles
632 set MomentumOutputArray momentum
633}
634
635##################
636# Scalar HT merger
637##################
638
639module Merger ScalarHT {
640# add InputArray InputArray
641 add InputArray RunPUPPI/PuppiParticles
642 set EnergyOutputArray energy
643}
644
645#################
646# Neutrino Filter
647#################
648
649module PdgCodeFilter NeutrinoFilter {
650
651 set InputArray Delphes/stableParticles
652 set OutputArray filteredParticles
653
654 set PTMin 0.0
655
656 add PdgCode {12}
657 add PdgCode {14}
658 add PdgCode {16}
659 add PdgCode {-12}
660 add PdgCode {-14}
661 add PdgCode {-16}
662
663}
664
665
666#####################
667# MC truth jet finder
668#####################
669
670module FastJetFinder GenJetFinder {
671 set InputArray NeutrinoFilter/filteredParticles
[740d430]672
[f7f33af]673 set OutputArray jets
674
675 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
676 set JetAlgorithm 6
677 set ParameterR 0.4
678
679 set JetPTMin 15.0
680}
681
682#########################
683# Gen Missing ET merger
684########################
685
686module Merger GenMissingET {
687
688# add InputArray InputArray
689 add InputArray NeutrinoFilter/filteredParticles
690 set MomentumOutputArray momentum
691}
692
693
694
695############
696# Jet finder
697############
698
699module FastJetFinder FastJetFinder {
700# set InputArray TowerMerger/towers
701 set InputArray RunPUPPI/PuppiParticles
702
703 set OutputArray jets
704
705 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
706 set JetAlgorithm 6
707 set ParameterR 0.4
708
709 set JetPTMin 15.0
710}
711
712##################
713# Jet Energy Scale
714##################
715
716module EnergyScale JetEnergyScale {
717 set InputArray FastJetFinder/jets
718 set OutputArray jets
719
720 # scale formula for jets
721 set ScaleFormula {1.00}
722}
723
724
[740d430]725#################
726# Photon filter
727#################
728
729module PdgCodeFilter PhotonFilter {
730 set InputArray ECal/eflowPhotons
731 set OutputArray photons
732 set Invert true
733 set PTMin 5.0
734 add PdgCode {22}
735}
736
737
[f7f33af]738####################
739# Photon isolation #
740####################
741
742module Isolation PhotonIsolation {
[740d430]743
[f7f33af]744 # particle for which calculate the isolation
[740d430]745 set CandidateInputArray PhotonFilter/photons
746
[f7f33af]747 # isolation collection
[740d430]748 set IsolationInputArray RunPUPPI/PuppiParticles
749
[f7f33af]750 # output array
751 set OutputArray photons
[740d430]752
753 # veto isolation cand. based on proximity to input cand.
754 set DeltaRMin 0.01
755 set UseMiniCone true
756
[f7f33af]757 # isolation cone
758 set DeltaRMax 0.3
[740d430]759
760 # minimum pT
[f7f33af]761 set PTMin 1.0
[740d430]762
[f7f33af]763 # iso ratio to cut
764 set PTRatioMax 9999.
[740d430]765
[f7f33af]766}
767
768
769
770#####################
771# Photon efficiency #
772#####################
773
774module Efficiency PhotonEfficiency {
[740d430]775
[f7f33af]776 ## input particles
[740d430]777 set InputArray PhotonIsolation/photons
[f7f33af]778 ## output particles
779 set OutputArray photons
780 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
781 # efficiency formula for photons
782 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
783 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
784 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
785 (abs(eta) > 4.0) * (0.00)}
786
787}
788
789
790######################
791# Electron isolation #
792######################
793
794module Isolation ElectronIsolation {
[740d430]795
[f7f33af]796 set CandidateInputArray ElectronFilter/electrons
[740d430]797
[f7f33af]798 # isolation collection
[740d430]799 set IsolationInputArray RunPUPPI/PuppiParticles
800 #set IsolationInputArray EFlowMerger/eflow
801
[f7f33af]802 set OutputArray electrons
[740d430]803
[f7f33af]804 set DeltaRMax 0.3
805 set PTMin 1.0
806 set PTRatioMax 9999.
807
808}
809
810
811
812#######################
813# Electron efficiency #
814#######################
815
816module Efficiency ElectronEfficiency {
[740d430]817
[f7f33af]818 set InputArray ElectronIsolation/electrons
819 set OutputArray electrons
[740d430]820
[f7f33af]821 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
822 # efficiency formula for electrons
823 set EfficiencyFormula {
824 (pt <= 4.0) * (0.00) + \
825 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
826 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
827 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
[740d430]828 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
829 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
830 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
831 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
[f7f33af]832 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
[740d430]833 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
834 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
835 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
[f7f33af]836 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
[740d430]837 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
838 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
839 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
840 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
[f7f33af]841 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
[740d430]842 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
843 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
844 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
845 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
[f7f33af]846 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
847 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
848 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
849 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
850 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
851 (abs(eta) > 4.0) * (0.00)
[740d430]852
[f7f33af]853 }
854}
855
856##################
857# Muon isolation #
858##################
859
860module Isolation MuonIsolation {
861 set CandidateInputArray MuonMomentumSmearing/muons
[740d430]862
[f7f33af]863 # isolation collection
[740d430]864 set IsolationInputArray RunPUPPI/PuppiParticles
865
[f7f33af]866 set OutputArray muons
[740d430]867
[f7f33af]868 set DeltaRMax 0.3
869 set PTMin 1.0
870 set PTRatioMax 9999.
871
872}
873
874
875##################
876# Muon Loose Id #
877##################
878
879module Efficiency MuonLooseIdEfficiency {
880 set InputArray MuonIsolation/muons
881 set OutputArray muons
882 # tracking + LooseID efficiency formula for muons
883 source muonLooseId.tcl
884
885}
886
887
888##################
889# Muon Tight Id #
890##################
891
892module Efficiency MuonTightIdEfficiency {
893 set InputArray MuonIsolation/muons
894 set OutputArray muons
895 # tracking + TightID efficiency formula for muons
896 source muonTightId.tcl
897}
898
899
900
901########################
902# Jet Flavor Association
903########################
904
905module JetFlavorAssociation JetFlavorAssociation {
906
907 set PartonInputArray Delphes/partons
908 set ParticleInputArray Delphes/allParticles
909 set ParticleLHEFInputArray Delphes/allParticlesLHEF
910 set JetInputArray JetEnergyScale/jets
911
912 set DeltaR 0.5
913 set PartonPTMin 10.0
914 set PartonEtaMax 4.0
915
916}
917
918
919#############
920# b-tagging #
921#############
922module BTagging BTaggingLoose {
923
924 set JetInputArray JetEnergyScale/jets
925
926 set BitNumber 0
927
928 source btagLoose.tcl
929}
930
931
932#############
933# b-tagging #
934#############
935module BTagging BTaggingMedium {
936
937 set JetInputArray JetEnergyScale/jets
938
939 set BitNumber 1
[740d430]940
[f7f33af]941 source btagMedium.tcl
942}
943
944
945#############
946# b-tagging #
947#############
948module BTagging BTaggingTight {
949
950 set JetInputArray JetEnergyScale/jets
951
952 set BitNumber 2
953
[740d430]954 source btagTight.tcl
[f7f33af]955}
956
957
958#############
959# tau-tagging
960#############
961
962
963module TauTagging TauTagging {
964 set ParticleInputArray Delphes/allParticles
965 set PartonInputArray Delphes/partons
966 set JetInputArray JetEnergyScale/jets
967
968 set DeltaR 0.5
969
970 set TauPTMin 20.0
971
972 set TauEtaMax 2.3
973
974 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
975
976 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)) + \
977 (abs(eta) > 2.3) * (0.000)
978 }
979 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) ) + \
980 (abs(eta) > 2.3) * (0.000)
[740d430]981 }
[f7f33af]982}
983
984
985###############################################################################################################
986# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
987###############################################################################################################
988
989module StatusPidFilter GenParticleFilter {
[740d430]990
[f7f33af]991 set InputArray Delphes/allParticles
992 set OutputArray filteredParticles
993 set PTMin 5.0
[740d430]994
[f7f33af]995}
996
997
998##################
999# ROOT tree writer
1000##################
1001
1002module TreeWriter TreeWriter {
1003# add Branch InputArray BranchName BranchClass
1004 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1005 add Branch PileUpMerger/vertices Vertex Vertex
[740d430]1006
[f7f33af]1007 add Branch GenJetFinder/jets GenJet Jet
1008 add Branch GenMissingET/momentum GenMissingET MissingET
1009
1010# add Branch HCal/eflowTracks EFlowTrack Track
1011# add Branch ECal/eflowPhotons EFlowPhoton Tower
1012# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
[740d430]1013
[f7f33af]1014 add Branch PhotonEfficiency/photons Photon Photon
1015 add Branch ElectronEfficiency/electrons Electron Electron
1016 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
1017 add Branch MuonTightIdEfficiency/muons MuonTight Muon
[740d430]1018
[f7f33af]1019 add Branch JetEnergyScale/jets Jet Jet
[740d430]1020
[f7f33af]1021 add Branch MissingET/momentum MissingET MissingET
[e921a28]1022 add Branch PuppiMissingET/momentum PuppiMissingET MissingET
[f7f33af]1023 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1024 add Branch ScalarHT/energy ScalarHT ScalarHT
1025}
1026
Note: See TracBrowser for help on using the repository browser.