Fork me on GitHub

source: git/examples/delphes_card_CMS_PileUp.tcl@ be2222c

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

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