analysis_INTT/plot/plotCluster.py

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, kRainBow, kBird, kHAlignRight, kVAlignBottom, kHAlignLeft, kVAlignTop, kThermometer, gROOT, gStyle
0008 import numpy
0009 import math
0010 import glob
0011 from plotUtil import *
0012
0013 # gROOT.LoadMacro('./sPHENIXStyle/sPhenixStyle.C')
0014 # gROOT.ProcessLine('SetsPhenixStyle()')
0015 gROOT.SetBatch(True)
0016
0017 if __name__ == '__main__':
0018     parser = OptionParser()
0019     parser.add_option("-f", "--histdir", dest="histdir", type="string", default='/sphenix/user/hjheng/TrackletAna/analysis_INTT/plot/hists/data_Run20869/Cluster', help="Input file name")
0020     parser.add_option("-d", "--plotdir", dest="plotdir", type="string", default='/sphenix/user/hjheng/TrackletAna/analysis_INTT/plot/RecoCluster/data_Run20869', help="Plot directory")
0021     parser.add_option("-s", "--isdata", dest="isdata", action="store_true", default=False, help="Is data")
0022     parser.add_option("-p", "--preliminary", dest="preliminary", action="store_true", default=False, help="Preliminary")
0023     (opt, args) = parser.parse_args()
0024     print('opt: {}'.format(opt))
0025
0026     histdir = opt.histdir
0027     plotdir = opt.plotdir
0028     isdata = opt.isdata
0029     preliminary = opt.preliminary
0030     os.makedirs(plotdir, exist_ok=True)
0031
0032     if os.path.isfile("{}/hists_merged.root".format(histdir)):
0033         # os.system("rm {}/hists_merged.root".format(histdir))
0034         os.system("hadd -f -j 20 {}/hists_merged.root {}/hists_*.root".format(histdir, histdir))
0035     else:
0036         os.system("hadd -f -j 20 {}/hists_merged.root {}/hists_*.root".format(histdir, histdir))
0037
0038     hM_ClusZ_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_all')
0039     hM_ClusZ_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_layer1')
0040     hM_ClusZ_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_layer2')
0041     hM_ClusEtaPV_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_all')
0042     hM_ClusEtaPV_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer1')
0043     hM_ClusEtaPV_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer2')
0044     hM_ClusEtaPV_all_ClusADCg35 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_all_ClusADCg35')
0045     hM_ClusEtaPV_layer1_ClusADCg35 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer1_ClusADCg35')
0046     hM_ClusEtaPV_layer2_ClusADCg35 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer2_ClusADCg35')
0047     hM_ClusEtaPV_all_weiClusADCoverPhiSize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_all_weiClusADCoverPhiSize')
0048     hM_ClusEtaPV_layer1_weiClusADCoverPhiSize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer1_weiClusADCoverPhiSize')
0049     hM_ClusEtaPV_layer2_weiClusADCoverPhiSize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer2_weiClusADCoverPhiSize')
0050     hM_ClusEtaPV_all_ClusADCg35_weiClusADCoverPhiSize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_all_ClusADCg35_weiClusADCoverPhiSize')
0051     hM_ClusEtaPV_layer1_ClusADCg35_weiClusADCoverPhiSize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer1_ClusADCg35_weiClusADCoverPhiSize')
0052     hM_ClusEtaPV_layer2_ClusADCg35_weiClusADCoverPhiSize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_layer2_ClusADCg35_weiClusADCoverPhiSize')
0053     hM_ClusPhiPV_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_all')
0054     hM_ClusPhiPV_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_layer1')
0055     hM_ClusPhiPV_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_layer2')
0056     hM_ClusPhiPV_ClusPhiSize43or46 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_ClusPhiSize43or46')
0057     hM_ClusADC_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusADC_all')
0058     hM_ClusADC_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusADC_layer1')
0059     hM_ClusADC_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusADC_layer2')
0060     # hM_ClusZSize_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_all')
0061     # hM_ClusZSize_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_layer1')
0062     # hM_ClusZSize_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_layer2')
0063     hM_ClusPhiSize_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_all')
0064     hM_ClusPhiSize_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_layer1')
0065     hM_ClusPhiSize_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_layer2')
0066
0067     hM_ClusX_ClusY_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusX_ClusY_all')
0068     hM_ClusX_ClusY_ClusPhiSize43or46 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusX_ClusY_ClusPhiSize43or46')
0069     hM_ClusX_ClusY_ClusPhiSize43 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusX_ClusY_ClusPhiSize43')
0070     hM_ClusX_ClusY_ClusPhiSize46 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusX_ClusY_ClusPhiSize46')
0071     hM_ClusX_ClusY_all_weiphisize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusX_ClusY_all_weiphisize')
0072     hM_ClusX_ClusY_all_weiclusadc = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusX_ClusY_all_weiclusadc')
0073     hM_ClusZ_ClusPhiPV_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_ClusPhiPV_all')
0074     hM_ClusZ_ClusPhiPV_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_ClusPhiPV_layer1')
0075     hM_ClusZ_ClusPhiPV_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_ClusPhiPV_layer2')
0076     hM_ClusZ_ClusPhiPV_all_coarse = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_ClusPhiPV_all_coarse')
0077     hM_ClusZ_ClusPhiPV_all_coarse_weiphisize = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_ClusPhiPV_all_coarse_weiphisize')
0078     hM_ClusZ_ClusPhiPV_all_coarse_weiclusadc = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZ_ClusPhiPV_all_coarse_weiclusadc')
0079     # hM_ClusEtaPV_ClusZSize_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusZSize_all')
0080     # hM_ClusEtaPV_ClusZSize_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusZSize_layer1')
0081     # hM_ClusEtaPV_ClusZSize_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusZSize_layer2')
0082     hM_ClusPhiPV_ClusPhiSize_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_ClusPhiSize_all')
0083     hM_ClusPhiPV_ClusPhiSize_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_ClusPhiSize_layer1')
0084     hM_ClusPhiPV_ClusPhiSize_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_ClusPhiSize_layer2')
0085     hM_ClusEtaPV_ClusADC_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADC_all')
0086     hM_ClusEtaPV_ClusADC_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADC_layer1')
0087     hM_ClusEtaPV_ClusADC_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADC_layer2')
0088     hM_ClusEtaPV_ClusADC_all_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADC_all_zoomin')
0089     hM_ClusEtaPV_ClusADC_layer1_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADC_layer1_zoomin')
0090     hM_ClusEtaPV_ClusADC_layer2_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADC_layer2_zoomin')
0091     hM_ClusEtaPV_ClusPhiSize_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusPhiSize_all')
0092     hM_ClusEtaPV_ClusPhiSize_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusPhiSize_layer1')
0093     hM_ClusEtaPV_ClusPhiSize_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusPhiSize_layer2')
0094     hM_ClusEtaPV_ClusPhiSize_all_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusPhiSize_all_zoomin')
0095     hM_ClusEtaPV_ClusPhiSize_layer1_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusPhiSize_layer1_zoomin')
0096     hM_ClusEtaPV_ClusPhiSize_layer2_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusPhiSize_layer2_zoomin')
0097     hM_ClusPhiPV_ClusADC_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_ClusADC_all')
0098     hM_ClusPhiPV_ClusADC_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_ClusADC_layer1')
0099     hM_ClusPhiPV_ClusADC_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiPV_ClusADC_layer2')
0100     # hM_ClusZSize_ClusPhiSize_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_ClusPhiSize_all')
0101     # hM_ClusZSize_ClusPhiSize_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_ClusPhiSize_layer1')
0102     # hM_ClusZSize_ClusPhiSize_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_ClusPhiSize_layer2')
0103     # hM_ClusZSize_ClusADC_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_ClusADC_all')
0104     # hM_ClusZSize_ClusADC_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_ClusADC_layer1')
0105     # hM_ClusZSize_ClusADC_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusZSize_ClusADC_layer2')
0106     hM_ClusPhiSize_ClusADC_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_ClusADC_all')
0107     hM_ClusPhiSize_ClusADC_layer1 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_ClusADC_layer1')
0108     hM_ClusPhiSize_ClusADC_layer2 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_ClusADC_layer2')
0109     hM_ClusPhiSize_ClusADC_all_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_ClusADC_all_zoomin')
0110     hM_ClusPhiSize_ClusADC_layer1_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_ClusADC_layer1_zoomin')
0111     hM_ClusPhiSize_ClusADC_layer2_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusPhiSize_ClusADC_layer2_zoomin')
0112     hM_ClusEtaPV_ClusADCoverClusPhiSize_all_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADCoverClusPhiSize_all_zoomin')
0113     hM_ClusEtaPV_ClusADCoverClusPhiSize_layer1_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADCoverClusPhiSize_layer1_zoomin')
0114     hM_ClusEtaPV_ClusADCoverClusPhiSize_layer2_zoomin = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADCoverClusPhiSize_layer2_zoomin')
0115     hM_ClusEtaPV_ClusADCoverClusPhiSize_all_ClusADCg35 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADCoverClusPhiSize_all_ClusADCg35')
0116     hM_ClusEtaPV_ClusADCoverClusPhiSize_layer1_ClusADCg35 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADCoverClusPhiSize_layer1_ClusADCg35')
0117     hM_ClusEtaPV_ClusADCoverClusPhiSize_layer2_ClusADCg35 = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_ClusEtaPV_ClusADCoverClusPhiSize_layer2_ClusADCg35')
0118
0119     hM_mutualdRcluster_all = GetHistogram('{}/hists_merged.root'.format(histdir), 'hM_mutualdRcluster_all')
0120
0121     # Draw_1Dhist(hist, IsData, norm1, logy, ymaxscale, XaxisName, Ytitle_unit, outname)
0122     Draw_1Dhist(hM_ClusZ_all, isdata, False, False, 1.3, 'Cluster Z [cm]', 'cm', '{}/ClusZ_all'.format(plotdir))
0123     Draw_1Dhist(hM_ClusZ_layer1, isdata, False, False, 1.3, 'Cluster Z [cm]', 'cm', '{}/ClusZ_layer1'.format(plotdir))
0124     Draw_1Dhist(hM_ClusZ_layer2, isdata, False, False, 1.3, 'Cluster Z [cm]', 'cm', '{}/ClusZ_layer2'.format(plotdir))
0125     Draw_1Dhist(hM_ClusEtaPV_all, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_all'.format(plotdir))
0126     Draw_1Dhist(hM_ClusEtaPV_layer1, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer1'.format(plotdir))
0127     Draw_1Dhist(hM_ClusEtaPV_layer2, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer2'.format(plotdir))
0128     Draw_1Dhist(hM_ClusEtaPV_all_ClusADCg35, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_all_ClusADCg35'.format(plotdir))
0129     Draw_1Dhist(hM_ClusEtaPV_layer1_ClusADCg35, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer1_ClusADCg35'.format(plotdir))
0130     Draw_1Dhist(hM_ClusEtaPV_layer2_ClusADCg35, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer2_ClusADCg35'.format(plotdir))
0131     Draw_1Dhist(hM_ClusEtaPV_all_weiClusADCoverPhiSize, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_all_weiClusADCoverPhiSize'.format(plotdir))
0132     Draw_1Dhist(hM_ClusEtaPV_layer1_weiClusADCoverPhiSize, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer1_weiClusADCoverPhiSize'.format(plotdir))
0133     Draw_1Dhist(hM_ClusEtaPV_layer2_weiClusADCoverPhiSize, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer2_weiClusADCoverPhiSize'.format(plotdir))
0134     Draw_1Dhist(hM_ClusEtaPV_all_ClusADCg35_weiClusADCoverPhiSize, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_all_ClusADCg35_weiClusADCoverPhiSize'.format(plotdir))
0135     Draw_1Dhist(hM_ClusEtaPV_layer1_ClusADCg35_weiClusADCoverPhiSize, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer1_ClusADCg35_weiClusADCoverPhiSize'.format(plotdir))
0136     Draw_1Dhist(hM_ClusEtaPV_layer2_ClusADCg35_weiClusADCoverPhiSize, isdata, False, False, 1.3, 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer2_ClusADCg35_weiClusADCoverPhiSize'.format(plotdir))
0137     Draw_1Dhist(hM_ClusPhiPV_all, isdata, False, False, 1.3, 'Cluster #phi (w.r.t PV)', '', '{}/ClusPhiPV_all'.format(plotdir))
0138     Draw_1Dhist(hM_ClusPhiPV_layer1, isdata, False, False, 1.3, 'Cluster #phi (w.r.t PV)', '', '{}/ClusPhiPV_layer1'.format(plotdir))
0139     Draw_1Dhist(hM_ClusPhiPV_layer2, isdata, False, False, 1.3, 'Cluster #phi (w.r.t PV)', '', '{}/ClusPhiPV_layer2'.format(plotdir))
0140     Draw_1Dhist(hM_ClusPhiPV_ClusPhiSize43or46, isdata, False, False, 1.3, 'Cluster #phi (w.r.t PV)', '', '{}/ClusPhiPV_ClusPhiSize43or46'.format(plotdir))
0141     Draw_1Dhist(hM_ClusADC_all, isdata, False, True, 100, 'Cluster ADC', '', '{}/ClusADC_all'.format(plotdir))
0142     Draw_1Dhist(hM_ClusADC_layer1, isdata, False, True, 100, 'Cluster ADC', '', '{}/ClusADC_layer1'.format(plotdir))
0143     Draw_1Dhist(hM_ClusADC_layer2, isdata, False, True, 100, 'Cluster ADC', '', '{}/ClusADC_layer2'.format(plotdir))
0144     # Draw_1Dhist(hM_ClusZSize_all, isdata, False, True, 100, 'Cluster Z size', '', '{}/ClusZSize_all'.format(plotdir))
0145     # Draw_1Dhist(hM_ClusZSize_layer1, isdata, False, True, 100, 'Cluster Z size', '', '{}/ClusZSize_layer1'.format(plotdir))
0146     # Draw_1Dhist(hM_ClusZSize_layer2, isdata, False, True, 100, 'Cluster Z size', '', '{}/ClusZSize_layer2'.format(plotdir))
0147     Draw_1Dhist(hM_ClusPhiSize_all, isdata, False, True, 100, 'Cluster #phi size', '', '{}/ClusPhiSize_all'.format(plotdir))
0148     Draw_1Dhist(hM_ClusPhiSize_layer1, isdata, False, True, 100, 'Cluster #phi size', '', '{}/ClusPhiSize_layer1'.format(plotdir))
0149     Draw_1Dhist(hM_ClusPhiSize_layer2, isdata, False, True, 100, 'Cluster #phi size', '', '{}/ClusPhiSize_layer2'.format(plotdir))
0150     hM_ClusEtaPV_all_avgclusadcoverphisize = hM_ClusEtaPV_all_weiClusADCoverPhiSize.Clone('hM_ClusEtaPV_all_avgclusadcoverphisize')
0151     hM_ClusEtaPV_layer1_avgclusadcoverphisize = hM_ClusEtaPV_layer1_weiClusADCoverPhiSize.Clone('hM_ClusEtaPV_layer1_avgclusadcoverphisize')
0152     hM_ClusEtaPV_layer2_avgclusadcoverphisize = hM_ClusEtaPV_layer2_weiClusADCoverPhiSize.Clone('hM_ClusEtaPV_layer2_avgclusadcoverphisize')
0153     hM_ClusEtaPV_all_avgclusadcoverphisize.Divide(hM_ClusEtaPV_all)
0154     hM_ClusEtaPV_layer1_avgclusadcoverphisize.Divide(hM_ClusEtaPV_layer1)
0155     hM_ClusEtaPV_layer2_avgclusadcoverphisize.Divide(hM_ClusEtaPV_layer2)
0156     print (hM_ClusEtaPV_all_avgclusadcoverphisize.GetBinContent(hM_ClusEtaPV_all_avgclusadcoverphisize.GetMaximumBin()))
0157     # reset the histogram range
0158     hM_ClusEtaPV_all_avgclusadcoverphisize.GetYaxis().SetRangeUser(0, hM_ClusEtaPV_all_avgclusadcoverphisize.GetBinContent(hM_ClusEtaPV_all_avgclusadcoverphisize.GetMaximumBin())*1.3)
0159     hM_ClusEtaPV_layer1_avgclusadcoverphisize.GetYaxis().SetRangeUser(0, hM_ClusEtaPV_layer1_avgclusadcoverphisize.GetBinContent(hM_ClusEtaPV_layer1_avgclusadcoverphisize.GetMaximumBin())*1.3)
0160     hM_ClusEtaPV_layer2_avgclusadcoverphisize.GetYaxis().SetRangeUser(0, hM_ClusEtaPV_layer2_avgclusadcoverphisize.GetBinContent(hM_ClusEtaPV_layer2_avgclusadcoverphisize.GetMaximumBin())*1.3)
0161     Draw_1Dhist(hM_ClusEtaPV_all_avgclusadcoverphisize, isdata, False, False, 1., 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_all_avgclusadcoverphisize'.format(plotdir))
0162     Draw_1Dhist(hM_ClusEtaPV_layer1_avgclusadcoverphisize, isdata, False, False, 1., 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer1_avgclusadcoverphisize'.format(plotdir))
0163     Draw_1Dhist(hM_ClusEtaPV_layer2_avgclusadcoverphisize, isdata, False, False, 1., 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer2_avgclusadcoverphisize'.format(plotdir))
0164     hM_ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize = hM_ClusEtaPV_all_ClusADCg35_weiClusADCoverPhiSize.Clone('hM_ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize')
0165     hM_ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize = hM_ClusEtaPV_layer1_ClusADCg35_weiClusADCoverPhiSize.Clone('hM_ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize')
0166     hM_ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize = hM_ClusEtaPV_layer2_ClusADCg35_weiClusADCoverPhiSize.Clone('hM_ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize')
0167     hM_ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize.Divide(hM_ClusEtaPV_all_ClusADCg35)
0168     hM_ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize.Divide(hM_ClusEtaPV_layer1_ClusADCg35)
0169     hM_ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize.Divide(hM_ClusEtaPV_layer2_ClusADCg35)
0170     hM_ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize.GetYaxis().SetRangeUser(0, hM_ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize.GetBinContent(hM_ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize.GetMaximumBin())*1.3)
0171     hM_ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize.GetYaxis().SetRangeUser(0, hM_ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize.GetBinContent(hM_ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize.GetMaximumBin())*1.3)
0172     hM_ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize.GetYaxis().SetRangeUser(0, hM_ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize.GetBinContent(hM_ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize.GetMaximumBin())*1.3)
0173     Draw_1Dhist(hM_ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize, isdata, False, False, 1., 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_all_ClusADCg35_avgclusadcoverphisize'.format(plotdir))
0174     Draw_1Dhist(hM_ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize, isdata, False, False, 1., 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer1_ClusADCg35_avgclusadcoverphisize'.format(plotdir))
0175     Draw_1Dhist(hM_ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize, isdata, False, False, 1., 'Cluster #eta (w.r.t PV)', '', '{}/ClusEtaPV_layer2_ClusADCg35_avgclusadcoverphisize'.format(plotdir))
0176
0177     Draw_1Dhist(hM_mutualdRcluster_all, isdata, False, False, 1.3, 'dR(cluster, cluster)', '', '{}/mutualdRcluster_all'.format(plotdir))
0178
0179     # Draw_2Dhist(hist, IsData, logz, norm1, rmargin, XaxisName, YaxisName, drawopt, outname)
0180     gStyle.SetPalette(kBird)
0181     # gStyle.SetPalette(kRainBow)
0182     Draw_2Dhist(hM_ClusX_ClusY_all, isdata, False, False, 0.15, 'Cluster X [cm]', 'Cluster Y [cm]', 'Entries', 'colz', '{}/ClusX_ClusY_all'.format(plotdir))
0183     Draw_2Dhist(hM_ClusX_ClusY_ClusPhiSize43or46, isdata, True, False, 0.15, 'Cluster X [cm]', 'Cluster Y [cm]', 'Entries', 'colz', '{}/ClusX_ClusY_ClusPhiSize43or46'.format(plotdir))
0184     Draw_2Dhist(hM_ClusX_ClusY_ClusPhiSize43, isdata, True, False, 0.15, 'Cluster X [cm]', 'Cluster Y [cm]', 'Entries', 'colz', '{}/ClusX_ClusY_ClusPhiSize43'.format(plotdir))
0185     Draw_2Dhist(hM_ClusX_ClusY_ClusPhiSize46, isdata, True, False, 0.15, 'Cluster X [cm]', 'Cluster Y [cm]', 'Entries', 'colz', '{}/ClusX_ClusY_ClusPhiSize46'.format(plotdir))
0186     Draw_2Dhist(hM_ClusX_ClusY_all_weiphisize, isdata, False, False, 0.15, 'Cluster X [cm]', 'Cluster Y [cm]','Entries', 'colz', '{}/ClusX_ClusY_all_weiphisize'.format(plotdir))
0187     hM_ClusX_ClusY_all_avgphisize = hM_ClusX_ClusY_all_weiphisize.Clone('hM_ClusX_ClusY_all_avgphisize')
0188     hM_ClusX_ClusY_all_avgclusadc = hM_ClusX_ClusY_all_weiclusadc.Clone('hM_ClusX_ClusY_all_avgclusadc')
0189     hM_ClusX_ClusY_all_avgphisize.Divide(hM_ClusX_ClusY_all)
0190     hM_ClusX_ClusY_all_avgclusadc.Divide(hM_ClusX_ClusY_all)
0191     Draw_2Dhist(hM_ClusX_ClusY_all_avgphisize, isdata, False, False, 0.15, 'Cluster X [cm]', 'Cluster Y [cm]','Average cluster #phi size', 'colz', '{}/ClusX_ClusY_all_avgphisize'.format(plotdir))
0192     Draw_2Dhist(hM_ClusX_ClusY_all_avgclusadc, isdata, False, False, 0.15, 'Cluster X [cm]', 'Cluster Y [cm]', 'Average cluster ADC', 'colz', '{}/ClusX_ClusY_all_avgclusadc'.format(plotdir))
0193
0194     Draw_2Dhist(hM_ClusZ_ClusPhiPV_all, isdata, False, False, 0.16, 'Cluster Z [cm]', 'Cluster #phi (w.r.t PV)', 'Entries', 'colz', '{}/ClusZ_ClusPhiPV_all'.format(plotdir))
0195     Draw_2Dhist(hM_ClusZ_ClusPhiPV_layer1, isdata, False, False, 0.16, 'Cluster Z [cm]', 'Cluster #phi (w.r.t PV)', 'Entries', 'colz', '{}/ClusZ_ClusPhiPV_layer1'.format(plotdir))
0196     Draw_2Dhist(hM_ClusZ_ClusPhiPV_layer2, isdata, False, False, 0.16, 'Cluster Z [cm]', 'Cluster #phi (w.r.t PV)', 'Entries', 'colz', '{}/ClusZ_ClusPhiPV_layer2'.format(plotdir))
0197     Draw_2Dhist(hM_ClusZ_ClusPhiPV_all_coarse, isdata, False, False, 0.16, 'Cluster Z [cm]', 'Cluster #phi (w.r.t PV)', 'Entries', 'colz', '{}/ClusZ_ClusPhiPV_all_coarse'.format(plotdir))
0198     hM_ClusZ_ClusPhiPV_all_coarse_avgphisize = hM_ClusZ_ClusPhiPV_all_coarse_weiphisize.Clone('hM_ClusZ_ClusPhiPV_all_coarse_avgphisize')
0199     hM_ClusZ_ClusPhiPV_all_coarse_avgclusadc = hM_ClusZ_ClusPhiPV_all_coarse_weiclusadc.Clone('hM_ClusZ_ClusPhiPV_all_coarse_avgclusadc')
0200     hM_ClusZ_ClusPhiPV_all_coarse_avgphisize.Divide(hM_ClusZ_ClusPhiPV_all_coarse)
0201     hM_ClusZ_ClusPhiPV_all_coarse_avgclusadc.Divide(hM_ClusZ_ClusPhiPV_all_coarse)
0202     Draw_2Dhist(hM_ClusZ_ClusPhiPV_all_coarse_avgphisize, isdata, False, False, 0.16, 'Cluster Z [cm]', 'Cluster #phi (w.r.t PV)', 'Average cluster #phi size', 'colz', '{}/ClusZ_ClusPhiPV_all_coarse_avgphisize'.format(plotdir))
0203     Draw_2Dhist(hM_ClusZ_ClusPhiPV_all_coarse_avgclusadc, isdata, False, False, 0.16, 'Cluster Z [cm]', 'Cluster #phi (w.r.t PV)', 'Average cluster ADC', 'colz', '{}/ClusZ_ClusPhiPV_all_coarse_avgclusadc'.format(plotdir))
0204     # Draw_2Dhist(hM_ClusEtaPV_ClusZSize_all, isdata, True, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster Z size', 'colz', '{}/ClusEtaPV_ClusZSize_all'.format(plotdir))
0205     # Draw_2Dhist(hM_ClusEtaPV_ClusZSize_layer1, isdata, True, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster Z size', 'colz', '{}/ClusEtaPV_ClusZSize_layer1'.format(plotdir))
0206     # Draw_2Dhist(hM_ClusEtaPV_ClusZSize_layer2, isdata, True, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster Z size', 'colz', '{}/ClusEtaPV_ClusZSize_layer2'.format(plotdir))
0207     Draw_2Dhist(hM_ClusPhiPV_ClusPhiSize_all, isdata, False, False, 0.16, 'Cluster #phi (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusPhiPV_ClusPhiSize_all'.format(plotdir))
0208     Draw_2Dhist(hM_ClusPhiPV_ClusPhiSize_layer1, isdata, False, False, 0.16, 'Cluster #phi (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusPhiPV_ClusPhiSize_layer1'.format(plotdir))
0209     Draw_2Dhist(hM_ClusPhiPV_ClusPhiSize_layer2, isdata, False, False, 0.16, 'Cluster #phi (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusPhiPV_ClusPhiSize_layer2'.format(plotdir))
0210     Draw_2Dhist(hM_ClusEtaPV_ClusADC_all, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusEtaPV_ClusADC_all'.format(plotdir))
0211     Draw_2Dhist(hM_ClusEtaPV_ClusADC_layer1, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusEtaPV_ClusADC_layer1'.format(plotdir))
0212     Draw_2Dhist(hM_ClusEtaPV_ClusADC_layer2, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusEtaPV_ClusADC_layer2'.format(plotdir))
0213     Draw_2Dhist(hM_ClusEtaPV_ClusADC_all_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusEtaPV_ClusADC_all_zoomin'.format(plotdir))
0214     Draw_2Dhist(hM_ClusEtaPV_ClusADC_layer1_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusEtaPV_ClusADC_layer1_zoomin'.format(plotdir))
0215     Draw_2Dhist(hM_ClusEtaPV_ClusADC_layer2_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusEtaPV_ClusADC_layer2_zoomin'.format(plotdir))
0216     Draw_2Dhist(hM_ClusEtaPV_ClusPhiSize_all, isdata, True, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusPhiSize_all'.format(plotdir))
0217     Draw_2Dhist(hM_ClusEtaPV_ClusPhiSize_layer1, isdata, True, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusPhiSize_layer1'.format(plotdir))
0218     Draw_2Dhist(hM_ClusEtaPV_ClusPhiSize_layer2, isdata, True, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusPhiSize_layer2'.format(plotdir))
0219     Draw_2Dhist(hM_ClusEtaPV_ClusPhiSize_all_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusPhiSize_all_zoomin'.format(plotdir))
0220     Draw_2Dhist(hM_ClusEtaPV_ClusPhiSize_layer1_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusPhiSize_layer1_zoomin'.format(plotdir))
0221     Draw_2Dhist(hM_ClusEtaPV_ClusPhiSize_layer2_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusPhiSize_layer2_zoomin'.format(plotdir))
0222     Draw_2Dhist(hM_ClusPhiPV_ClusADC_all, isdata, False, False, 0.16, 'Cluster #phi (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiPV_ClusADC_all'.format(plotdir))
0223     Draw_2Dhist(hM_ClusPhiPV_ClusADC_layer1, isdata, False, False, 0.16, 'Cluster #phi (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiPV_ClusADC_layer1'.format(plotdir))
0224     Draw_2Dhist(hM_ClusPhiPV_ClusADC_layer2, isdata, False, False, 0.16, 'Cluster #phi (w.r.t PV)', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiPV_ClusADC_layer2'.format(plotdir))
0225     # Draw_2Dhist(hM_ClusZSize_ClusPhiSize_all, isdata, True, False, 0.16, 'Cluster Z size', 'Cluster #phi size', 'colz', '{}/ClusZSize_ClusPhiSize_all'.format(plotdir))
0226     # Draw_2Dhist(hM_ClusZSize_ClusPhiSize_layer1, isdata, True, False, 0.16, 'Cluster Z size', 'Cluster #phi size', 'colz', '{}/ClusZSize_ClusPhiSize_layer1'.format(plotdir))
0227     # Draw_2Dhist(hM_ClusZSize_ClusPhiSize_layer2, isdata, True, False, 0.16, 'Cluster Z size', 'Cluster #phi size', 'colz', '{}/ClusZSize_ClusPhiSize_layer2'.format(plotdir))
0228     # Draw_2Dhist(hM_ClusZSize_ClusADC_all, isdata, True, False, 0.16, 'Cluster Z size', 'Cluster #phi size', 'colz', '{}/ClusZSize_ClusADC_all'.format(plotdir))
0229     # Draw_2Dhist(hM_ClusZSize_ClusADC_layer1, isdata, True, False, 0.16, 'Cluster Z size', 'Cluster ADC', 'colz', '{}/ClusZSize_ClusADC_layer1'.format(plotdir))
0230     # Draw_2Dhist(hM_ClusZSize_ClusADC_layer2, isdata, True, False, 0.16, 'Cluster Z size', 'Cluster ADC', 'colz', '{}/ClusZSize_ClusADC_layer2'.format(plotdir))
0231     Draw_2Dhist(hM_ClusPhiSize_ClusADC_all, isdata, True, False, 0.16, 'Cluster #phi size', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiSize_ClusADC_all'.format(plotdir))
0232     Draw_2Dhist(hM_ClusPhiSize_ClusADC_layer1, isdata, True, False, 0.16, 'Cluster #phi size', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiSize_ClusADC_layer1'.format(plotdir))
0233     Draw_2Dhist(hM_ClusPhiSize_ClusADC_layer2, isdata, True, False, 0.16, 'Cluster #phi size', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiSize_ClusADC_layer2'.format(plotdir))
0234     Draw_2Dhist(hM_ClusPhiSize_ClusADC_all_zoomin, isdata, False, False, 0.16, 'Cluster #phi size', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiSize_ClusADC_all_zoomin'.format(plotdir))
0235     Draw_2Dhist(hM_ClusPhiSize_ClusADC_layer1_zoomin, isdata, False, False, 0.16, 'Cluster #phi size', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiSize_ClusADC_layer1_zoomin'.format(plotdir))
0236     Draw_2Dhist(hM_ClusPhiSize_ClusADC_layer2_zoomin, isdata, False, False, 0.16, 'Cluster #phi size', 'Cluster ADC', 'Entries', 'colz', '{}/ClusPhiSize_ClusADC_layer2_zoomin'.format(plotdir))
0237     Draw_2Dhist(hM_ClusEtaPV_ClusADCoverClusPhiSize_all_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC/Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusADCoverClusPhiSize_all_zoomin'.format(plotdir))
0238     Draw_2Dhist(hM_ClusEtaPV_ClusADCoverClusPhiSize_layer1_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC/Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusADCoverClusPhiSize_layer1_zoomin'.format(plotdir))
0239     Draw_2Dhist(hM_ClusEtaPV_ClusADCoverClusPhiSize_layer2_zoomin, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC/Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusADCoverClusPhiSize_layer2_zoomin'.format(plotdir))
0240     Draw_2Dhist(hM_ClusEtaPV_ClusADCoverClusPhiSize_all_ClusADCg35, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC/Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusADCoverClusPhiSize_all_ClusADCg35'.format(plotdir))
0241     Draw_2Dhist(hM_ClusEtaPV_ClusADCoverClusPhiSize_layer1_ClusADCg35, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC/Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusADCoverClusPhiSize_layer1_ClusADCg35'.format(plotdir))
0242     Draw_2Dhist(hM_ClusEtaPV_ClusADCoverClusPhiSize_layer2_ClusADCg35, isdata, False, False, 0.16, 'Cluster #eta (w.r.t PV)', 'Cluster ADC/Cluster #phi size', 'Entries', 'colz', '{}/ClusEtaPV_ClusADCoverClusPhiSize_layer2_ClusADCg35'.format(plotdir))