Fork me on GitHub

Ticket #287: delphes_card_Snowmass_NoPileUp_v2.tcl

File delphes_card_Snowmass_NoPileUp_v2.tcl, 22.1 KB (added by Michele Selvaggi, 10 years ago)
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7 StatusPid
8
9 ChargedHadronTrackingEfficiency
10 ElectronTrackingEfficiency
11 MuonTrackingEfficiency
12
13 ChargedHadronMomentumSmearing
14 ElectronEnergySmearing
15 MuonMomentumSmearing
16
17 TrackMerger
18 Calorimeter
19 TowerMerger
20 EFlowMerger
21
22 PhotonEfficiency
23 PhotonIsolation
24
25 ElectronEfficiency
26 ElectronIsolation
27
28 MuonEfficiency
29 MuonIsolation
30
31 MissingET
32
33 GenJetFinder
34 FastJetFinder
35 CAJetFinder
36
37 UniqueObjectFinderGJ
38 UniqueObjectFinderEJ
39 UniqueObjectFinderMJ
40
41 ConstituentFilter
42
43 BTagging
44 BTaggingLoose
45 TauTagging
46
47 ScalarHT
48
49 TreeWriter
50}
51
52#################################
53# Propagate particles in cylinder
54#################################
55
56module ParticlePropagator ParticlePropagator {
57 set InputArray Delphes/stableParticles
58
59 set OutputArray stableParticles
60 set ChargedHadronOutputArray chargedHadrons
61 set ElectronOutputArray electrons
62 set MuonOutputArray muons
63
64 # radius of the magnetic field coverage, in m
65 set Radius 1.29
66 # half-length of the magnetic field coverage, in m
67 set HalfLength 3.00
68
69 # magnetic field
70 set Bz 3.8
71}
72
73####################################
74# StatusPidFilter
75# This module removes all generated particles
76# except electrons, muons, taus, and status==3
77####################################
78
79module StatusPidFilter StatusPid {
80# set InputArray Delphes/stableParticles
81 set InputArray Delphes/allParticles
82 set OutputArray filteredParticles
83
84 set PTMin 0.5
85}
86
87
88####################################
89# Charged hadron tracking efficiency
90####################################
91
92module Efficiency ChargedHadronTrackingEfficiency {
93 set InputArray ParticlePropagator/chargedHadrons
94 set OutputArray chargedHadrons
95
96 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
97 # tracking efficiency formula for charged hadrons
98 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
99 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.85) + \
100 (abs(eta) <= 1.5) * (pt > 1.0) * (0.97) + \
101 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.85) + \
102 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90) + \
103 (abs(eta) > 2.5) * (0.00)}
104
105}
106
107##############################
108# Electron tracking efficiency
109##############################
110
111module Efficiency ElectronTrackingEfficiency {
112 set InputArray ParticlePropagator/electrons
113 set OutputArray electrons
114
115 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
116 # tracking efficiency formula for electrons
117 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
118 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.85) + \
119 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.97) + \
120 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
121 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.85) + \
122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.90) + \
123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.95) + \
124 (abs(eta) > 2.5) * (0.00)}
125
126}
127
128##########################
129# Muon tracking efficiency
130##########################
131
132module Efficiency MuonTrackingEfficiency {
133 set InputArray ParticlePropagator/muons
134 set OutputArray muons
135
136 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
137
138 # tracking efficiency formula for muons
139 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
140 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.998) + \
141 (abs(eta) <= 1.5) * (pt > 1.0) * (0.9998) + \
142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.98) + \
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
144 (abs(eta) > 2.5) * (0.00)}
145}
146
147########################################
148# Momentum resolution for charged tracks
149########################################
150
151module MomentumSmearing ChargedHadronMomentumSmearing {
152 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
153 set OutputArray chargedHadrons
154
155 # set ResolutionFormula {resolution formula as a function of eta and pt}
156
157 # resolution formula for charged hadrons
158 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
159 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
160 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
161 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
162 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
163 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
164 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
166
167}
168
169#################################
170# Energy resolution for electrons
171#################################
172
173module EnergySmearing ElectronEnergySmearing {
174 set InputArray ElectronTrackingEfficiency/electrons
175 set OutputArray electrons
176
177 # set ResolutionFormula {resolution formula as a function of eta and energy}
178 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
179 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
180 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
181 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.97^2)}
182
183}
184
185###############################
186# Momentum resolution for muons
187###############################
188
189module MomentumSmearing MuonMomentumSmearing {
190 set InputArray MuonTrackingEfficiency/muons
191 set OutputArray muons
192
193 # set ResolutionFormula {resolution formula as a function of eta and pt}
194 # resolution formula for muons
195 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
196 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
197 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
198 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
199 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
200 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
201 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
202 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
203
204}
205
206##############
207# Track merger
208##############
209
210module Merger TrackMerger {
211# add InputArray InputArray
212 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
213 add InputArray ElectronEnergySmearing/electrons
214 set OutputArray tracks
215}
216
217#############
218# Calorimeter
219#############
220
221module Calorimeter Calorimeter {
222 set ParticleInputArray ParticlePropagator/stableParticles
223 set TrackInputArray TrackMerger/tracks
224
225 set TowerOutputArray towers
226 set PhotonOutputArray photons
227
228 set EFlowTrackOutputArray eflowTracks
229 set EFlowPhotonOutputArray eflowPhotons
230 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
231
232 set pi [expr {acos(-1)}]
233
234 # lists of the edges of each tower in eta and phi
235 # each list starts with the lower edge of the first tower
236 # the list ends with the higher edged of the last tower
237
238 # 5 degrees towers
239 set PhiBins {}
240 for {set i -36} {$i <= 36} {incr i} {
241 add PhiBins [expr {$i * $pi/36.0}]
242 }
243 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} {
244 add EtaPhiBins $eta $PhiBins
245 }
246
247 # 10 degrees towers
248 set PhiBins {}
249 for {set i -18} {$i <= 18} {incr i} {
250 add PhiBins [expr {$i * $pi/18.0}]
251 }
252 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} {
253 add EtaPhiBins $eta $PhiBins
254 }
255
256 # 20 degrees towers
257 set PhiBins {}
258 for {set i -9} {$i <= 9} {incr i} {
259 add PhiBins [expr {$i * $pi/9.0}]
260 }
261 foreach eta {-5 -4.7 -4.525 4.7 5} {
262 add EtaPhiBins $eta $PhiBins
263 }
264
265 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
266 add EnergyFraction {0} {0.0 1.0}
267 # energy fractions for e, gamma and pi0
268 add EnergyFraction {11} {1.0 0.0}
269 add EnergyFraction {22} {1.0 0.0}
270 add EnergyFraction {111} {1.0 0.0}
271 # energy fractions for muon, neutrinos and neutralinos
272 add EnergyFraction {12} {0.0 0.0}
273 add EnergyFraction {13} {0.0 0.0}
274 add EnergyFraction {14} {0.0 0.0}
275 add EnergyFraction {16} {0.0 0.0}
276 add EnergyFraction {1000022} {0.0 0.0}
277 add EnergyFraction {1000023} {0.0 0.0}
278 add EnergyFraction {1000025} {0.0 0.0}
279 add EnergyFraction {1000035} {0.0 0.0}
280 add EnergyFraction {1000045} {0.0 0.0}
281 # energy fractions for K0short and Lambda
282 add EnergyFraction {310} {0.3 0.7}
283 add EnergyFraction {3122} {0.3 0.7}
284
285 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
286 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
287 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.97^2)}
288
289 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
290 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
291 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.050^2 + energy*0.706^2) + \
292 (abs(eta) > 3.0 && abs(eta) <= 4.9) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
293}
294
295####################
296# Tower Merger (in case not using e-flow algorithm)
297####################
298
299module Merger TowerMerger {
300# add InputArray InputArray
301 add InputArray Calorimeter/eflowPhotons
302 add InputArray Calorimeter/eflowNeutralHadrons
303 set OutputArray towers
304}
305
306
307
308####################
309# Energy flow merger
310####################
311
312module Merger EFlowMerger {
313# add InputArray InputArray
314 add InputArray Calorimeter/eflowTracks
315 add InputArray TowerMerger/towers
316 add InputArray MuonMomentumSmearing/muons
317 set OutputArray eflow
318}
319
320###################
321# Photon efficiency
322###################
323
324module Efficiency PhotonEfficiency {
325 set InputArray Calorimeter/photons
326 set OutputArray photons
327
328 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
329 # efficiency formula for photons
330 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
331 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
332 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.9624) + \
333 (abs(eta) > 2.5) * (0.00)}
334}
335
336##################
337# Photon isolation
338##################
339
340module Isolation PhotonIsolation {
341 set CandidateInputArray PhotonEfficiency/photons
342 set IsolationInputArray EFlowMerger/eflow
343
344 set OutputArray photons
345
346 set DeltaRMax 0.3
347
348 set PTMin 1.0
349
350 set PTRatioMax 0.1
351}
352
353#####################
354# Electron efficiency
355#####################
356
357module Efficiency ElectronEfficiency {
358 set InputArray ElectronEnergySmearing/electrons
359 set OutputArray electrons
360
361 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
362 # efficiency formula for electrons
363 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
364 (abs(eta) <= 1.5) * (pt > 10.0) * (0.98) + \
365 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
366 (abs(eta) > 2.5) * (0.00)}
367}
368
369####################
370# Electron isolation
371####################
372
373module Isolation ElectronIsolation {
374 set CandidateInputArray ElectronEfficiency/electrons
375 set IsolationInputArray EFlowMerger/eflow
376
377 set OutputArray electrons
378
379 set DeltaRMax 0.3
380
381 set PTMin 0.5
382
383 set PTRatioMax 0.1
384}
385
386#################
387# Muon efficiency
388#################
389
390module Efficiency MuonEfficiency {
391 set InputArray MuonMomentumSmearing/muons
392 set OutputArray muons
393
394 # set EfficiencyFormula {efficiency as a function of eta and pt}
395 # efficiency formula for muons
396 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
397 (abs(eta) <= 1.5) * (pt > 10.0) * (0.99) + \
398 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.97) + \
399 (abs(eta) > 2.4) * (0.00)}
400
401}
402
403################
404# Muon isolation
405################
406
407module Isolation MuonIsolation {
408 set CandidateInputArray MuonEfficiency/muons
409 set IsolationInputArray EFlowMerger/eflow
410
411 set OutputArray muons
412
413 set DeltaRMax 0.3
414
415 set PTMin 0.5
416
417 set PTRatioMax 0.1
418}
419
420###################
421# Missing ET merger
422###################
423
424module Merger MissingET {
425# add InputArray InputArray
426 add InputArray EFlowMerger/eflow
427 set MomentumOutputArray momentum
428}
429
430##################
431# Scalar HT merger
432##################
433
434module Merger ScalarHT {
435# add InputArray InputArray
436 add InputArray UniqueObjectFinderMJ/jets
437 add InputArray UniqueObjectFinderEJ/electrons
438 add InputArray UniqueObjectFinderGJ/photons
439 add InputArray UniqueObjectFinderMJ/muons
440 set EnergyOutputArray energy
441}
442
443#####################
444# MC truth jet finder
445#####################
446
447module FastJetFinder GenJetFinder {
448 set InputArray Delphes/stableParticles
449
450 set OutputArray jets
451
452 # 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 ConeRadius 0.7
457 set SeedThreshold 1.0
458 set ConeAreaFraction 1.0
459 set AdjacencyCut 2
460 set OverlapThreshold 0.75
461
462 set MaxIterations 100
463 set MaxPairSize 2
464 set Iratch 1
465
466 set JetPTMin 20.0
467}
468
469############
470# Jet finder
471############
472
473module FastJetFinder FastJetFinder {
474# set InputArray TowerMerger/towers
475 set InputArray EFlowMerger/eflow
476
477 set OutputArray jets
478
479 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
480 set JetAlgorithm 6
481 set ParameterR 0.5
482
483 set ConeRadius 0.5
484 set SeedThreshold 1.0
485 set ConeAreaFraction 1.0
486 set AdjacencyCut 2
487 set OverlapThreshold 0.75
488
489 set MaxIterations 100
490 set MaxPairSize 2
491 set Iratch 1
492
493 set JetPTMin 20.0
494}
495
496
497############
498# Cambridge-Aachen Jet finder
499############
500
501module FastJetFinder CAJetFinder {
502# set InputArray TowerMerger/towers
503 set InputArray EFlowMerger/eflow
504 set OutputArray jets
505 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
506 set JetAlgorithm 5
507 set ParameterR 0.8
508 # 200 GeV needed for boosted W bosons, 300 GeV is safe for boosted tops
509 set JetPTMin 200.0
510}
511
512
513####################
514# Constituent filter
515####################
516
517module ConstituentFilter ConstituentFilter {
518
519# add JetInputArray InputArray
520# add JetInputArray GenJetFinder/jets
521 add JetInputArray CAJetFinder/jets
522
523# add ConstituentInputArray InputArray OutputArray
524# add ConstituentInputArray Delphes/stableParticles stableParticles
525 add ConstituentInputArray Calorimeter/eflowTracks eflowTracks
526 add ConstituentInputArray Calorimeter/eflowPhotons eflowPhotons
527 add ConstituentInputArray Calorimeter/eflowNeutralHadrons eflowNeutralHadrons
528 add ConstituentInputArray MuonMomentumSmearing/muons muons
529}
530
531###########
532# b-tagging
533###########
534
535module BTagging BTagging {
536 set PartonInputArray Delphes/partons
537 set JetInputArray FastJetFinder/jets
538
539 set BitNumber 0
540 set DeltaR 0.5
541 set PartonPTMin 1.0
542 set PartonEtaMax 2.5
543
544 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
545 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
546 # gluon's PDG code has the lowest priority
547
548 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
549 # default efficiency formula (misidentification rate)
550 add EfficiencyFormula {0} {0.001}
551
552 # efficiency formula for c-jets (misidentification rate)
553 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
554 (abs(eta) <= 1.2) * (pt > 15.0) * (0.1873*tanh(pt*0.0183 - 0.2196)) + \
555 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1898*tanh(pt*0.00997 - 0.143)) + \
556 (abs(eta) > 2.5) * (0.000)}
557
558 # efficiency formula for b-jets
559 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
560 (abs(eta) <= 1.2) * (pt > 15.0) * (0.7*tanh(pt*0.01317 - 0.062)) + \
561 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.6*tanh(pt*0.0105 - 0.101)) + \
562 (abs(eta) > 2.5) * (0.000)}
563}
564
565#### BTaggingLoose ###################
566
567module BTagging BTaggingLoose {
568 set PartonInputArray Delphes/partons
569 set JetInputArray FastJetFinder/jets
570
571 set BitNumber 1
572 set DeltaR 0.5
573 set PartonPTMin 1.0
574 set PartonEtaMax 2.5
575
576 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
577 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
578 # gluon's PDG code has the lowest priority
579
580 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
581 # default efficiency formula (misidentification rate)
582 add EfficiencyFormula {0} {0.03}
583
584 # efficiency formula for c-jets (misidentification rate)
585 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
586 (abs(eta) <= 1.2) * (pt > 15.0) * (0.29*tanh(pt*0.0183 - 0.2196)) + \
587 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.29*tanh(pt*0.00997 - 0.143)) + \
588 (abs(eta) > 2.5) * (0.000)}
589
590 # efficiency formula for b-jets
591 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
592 (abs(eta) <= 1.2) * (pt > 15.0) * (0.75*tanh(pt*0.01317 - 0.062)) + \
593 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.69*tanh(pt*0.0105 - 0.101)) + \
594 (abs(eta) > 2.5) * (0.000)}
595}
596
597
598module TauTagging TauTagging {
599 set ParticleInputArray Delphes/allParticles
600 set PartonInputArray Delphes/partons
601 set JetInputArray FastJetFinder/jets
602
603 set DeltaR 0.5
604
605 set TauPTMin 1.0
606
607 set TauEtaMax 2.5
608
609 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
610
611 # default efficiency formula (misidentification rate)
612 add EfficiencyFormula {0} {0.004}
613 # efficiency formula for tau-jets
614 add EfficiencyFormula {15} {0.65}
615}
616
617#####################################################
618# Find uniquely identified photons/electrons/tau/jets
619#####################################################
620
621#module UniqueObjectFinder UniqueObjectFinder {
622# earlier arrays take precedence over later ones
623# add InputArray InputArray OutputArray
624# add InputArray PhotonIsolation/photons photons
625# add InputArray ElectronIsolation/electrons electrons
626# add InputArray FastJetFinder/jets jets
627#}
628
629module UniqueObjectFinder UniqueObjectFinderGJ {
630 add InputArray PhotonIsolation/photons photons
631 add InputArray FastJetFinder/jets jets
632}
633
634module UniqueObjectFinder UniqueObjectFinderEJ {
635 add InputArray ElectronIsolation/electrons electrons
636 add InputArray UniqueObjectFinderGJ/jets jets
637}
638
639module UniqueObjectFinder UniqueObjectFinderMJ {
640 add InputArray MuonIsolation/muons muons
641 add InputArray UniqueObjectFinderEJ/jets jets
642}
643
644
645##################
646# ROOT tree writer
647##################
648
649module TreeWriter TreeWriter {
650# add Branch InputArray BranchName BranchClass
651# add Branch Delphes/allParticles Particle GenParticle
652 add Branch StatusPid/filteredParticles Particle GenParticle
653# add Branch TrackMerger/tracks Track Track
654# add Branch Calorimeter/towers Tower Tower
655 add Branch ConstituentFilter/eflowTracks EFlowTrack Track
656 add Branch ConstituentFilter/eflowPhotons EFlowPhoton Tower
657 add Branch ConstituentFilter/eflowNeutralHadrons EFlowNeutralHadron Tower
658 add Branch ConstituentFilter/muons EFlowMuon Muon
659 add Branch GenJetFinder/jets GenJet Jet
660 add Branch CAJetFinder/jets CAJet Jet
661 add Branch UniqueObjectFinderMJ/jets Jet Jet
662 add Branch UniqueObjectFinderEJ/electrons Electron Electron
663 add Branch UniqueObjectFinderGJ/photons Photon Photon
664 add Branch UniqueObjectFinderMJ/muons Muon Muon
665 add Branch MissingET/momentum MissingET MissingET
666 add Branch ScalarHT/energy ScalarHT ScalarHT
667}
668