Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ 532016b

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 532016b was a2983ec, checked in by Pavel Demin <pavel.demin@…>, 10 years ago

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