Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ 914a1f5

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

fix rapidity ranges for GridMedianBackgroundEstimator
  • Property mode set to 100644
File size: 23.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 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 ECalEnergyMin 0.5
258 set HCalEnergyMin 1.0
259
260 set ECalEnergySignificanceMin 1.0
261 set HCalEnergySignificanceMin 1.0
262
263 set EFlowTrackOutputArray eflowTracks
264 set EFlowPhotonOutputArray eflowPhotons
265 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
266
267 set pi [expr {acos(-1)}]
268
269 # lists of the edges of each tower in eta and phi
270 # each list starts with the lower edge of the first tower
271 # the list ends with the higher edged of the last tower
272
273 # 5 degrees towers
274 set PhiBins {}
275 for {set i -36} {$i <= 36} {incr i} {
276 add PhiBins [expr {$i * $pi/36.0}]
277 }
278 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} {
279 add EtaPhiBins $eta $PhiBins
280 }
281
282 # 10 degrees towers
283 set PhiBins {}
284 for {set i -18} {$i <= 18} {incr i} {
285 add PhiBins [expr {$i * $pi/18.0}]
286 }
287 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} {
288 add EtaPhiBins $eta $PhiBins
289 }
290
291 # 20 degrees towers
292 set PhiBins {}
293 for {set i -9} {$i <= 9} {incr i} {
294 add PhiBins [expr {$i * $pi/9.0}]
295 }
296 foreach eta {-5 -4.7 -4.525 4.7 5} {
297 add EtaPhiBins $eta $PhiBins
298 }
299
300 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
301 add EnergyFraction {0} {0.0 1.0}
302 # energy fractions for e, gamma and pi0
303 add EnergyFraction {11} {1.0 0.0}
304 add EnergyFraction {22} {1.0 0.0}
305 add EnergyFraction {111} {1.0 0.0}
306 # energy fractions for muon, neutrinos and neutralinos
307 add EnergyFraction {12} {0.0 0.0}
308 add EnergyFraction {13} {0.0 0.0}
309 add EnergyFraction {14} {0.0 0.0}
310 add EnergyFraction {16} {0.0 0.0}
311 add EnergyFraction {1000022} {0.0 0.0}
312 add EnergyFraction {1000023} {0.0 0.0}
313 add EnergyFraction {1000025} {0.0 0.0}
314 add EnergyFraction {1000035} {0.0 0.0}
315 add EnergyFraction {1000045} {0.0 0.0}
316 # energy fractions for K0short and Lambda
317 add EnergyFraction {310} {0.3 0.7}
318 add EnergyFraction {3122} {0.3 0.7}
319
320 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
321 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
322 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
323
324 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
325 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
326 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
327}
328
329##########################
330# Track pile-up subtractor
331##########################
332
333module TrackPileUpSubtractor TrackPileUpSubtractor {
334# add InputArray InputArray OutputArray
335 add InputArray Calorimeter/eflowTracks eflowTracks
336 add InputArray ElectronEnergySmearing/electrons electrons
337 add InputArray MuonMomentumSmearing/muons muons
338
339 set VertexInputArray PileUpMerger/vertices
340 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
341 # Z vertex resolution in m
342 set ZVertexResolution 0.0001
343}
344
345####################
346# Neutral tower merger
347####################
348
349module Merger NeutralTowerMerger {
350# add InputArray InputArray
351 add InputArray Calorimeter/eflowPhotons
352 add InputArray Calorimeter/eflowNeutralHadrons
353 set OutputArray eflowTowers
354}
355
356##################################
357# Energy flow merger (all tracks)
358##################################
359
360module Merger EFlowMergerAllTracks {
361# add InputArray InputArray
362 add InputArray TrackMerger/tracks
363 add InputArray Calorimeter/eflowPhotons
364 add InputArray Calorimeter/eflowNeutralHadrons
365 set OutputArray eflow
366}
367
368
369####################
370# Energy flow merger
371####################
372
373module Merger EFlowMerger {
374# add InputArray InputArray
375 add InputArray TrackPileUpSubtractor/eflowTracks
376 add InputArray Calorimeter/eflowPhotons
377 add InputArray Calorimeter/eflowNeutralHadrons
378 set OutputArray eflow
379}
380
381
382
383#############
384# Rho pile-up
385#############
386
387module FastJetGridMedianEstimator Rho {
388
389 set InputArray EFlowMerger/eflow
390 set RhoOutputArray rho
391
392 # add GridRange rapmin rapmax drap dphi
393 # rapmin - the minimum rapidity extent of the grid
394 # rapmax - the maximum rapidity extent of the grid
395 # drap - the grid spacing in rapidity
396 # dphi - the grid spacing in azimuth
397
398 add GridRange -5.0 -2.5 1.0 1.0
399 add GridRange -2.5 2.5 1.0 1.0
400 add GridRange 2.5 5.0 1.0 1.0
401
402}
403
404#####################
405# Neutrino Filter
406#####################
407
408module PdgCodeFilter NeutrinoFilter {
409
410 set InputArray Delphes/stableParticles
411 set OutputArray filteredParticles
412
413 set PTMin 0.0
414
415 add PdgCode {12}
416 add PdgCode {14}
417 add PdgCode {16}
418 add PdgCode {-12}
419 add PdgCode {-14}
420 add PdgCode {-16}
421
422}
423
424
425
426#####################
427# MC truth jet finder
428#####################
429
430module FastJetFinder GenJetFinder {
431 set InputArray NeutrinoFilter/filteredParticles
432
433 set OutputArray jets
434
435 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
436 set JetAlgorithm 6
437 set ParameterR 0.5
438
439 set JetPTMin 20.0
440}
441
442############
443# Jet finder
444############
445
446module FastJetFinder FastJetFinder {
447# set InputArray Calorimeter/towers
448 set InputArray EFlowMerger/eflow
449
450 set OutputArray jets
451
452 # 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
453 set AreaAlgorithm 5
454
455 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
456 set JetAlgorithm 6
457 set ParameterR 0.5
458
459 set JetPTMin 20.0
460}
461
462###########################
463# Jet Pile-Up ID
464###########################
465
466module PileUpJetID PileUpJetID {
467 set JetInputArray FastJetFinder/jets
468 set TrackInputArray Calorimeter/eflowTracks
469 set NeutralInputArray NeutralTowerMerger/eflowTowers
470
471 set VertexInputArray PileUpMerger/vertices
472 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
473 # Z vertex resolution in m
474 set ZVertexResolution 0.0001
475
476 set OutputArray jets
477
478 set UseConstituents 0
479 set ParameterR 0.5
480
481 set JetPTMin 20.0
482}
483
484###########################
485# Jet Pile-Up Subtraction
486###########################
487
488module JetPileUpSubtractor JetPileUpSubtractor {
489 set JetInputArray PileUpJetID/jets
490 set RhoInputArray Rho/rho
491
492 set OutputArray jets
493
494 set JetPTMin 20.0
495}
496
497##################
498# Jet Energy Scale
499##################
500
501module EnergyScale JetEnergyScale {
502 set InputArray JetPileUpSubtractor/jets
503 set OutputArray jets
504
505 # scale formula for jets
506 set ScaleFormula {1.0}
507}
508
509###################
510# Photon efficiency
511###################
512
513module Efficiency PhotonEfficiency {
514 set InputArray Calorimeter/eflowPhotons
515 set OutputArray photons
516
517 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
518
519 # efficiency formula for photons
520 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
521 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
522 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
523 (abs(eta) > 2.5) * (0.00)}
524}
525
526
527##################
528# Photon isolation
529##################
530
531module Isolation PhotonIsolation {
532 set CandidateInputArray PhotonEfficiency/photons
533 set IsolationInputArray EFlowMerger/eflow
534 set RhoInputArray Rho/rho
535
536 set OutputArray photons
537
538 set DeltaRMax 0.5
539
540 set PTMin 0.5
541
542 set PTRatioMax 0.1
543}
544
545#####################
546# Electron efficiency
547#####################
548
549module Efficiency ElectronEfficiency {
550 set InputArray TrackPileUpSubtractor/electrons
551 set OutputArray electrons
552
553 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
554
555 # efficiency formula for electrons
556 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
557 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
558 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
559 (abs(eta) > 2.5) * (0.00)}
560}
561
562####################
563# Electron isolation
564####################
565
566module Isolation ElectronIsolation {
567 set CandidateInputArray ElectronEfficiency/electrons
568 set IsolationInputArray EFlowMerger/eflow
569 set RhoInputArray Rho/rho
570
571 set OutputArray electrons
572
573 set DeltaRMax 0.5
574
575 set PTMin 0.5
576
577 set PTRatioMax 0.1
578}
579
580#################
581# Muon efficiency
582#################
583
584module Efficiency MuonEfficiency {
585 set InputArray TrackPileUpSubtractor/muons
586 set OutputArray muons
587
588 # set EfficiencyFormula {efficiency as a function of eta and pt}
589
590 # efficiency formula for muons
591 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
592 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
593 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
594 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
595 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
596 (abs(eta) > 2.4) * (0.00)}
597}
598
599################
600# Muon isolation
601################
602
603module Isolation MuonIsolation {
604 set CandidateInputArray MuonEfficiency/muons
605 set IsolationInputArray EFlowMerger/eflow
606 set RhoInputArray Rho/rho
607
608 set OutputArray muons
609
610 set DeltaRMax 0.5
611
612 set PTMin 0.5
613
614 set PTRatioMax 0.1
615}
616
617###################
618# Missing ET merger
619###################
620
621module Merger MissingET {
622# add InputArray InputArray
623 add InputArray EFlowMergerAllTracks/eflow
624 set MomentumOutputArray momentum
625}
626
627
628
629##################
630# Scalar HT merger
631##################
632
633module Merger ScalarHT {
634# add InputArray InputArray
635 add InputArray UniqueObjectFinder/jets
636 add InputArray UniqueObjectFinder/electrons
637 add InputArray UniqueObjectFinder/photons
638 add InputArray UniqueObjectFinder/muons
639 set EnergyOutputArray energy
640}
641
642###########
643# b-tagging
644###########
645
646module BTagging BTagging {
647 set PartonInputArray Delphes/partons
648 set JetInputArray JetEnergyScale/jets
649
650 set BitNumber 0
651
652 set DeltaR 0.5
653
654 set PartonPTMin 1.0
655
656 set PartonEtaMax 2.5
657
658 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
659 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
660 # gluon's PDG code has the lowest priority
661
662 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
663 # default efficiency formula (misidentification rate)
664 add EfficiencyFormula {0} {0.001}
665
666 # efficiency formula for c-jets (misidentification rate)
667 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
668 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
669 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
670 (abs(eta) > 2.5) * (0.000)}
671
672 # efficiency formula for b-jets
673 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
674 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
675 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
676 (abs(eta) > 2.5) * (0.000)}
677}
678
679module TauTagging TauTagging {
680 set ParticleInputArray Delphes/allParticles
681 set PartonInputArray Delphes/partons
682 set JetInputArray JetEnergyScale/jets
683
684 set DeltaR 0.5
685
686 set TauPTMin 1.0
687
688 set TauEtaMax 2.5
689
690 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
691
692 # default efficiency formula (misidentification rate)
693 add EfficiencyFormula {0} {0.001}
694 # efficiency formula for tau-jets
695 add EfficiencyFormula {15} {0.4}
696}
697
698#####################################################
699# Find uniquely identified photons/electrons/tau/jets
700#####################################################
701
702module UniqueObjectFinder UniqueObjectFinder {
703# earlier arrays take precedence over later ones
704# add InputArray InputArray OutputArray
705 add InputArray PhotonIsolation/photons photons
706 add InputArray ElectronIsolation/electrons electrons
707 add InputArray MuonIsolation/muons muons
708 add InputArray JetEnergyScale/jets jets
709}
710
711##################
712# ROOT tree writer
713##################
714
715# tracks, towers and eflow objects are not stored by default in the output.
716# if needed (for jet constituent or other studies), uncomment the relevant
717# "add Branch ..." lines.
718
719module TreeWriter TreeWriter {
720# add Branch InputArray BranchName BranchClass
721 add Branch Delphes/allParticles Particle GenParticle
722
723# add Branch TrackMerger/tracks Track Track
724# add Branch Calorimeter/towers Tower Tower
725
726# add Branch Calorimeter/eflowTracks EFlowTrack Track
727# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
728# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
729
730 add Branch GenJetFinder/jets GenJet Jet
731 add Branch UniqueObjectFinder/jets Jet Jet
732 add Branch UniqueObjectFinder/electrons Electron Electron
733 add Branch UniqueObjectFinder/photons Photon Photon
734 add Branch UniqueObjectFinder/muons Muon Muon
735 add Branch MissingET/momentum MissingET MissingET
736 add Branch ScalarHT/energy ScalarHT ScalarHT
737 add Branch Rho/rho Rho Rho
738 add Branch PileUpMerger/vertices Vertex Vertex
739}
740
Note: See TracBrowser for help on using the repository browser.