#include "TCanvas.h"
#include "TStyle.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TGaxis.h"
#include "TRandom.h"
#include "TFile.h"
#include "TTree.h"
#include "TPaveStats.h"
#include "TPaveLabel.h"
#include "TPaletteAxis.h"
#include "TLegend.h"
#include "TBranch.h"
#include "TSystem.h"
#include "TLorentzVector.h"
#include "TVector3.h"
#include "TEventList.h"
#include "TGraph.h"
#include "THashList.h"
#include "TMultiGraph.h"
#include "TLine.h"
#include "TText.h"
#include "TLatex.h"
#include "TTreeReader.h"
#include "TTreeReaderValue.h"
#include "TStopwatch.h"
#include "TExec.h"
#include <vector>
#include <set>
#include <string>
#include <algorithm>
#include <cmath>
using namespace std;
void dyStats(){
    TCanvas *cc = new TCanvas("cc","cc",800,600);
    gStyle->SetOptDate(0);
    gStyle->SetLineWidth(1);
    gStyle->SetLabelSize(0.06,"XY");
    gStyle->SetStatFontSize(0.08);
    gStyle->SetTitleFontSize(0.08);
    gStyle->SetStatBorderSize(0.0);
    gStyle->SetPadBorderSize(0);
    gStyle->SetPadBottomMargin(0.15);
    gStyle->SetPadTopMargin(0.07);
    gStyle->SetPadLeftMargin(0.10);
    gStyle->SetPadRightMargin(0.10);
    gStyle->SetStatStyle(0);
    gStyle->SetStatX(0.9);
    gStyle->SetStatY(0.92);
    gStyle->SetStatW(0.4);
    //gStyle->SetOptStat(11);
    gStyle->SetTitleOffset(0.8,"X");
    gStyle->SetTitleOffset(0.8,"Y");
    gStyle->SetTitleSize(0.08,"XY");
    gStyle->SetTitleAlign(23);
    gStyle->SetTitleY(1.01);
    gStyle->SetCanvasPreferGL(kTRUE);
    gPad->SetFillStyle(0);
    gPad->SetFrameFillStyle(0);
    //TGaxis::SetMaxDigits(3);
    //gStyle->SetPalette(kFruitPunch);
    //TColor::InvertPalette();
    cc->Print("dyStats.pdf[");
    TFile *f;
    vector<const char*> eventtype,cut;
    cut.push_back("");
    eventtype.push_back("DY"); 
    //eventtype.push_back("QCD"); 
    TPaveLabel *title=0; TPad *ct=0;
    TPaveStats *st=0; TH1F *hist=0;
    //TExec ex1("ex1","TGaxis::SetMaxDigits(10);");
    //ex1.Draw();
    for(vector<const char*>::iterator ccut=cut.begin();ccut!=cut.end();ccut++){
        cc->Clear();
        cc->cd();
        gPad->SetFillStyle(0);
        title = new TPaveLabel(.10,.95,.90,.99,Form(""));
        title->Draw();
        ct= new TPad("ct","ct",0.0,0.0,1.0,0.95);
        ct->Draw();
        for(vector<const char*>::iterator it=eventtype.begin();it!=eventtype.end();it++){
            cc->cd();
            gPad->SetFillStyle(0);
            title->DrawPaveLabel(.10,.95,.90,.99,Form("%s cross section tag2017 vertex_z(0,30)",*it));
            ct->cd();
            gPad->SetFillStyle(0);
            f = new TFile(Form("./dyStats.root"));
            TTree* cs = (TTree*)f->Get("dyStats");
            cs->Draw("crossSection>>cross(100,1.4e-5,1.57e-5)");
            //cs->Draw("crossSection>>cross");
            hist=(TH1F*)f->Get("cross");
            hist->SetNameTitle(Form("%s cross section",*it),Form("%s cross section",*it));
            hist->SetXTitle("cross section(mb)");
            //TGaxis* gx = (TGaxis*)hist->GetXaxis();
            //gx->SetMaxDigits(2);
            cc->Print("dyStats.pdf");
        }
    }
    cc->Print("dyStats.pdf]");
}