Fork me on GitHub

source: git/examples/delphes_card_CMS_PileUp.tcl@ b25d4cf

ImprovedOutputFile Timing dual_readout llp
Last change on this file since b25d4cf was 1c8d9db, checked in by Michele <michele.selvaggi@…>, 10 years ago

added new parameters to Calorimeter module
  • Property mode set to 100644
File size: 23.2 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 EcalTowerMinEnergy 0.5
258 set HcalTowerMinEnergy 1.0
259
260 set EcalTowerMinSignificance 1.0
261 set HcalTowerMinSignificance 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# Rho pile-up
384#############
385
386module FastJetFinder Rho {
387# set InputArray Calorimeter/towers
388 set InputArray EFlowMerger/eflow
389
390 set ComputeRho true
391 set RhoOutputArray rho
392
393 # 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
394 set AreaAlgorithm 5
395
396 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
397 set JetAlgorithm 4
398 set ParameterR 0.6
399 set GhostEtaMax 5.0
400
401 add RhoEtaRange 0.0 2.5
402 add RhoEtaRange 2.5 5.0
403
404 set JetPTMin 0.0
405}
406
407#####################
408# Neutrino Filter
409#####################
410
411module PdgCodeFilter NeutrinoFilter {
412
413 set InputArray Delphes/stableParticles
414 set OutputArray filteredParticles
415
416 set PTMin 0.0
417
418 add PdgCode {12}
419 add PdgCode {14}
420 add PdgCode {16}
421 add PdgCode {-12}
422 add PdgCode {-14}
423 add PdgCode {-16}
424
425}
426
427
428
429#####################
430# MC truth jet finder
431#####################
432
433module FastJetFinder GenJetFinder {
434 set InputArray NeutrinoFilter/filteredParticles
435
436 set OutputArray jets
437
438 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
439 set JetAlgorithm 6
440 set ParameterR 0.5
441
442 set JetPTMin 20.0
443}
444
445############
446# Jet finder
447############
448
449module FastJetFinder FastJetFinder {
450# set InputArray Calorimeter/towers
451 set InputArray EFlowMerger/eflow
452
453 set OutputArray jets
454
455 # 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
456 set AreaAlgorithm 5
457
458 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
459 set JetAlgorithm 6
460 set ParameterR 0.5
461
462 set JetPTMin 20.0
463}
464
465###########################
466# Jet Pile-Up ID
467###########################
468
469module PileUpJetID PileUpJetID {
470 set JetInputArray FastJetFinder/jets
471 set TrackInputArray Calorimeter/eflowTracks
472 set NeutralInputArray NeutralTowerMerger/eflowTowers
473
474 set VertexInputArray PileUpMerger/vertices
475 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
476 # Z vertex resolution in m
477 set ZVertexResolution 0.0001
478
479 set OutputArray jets
480
481 set UseConstituents 0
482 set ParameterR 0.5
483
484 set JetPTMin 20.0
485}
486
487###########################
488# Jet Pile-Up Subtraction
489###########################
490
491module JetPileUpSubtractor JetPileUpSubtractor {
492 set JetInputArray PileUpJetID/jets
493 set RhoInputArray Rho/rho
494
495 set OutputArray jets
496
497 set JetPTMin 20.0
498}
499
500##################
501# Jet Energy Scale
502##################
503
504module EnergyScale JetEnergyScale {
505 set InputArray JetPileUpSubtractor/jets
506 set OutputArray jets
507
508 # scale formula for jets
509 set ScaleFormula {1.0}
510}
511
512###################
513# Photon efficiency
514###################
515
516module Efficiency PhotonEfficiency {
517 set InputArray Calorimeter/eflowPhotons
518 set OutputArray photons
519
520 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
521
522 # efficiency formula for photons
523 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
524 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
525 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
526 (abs(eta) > 2.5) * (0.00)}
527}
528
529
530##################
531# Photon isolation
532##################
533
534module Isolation PhotonIsolation {
535 set CandidateInputArray PhotonEfficiency/photons
536 set IsolationInputArray EFlowMerger/eflow
537 set RhoInputArray Rho/rho
538
539 set OutputArray photons
540
541 set DeltaRMax 0.5
542
543 set PTMin 0.5
544
545 set PTRatioMax 0.1
546}
547
548#####################
549# Electron efficiency
550#####################
551
552module Efficiency ElectronEfficiency {
553 set InputArray TrackPileUpSubtractor/electrons
554 set OutputArray electrons
555
556 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
557
558 # efficiency formula for electrons
559 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
560 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
561 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
562 (abs(eta) > 2.5) * (0.00)}
563}
564
565####################
566# Electron isolation
567####################
568
569module Isolation ElectronIsolation {
570 set CandidateInputArray ElectronEfficiency/electrons
571 set IsolationInputArray EFlowMerger/eflow
572 set RhoInputArray Rho/rho
573
574 set OutputArray electrons
575
576 set DeltaRMax 0.5
577
578 set PTMin 0.5
579
580 set PTRatioMax 0.1
581}
582
583#################
584# Muon efficiency
585#################
586
587module Efficiency MuonEfficiency {
588 set InputArray TrackPileUpSubtractor/muons
589 set OutputArray muons
590
591 # set EfficiencyFormula {efficiency as a function of eta and pt}
592
593 # efficiency formula for muons
594 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
595 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
596 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
597 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
598 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) + \
599 (abs(eta) > 2.4) * (0.00)}
600}
601
602################
603# Muon isolation
604################
605
606module Isolation MuonIsolation {
607 set CandidateInputArray MuonEfficiency/muons
608 set IsolationInputArray EFlowMerger/eflow
609 set RhoInputArray Rho/rho
610
611 set OutputArray muons
612
613 set DeltaRMax 0.5
614
615 set PTMin 0.5
616
617 set PTRatioMax 0.1
618}
619
620###################
621# Missing ET merger
622###################
623
624module Merger MissingET {
625# add InputArray InputArray
626 add InputArray EFlowMergerAllTracks/eflow
627 set MomentumOutputArray momentum
628}
629
630
631
632##################
633# Scalar HT merger
634##################
635
636module Merger ScalarHT {
637# add InputArray InputArray
638 add InputArray UniqueObjectFinder/jets
639 add InputArray UniqueObjectFinder/electrons
640 add InputArray UniqueObjectFinder/photons
641 add InputArray UniqueObjectFinder/muons
642 set EnergyOutputArray energy
643}
644
645###########
646# b-tagging
647###########
648
649module BTagging BTagging {
650 set PartonInputArray Delphes/partons
651 set JetInputArray JetEnergyScale/jets
652
653 set BitNumber 0
654
655 set DeltaR 0.5
656
657 set PartonPTMin 1.0
658
659 set PartonEtaMax 2.5
660
661 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
662 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
663 # gluon's PDG code has the lowest priority
664
665 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
666 # default efficiency formula (misidentification rate)
667 add EfficiencyFormula {0} {0.001}
668
669 # efficiency formula for c-jets (misidentification rate)
670 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
671 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
672 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
673 (abs(eta) > 2.5) * (0.000)}
674
675 # efficiency formula for b-jets
676 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
677 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
678 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
679 (abs(eta) > 2.5) * (0.000)}
680}
681
682module TauTagging TauTagging {
683 set ParticleInputArray Delphes/allParticles
684 set PartonInputArray Delphes/partons
685 set JetInputArray JetEnergyScale/jets
686
687 set DeltaR 0.5
688
689 set TauPTMin 1.0
690
691 set TauEtaMax 2.5
692
693 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
694
695 # default efficiency formula (misidentification rate)
696 add EfficiencyFormula {0} {0.001}
697 # efficiency formula for tau-jets
698 add EfficiencyFormula {15} {0.4}
699}
700
701#####################################################
702# Find uniquely identified photons/electrons/tau/jets
703#####################################################
704
705module UniqueObjectFinder UniqueObjectFinder {
706# earlier arrays take precedence over later ones
707# add InputArray InputArray OutputArray
708 add InputArray PhotonIsolation/photons photons
709 add InputArray ElectronIsolation/electrons electrons
710 add InputArray MuonIsolation/muons muons
711 add InputArray JetEnergyScale/jets jets
712}
713
714##################
715# ROOT tree writer
716##################
717
718# tracks, towers and eflow objects are not stored by default in the output.
719# if needed (for jet constituent or other studies), uncomment the relevant
720# "add Branch ..." lines.
721
722module TreeWriter TreeWriter {
723# add Branch InputArray BranchName BranchClass
724 add Branch Delphes/allParticles Particle GenParticle
725
726# add Branch TrackMerger/tracks Track Track
727# add Branch Calorimeter/towers Tower Tower
728
729# add Branch Calorimeter/eflowTracks EFlowTrack Track
730# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
731# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
732
733 add Branch GenJetFinder/jets GenJet Jet
734 add Branch UniqueObjectFinder/jets Jet Jet
735 add Branch UniqueObjectFinder/electrons Electron Electron
736 add Branch UniqueObjectFinder/photons Photon Photon
737 add Branch UniqueObjectFinder/muons Muon Muon
738 add Branch MissingET/momentum MissingET MissingET
739 add Branch ScalarHT/energy ScalarHT ScalarHT
740 add Branch Rho/rho Rho Rho
741 add Branch PileUpMerger/vertices Vertex Vertex
742}
743
Note: See TracBrowser for help on using the repository browser.