void EP_ana()
{
    gROOT->SetStyle("Plain");
    gStyle->SetOptStat(0);
    TFile* f = TFile::Open("V2_RP_RAW.root");

    TH1F* hEvtVsRP[2];
    hEvtVsRP[1] = (TH1F*)f->Get("compEvtPlane/evtplane_neg_versus_RP_handcorr");

    hEvtVsRP[0]->Draw();
    hEvtVsRP[0]->GetXaxis()->SetTitle("#Delta#phi");
    hEvtVsRP[0]->GetYaxis()->SetTitle("# of Events");
    hEvtVsRP[0]->GetYaxis()->SetLabelSize(0.03);
    hEvtVsRP[0]->SetTitle("#Psi^{+}_{2}- #Psi_{rp}");

    TCanvas* c1 = new TCanvas;
    hEvtVsRP[1]->Draw();
    hEvtVsRP[1]->GetXaxis()->SetTitle("#Delta#phi");
    hEvtVsRP[1]->GetYaxis()->SetTitle("# of Events");
    hEvtVsRP[1]->GetYaxis()->SetLabelSize(0.03);

   TCanvas* c1 = new TCanvas;
    ResolutionHandmade->SetLineColor(kRed);
    ResolutionHandmade->SetMarkerStyle(29);
    ResolutionHandmade->SetMarkerColor(kRed);
    ResolutionHandmade->SetTitle("Profile of sub-event Resolution");


    
    ResolutionHandcorr->SetLineColor(kBlue);
    ResolutionHandcorr->SetMarkerStyle(20);
    ResolutionHandcorr->SetMarkerColor(kBlue);
    ResolutionHandcorr->SetMarkerSize(1);
    ResolutionHandcorr->SetTitle("Profile of corrected sub-event Resolution");

    ResolutionHandcorr->GetYaxis()->SetTitle("#sqrt{ <cos(2(#Psi^{+}_{2} - #Psi^{-}_{2}))> }");
    ResolutionHandcorr->GetXaxis()->SetTitle("Centrality (%)");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(1,"10");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(2,"20");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(3,"30");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(4,"40");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(5,"50");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(6,"60");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(7,"70");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(8,"80");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(9,"90");
    ResolutionHandcorr->GetXaxis()->SetBinLabel(10,"100");


}

-- GunmoRyu - 04 Oct 2011