Fork me on GitHub

source: git/cards/delphes_card_CMS_PileUp.tcl@ 1343a07

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 1343a07 was 6153fb0, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

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