Fork me on GitHub

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