Fork me on GitHub

source: git/examples/delphes_card_CMS_PileUp.tcl@ 4b9a2dc

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 4b9a2dc was 492aef9, checked in by Michele <michele.selvaggi@…>, 10 years ago

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