MA5SandBox: delphesMA5tune_card_ATLAS_05.tcl

File delphesMA5tune_card_ATLAS_05.tcl, 20.7 KB (added by Benjamin Fuks, 6 years ago)
Line 
1#######################################
2# MODIF OF STANDARD DELFES CARD #
3#######################################
4
5# => Medium electron 2013 efficiency implemented
6# => 2013 efficiency for muons implemented
7# => Efficiency for b-tag, c-jets, light jets taken
8# from ATLAS 0-leptons + 2b-jets (ATLAS SUSY 13 05)
9
10
11#######################################
12# Order of execution of various modules
13#######################################
14
15set ExecutionPath {
16 ParticlePropagator
17
18 ChargedHadronTrackingEfficiency
19 ElectronTrackingEfficiency
20 MuonTrackingEfficiency
21
22 ChargedHadronMomentumSmearing
23 ElectronEnergySmearing
24 MuonMomentumSmearing
25
26 TrackMerger
27 TrackIsolationCalculation
28 Calorimeter
29 EFlowMerger
30
31 PhotonEfficiency
32 PhotonIsolationCalculation
33 PhotonIsolation
34
35 ElectronEfficiency
36 ElectronIsolationCalculation
37 ElectronIsolation
38
39 MuonEfficiency
40 MuonIsolationCalculation
41 MuonIsolation
42
43 MissingET
44
45 GenJetFinder
46 FastJetFinder
47
48 JetEnergyScale
49
50 BTagging
51 TauTagging
52
53 UniqueObjectFinder
54
55 ScalarHT
56
57 TreeWriter
58}
59
60#MA5 BEGIN
61#################################
62# Isolation Calculation
63#################################
64module IsolationCalculation PhotonIsolationCalculation {
65 set CandidateInputArray PhotonEfficiency/photons
66 set IsolationInputArray EFlowMerger/eflow
67
68 set OutputArray DelphesMA5tunePhotons
69 set PTMin 0.5
70}
71
72module IsolationCalculation ElectronIsolationCalculation {
73 set CandidateInputArray ElectronEfficiency/electrons
74 set IsolationInputArray TrackMerger/tracks
75
76 set OutputArray DelphesMA5tuneElectrons
77 set PTMin 0.5
78}
79
80module IsolationCalculation MuonIsolationCalculation {
81 set CandidateInputArray MuonEfficiency/muons
82 set IsolationInputArray TrackMerger/tracks
83
84 set OutputArray DelphesMA5tuneMuons
85 set PTMin 0.5
86}
87module IsolationCalculation TrackIsolationCalculation {
88 set CandidateInputArray TrackMerger/tracks
89 set IsolationInputArray TrackMerger/tracks
90
91 set OutputArray DelphesMA5tuneTracks
92 set PTMin 0.5
93}
94#MA5 END
95
96#################################
97# Propagate particles in cylinder
98#################################
99
100module ParticlePropagator ParticlePropagator {
101 set InputArray Delphes/stableParticles
102
103 set OutputArray stableParticles
104 set ChargedHadronOutputArray chargedHadrons
105 set ElectronOutputArray electrons
106 set MuonOutputArray muons
107
108 # radius of the magnetic field coverage, in m
109 set Radius 1.15
110 # half-length of the magnetic field coverage, in m
111 set HalfLength 3.51
112
113 # magnetic field
114 set Bz 2.0
115}
116
117####################################
118# Charged hadron tracking efficiency
119####################################
120
121module Efficiency ChargedHadronTrackingEfficiency {
122 set InputArray ParticlePropagator/chargedHadrons
123 set OutputArray chargedHadrons
124
125 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
126
127 # tracking efficiency formula for charged hadrons
128 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
129 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
130 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
131 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
132 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
133 (abs(eta) > 2.5) * (0.00)}
134}
135
136##############################
137# Electron tracking efficiency
138##############################
139
140module Efficiency ElectronTrackingEfficiency {
141 set InputArray ParticlePropagator/electrons
142 set OutputArray electrons
143
144 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
145
146 # tracking efficiency formula for electrons
147 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
148 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
149 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
150 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
151 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
152 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
153 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
154 (abs(eta) > 2.5) * (0.00)}
155}
156
157##########################
158# Muon tracking efficiency
159##########################
160
161module Efficiency MuonTrackingEfficiency {
162 set InputArray ParticlePropagator/muons
163 set OutputArray muons
164
165 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
166
167 # tracking efficiency formula for muons
168 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
169 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
170 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
171 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
172 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
173 (abs(eta) > 2.5) * (0.00)}
174}
175
176########################################
177# Momentum resolution for charged tracks
178########################################
179
180module MomentumSmearing ChargedHadronMomentumSmearing {
181 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
182 set OutputArray chargedHadrons
183
184 # set ResolutionFormula {resolution formula as a function of eta and pt}
185
186 # resolution formula for charged hadrons
187 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
188 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
189 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
190 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
191 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
192 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
193 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
194 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
195}
196
197#################################
198# Energy resolution for electrons
199#################################
200
201module EnergySmearing ElectronEnergySmearing {
202 set InputArray ElectronTrackingEfficiency/electrons
203 set OutputArray electrons
204
205 # set ResolutionFormula {resolution formula as a function of eta and energy}
206
207 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
208 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
209 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
210 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
211
212}
213
214###############################
215# Momentum resolution for muons
216###############################
217
218module MomentumSmearing MuonMomentumSmearing {
219 set InputArray MuonTrackingEfficiency/muons
220 set OutputArray muons
221
222 # set ResolutionFormula {resolution formula as a function of eta and pt}
223
224 # resolution formula for muons
225 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
226 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \
227 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
228 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \
229 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \
230 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \
231 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
232 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)}
233}
234
235##############
236# Track merger
237##############
238
239module Merger TrackMerger {
240# add InputArray InputArray
241 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
242 add InputArray ElectronEnergySmearing/electrons
243 add InputArray MuonMomentumSmearing/muons
244 set OutputArray tracks
245}
246
247#############
248# Calorimeter
249#############
250
251module Calorimeter Calorimeter {
252 set ParticleInputArray ParticlePropagator/stableParticles
253 set TrackInputArray TrackMerger/tracks
254
255 set TowerOutputArray towers
256 set PhotonOutputArray photons
257
258 set EFlowTrackOutputArray eflowTracks
259 set EFlowPhotonOutputArray eflowPhotons
260 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
261
262 set pi [expr {acos(-1)}]
263
264 # lists of the edges of each tower in eta and phi
265 # each list starts with the lower edge of the first tower
266 # the list ends with the higher edged of the last tower
267
268 # 10 degrees towers
269 set PhiBins {}
270 for {set i -18} {$i <= 18} {incr i} {
271 add PhiBins [expr {$i * $pi/18.0}]
272 }
273 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} {
274 add EtaPhiBins $eta $PhiBins
275 }
276
277 # 20 degrees towers
278 set PhiBins {}
279 for {set i -9} {$i <= 9} {incr i} {
280 add PhiBins [expr {$i * $pi/9.0}]
281 }
282 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} {
283 add EtaPhiBins $eta $PhiBins
284 }
285
286 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
287 add EnergyFraction {0} {0.0 1.0}
288 # energy fractions for e, gamma and pi0
289 add EnergyFraction {11} {1.0 0.0}
290 add EnergyFraction {22} {1.0 0.0}
291 add EnergyFraction {111} {1.0 0.0}
292 # energy fractions for muon, neutrinos and neutralinos
293 add EnergyFraction {12} {0.0 0.0}
294 add EnergyFraction {13} {0.0 0.0}
295 add EnergyFraction {14} {0.0 0.0}
296 add EnergyFraction {16} {0.0 0.0}
297 add EnergyFraction {1000022} {0.0 0.0}
298 add EnergyFraction {1000023} {0.0 0.0}
299 add EnergyFraction {1000025} {0.0 0.0}
300 add EnergyFraction {1000035} {0.0 0.0}
301 add EnergyFraction {1000045} {0.0 0.0}
302 # energy fractions for K0short and Lambda
303 add EnergyFraction {310} {0.3 0.7}
304 add EnergyFraction {3122} {0.3 0.7}
305
306 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
307 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
308 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
309 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
310 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
311 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
312
313 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
314 # http://arxiv.org/pdf/hep-ex/0004009v1
315 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
316 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
317 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
318 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
319}
320
321####################
322# Energy flow merger
323####################
324
325module Merger EFlowMerger {
326# add InputArray InputArray
327 add InputArray Calorimeter/eflowTracks
328 add InputArray Calorimeter/eflowPhotons
329 add InputArray Calorimeter/eflowNeutralHadrons
330 set OutputArray eflow
331}
332
333###################
334# Photon efficiency
335###################
336
337module Efficiency PhotonEfficiency {
338 set InputArray Calorimeter/photons
339 set OutputArray photons
340
341 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
342
343 # efficiency formula for photons
344 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
345 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
346 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
347 (abs(eta) > 2.5) * (0.00)}
348}
349
350##################
351# Photon isolation
352##################
353
354module Isolation PhotonIsolation {
355 set CandidateInputArray PhotonEfficiency/photons
356 set IsolationInputArray EFlowMerger/eflow
357
358 set OutputArray photons
359
360 set DeltaRMax 0.5
361
362 set PTMin 0.5
363
364 set PTRatioMax 0.1
365}
366
367#####################
368# Electron efficiency
369#####################
370
371module Efficiency ElectronEfficiency {
372 set InputArray ElectronEnergySmearing/electrons
373 set OutputArray electrons
374
375 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
376
377 # efficiency formula for "medium" electrons
378 set EfficiencyFormula { (pt <= 7.0) * (0.000) + \
379 (abs(eta) <= 2.47) * (pt > 7.0 && pt <= 25.0) * (0.800) + \
380 (abs(eta) <= 2.47) * (pt > 25.0 && pt <= 50.0) * (0.670 + pt*4.8e-3) + \
381 (abs(eta) <= 2.47) * (pt > 50.0 && pt <= 75.0) * (0.910) + \
382 (abs(eta) <= 2.47) * (pt > 75.0 && pt <= 90.0) * (0.810 + pt*1.5e-3) + \
383 (abs(eta) <= 2.47) * (pt > 90) * (0.945) + \
384 (abs(eta) > 2.47) * (0.000)}
385}
386
387# }
388# set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
389# (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
390# (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
391# (abs(eta) > 2.5) * (0.00)}
392####################
393# Electron isolation
394####################
395
396module Isolation ElectronIsolation {
397 set CandidateInputArray ElectronEfficiency/electrons
398 set IsolationInputArray EFlowMerger/eflow
399
400 set OutputArray electrons
401
402 set DeltaRMax 0.5
403
404 set PTMin 0.5
405
406 set PTRatioMax 0.1
407}
408
409#################
410# Muon efficiency
411#################
412
413module Efficiency MuonEfficiency {
414 set InputArray MuonMomentumSmearing/muons
415 set OutputArray muons
416
417 # set EfficiencyFormula {efficiency as a function of eta and pt}
418 set EfficiencyFormula {
419 (pt <= 6.0) * (0.00) + \
420 (abs(eta) <= 0.1) * (pt > 6.0) * (0.65) + \
421 (abs(eta) > 0.1 && abs(eta) <=2.5) * (pt > 6.0) * (0.99) + \
422 (abs(eta) > 2.5) * (0.00)}
423}
424# default efficiency for muons # efficiency formula for muons
425# set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
426# (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
427# (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
428# (abs(eta) > 2.7) * (0.00)}
429# }
430
431################
432# Muon isolation
433################
434
435module Isolation MuonIsolation {
436 set CandidateInputArray MuonEfficiency/muons
437 set IsolationInputArray EFlowMerger/eflow
438
439 set OutputArray muons
440
441 set DeltaRMax 0.5
442
443 set PTMin 0.5
444
445 set PTRatioMax 0.1
446}
447
448###################
449# Missing ET merger
450###################
451
452module Merger MissingET {
453# add InputArray InputArray
454 add InputArray EFlowMerger/eflow
455 set MomentumOutputArray momentum
456}
457
458##################
459# Scalar HT merger
460##################
461
462module Merger ScalarHT {
463# add InputArray InputArray
464 add InputArray UniqueObjectFinder/jets
465 add InputArray UniqueObjectFinder/electrons
466 add InputArray UniqueObjectFinder/photons
467 add InputArray UniqueObjectFinder/muons
468 set EnergyOutputArray energy
469}
470
471#####################
472# MC truth jet finder
473#####################
474
475module FastJetFinder GenJetFinder {
476 set InputArray Delphes/stableParticles
477
478 set OutputArray jets
479
480 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
481 set JetAlgorithm 6
482 set ParameterR 0.4
483
484 set JetPTMin 20.0
485}
486
487############
488# Jet finder
489############
490
491module FastJetFinder FastJetFinder {
492 set InputArray Calorimeter/towers
493
494 set OutputArray jets
495
496 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
497 set JetAlgorithm 6
498 set ParameterR 0.4
499
500 set JetPTMin 20.0
501}
502
503##################
504# Jet Energy Scale
505##################
506
507module EnergyScale JetEnergyScale {
508 set InputArray FastJetFinder/jets
509 set OutputArray jets
510
511 # scale formula for jets
512 set ScaleFormula {1.08}
513}
514
515###########
516# b-tagging
517###########
518
519module BTagging BTagging {
520 set PartonInputArray Delphes/partons
521 set JetInputArray JetEnergyScale/jets
522
523 set BitNumber 0
524
525 set DeltaR 0.5
526
527 set PartonPTMin 1.0
528
529 set PartonEtaMax 2.5
530
531 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
532 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
533 # gluon's PDG code has the lowest priority
534
535 # default efficiency formula (misidentification rate)
536 add EfficiencyFormula {0} { (pt <= 20.0) * (0.00) + \
537 (abs(eta) <= 2.5) * (pt >= 20.0) * (-2.657e-06*pt + 9.11e-09*pt^2 - 6.19e-3*log(pt)/pt+2.02e-3) + \
538 (abs(eta) >= 2.5) * (0.00)}
539
540 # efficiency formula for c-jets (misidentification rate)
541 add EfficiencyFormula {4} { (pt <= 20.0) * (0.00) + \
542 (abs(eta) <= 2.5) * (pt >= 20.0) * (-4.533e-04*pt + 2.595e-07*pt^2 - 1.243*log(pt)/pt+0.249) + \
543 (abs(eta) >= 2.5) * (0.00)}
544
545 # efficiency formula for b-jets
546 add EfficiencyFormula {5} { (pt <= 20.0) * (0.00) + \
547 (abs(eta) <= 2.5) * (pt >= 20) * (-1.28247e-03*pt + 1.3885e-06*pt^2 - 3.56999*log(pt)/pt+0.817471) + \
548 (abs(eta) >= 2.5) * (0.00)}
549}
550
551module TauTagging TauTagging {
552 set ParticleInputArray Delphes/allParticles
553 set PartonInputArray Delphes/partons
554 set JetInputArray JetEnergyScale/jets
555
556 set DeltaR 0.5
557
558 set TauPTMin 1.0
559
560 set TauEtaMax 2.5
561
562 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
563
564 # default efficiency formula (misidentification rate)
565 add EfficiencyFormula {0} {0.001}
566 # efficiency formula for tau-jets
567 add EfficiencyFormula {15} {0.4}
568}
569
570
571#####################################################
572# Find uniquely identified photons/electrons/tau/jets
573#####################################################
574
575module UniqueObjectFinder UniqueObjectFinder {
576# earlier arrays take precedence over later ones
577# add InputArray InputArray OutputArray
578 add InputArray PhotonIsolation/photons photons
579 add InputArray ElectronIsolation/electrons electrons
580 add InputArray MuonIsolation/muons muons
581 add InputArray JetEnergyScale/jets jets
582}
583
584##################
585# ROOT tree writer
586##################
587
588# tracks, towers and eflow objects are not stored by default in the output.
589# if needed (for jet constituent or other studies), uncomment the relevant
590# "add Branch ..." lines.
591
592# MA5 BEGIN (TO TUNE)
593
594module TreeWriter TreeWriter {
595# add Branch InputArray BranchName BranchClass
596 add Branch Delphes/allParticles Particle GenParticle
597
598# add Branch TrackMerger/tracks Track Track
599# add Branch Calorimeter/towers Tower Tower
600
601# add Branch Calorimeter/eflowTracks EFlowTrack Track
602# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
603# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
604
605 add Branch GenJetFinder/jets GenJet Jet
606 add Branch UniqueObjectFinder/jets Jet Jet
607# add Branch UniqueObjectFinder/electrons Electron Electron
608# add Branch UniqueObjectFinder/photons Photon Photon
609# add Branch UniqueObjectFinder/muons Muon Muon
610 add Branch MissingET/momentum MissingET MissingET
611 add Branch ScalarHT/energy ScalarHT ScalarHT
612
613 add Branch ElectronIsolationCalculation/DelphesMA5tuneElectrons DelphesMA5tuneElectron Electron
614 add Branch MuonIsolationCalculation/DelphesMA5tuneMuons DelphesMA5tuneMuon Muon
615 add Branch PhotonIsolationCalculation/DelphesMA5tunePhotons DelphesMA5tunePhoton Photon
616 add Branch TrackIsolationCalculation/DelphesMA5tuneTracks DelphesMA5tuneTrack Track
617
618}
619
620# MA5 END (TO TUNE)