PublicAnalysisDatabase: delphesMA5tune_card_CMS_sus14004.tcl

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