Fork me on GitHub

source: git/examples/delphes_card_CMS_PileUp.tcl@ 3d963d5

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 3d963d5 was 6fb1a5d, checked in by Michele <michele.selvaggi@…>, 10 years ago

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