Back to home page

sPhenix code displayed by LXR

 
 

    


File indexing completed on 2025-08-05 08:11:19

0001 #! /usr/bin/env python
0002 from optparse import OptionParser
0003 import sys
0004 import os
0005 import datetime
0006 from array import *
0007 from ROOT import TH1F, TH2F, TCanvas, TFile, TLegend, TColor, TLine, gROOT, gPad, gSystem, kBlack, kRed, kBlue, kGreen, kOrange, kMagenta, kCyan, kYellow, kViolet, kAzure, kTeal, kSpring, kPink, kGray, kWhite, kDashed, kSolid, kDashDotted, kDashDotted, kDotted, kTRUE, kFALSE
0008 import numpy
0009 import math
0010 import glob
0011 from plotUtil import *
0012 
0013 gROOT.SetBatch(True)
0014 
0015 def colorset_alt(i):
0016     if i == 1:
0017         return ['#f2777a']
0018     elif i == 2:
0019         return ['#f2777a', '#6699cc']
0020     elif i == 3:
0021         return ['#f2777a', '#6699cc', '#99cc99']
0022     elif i == 4:
0023         return ['#f2777a', '#6699cc', '#99cc99', '#9999cc']
0024     elif i == 5:
0025         return ['#f2777a', '#6699cc', '#99cc99', '#9999cc', '#e99960']
0026     elif i == 6:
0027         return ['#f2777a', '#6699cc', '#99cc99', '#9999cc', '#e99960', '#FFC0CB']
0028     elif i == 7:
0029         return ['#f2777a', '#6699cc', '#99cc99', '#9999cc', '#e99960', '#FFC0CB', '#ffcc66']
0030     elif i == 8:
0031         return ['#f2777a', '#6699cc', '#99cc99', '#9999cc', '#e99960', '#FFC0CB', '#ffcc66', '#7FE9DE']
0032     else:
0033         print ("Attempt to use more than 5 colors")
0034         return ['#f2777a', '#6699cc', '#99cc99', '#9999cc', '#e99960', '#FFC0CB', '#ffcc66', '#7FE9DE']
0035 
0036 def Draw_1Dhist_datasimcomp(hdata, hsims, gpadmargin, norm, logy, ymaxscale, XaxisName, Ytitle_unit, prelim, simlegtex, evtseltexts, outname):
0037     hsimcolor = colorset_alt(len(hsims))
0038 
0039     hdata.Sumw2()
0040     for hsim in hsims:
0041         hsim.Sumw2()
0042 
0043     binwidth = hdata.GetXaxis().GetBinWidth(1)
0044 
0045     if norm == 'unity':
0046         hdata.Scale(1. / hdata.Integral(-1, -1))
0047         for hsim in hsims:
0048             hsim.Scale(1. / hsim.Integral(-1, -1))
0049     elif norm == 'data':
0050         for hsim in hsims:
0051             hsim.Scale(hdata.Integral(-1, -1) / hsim.Integral(-1, -1))
0052     else:
0053         if norm != 'none':
0054             print('Invalid normalization option: {}'.format(norm))
0055             sys.exit(1)
0056     
0057     # Get the maximum bin content 
0058     maxbincontent = max(hdata.GetMaximum(), hsim.GetMaximum())
0059     minbincontent = min(hdata.GetMinimum(0), hsim.GetMinimum(0))
0060 
0061     c = TCanvas('c', 'c', 800, 700)
0062     pad1 = TPad( 'pad1', ' ', 0, 0.3, 1, 1)
0063     pad2 = TPad( 'pad2', ' ', 0, 0, 1, 0.3)
0064     pad1.SetRightMargin(gpadmargin[0])
0065     pad1.SetTopMargin(gpadmargin[1])
0066     pad1.SetLeftMargin(gpadmargin[2])
0067     pad1.SetBottomMargin(0.0)
0068     pad1.Draw()
0069     pad2.SetRightMargin(gpadmargin[0])
0070     pad2.SetTopMargin(0.0)
0071     pad2.SetLeftMargin(gpadmargin[2])
0072     pad2.SetBottomMargin(gpadmargin[3])
0073     pad2.Draw() # Draw the TPad on the TCanvas before plotting something on TPad
0074     # cd to the pad1
0075     pad1.cd()
0076     if logy:
0077         pad1.SetLogy()
0078     
0079     for i, hsim in enumerate(hsims):
0080         if i == 0:
0081             if norm == 'unity' or norm == 'data':
0082                 if Ytitle_unit == '':
0083                     hsim.GetYaxis().SetTitle(
0084                         'Normalized entries / ({:g})'.format(binwidth))
0085                 else:
0086                     hsim.GetYaxis().SetTitle(
0087                         'Normalized entries / ({:g} {unit})'.format(binwidth, unit=Ytitle_unit))
0088             else:
0089                 if Ytitle_unit == '':
0090                     hsim.GetYaxis().SetTitle('Entries / ({:g})'.format(binwidth))
0091                 else:
0092                     hsim.GetYaxis().SetTitle(
0093                         'Entries / ({:g} {unit})'.format(binwidth, unit=Ytitle_unit))
0094 
0095             if logy:
0096                 hsim.GetYaxis().SetRangeUser(minbincontent*0.5, maxbincontent * ymaxscale)
0097             else:
0098                 hsim.GetYaxis().SetRangeUser(1E-3, (hsim.GetMaximum()) * ymaxscale)
0099 
0100             # hsim.GetXaxis().SetTitle(XaxisName)
0101             # hsim.GetXaxis().SetTitleOffset(1.1)
0102             hsim.GetXaxis().SetLabelOffset(999)
0103             hsim.GetXaxis().SetLabelSize(0)
0104             hsim.GetYaxis().SetTitleOffset(1.35)
0105             hsim.SetLineColor(TColor.GetColor(hsimcolor[i]))
0106             hsim.SetLineWidth(2)
0107             hsim.SetMarkerSize(0)
0108             hsim.Draw('histe')
0109         else:
0110             hsim.SetLineColor(TColor.GetColor(hsimcolor[i]))
0111             hsim.SetLineWidth(2)
0112             hsim.SetMarkerSize(0)
0113             hsim.Draw('histe same')
0114 
0115     hdata.SetMarkerStyle(20)
0116     hdata.SetMarkerSize(1)
0117     hdata.SetMarkerColor(1)
0118     hdata.SetLineColor(1)
0119     hdata.SetLineWidth(2)
0120     hdata.Draw('same PE1')
0121     shift = 0.45 if prelim else 0.75
0122     legylow = 0.2 + 0.09 * (len(hsims))
0123     leg = TLegend(gPad.GetLeftMargin(), 1-gPad.GetTopMargin()+0.03, 1-gPad.GetRightMargin(), 0.98)
0124     leg.SetNColumns(2)
0125     leg.SetTextSize(0.05)
0126     leg.SetFillStyle(0)
0127     leg.AddEntry(hdata, 'Data', "PE1")
0128     for i, lt in enumerate(simlegtex):
0129         leg.AddEntry(hsims[i], lt, "l")
0130     leg.Draw('same')
0131     
0132     prelimtext = 'Preliminary' if prelim else 'Internal'
0133     sphenixleg = TLegend(gPad.GetLeftMargin()+0.03, 1-gPad.GetTopMargin()-0.18, gPad.GetLeftMargin()+0.15, 1-gPad.GetTopMargin()-0.05)
0134     sphenixleg.SetTextSize(0.05)
0135     sphenixleg.SetFillStyle(0)
0136     sphenixleg.AddEntry('', '#it{#bf{sPHENIX}}'+prelimtext, '')
0137     sphenixleg.AddEntry('', 'Au+Au #sqrt{s_{NN}}=200 GeV', '')
0138     sphenixleg.Draw('same')
0139     
0140     # event selection text
0141     legylow_evtselshift = 0.06 * len(evtseltexts)
0142     leg2 = TLegend(gPad.GetLeftMargin()+0.03, 1-gPad.GetTopMargin()-0.18-legylow_evtselshift, gPad.GetLeftMargin()+0.15, 1-gPad.GetTopMargin()-0.18)
0143     leg2.SetTextSize(0.045)
0144     leg2.SetFillStyle(0)
0145     for i, evtseltext in enumerate(evtseltexts):
0146         leg2.AddEntry('', evtseltext, '')
0147     leg2.Draw('same')
0148     c.RedrawAxis()
0149     
0150     c.Update()
0151     # cd to the pad2
0152     textscale = 2.7
0153     pad2.cd()
0154     pad2.SetLogy(0)
0155     pad2.SetGridy(0)
0156     # take the ratio = data/simulation
0157     l_hM_ratio = []
0158     for i, hsim in enumerate(hsims):
0159         hM_ratio = hdata.Clone('hM_ratio'+str(i))
0160         hM_ratio.Divide(hsim)
0161         hM_ratio.SetMarkerStyle(20)
0162         hM_ratio.SetMarkerSize(0.7)
0163         hM_ratio.SetMarkerColor(TColor.GetColor(hsimcolor[i]))
0164         hM_ratio.SetLineColor(TColor.GetColor(hsimcolor[i]))
0165         hM_ratio.SetLineWidth(2)
0166         if i == 0:
0167             # hM_ratio.GetYaxis().SetRangeUser(0.001, 1.999)
0168             hM_ratio.GetYaxis().SetRangeUser(0.501, 1.499)
0169             hM_ratio.GetYaxis().SetNdivisions(505)
0170             hM_ratio.GetYaxis().SetTitle('Data/Sim.')
0171             hM_ratio.GetYaxis().SetTitleOffset(0.5)
0172             hM_ratio.GetYaxis().SetTickSize(TickSize)
0173             hM_ratio.GetYaxis().SetTitleSize(AxisTitleSize*textscale)
0174             hM_ratio.GetYaxis().SetLabelSize(AxisLabelSize*textscale)
0175             hM_ratio.GetXaxis().SetTickSize(TickSize*textscale)
0176             hM_ratio.GetXaxis().SetTitleSize(AxisTitleSize*textscale)
0177             hM_ratio.GetXaxis().SetLabelSize(AxisLabelSize*textscale)
0178             hM_ratio.GetXaxis().SetTitleOffset(1.1)
0179             hM_ratio.GetXaxis().SetTitle(XaxisName)
0180             hM_ratio.Draw('PE1X0')
0181         else:
0182             hM_ratio.Draw('PE1X0 same')
0183             
0184         l_hM_ratio.append(hM_ratio)
0185     
0186     # Draw a line at 1
0187     line = TLine(hM_ratio.GetXaxis().GetXmin(), 1, hM_ratio.GetXaxis().GetXmax(), 1)
0188     line.SetLineColor(kBlack)
0189     line.SetLineStyle(2)
0190     line.Draw()
0191     c.Draw()
0192     c.SaveAs(outname+'.pdf')
0193     c.SaveAs(outname+'.png')
0194     if(c):
0195         c.Close()
0196         gSystem.ProcessEvents()
0197         del c
0198         c = 0
0199         
0200 
0201 # comparae data and (one of the) simulation
0202 def Draw_2Dhist_datasimcomp(hdata, hsim, logz, norm, rmargin, XaxisName, YaxisName, outname):
0203     if norm == 'unity':
0204         hdata.Scale(1. / hdata.Integral(-1, -1, -1, -1))
0205         hsim.Scale(1. / hsim.Integral(-1, -1, -1, -1))
0206     elif norm == 'data':
0207         hsim.Scale(hdata.Integral(-1, -1, -1, -1) / hsim.Integral(-1, -1, -1, -1))
0208     else:
0209         if norm != 'none':
0210             print('Invalid normalization option: {}'.format(norm))
0211             sys.exit(1)
0212     
0213     c = TCanvas('c', 'c', 800, 700)
0214     if logz:
0215         c.SetLogz()
0216     c.cd()
0217     gPad.SetRightMargin(rmargin)
0218     gPad.SetTopMargin(TopMargin)
0219     gPad.SetLeftMargin(LeftMargin)
0220     gPad.SetBottomMargin(BottomMargin)
0221     hdata.GetXaxis().SetTitle(XaxisName)
0222     hdata.GetYaxis().SetTitle(YaxisName)
0223     hdata.GetXaxis().SetTickSize(TickSize)
0224     hdata.GetYaxis().SetTickSize(TickSize)
0225     hdata.GetXaxis().SetTitleSize(AxisTitleSize)
0226     hdata.GetYaxis().SetTitleSize(AxisTitleSize)
0227     hdata.GetXaxis().SetLabelSize(AxisLabelSize)
0228     hdata.GetYaxis().SetLabelSize(AxisLabelSize)
0229     hdata.GetXaxis().SetTitleOffset(1.1)
0230     hdata.GetYaxis().SetTitleOffset(1.3)
0231     hdata.GetZaxis().SetLabelSize(AxisLabelSize)
0232     hdata.SetLineColor(kBlack)
0233     hdata.SetLineWidth(1)
0234     hdata.Draw('box')
0235     hsim.SetLineColorAlpha(kRed, 0.5)
0236     hsim.SetLineWidth(1)
0237     hsim.Draw('box same')
0238     c.RedrawAxis()
0239     c.Draw()
0240     c.SaveAs(outname+'.pdf')
0241     c.SaveAs(outname+'.png')
0242     if(c):
0243         c.Close()
0244         gSystem.ProcessEvents()
0245         del c
0246         c = 0
0247 
0248 
0249 if __name__ == '__main__':
0250     parser = OptionParser(usage='usage: %prog ver [options -h]')
0251     parser.add_option('-d', '--datahistdir', dest='datahistdir', type='string', default='/sphenix/user/hjheng/TrackletAna/analysis_INTT/plot/hists/data_run20869/', help='Histogram file name (data)')
0252     parser.add_option('-s', '--simhistdir', action='append', dest='simhistdir', type='string', help='Histogram file name (simulation). Example: /sphenix/user/hjheng/TrackletAna/analysis_INTT/plot/hists/ana382_zvtx-20cm_dummyAlignParams/Hists_RecoTracklets_merged.root')
0253     parser.add_option('-l', '--simlegtext', action='append', dest='simlegtext', type='string', help='Legend text for simulation. Example: HIJING/EPOS/AMPT)')
0254     parser.add_option('-p', '--plotdir', dest='plotdir', type='string', default='ana382_zvtx-20cm_dummyAlignParams', help='Plot directory')
0255 
0256     (opt, args) = parser.parse_args()
0257 
0258     print('opt: {}'.format(opt))
0259 
0260     datahistdir = opt.datahistdir
0261     simhistdir = opt.simhistdir
0262     simlegtext = opt.simlegtext
0263     plotdir = opt.plotdir
0264 
0265     if os.path.isfile("{}/hists_merged.root".format(datahistdir)):
0266         os.system("rm {}/hists_merged.root".format(datahistdir))
0267         os.system("hadd -f -j 20 {}/hists_merged.root {}/hists_0*.root".format(datahistdir, datahistdir))
0268     else:
0269         os.system("hadd -f -j 20 {}/hists_merged.root {}/hists_0*.root".format(datahistdir, datahistdir))
0270     
0271     for simhistd in simhistdir:
0272         if os.path.isfile("{}/hists_merged.root".format(simhistd)):
0273             os.system("rm {}/hists_merged.root".format(simhistd))
0274             os.system("hadd -f -j 20 {}/hists_merged.root {}/hists_0*.root".format(simhistd, simhistd))
0275         else:
0276             os.system("hadd -f -j 20 {}/hists_merged.root {}/hists_0*.root".format(simhistd, simhistd))
0277 
0278     os.makedirs('./DataSimComp/{}'.format(plotdir), exist_ok=True)
0279 
0280     hM_NClusLayer1_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_NClusLayer1')
0281     hM_NClusLayer2_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_NClusLayer2')
0282     hM_NClusLayer1_zvtxwei_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_NClusLayer1_zvtxwei')
0283     hM_NTklclusLayer1_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_NTklclusLayer1')
0284     hM_NRecotkl_Raw_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_NRecotkl_Raw')
0285     hM_dEta_reco_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco') 
0286     hM_dEta_reco_altrange_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_altrange')
0287     hM_dEta_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20')
0288     hM_dEta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20')
0289     hM_dEta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20')
0290     hM_dEta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20')
0291     hM_dEta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20')
0292     hM_dEta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20')
0293     hM_dEta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20')
0294     hM_dEta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20')
0295     hM_dEta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20')
0296     hM_dEta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20')
0297     hM_dEta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20')
0298     hM_dEta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20')
0299     hM_dEta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20')
0300     hM_dEta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20')
0301     hM_dEta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20')
0302     hM_dEta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20')
0303     hM_dEta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20')
0304     hM_dEta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20')
0305     hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20')
0306     hM_dPhi_reco_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco')
0307     hM_dPhi_reco_altrange_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_altrange')
0308     hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20')
0309     hM_dPhi_reco_Centrality_0to3_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_0to3')
0310     hM_dPhi_reco_Centrality_3to6_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_3to6')
0311     hM_dPhi_reco_Centrality_6to10_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_6to10')
0312     hM_dPhi_reco_Centrality_10to15_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_10to15')
0313     hM_dPhi_reco_Centrality_15to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_15to20')
0314     hM_dPhi_reco_Centrality_20to25_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_20to25')
0315     hM_dPhi_reco_Centrality_25to30_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_25to30')
0316     hM_dPhi_reco_Centrality_30to35_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_30to35')
0317     hM_dPhi_reco_Centrality_35to40_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_35to40')
0318     hM_dPhi_reco_Centrality_40to45_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_40to45')
0319     hM_dPhi_reco_Centrality_45to50_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_45to50')
0320     hM_dPhi_reco_Centrality_50to55_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_50to55')
0321     hM_dPhi_reco_Centrality_55to60_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_55to60')
0322     hM_dPhi_reco_Centrality_60to65_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_60to65')
0323     hM_dPhi_reco_Centrality_65to70_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_65to70')
0324     hM_dPhi_reco_Centrality_70to80_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_70to80')
0325     hM_dPhi_reco_Centrality_80to90_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_80to90')
0326     hM_dPhi_reco_Centrality_90to100_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_90to100')
0327     hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75')
0328     hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75')
0329     hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75')
0330     hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75')
0331     hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75')
0332     hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75')
0333     hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75')
0334     hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75')
0335     hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75')
0336     hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75')
0337     hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75')
0338     hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75')
0339     hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75')
0340     hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75')
0341     hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75')
0342     hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75')
0343     hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75')
0344     hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75')
0345     hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20')
0346     hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20')
0347     hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20')
0348     hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20')
0349     hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20')
0350     hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20')
0351     hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20')
0352     hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20')
0353     hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20')
0354     hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20')
0355     hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20')
0356     hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20')
0357     hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20')
0358     hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20')
0359     hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20')
0360     hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20')
0361     hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20')
0362     hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20')
0363     hM_dR_reco_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dR_reco')
0364     hM_dR_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_dR_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20')
0365     hM_Eta_reco_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco')
0366     hM_Eta_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20')
0367     hM_Eta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20')
0368     hM_Eta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20')
0369     hM_Eta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20')
0370     hM_Eta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20')
0371     hM_Eta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20')
0372     hM_Eta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20')
0373     hM_Eta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20')
0374     hM_Eta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20')
0375     hM_Eta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20')
0376     hM_Eta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20')
0377     hM_Eta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20')
0378     hM_Eta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20')
0379     hM_Eta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20')
0380     hM_Eta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20')
0381     hM_Eta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20')
0382     hM_Eta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20')
0383     hM_Eta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20')
0384     hM_Eta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Eta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20')
0385     hM_Phi_reco_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Phi_reco')
0386     hM_Phi_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_Phi_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20') 
0387     hM_RecoPVz_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_RecoPVz')
0388     hM_RecoPVz_MBDAsymLe0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_RecoPVz_MBDAsymLe0p75_VtxZm20to20')
0389     hM_MBDChargeAsymm_Le0p75_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_MBDChargeAsymm_Le0p75')
0390     hM_MBDChargeAsymm_Le0p75_VtxZm20to20_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_MBDChargeAsymm_Le0p75_VtxZm20to20')
0391     hM_clusphi_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_clusphi')
0392     hM_clusphi_zvtxwei_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_clusphi_zvtxwei')
0393     hM_cluseta_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_cluseta')
0394     hM_cluseta_zvtxwei_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_cluseta_zvtxwei')
0395     hM_clusphisize_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_clusphisize')
0396     hM_clusadc_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_clusadc')
0397     hM_cluseta_clusphisize_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_cluseta_clusphisize')
0398     hM_clusphi_clusphisize_data = GetHistogram("{}/hists_merged.root".format(datahistdir), 'hM_clusphi_clusphisize')
0399 
0400     l_hM_NClusLayer1_sim = []
0401     l_hM_NClusLayer2_sim = []
0402     l_hM_NClusLayer1_zvtxwei_sim = []
0403     l_hM_NTklclusLayer1_sim = []
0404     l_hM_NRecotkl_Raw_sim = []
0405     l_hM_dEta_reco_sim = []
0406     l_hM_dEta_reco_altrange_sim = []
0407     l_hM_dEta_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0408     l_hM_dEta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim = []
0409     l_hM_dEta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim = []
0410     l_hM_dEta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim = []
0411     l_hM_dEta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim = []
0412     l_hM_dEta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim = []
0413     l_hM_dEta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim = []
0414     l_hM_dEta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim = []
0415     l_hM_dEta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim = []
0416     l_hM_dEta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim = []
0417     l_hM_dEta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim = []
0418     l_hM_dEta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim = []
0419     l_hM_dEta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim = []
0420     l_hM_dEta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim = []
0421     l_hM_dEta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim = []
0422     l_hM_dEta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0423     l_hM_dEta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim = []
0424     l_hM_dEta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim = []
0425     l_hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim = []
0426     l_hM_dPhi_reco_sim = []
0427     l_hM_dPhi_reco_altrange_sim = []
0428     l_hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0429     l_hM_dPhi_reco_Centrality_0to3_sim = []
0430     l_hM_dPhi_reco_Centrality_3to6_sim = []
0431     l_hM_dPhi_reco_Centrality_6to10_sim = []
0432     l_hM_dPhi_reco_Centrality_10to15_sim = []
0433     l_hM_dPhi_reco_Centrality_15to20_sim = []
0434     l_hM_dPhi_reco_Centrality_20to25_sim = []
0435     l_hM_dPhi_reco_Centrality_25to30_sim = []
0436     l_hM_dPhi_reco_Centrality_30to35_sim = []
0437     l_hM_dPhi_reco_Centrality_35to40_sim = []
0438     l_hM_dPhi_reco_Centrality_40to45_sim = []
0439     l_hM_dPhi_reco_Centrality_45to50_sim = []
0440     l_hM_dPhi_reco_Centrality_50to55_sim = []
0441     l_hM_dPhi_reco_Centrality_55to60_sim = []
0442     l_hM_dPhi_reco_Centrality_60to65_sim = []
0443     l_hM_dPhi_reco_Centrality_65to70_sim = []
0444     l_hM_dPhi_reco_Centrality_70to80_sim = []
0445     l_hM_dPhi_reco_Centrality_80to90_sim = []
0446     l_hM_dPhi_reco_Centrality_90to100_sim = []
0447     l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_sim = []
0448     l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_sim = []
0449     l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_sim = []
0450     l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_sim = []
0451     l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_sim = []
0452     l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_sim = []
0453     l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_sim = []
0454     l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_sim = []
0455     l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_sim = []
0456     l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_sim = []
0457     l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_sim = []
0458     l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_sim = []
0459     l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_sim = []
0460     l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_sim = []
0461     l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_sim = []
0462     l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_sim = []
0463     l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_sim = []
0464     l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_sim = []
0465     l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim = []
0466     l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim = []
0467     l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim = []
0468     l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim = []
0469     l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim = []
0470     l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim = []
0471     l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim = []
0472     l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim = []
0473     l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim = []
0474     l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim = []
0475     l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim = []
0476     l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim = []
0477     l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim = []
0478     l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim = []
0479     l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0480     l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim = []
0481     l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim = []
0482     l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim = []
0483     l_hM_dR_reco_sim = []
0484     l_hM_dR_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0485     l_hM_Eta_reco_sim = []
0486     l_hM_Eta_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0487     l_hM_Eta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim = []
0488     l_hM_Eta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim = []
0489     l_hM_Eta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim = []
0490     l_hM_Eta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim = []
0491     l_hM_Eta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim = []
0492     l_hM_Eta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim = []
0493     l_hM_Eta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim = []
0494     l_hM_Eta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim = []
0495     l_hM_Eta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim = []
0496     l_hM_Eta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim = []
0497     l_hM_Eta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim = []
0498     l_hM_Eta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim = []
0499     l_hM_Eta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim = []
0500     l_hM_Eta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim = []
0501     l_hM_Eta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0502     l_hM_Eta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim = []
0503     l_hM_Eta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim = []
0504     l_hM_Eta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim = []
0505     l_hM_Phi_reco_sim = []
0506     l_hM_Phi_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim = []
0507     l_hM_RecoPVz_sim = []
0508     l_hM_RecoPVz_MBDAsymLe0p75_VtxZm20to20_sim = []
0509     l_hM_MBDChargeAsymm_Le0p75_sim = []
0510     l_hM_MBDChargeAsymm_Le0p75_VtxZm20to20_sim = []
0511     l_hM_clusphi_sim = []
0512     l_hM_clusphi_zvtxwei_sim = []
0513     l_hM_cluseta_sim = []
0514     l_hM_cluseta_zvtxwei_sim = []
0515     l_hM_clusphisize_sim = []
0516     l_hM_clusadc_sim = []
0517     l_hM_cluseta_clusphisize_sim = []
0518     l_hM_clusphi_clusphisize_sim = []
0519     for i, simhistd in enumerate(simhistdir):
0520         l_hM_NClusLayer1_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_NClusLayer1'))
0521         l_hM_NClusLayer2_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_NClusLayer2'))
0522         l_hM_NClusLayer1_zvtxwei_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_NClusLayer1_zvtxwei'))
0523         l_hM_NTklclusLayer1_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_NTklclusLayer1'))
0524         l_hM_NRecotkl_Raw_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_NRecotkl_Raw'))
0525         l_hM_dEta_reco_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco'))
0526         l_hM_dEta_reco_altrange_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_altrange'))
0527         l_hM_dEta_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'))
0528         l_hM_dEta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20'))
0529         l_hM_dEta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20'))
0530         l_hM_dEta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20'))
0531         l_hM_dEta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20'))
0532         l_hM_dEta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20'))
0533         l_hM_dEta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20'))
0534         l_hM_dEta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20'))
0535         l_hM_dEta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20'))
0536         l_hM_dEta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20'))
0537         l_hM_dEta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20'))
0538         l_hM_dEta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20'))
0539         l_hM_dEta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20'))
0540         l_hM_dEta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20'))
0541         l_hM_dEta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20'))
0542         l_hM_dEta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20'))
0543         l_hM_dEta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20'))
0544         l_hM_dEta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20'))
0545         l_hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'))
0546         l_hM_dPhi_reco_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco'))
0547         l_hM_dPhi_reco_altrange_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_altrange'))
0548         l_hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'))
0549         l_hM_dPhi_reco_Centrality_0to3_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_0to3'))
0550         l_hM_dPhi_reco_Centrality_3to6_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_3to6'))
0551         l_hM_dPhi_reco_Centrality_6to10_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_6to10'))
0552         l_hM_dPhi_reco_Centrality_10to15_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_10to15'))
0553         l_hM_dPhi_reco_Centrality_15to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_15to20'))
0554         l_hM_dPhi_reco_Centrality_20to25_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_20to25'))
0555         l_hM_dPhi_reco_Centrality_25to30_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_25to30'))
0556         l_hM_dPhi_reco_Centrality_30to35_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_30to35'))
0557         l_hM_dPhi_reco_Centrality_35to40_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_35to40'))
0558         l_hM_dPhi_reco_Centrality_40to45_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_40to45'))
0559         l_hM_dPhi_reco_Centrality_45to50_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_45to50'))
0560         l_hM_dPhi_reco_Centrality_50to55_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_50to55'))
0561         l_hM_dPhi_reco_Centrality_55to60_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_55to60'))
0562         l_hM_dPhi_reco_Centrality_60to65_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_60to65'))
0563         l_hM_dPhi_reco_Centrality_65to70_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_65to70'))
0564         l_hM_dPhi_reco_Centrality_70to80_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_70to80'))
0565         l_hM_dPhi_reco_Centrality_80to90_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_80to90'))
0566         l_hM_dPhi_reco_Centrality_90to100_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_90to100'))
0567         l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75'))
0568         l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75'))
0569         l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75'))
0570         l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75'))
0571         l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75'))
0572         l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75'))
0573         l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75'))
0574         l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75'))
0575         l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75'))
0576         l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75'))
0577         l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75'))
0578         l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75'))
0579         l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75'))
0580         l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75'))
0581         l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75'))
0582         l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75'))
0583         l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75'))
0584         l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75'))
0585         l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20'))
0586         l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20'))
0587         l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20'))
0588         l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20'))
0589         l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20'))
0590         l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20'))
0591         l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20'))
0592         l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20'))
0593         l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20'))
0594         l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20'))
0595         l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20'))
0596         l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20'))
0597         l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20'))
0598         l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20'))
0599         l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20'))
0600         l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20'))
0601         l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20'))
0602         l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'))
0603         l_hM_dR_reco_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dR_reco'))
0604         l_hM_dR_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_dR_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'))
0605         l_hM_Eta_reco_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco'))
0606         l_hM_Eta_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'))
0607         l_hM_Eta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20'))
0608         l_hM_Eta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20'))
0609         l_hM_Eta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20'))
0610         l_hM_Eta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20'))
0611         l_hM_Eta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20'))
0612         l_hM_Eta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20'))
0613         l_hM_Eta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20'))
0614         l_hM_Eta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20'))
0615         l_hM_Eta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20'))
0616         l_hM_Eta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20'))
0617         l_hM_Eta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20'))
0618         l_hM_Eta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20'))
0619         l_hM_Eta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20'))
0620         l_hM_Eta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20'))
0621         l_hM_Eta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20'))
0622         l_hM_Eta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20'))
0623         l_hM_Eta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20'))
0624         l_hM_Eta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Eta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'))
0625         l_hM_Phi_reco_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Phi_reco'))
0626         l_hM_Phi_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_Phi_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'))
0627         l_hM_RecoPVz_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_RecoPVz'))
0628         l_hM_RecoPVz_MBDAsymLe0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_RecoPVz_MBDAsymLe0p75_VtxZm20to20'))
0629         l_hM_MBDChargeAsymm_Le0p75_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_MBDChargeAsymm_Le0p75'))
0630         l_hM_MBDChargeAsymm_Le0p75_VtxZm20to20_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_MBDChargeAsymm_Le0p75_VtxZm20to20'))
0631         l_hM_clusphi_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_clusphi'))
0632         l_hM_clusphi_zvtxwei_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_clusphi_zvtxwei'))
0633         l_hM_cluseta_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_cluseta'))
0634         l_hM_cluseta_zvtxwei_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_cluseta_zvtxwei'))
0635         l_hM_clusphisize_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_clusphisize'))
0636         l_hM_clusadc_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_clusadc'))
0637         l_hM_cluseta_clusphisize_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_cluseta_clusphisize'))
0638         l_hM_clusphi_clusphisize_sim.append(GetHistogram("{}/hists_merged.root".format(simhistd), 'hM_clusphi_clusphisize'))
0639         
0640     str_clussel = 'Cluster ADC>35'
0641     str_zvtxsel = '|INTT&MBD Z_{vtx}|#leq10cm'
0642     str_evtsel = 'Trigger bit 10 (MBD N&S#geq 2), Is MinBias'
0643     str_globsel = ''
0644 
0645     padmargin = [0.08,0.2,0.15,0.32]
0646 
0647     # Draw_1Dhist_datasimcomp(hdata, hsims, gpadmargin, norm, logy, ymaxscale, XaxisName, Ytitle_unit, prelim, simlegtex, evtseltexts, evtseltextpos, outname)
0648     hM_NClusLayer1_data.GetXaxis().SetMaxDigits(2)
0649     for hM_NClusLayer1_sim in l_hM_NClusLayer1_sim:
0650         hM_NClusLayer1_sim.GetXaxis().SetMaxDigits(2)
0651     Draw_1Dhist_datasimcomp(hM_NClusLayer1_data, l_hM_NClusLayer1_sim, padmargin, 'data', True, 500, 'Number of clusters (inner)', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/NClusLayer1'.format(plotdir))
0652     
0653     hM_NClusLayer2_data.GetXaxis().SetMaxDigits(2)
0654     for hM_NClusLayer2_sim in l_hM_NClusLayer2_sim:
0655         hM_NClusLayer2_sim.GetXaxis().SetMaxDigits(2)
0656     Draw_1Dhist_datasimcomp(hM_NClusLayer2_data, l_hM_NClusLayer2_sim, padmargin, 'data', True, 500, 'Number of clusters (outer)', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/NClusLayer2'.format(plotdir))
0657     
0658     hM_NClusLayer1_zvtxwei_data.GetXaxis().SetMaxDigits(2)
0659     for hM_NClusLayer1_zvtxwei_sim in l_hM_NClusLayer1_zvtxwei_sim:
0660         hM_NClusLayer1_zvtxwei_sim.GetXaxis().SetMaxDigits(2)
0661     Draw_1Dhist_datasimcomp(hM_NClusLayer1_zvtxwei_data, l_hM_NClusLayer1_zvtxwei_sim, padmargin, 'data', True, 500, 'Number of clusters (inner)', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/NClusLayer1_zvtxwei'.format(plotdir))
0662     
0663     hM_NTklclusLayer1_data.GetXaxis().SetMaxDigits(2)
0664     for hM_NTklclusLayer1_sim in l_hM_NTklclusLayer1_sim:
0665         hM_NTklclusLayer1_sim.GetXaxis().SetMaxDigits(2)
0666     Draw_1Dhist_datasimcomp(hM_NTklclusLayer1_data, l_hM_NTklclusLayer1_sim, padmargin, 'data', True, 500, 'Number of tracklet clusters (inner)', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/NTklClusLayer1'.format(plotdir))
0667     
0668     hM_NRecotkl_Raw_data.GetXaxis().SetMaxDigits(2)
0669     for hM_NRecotkl_Raw_sim in l_hM_NRecotkl_Raw_sim:
0670         hM_NRecotkl_Raw_sim.GetXaxis().SetMaxDigits(2)
0671     Draw_1Dhist_datasimcomp(hM_NRecotkl_Raw_data, l_hM_NRecotkl_Raw_sim, padmargin, 'data', True, 500, 'Number of reco-tracklets', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/NRecoTracklets_Raw'.format(plotdir))
0672     
0673     Draw_1Dhist_datasimcomp(hM_RecoPVz_data, l_hM_RecoPVz_sim, padmargin, 'data', False, 1.7, 'INTT Z_{vtx} [cm]', 'cm', False, simlegtext, ['{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoPVz'.format(plotdir))
0674     Draw_1Dhist_datasimcomp(hM_RecoPVz_MBDAsymLe0p75_VtxZm20to20_data, l_hM_RecoPVz_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 1.7, 'INTT Z_{vtx} [cm]', 'cm', False, simlegtext, ['{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))+', |Asymm._{MBD}|#leq0.75'], './DataSimComp/{}/RecoPVz_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0675     Draw_1Dhist_datasimcomp(hM_MBDChargeAsymm_Le0p75_data, l_hM_MBDChargeAsymm_Le0p75_sim, padmargin, 'data', True, 15, 'MBD charge asymmetry', '', False, simlegtext, ['|Asymm._{MBD}|#leq0.75'], './DataSimComp/{}/MBDChargeAsymm'.format(plotdir))
0676     Draw_1Dhist_datasimcomp(hM_MBDChargeAsymm_Le0p75_VtxZm20to20_data, l_hM_MBDChargeAsymm_Le0p75_VtxZm20to20_sim, padmargin, 'data', True, 15, 'MBD charge asymmetry', '', False, simlegtext, ['{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))+', |Asymm._{MBD}|#leq0.75'], './DataSimComp/{}/MBDChargeAsymm_VtxZm20to20'.format(plotdir))
0677     
0678     Draw_1Dhist_datasimcomp(hM_clusphi_data, l_hM_clusphi_sim, padmargin, 'data', False, 1.8, 'Cluster #phi', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/Cluster_Phi'.format(plotdir))
0679     Draw_1Dhist_datasimcomp(hM_clusphi_zvtxwei_data, l_hM_clusphi_zvtxwei_sim, padmargin, 'data', False, 1.8, 'Cluster #phi', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/Cluster_Phi_zvtxwei'.format(plotdir))
0680     Draw_1Dhist_datasimcomp(hM_cluseta_data, l_hM_cluseta_sim, padmargin, 'data', False, 1.8, 'Cluster #eta', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/Cluster_Eta'.format(plotdir))
0681     Draw_1Dhist_datasimcomp(hM_cluseta_zvtxwei_data, l_hM_cluseta_zvtxwei_sim, padmargin, 'data', False, 1.8, 'Cluster #eta', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/Cluster_Eta_zvtxwei'.format(plotdir))
0682     Draw_1Dhist_datasimcomp(hM_clusphisize_data, l_hM_clusphisize_sim, padmargin, 'data', True, 100, 'Cluster #phi size', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/Cluster_PhiSize'.format(plotdir))
0683     hM_clusadc_data.GetXaxis().SetMaxDigits(3)
0684     hM_clusadc_data.GetXaxis().SetNdivisions(-10)
0685     for hM_clusadc_sim in l_hM_clusadc_sim:
0686         hM_clusadc_sim.GetXaxis().SetMaxDigits(3)
0687         hM_clusadc_sim.GetXaxis().SetNdivisions(-10)
0688     Draw_1Dhist_datasimcomp(hM_clusadc_data, l_hM_clusadc_sim, padmargin, 'data', True, 100, 'Cluster ADC', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/Cluster_ADC'.format(plotdir))
0689 
0690     Draw_1Dhist_datasimcomp(hM_dEta_reco_data, l_hM_dEta_reco_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dEta'.format(plotdir))
0691     Draw_1Dhist_datasimcomp(hM_dEta_reco_altrange_data, l_hM_dEta_reco_altrange_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dEta_altrange'.format(plotdir))
0692     Draw_1Dhist_datasimcomp(hM_dEta_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 0-70%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0693     # Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 800, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 0-5%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0694     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 0-3%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0695     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 3-6%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0696     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 6-10%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0697     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 10-15%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0698     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 15-20%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0699     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 20-25%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0700     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 25-30%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0701     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 30-35%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0702     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 35-40%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0703     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 40-45%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0704     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 45-50%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0705     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 50-55%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0706     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 55-60%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0707     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 60-65%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0708     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 65-70%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0709     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 70-80%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0710     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 80-90%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0711     Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0712     # Draw_1Dhist_datasimcomp(hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dEta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#eta', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dEta_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0713     Draw_1Dhist_datasimcomp(hM_dPhi_reco_data, l_hM_dPhi_reco_sim, padmargin, 'data', True, 750, 'Reco-tracklet #Delta#phi', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi'.format(plotdir))
0714     
0715     hM_dPhi_reco_altrange_data.GetXaxis().SetMaxDigits(2)
0716     hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0717     hM_dPhi_reco_Centrality_0to3_data.GetXaxis().SetMaxDigits(2)
0718     hM_dPhi_reco_Centrality_3to6_data.GetXaxis().SetMaxDigits(2)
0719     hM_dPhi_reco_Centrality_6to10_data.GetXaxis().SetMaxDigits(2)
0720     hM_dPhi_reco_Centrality_10to15_data.GetXaxis().SetMaxDigits(2)
0721     hM_dPhi_reco_Centrality_15to20_data.GetXaxis().SetMaxDigits(2)
0722     hM_dPhi_reco_Centrality_20to25_data.GetXaxis().SetMaxDigits(2)
0723     hM_dPhi_reco_Centrality_25to30_data.GetXaxis().SetMaxDigits(2)
0724     hM_dPhi_reco_Centrality_30to35_data.GetXaxis().SetMaxDigits(2)
0725     hM_dPhi_reco_Centrality_35to40_data.GetXaxis().SetMaxDigits(2)
0726     hM_dPhi_reco_Centrality_40to45_data.GetXaxis().SetMaxDigits(2)
0727     hM_dPhi_reco_Centrality_45to50_data.GetXaxis().SetMaxDigits(2)
0728     hM_dPhi_reco_Centrality_50to55_data.GetXaxis().SetMaxDigits(2)
0729     hM_dPhi_reco_Centrality_55to60_data.GetXaxis().SetMaxDigits(2)
0730     hM_dPhi_reco_Centrality_60to65_data.GetXaxis().SetMaxDigits(2)
0731     hM_dPhi_reco_Centrality_65to70_data.GetXaxis().SetMaxDigits(2)
0732     hM_dPhi_reco_Centrality_70to80_data.GetXaxis().SetMaxDigits(2)
0733     hM_dPhi_reco_Centrality_80to90_data.GetXaxis().SetMaxDigits(2)
0734     hM_dPhi_reco_Centrality_90to100_data.GetXaxis().SetMaxDigits(2)
0735     hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0736     hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0737     hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0738     hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0739     hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0740     hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0741     hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0742     hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0743     hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0744     hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0745     hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0746     hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0747     hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0748     hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0749     hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0750     hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0751     hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0752     hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_data.GetXaxis().SetMaxDigits(2)
0753     hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0754     hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0755     hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0756     hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0757     hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0758     hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0759     hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0760     hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0761     hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0762     hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0763     hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0764     hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0765     hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0766     hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0767     hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0768     hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0769     hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0770     hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data.GetXaxis().SetMaxDigits(2)
0771     for i in range(len(l_hM_dPhi_reco_altrange_sim)):
0772         l_hM_dPhi_reco_altrange_sim[i].GetXaxis().SetMaxDigits(2)
0773         l_hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0774         l_hM_dPhi_reco_Centrality_0to3_sim[i].GetXaxis().SetMaxDigits(2)
0775         l_hM_dPhi_reco_Centrality_3to6_sim[i].GetXaxis().SetMaxDigits(2)
0776         l_hM_dPhi_reco_Centrality_6to10_sim[i].GetXaxis().SetMaxDigits(2)
0777         l_hM_dPhi_reco_Centrality_10to15_sim[i].GetXaxis().SetMaxDigits(2)
0778         l_hM_dPhi_reco_Centrality_15to20_sim[i].GetXaxis().SetMaxDigits(2)
0779         l_hM_dPhi_reco_Centrality_20to25_sim[i].GetXaxis().SetMaxDigits(2)
0780         l_hM_dPhi_reco_Centrality_25to30_sim[i].GetXaxis().SetMaxDigits(2)
0781         l_hM_dPhi_reco_Centrality_30to35_sim[i].GetXaxis().SetMaxDigits(2)
0782         l_hM_dPhi_reco_Centrality_35to40_sim[i].GetXaxis().SetMaxDigits(2)
0783         l_hM_dPhi_reco_Centrality_40to45_sim[i].GetXaxis().SetMaxDigits(2)
0784         l_hM_dPhi_reco_Centrality_45to50_sim[i].GetXaxis().SetMaxDigits(2)
0785         l_hM_dPhi_reco_Centrality_50to55_sim[i].GetXaxis().SetMaxDigits(2)
0786         l_hM_dPhi_reco_Centrality_55to60_sim[i].GetXaxis().SetMaxDigits(2)
0787         l_hM_dPhi_reco_Centrality_60to65_sim[i].GetXaxis().SetMaxDigits(2)
0788         l_hM_dPhi_reco_Centrality_65to70_sim[i].GetXaxis().SetMaxDigits(2)
0789         l_hM_dPhi_reco_Centrality_70to80_sim[i].GetXaxis().SetMaxDigits(2)
0790         l_hM_dPhi_reco_Centrality_80to90_sim[i].GetXaxis().SetMaxDigits(2)
0791         l_hM_dPhi_reco_Centrality_90to100_sim[i].GetXaxis().SetMaxDigits(2)
0792         l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0793         l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0794         l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0795         l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0796         l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0797         l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0798         l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0799         l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0800         l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0801         l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0802         l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0803         l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0804         l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0805         l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0806         l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0807         l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0808         l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0809         l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_sim[i].GetXaxis().SetMaxDigits(2)
0810         l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0811         l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0812         l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0813         l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0814         l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0815         l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0816         l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0817         l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0818         l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0819         l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0820         l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0821         l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0822         l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0823         l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0824         l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0825         l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0826         l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0827         l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim[i].GetXaxis().SetMaxDigits(2)
0828         
0829     Draw_1Dhist_datasimcomp(hM_dPhi_reco_altrange_data, l_hM_dPhi_reco_altrange_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_altrange'.format(plotdir))
0830     Draw_1Dhist_datasimcomp(hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 0-70%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dPhi_altrange_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0831     # Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_0to5_data, l_hM_dPhi_reco_Centrality_0to5_sim, padmargin, 'data', True, 100, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 0-5%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_0to5'.format(plotdir))
0832     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_0to3_data, l_hM_dPhi_reco_Centrality_0to3_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 0-3%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_0to3'.format(plotdir))
0833     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_3to6_data, l_hM_dPhi_reco_Centrality_3to6_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 3-6%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_3to6'.format(plotdir))
0834     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_6to10_data, l_hM_dPhi_reco_Centrality_6to10_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 6-10%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_6to10'.format(plotdir))
0835     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_10to15_data, l_hM_dPhi_reco_Centrality_10to15_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 10-15%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_10to15'.format(plotdir))
0836     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_15to20_data, l_hM_dPhi_reco_Centrality_15to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 15-20%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_15to20'.format(plotdir))
0837     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_20to25_data, l_hM_dPhi_reco_Centrality_20to25_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 20-25%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_20to25'.format(plotdir))
0838     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_25to30_data, l_hM_dPhi_reco_Centrality_25to30_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 25-30%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_25to30'.format(plotdir))
0839     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_30to35_data, l_hM_dPhi_reco_Centrality_30to35_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 30-35%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_30to35'.format(plotdir))
0840     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_35to40_data, l_hM_dPhi_reco_Centrality_35to40_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 35-40%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_35to40'.format(plotdir))
0841     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_40to45_data, l_hM_dPhi_reco_Centrality_40to45_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 40-45%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_40to45'.format(plotdir))
0842     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_45to50_data, l_hM_dPhi_reco_Centrality_45to50_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 45-50%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_45to50'.format(plotdir))
0843     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_50to55_data, l_hM_dPhi_reco_Centrality_50to55_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 50-55%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_50to55'.format(plotdir))
0844     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_55to60_data, l_hM_dPhi_reco_Centrality_55to60_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 55-60%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_55to60'.format(plotdir))
0845     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_60to65_data, l_hM_dPhi_reco_Centrality_60to65_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 60-65%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_60to65'.format(plotdir))
0846     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_65to70_data, l_hM_dPhi_reco_Centrality_65to70_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 65-70%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_65to70'.format(plotdir))
0847     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_70to80_data, l_hM_dPhi_reco_Centrality_70to80_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 70-80%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_70to80'.format(plotdir))
0848     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_80to90_data, l_hM_dPhi_reco_Centrality_80to90_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 80-90%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_80to90'.format(plotdir))
0849     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_90to100_data, l_hM_dPhi_reco_Centrality_90to100_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_90to100'.format(plotdir))
0850     # Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_90to100_data, l_hM_dPhi_reco_Centrality_90to100_sim, padmargin, 'data', True, 100, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_90to100'.format(plotdir))
0851     # Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_0to5_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_0to5_MBDAsymLe0p75_sim, padmargin, 'data', True, 100, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 0-5%, |Asymm._{MBD}|#leq0.75','|Asymm._{MBD}|#leq0.75}'], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_0to5_MBDAsymLe0p75'.format(plotdir))
0852     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 0-3%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_0to3_MBDAsymLe0p75'.format(plotdir))
0853     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 3-6%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_3to6_MBDAsymLe0p75'.format(plotdir))
0854     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 6-10%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_6to10_MBDAsymLe0p75'.format(plotdir))
0855     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 10-15%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_10to15_MBDAsymLe0p75'.format(plotdir))
0856     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 15-20%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_15to20_MBDAsymLe0p75'.format(plotdir))
0857     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 20-25%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_20to25_MBDAsymLe0p75'.format(plotdir))
0858     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 25-30%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_25to30_MBDAsymLe0p75'.format(plotdir))
0859     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 30-35%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_30to35_MBDAsymLe0p75'.format(plotdir))
0860     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 35-40%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_35to40_MBDAsymLe0p75'.format(plotdir))
0861     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 40-45%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_40to45_MBDAsymLe0p75'.format(plotdir))
0862     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 45-50%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_45to50_MBDAsymLe0p75'.format(plotdir))
0863     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 50-55%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_50to55_MBDAsymLe0p75'.format(plotdir))
0864     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 55-60%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_55to60_MBDAsymLe0p75'.format(plotdir))
0865     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 60-65%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_60to65_MBDAsymLe0p75'.format(plotdir))
0866     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 65-70%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_65to70_MBDAsymLe0p75'.format(plotdir))
0867     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 70-80%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_70to80_MBDAsymLe0p75'.format(plotdir))
0868     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 80-90%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_80to90_MBDAsymLe0p75'.format(plotdir))
0869     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_90to100_MBDAsymLe0p75'.format(plotdir))
0870     # Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_data, l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_sim, padmargin, 'data', True, 100, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_90to100_MBDAsymLe0p75'.format(plotdir))
0871     # Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 100, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 0-5%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0872     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 0-3%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0873     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 3-6%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0874     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 6-10%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0875     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 10-15%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0876     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 15-20%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0877     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 20-25%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0878     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 25-30%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0879     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 30-35%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0880     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 35-40%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0881     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 40-45%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0882     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 45-50%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0883     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 50-55%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0884     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 55-60%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0885     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 60-65%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0886     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 65-70%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0887     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 70-80%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0888     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 80-90%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0889     Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 1E3, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0890     # Draw_1Dhist_datasimcomp(hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dPhi_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 100, 'Reco-tracklet #Delta#phi', '', False, simlegtext, ['Centrality 90-10%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dPhi_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0891     
0892     Draw_1Dhist_datasimcomp(hM_dR_reco_data, l_hM_dR_reco_sim, padmargin, 'data', True, 30, 'Reco-tracklet #DeltaR', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_dR'.format(plotdir))
0893     Draw_1Dhist_datasimcomp(hM_dR_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_dR_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', True, 100, 'Reco-tracklet #DeltaR', '', False, simlegtext, ['Centrality 0-70%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_dR_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0894     Draw_1Dhist_datasimcomp(hM_Eta_reco_data, l_hM_Eta_reco_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta'.format(plotdir))
0895     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 0-70%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_Eta_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0896     # Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 0-5%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_Eta_Centrality_0to5_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0897     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 0-3%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_0to3_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0898     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 3-6%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_3to6_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0899     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 6-10%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_6to10_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0900     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 10-15%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_10to15_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0901     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 15-20%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_15to20_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0902     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 20-25%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_20to25_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0903     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 25-30%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_25to30_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0904     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 30-35%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_30to35_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0905     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 35-40%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_35to40_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0906     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 40-45%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_40to45_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0907     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 45-50%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_45to50_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0908     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 50-55%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_50to55_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0909     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 55-60%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_55to60_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0910     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 60-65%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_60to65_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0911     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 65-70%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_65to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0912     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 70-80%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_70to80_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0913     Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 80-90%', str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Eta_Centrality_80to90_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0914     # Draw_1Dhist_datasimcomp(hM_Eta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Eta_reco_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2, 'Reco-tracklet #eta', '', False, simlegtext, ['Centrality 90-100%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_Eta_Centrality_90to100_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0915     
0916     for hM_Phi_reco_sim in l_hM_Phi_reco_sim:
0917         hM_Phi_reco_sim.GetYaxis().SetMaxDigits(2)
0918     Draw_1Dhist_datasimcomp(hM_Phi_reco_data, l_hM_Phi_reco_sim, padmargin, 'data', False, 1.8, 'Reco-tracklet #phi', '', False, simlegtext, [str_evtsel, '{}, {}'.format(str_zvtxsel,str_clussel)], './DataSimComp/{}/RecoTracklet_Phi'.format(plotdir))
0919     Draw_1Dhist_datasimcomp(hM_Phi_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_data, l_hM_Phi_reco_Centrality0to70_MBDAsymLe0p75_VtxZm20to20_sim, padmargin, 'data', False, 2.0, 'Reco-tracklet #phi', '', False, simlegtext, ['Centrality 0-70%', str_evtsel, '{}'.format('{}, {}'.format(str_zvtxsel,str_clussel))], './DataSimComp/{}/RecoTracklet_Phi_Centrality0to70_MBDAsymLe0p75_VtxZm20to20'.format(plotdir))
0920     
0921     # Draw_2Dhist_datasimcomp(hdata, hsim, logz, norm, rmargin, XaxisName, YaxisName, outname)
0922     Draw_2Dhist_datasimcomp(hM_cluseta_clusphisize_data, l_hM_cluseta_clusphisize_sim[0], False, 'data', 0.1, 'Cluster #eta', 'Cluster #phi size', './DataSimComp/{}/ClusEta_ClusPhiSize'.format(plotdir))
0923     Draw_2Dhist_datasimcomp(hM_clusphi_clusphisize_data, l_hM_clusphi_clusphisize_sim[0], False, 'data', 0.1, 'Cluster #phi', 'Cluster #phi size', './DataSimComp/{}/ClusPhi_ClusPhiSize'.format(plotdir))
0924 
0925