Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ b195ba1

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

replace calo by simple calo in CMS cards
  • Property mode set to 100644
File size: 24.1 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6
7 PileUpMerger
8 ParticlePropagator
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronMomentumSmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19
20 ECal
21 HCal
22
23 ElectronFilter
24 TrackPileUpSubtractor
25 NeutralTowerMerger
26 EFlowMergerAllTracks
27 EFlowMerger
28
29 NeutrinoFilter
30 GenJetFinder
31 GenMissingET
32
33 Rho
34 FastJetFinder
35 PileUpJetID
36 JetPileUpSubtractor
37
38 JetEnergyScale
39
40 PhotonEfficiency
41 PhotonIsolation
42
43 ElectronEfficiency
44 ElectronIsolation
45
46 MuonEfficiency
47 MuonIsolation
48
49 MissingET
50
51 JetFlavorAssociation
52
53 BTagging
54 TauTagging
55
56 UniqueObjectFinder
57
58 ScalarHT
59
60 TreeWriter
61}
62
63###############
64# PileUp Merger
65###############
66
67module PileUpMerger PileUpMerger {
68 set InputArray Delphes/stableParticles
69
70 set ParticleOutputArray stableParticles
71 set VertexOutputArray vertices
72
73 # pre-generated minbias input file
74 set PileUpFile MinBias.pileup
75
76 # average expected pile up
77 set MeanPileUp 50
78
79 # maximum spread in the beam direction in m
80 set ZVertexSpread 0.25
81
82 # maximum spread in time in s
83 set TVertexSpread 800E-12
84
85 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
86 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
87
88
89}
90
91#################################
92# Propagate particles in cylinder
93#################################
94
95module ParticlePropagator ParticlePropagator {
96 set InputArray PileUpMerger/stableParticles
97
98 set OutputArray stableParticles
99 set ChargedHadronOutputArray chargedHadrons
100 set ElectronOutputArray electrons
101 set MuonOutputArray muons
102
103 # radius of the magnetic field coverage, in m
104 set Radius 1.29
105 # half-length of the magnetic field coverage, in m
106 set HalfLength 3.00
107
108 # magnetic field
109 set Bz 3.8
110}
111
112####################################
113# Charged hadron tracking efficiency
114####################################
115
116module Efficiency ChargedHadronTrackingEfficiency {
117 set InputArray ParticlePropagator/chargedHadrons
118 set OutputArray chargedHadrons
119
120 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
121
122 # tracking efficiency formula for charged hadrons
123 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
124 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
125 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
126 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
127 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
128 (abs(eta) > 2.5) * (0.00)}
129}
130
131##############################
132# Electron tracking efficiency
133##############################
134
135module Efficiency ElectronTrackingEfficiency {
136 set InputArray ParticlePropagator/electrons
137 set OutputArray electrons
138
139 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
140
141 # tracking efficiency formula for electrons
142 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
143 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
144 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
145 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
148 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
149 (abs(eta) > 2.5) * (0.00)}
150}
151
152##########################
153# Muon tracking efficiency
154##########################
155
156module Efficiency MuonTrackingEfficiency {
157 set InputArray ParticlePropagator/muons
158 set OutputArray muons
159
160 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
161
162 # tracking efficiency formula for muons
163 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
164 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
165 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
166 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
167
168 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
169 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
170 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
171 (abs(eta) > 2.5) * (0.00)}
172}
173
174########################################
175# Momentum resolution for charged tracks
176########################################
177
178module MomentumSmearing ChargedHadronMomentumSmearing {
179 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
180 set OutputArray chargedHadrons
181
182 # set ResolutionFormula {resolution formula as a function of eta and pt}
183
184 # resolution formula for charged hadrons
185 # based on arXiv:1405.6569
186 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
187 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
188 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
189}
190
191###################################
192# Momentum resolution for electrons
193###################################
194
195module MomentumSmearing ElectronMomentumSmearing {
196 set InputArray ElectronTrackingEfficiency/electrons
197 set OutputArray electrons
198
199 # set ResolutionFormula {resolution formula as a function of eta and energy}
200
201 # resolution formula for electrons
202 # based on arXiv:1405.6569
203 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
204 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
205 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
206}
207
208###############################
209# Momentum resolution for muons
210###############################
211
212module MomentumSmearing MuonMomentumSmearing {
213 set InputArray MuonTrackingEfficiency/muons
214 set OutputArray muons
215
216 # set ResolutionFormula {resolution formula as a function of eta and pt}
217
218 # resolution formula for muons
219 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
220 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
221 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
222}
223
224##############
225# Track merger
226##############
227
228module Merger TrackMerger {
229# add InputArray InputArray
230 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
231 add InputArray ElectronMomentumSmearing/electrons
232 add InputArray MuonMomentumSmearing/muons
233 set OutputArray tracks
234}
235
236#############
237# ECAL
238#############
239
240module SimpleCalorimeter ECal {
241 set ParticleInputArray ParticlePropagator/stableParticles
242 set TrackInputArray TrackMerger/tracks
243
244 set TowerOutputArray ecalTowers
245 set EFlowTrackOutputArray eflowTracks
246 set EFlowTowerOutputArray eflowPhotons
247
248 set IsEcal true
249
250 set EnergyMin 0.5
251 set EnergySignificanceMin 2.0
252
253 set SmearTowerCenter true
254
255 set pi [expr {acos(-1)}]
256
257 # lists of the edges of each tower in eta and phi
258 # each list starts with the lower edge of the first tower
259 # the list ends with the higher edged of the last tower
260
261 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
262
263 set PhiBins {}
264 for {set i -180} {$i <= 180} {incr i} {
265 add PhiBins [expr {$i * $pi/180.0}]
266 }
267
268 # 0.02 unit in eta up to eta = 1.5 (barrel)
269 for {set i -85} {$i <= 86} {incr i} {
270 set eta [expr {$i * 0.0174}]
271 add EtaPhiBins $eta $PhiBins
272 }
273
274 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
275
276 set PhiBins {}
277 for {set i -180} {$i <= 180} {incr i} {
278 add PhiBins [expr {$i * $pi/180.0}]
279 }
280
281 # 0.02 unit in eta up to eta = 3
282 for {set i 1} {$i <= 84} {incr i} {
283 set eta [expr { -2.958 + $i * 0.0174}]
284 add EtaPhiBins $eta $PhiBins
285 }
286
287 for {set i 1} {$i <= 84} {incr i} {
288 set eta [expr { 1.4964 + $i * 0.0174}]
289 add EtaPhiBins $eta $PhiBins
290 }
291
292 # take present CMS granularity for HF
293
294 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
295 set PhiBins {}
296 for {set i -18} {$i <= 18} {incr i} {
297 add PhiBins [expr {$i * $pi/18.0}]
298 }
299
300 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} {
301 add EtaPhiBins $eta $PhiBins
302 }
303
304
305 add EnergyFraction {0} {0.0}
306 # energy fractions for e, gamma and pi0
307 add EnergyFraction {11} {1.0}
308 add EnergyFraction {22} {1.0}
309 add EnergyFraction {111} {1.0}
310 # energy fractions for muon, neutrinos and neutralinos
311 add EnergyFraction {12} {0.0}
312 add EnergyFraction {13} {0.0}
313 add EnergyFraction {14} {0.0}
314 add EnergyFraction {16} {0.0}
315 add EnergyFraction {1000022} {0.0}
316 add EnergyFraction {1000023} {0.0}
317 add EnergyFraction {1000025} {0.0}
318 add EnergyFraction {1000035} {0.0}
319 add EnergyFraction {1000045} {0.0}
320 # energy fractions for K0short and Lambda
321 add EnergyFraction {310} {0.3}
322 add EnergyFraction {3122} {0.3}
323
324 # set ResolutionFormula {resolution formula as a function of eta and energy}
325
326 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
327
328 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
329 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
330 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
331 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
332 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
333
334}
335
336
337#############
338# HCAL
339#############
340
341module SimpleCalorimeter HCal {
342 set ParticleInputArray ParticlePropagator/stableParticles
343 set TrackInputArray ECal/eflowTracks
344
345 set TowerOutputArray hcalTowers
346 set EFlowTrackOutputArray eflowTracks
347 set EFlowTowerOutputArray eflowNeutralHadrons
348
349 set IsEcal false
350
351 set EnergyMin 1.0
352 set EnergySignificanceMin 2.0
353
354 set SmearTowerCenter true
355
356 set pi [expr {acos(-1)}]
357
358 # lists of the edges of each tower in eta and phi
359 # each list starts with the lower edge of the first tower
360 # the list ends with the higher edged of the last tower
361
362 # 5 degrees towers
363 set PhiBins {}
364 for {set i -36} {$i <= 36} {incr i} {
365 add PhiBins [expr {$i * $pi/36.0}]
366 }
367 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.653} {
368 add EtaPhiBins $eta $PhiBins
369 }
370
371 # 10 degrees towers
372 set PhiBins {}
373 for {set i -18} {$i <= 18} {incr i} {
374 add PhiBins [expr {$i * $pi/18.0}]
375 }
376 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
377 add EtaPhiBins $eta $PhiBins
378 }
379
380 # 20 degrees towers
381 set PhiBins {}
382 for {set i -9} {$i <= 9} {incr i} {
383 add PhiBins [expr {$i * $pi/9.0}]
384 }
385 foreach eta {-5 -4.7 -4.525 4.7 5} {
386 add EtaPhiBins $eta $PhiBins
387 }
388
389 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
390 add EnergyFraction {0} {1.0}
391 # energy fractions for e, gamma and pi0
392 add EnergyFraction {11} {0.0}
393 add EnergyFraction {22} {0.0}
394 add EnergyFraction {111} {0.0}
395 # energy fractions for muon, neutrinos and neutralinos
396 add EnergyFraction {12} {0.0}
397 add EnergyFraction {13} {0.0}
398 add EnergyFraction {14} {0.0}
399 add EnergyFraction {16} {0.0}
400 add EnergyFraction {1000022} {0.0}
401 add EnergyFraction {1000023} {0.0}
402 add EnergyFraction {1000025} {0.0}
403 add EnergyFraction {1000035} {0.0}
404 add EnergyFraction {1000045} {0.0}
405 # energy fractions for K0short and Lambda
406 add EnergyFraction {310} {0.7}
407 add EnergyFraction {3122} {0.7}
408
409 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
410 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
411 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
412
413}
414
415#################
416# Electron filter
417#################
418
419module PdgCodeFilter ElectronFilter {
420 set InputArray HCal/eflowTracks
421 set OutputArray electrons
422 set Invert true
423 add PdgCode {11}
424 add PdgCode {-11}
425}
426
427###################################################
428# Tower Merger (in case not using e-flow algorithm)
429###################################################
430
431module Merger Calorimeter {
432# add InputArray InputArray
433 add InputArray ECal/ecalTowers
434 add InputArray HCal/hcalTowers
435 set OutputArray towers
436}
437
438
439##########################
440# Track pile-up subtractor
441##########################
442
443module TrackPileUpSubtractor TrackPileUpSubtractor {
444# add InputArray InputArray OutputArray
445 add InputArray HCal/eflowTracks eflowTracks
446 add InputArray ElectronFilter/electrons electrons
447 add InputArray MuonMomentumSmearing/muons muons
448
449 set VertexInputArray PileUpMerger/vertices
450 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
451 # Z vertex resolution in m
452 set ZVertexResolution 0.0001
453}
454
455
456####################
457# Neutral Tower merger
458####################
459
460module Merger NeutralTowerMerger {
461# add InputArray InputArray
462 add InputArray ECal/eflowPhotons
463 add InputArray HCal/eflowNeutralHadrons
464 set OutputArray towers
465}
466
467
468####################
469# Energy flow merger
470####################
471
472module Merger EFlowMergerAllTracks {
473# add InputArray InputArray
474 add InputArray HCal/eflowTracks
475 add InputArray ECal/eflowPhotons
476 add InputArray HCal/eflowNeutralHadrons
477 set OutputArray eflow
478}
479
480
481
482
483####################
484# Energy flow merger
485####################
486
487module Merger EFlowMerger {
488# add InputArray InputArray
489 add InputArray TrackPileUpSubtractor/eflowTracks
490 add InputArray ECal/eflowPhotons
491 add InputArray HCal/eflowNeutralHadrons
492 set OutputArray eflow
493}
494
495#############
496# Rho pile-up
497#############
498
499module FastJetGridMedianEstimator Rho {
500
501 set InputArray EFlowMerger/eflow
502 set RhoOutputArray rho
503
504 # add GridRange rapmin rapmax drap dphi
505 # rapmin - the minimum rapidity extent of the grid
506 # rapmax - the maximum rapidity extent of the grid
507 # drap - the grid spacing in rapidity
508 # dphi - the grid spacing in azimuth
509
510 add GridRange -5.0 -2.5 1.0 1.0
511 add GridRange -2.5 2.5 1.0 1.0
512 add GridRange 2.5 5.0 1.0 1.0
513
514}
515
516#####################
517# Neutrino Filter
518#####################
519
520module PdgCodeFilter NeutrinoFilter {
521
522 set InputArray Delphes/stableParticles
523 set OutputArray filteredParticles
524
525 set PTMin 0.0
526
527 add PdgCode {12}
528 add PdgCode {14}
529 add PdgCode {16}
530 add PdgCode {-12}
531 add PdgCode {-14}
532 add PdgCode {-16}
533
534}
535
536
537
538#####################
539# MC truth jet finder
540#####################
541
542module FastJetFinder GenJetFinder {
543 set InputArray NeutrinoFilter/filteredParticles
544
545 set OutputArray jets
546
547 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
548 set JetAlgorithm 6
549 set ParameterR 0.5
550
551 set JetPTMin 20.0
552}
553
554#########################
555# Gen Missing ET merger
556########################
557
558module Merger GenMissingET {
559# add InputArray InputArray
560 add InputArray NeutrinoFilter/filteredParticles
561 set MomentumOutputArray momentum
562}
563
564############
565# Jet finder
566############
567
568module FastJetFinder FastJetFinder {
569# set InputArray Calorimeter/towers
570 set InputArray EFlowMerger/eflow
571
572 set OutputArray jets
573
574 # 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
575 set AreaAlgorithm 5
576
577 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
578 set JetAlgorithm 6
579 set ParameterR 0.5
580
581 set JetPTMin 20.0
582}
583
584###########################
585# Jet Pile-Up ID
586###########################
587
588module PileUpJetID PileUpJetID {
589 set JetInputArray FastJetFinder/jets
590 set TrackInputArray HCal/eflowTracks
591 set NeutralInputArray NeutralTowerMerger/towers
592
593 set VertexInputArray PileUpMerger/vertices
594 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
595 # Z vertex resolution in m
596 set ZVertexResolution 0.0001
597
598 set OutputArray jets
599
600 set UseConstituents 0
601 set ParameterR 0.5
602
603 set JetPTMin 20.0
604}
605
606###########################
607# Jet Pile-Up Subtraction
608###########################
609
610module JetPileUpSubtractor JetPileUpSubtractor {
611 set JetInputArray PileUpJetID/jets
612 set RhoInputArray Rho/rho
613
614 set OutputArray jets
615
616 set JetPTMin 20.0
617}
618
619##################
620# Jet Energy Scale
621##################
622
623module EnergyScale JetEnergyScale {
624 set InputArray JetPileUpSubtractor/jets
625 set OutputArray jets
626
627 # scale formula for jets
628 set ScaleFormula {1.0}
629}
630
631###################
632# Photon efficiency
633###################
634
635module Efficiency PhotonEfficiency {
636 set InputArray ECal/eflowPhotons
637 set OutputArray photons
638
639 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
640
641 # efficiency formula for photons
642 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
643 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
644 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
645 (abs(eta) > 2.5) * (0.00)}
646}
647
648
649##################
650# Photon isolation
651##################
652
653module Isolation PhotonIsolation {
654 set CandidateInputArray PhotonEfficiency/photons
655 set IsolationInputArray EFlowMergerAllTracks/eflow
656 set RhoInputArray Rho/rho
657
658 set OutputArray photons
659
660 set DeltaRMax 0.5
661
662 set PTMin 0.5
663
664 set PTRatioMax 0.12
665}
666
667#####################
668# Electron efficiency
669#####################
670
671module Efficiency ElectronEfficiency {
672 set InputArray TrackPileUpSubtractor/electrons
673 set OutputArray electrons
674
675 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
676
677 # efficiency formula for electrons
678 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
679 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
680 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
681 (abs(eta) > 2.5) * (0.00)}
682}
683
684####################
685# Electron isolation
686####################
687
688module Isolation ElectronIsolation {
689 set CandidateInputArray ElectronEfficiency/electrons
690 set IsolationInputArray EFlowMergerAllTracks/eflow
691 set RhoInputArray Rho/rho
692
693 set OutputArray electrons
694
695 set DeltaRMax 0.5
696
697 set PTMin 0.5
698
699 set PTRatioMax 0.12
700}
701
702#################
703# Muon efficiency
704#################
705
706module Efficiency MuonEfficiency {
707 set InputArray TrackPileUpSubtractor/muons
708 set OutputArray muons
709
710 # set EfficiencyFormula {efficiency as a function of eta and pt}
711
712 # efficiency formula for muons
713 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
714 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
715 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
716 (abs(eta) > 2.4) * (0.00)}
717
718}
719
720################
721# Muon isolation
722################
723
724module Isolation MuonIsolation {
725 set CandidateInputArray MuonEfficiency/muons
726 set IsolationInputArray EFlowMergerAllTracks/eflow
727 set RhoInputArray Rho/rho
728
729 set OutputArray muons
730
731 set DeltaRMax 0.5
732
733 set PTMin 0.5
734
735 set PTRatioMax 0.25
736}
737
738###################
739# Missing ET merger
740###################
741
742module Merger MissingET {
743# add InputArray InputArray
744 add InputArray EFlowMergerAllTracks/eflow
745 set MomentumOutputArray momentum
746}
747
748
749
750##################
751# Scalar HT merger
752##################
753
754module Merger ScalarHT {
755# add InputArray InputArray
756 add InputArray UniqueObjectFinder/jets
757 add InputArray UniqueObjectFinder/electrons
758 add InputArray UniqueObjectFinder/photons
759 add InputArray UniqueObjectFinder/muons
760 set EnergyOutputArray energy
761}
762
763########################
764# Jet Flavor Association
765########################
766
767module JetFlavorAssociation JetFlavorAssociation {
768
769 set PartonInputArray Delphes/partons
770 set ParticleInputArray Delphes/allParticles
771 set ParticleLHEFInputArray Delphes/allParticlesLHEF
772 set JetInputArray JetEnergyScale/jets
773
774 set DeltaR 0.5
775 set PartonPTMin 1.0
776 set PartonEtaMax 2.5
777
778}
779
780###########
781# b-tagging
782###########
783
784module BTagging BTagging {
785 set JetInputArray JetEnergyScale/jets
786
787 set BitNumber 0
788
789 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
790 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
791 # gluon's PDG code has the lowest priority
792
793 # based on arXiv:1211.4462
794
795 # default efficiency formula (misidentification rate)
796 add EfficiencyFormula {0} {0.01+0.000038*pt}
797
798 # efficiency formula for c-jets (misidentification rate)
799 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
800
801 # efficiency formula for b-jets
802 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
803}
804
805#############
806# tau-tagging
807#############
808
809module TauTagging TauTagging {
810 set ParticleInputArray Delphes/allParticles
811 set PartonInputArray Delphes/partons
812 set JetInputArray JetEnergyScale/jets
813
814 set DeltaR 0.5
815
816 set TauPTMin 1.0
817
818 set TauEtaMax 2.5
819
820 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
821
822 # default efficiency formula (misidentification rate)
823 add EfficiencyFormula {0} {0.01}
824 # efficiency formula for tau-jets
825 add EfficiencyFormula {15} {0.6}
826}
827
828#####################################################
829# Find uniquely identified photons/electrons/tau/jets
830#####################################################
831
832module UniqueObjectFinder UniqueObjectFinder {
833# earlier arrays take precedence over later ones
834# add InputArray InputArray OutputArray
835 add InputArray PhotonIsolation/photons photons
836 add InputArray ElectronIsolation/electrons electrons
837 add InputArray MuonIsolation/muons muons
838 add InputArray JetEnergyScale/jets jets
839}
840
841##################
842# ROOT tree writer
843##################
844
845# tracks, towers and eflow objects are not stored by default in the output.
846# if needed (for jet constituent or other studies), uncomment the relevant
847# "add Branch ..." lines.
848
849module TreeWriter TreeWriter {
850# add Branch InputArray BranchName BranchClass
851 add Branch Delphes/allParticles Particle GenParticle
852
853# add Branch TrackMerger/tracks Track Track
854# add Branch Calorimeter/towers Tower Tower
855
856# add Branch Calorimeter/eflowTracks EFlowTrack Track
857# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
858# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
859
860 add Branch GenJetFinder/jets GenJet Jet
861 add Branch GenMissingET/momentum GenMissingET MissingET
862
863 add Branch UniqueObjectFinder/jets Jet Jet
864 add Branch UniqueObjectFinder/electrons Electron Electron
865 add Branch UniqueObjectFinder/photons Photon Photon
866 add Branch UniqueObjectFinder/muons Muon Muon
867 add Branch MissingET/momentum MissingET MissingET
868 add Branch ScalarHT/energy ScalarHT ScalarHT
869 add Branch Rho/rho Rho Rho
870 add Branch PileUpMerger/vertices Vertex Vertex
871}
Note: See TracBrowser for help on using the repository browser.