bdt_com.C

Go to the documentation of this file.

{
  const int nhists=7;
  TString hists[nhists]={"efrac_2sig","emaxfrac_6plane","Efrac_sumplane9","efrac_plane2","efrac_plane3","efrac_plane4","BDT7ratio"};
  TString hists_name[nhists]={"Fraction of energy in 2#sigma road","Maximum fraction of energy in 6-continuous planes","Fraction of energy in frist 10 planes","Fraction of energy in 2nd plane","Fraction of energy in 3rd plane","Fraction energy in 4th plane","BDT7 output"};

  double minXbin[nhists]={0.9,0,0,0,0,0,-0.3};
  double maxXbin[nhists]={1.01,0.2,0.2,0.025,0.025,0.025,0.3};

  int nxbins[nhists]={120,20,40,200,200,200,50};
  double lxbins[nhists]={0,0,0,0,0,0,-1};
  double hxbins[nhists]={1.2,0.2,0.2,0.2,0.2,0.2,1};

  TFile *fd = new TFile("data_results.root");
  TFile *fmc = new TFile("makeplots_LID02/data_results.root");

  gStyle->SetOptStat(0);

  for( int ih=0; ih<nhists; ++ih ){

     TH1F *hmc = (TH1F*)fmc->Get(hists[ih]);
     TH1F *hd = (TH1F*)fd->Get(hists[ih]);

     TString ratio_name="ratio_"+nhists[ih];
     TH1F *hratio=new TH1F(ratio_name,"",nxbins[ih],lxbins[ih],hxbins[ih]);
     hratio->Sumw2();
     for( int ibin=1; ibin<hratio->GetNbinsX()+1; ++ibin ){
       double bin_data=hd->GetBinContent(ibin);
       double err_data=hd->GetBinError(ibin);
       double bin_mc=hmc->GetBinContent(ibin);
       double err_mc=hmc->GetBinError(ibin);
       double bin_con=0.;
       double bin_err=0.;
       if( bin_data>0. && bin_mc>0. ){
         bin_con=bin_data/bin_mc;
         bin_err = sqrt((err_mc/bin_mc)**2+(err_data/bin_data)**2);
       }
       hratio->SetBinContent(ibin,bin_con);
       hratio->SetBinError(ibin,bin_con*bin_err);
     }


     //create canvas
     TString can_name="variable_"+hists[ih];
     TCanvas *c1 = new TCanvas(can_name,"",700,600);
     c1->cd();

     pad1 = new TPad("pad1", "",0,0,0.99,0.99);
     pad1->Draw();
     pad1->cd();
     pad1->Range(-1.081083,0.6027835,104.4476,1.087781);
     pad1->SetFillColor(0);
     pad1->SetBorderMode(0);
     pad1->SetBorderSize(2);
     pad1->SetLeftMargin(0.1523859);
     pad1->SetRightMargin(0.04214559);
     pad1->SetTopMargin(0.07790022);
     pad1->SetBottomMargin(0.3035405);
     pad1->SetFrameBorderMode(0);
     pad1->SetFrameBorderMode(0);

     //2D hist template
     TString h2f_name="h2f_"+hists[ih];
     TH2F *h2f = new TH2F(h2f_name,"Rock muons",100,minXbin[ih],maxXbin[ih],100,0.5,8000);
     //h2f->SetXTitle(hists_name[ih]);
     //h2f->SetYTitle(ytitles[ih]);
     //h2f->GetYaxis()->SetTitleOffset(2.0);
     h2f->Draw();

     hmc->SetLineColor(4);
     hmc->SetMarkerColor(4);
     hmc->SetMarkerStyle(24);
     hmc->Draw("E0 same");

     hd->SetMarkerStyle(8);
     hd->SetLineColor(1);
     hd->Draw("E0 same");

     TLegend *leg = new TLegend(0.6,0.6,0.89,0.89);
     leg->SetLineColor(0);
     leg->SetFillColor(0);
     leg->AddEntry(hd,"new","p");
     leg->AddEntry(hmc,"old","p");
     leg->Draw("same");

    //TString gifname="gif/"+can_name+".gif";
    //c1->Print(gifname);   

     pad2 = new TPad("pad2", "",0,0,0.99,0.3);
  pad2->Draw();
  pad2->cd();
  pad2->Range(-1.081081,0.6526316,103.3923,1.202863);
  pad2->SetFillColor(0);
  pad2->SetBorderMode(0);
  pad2->SetBorderSize(2);
  pad2->SetLeftMargin(0.1539252);
  pad2->SetRightMargin(0.0324703);
  pad2->SetTopMargin(0.005203551);
  pad2->SetBottomMargin(0.2678298);
  pad2->SetFrameBorderMode(0);
  pad2->SetFrameBorderMode(0);
  pad2->SetGridy();

     TString h2f2_name="h2f2_"+hists[ih];
     TH2F *h2f2 = new TH2F(h2f2_name,"",100,minXbin[ih],maxXbin[ih],100,0.01,1.99);
     h2f2->SetXTitle(hists_name[ih]);
     h2f2->SetYTitle("data/MC");
     h2f2->GetXaxis()->SetLabelSize(0.1);
     h2f2->GetXaxis()->SetTitleSize(0.13);
     h2f2->GetXaxis()->SetTitleOffset(0.85);
     h2f2->GetYaxis()->SetLabelSize(0.12);
     h2f2->GetYaxis()->SetTitleSize(0.13);
     h2f2->GetYaxis()->SetTitleOffset(0.5);
     h2f2->GetYaxis()->SetNdivisions(011);
     h2f2->Draw();

     hratio->SetMarkerStyle(24);
     hratio->SetLineColor(2);
     hratio->SetMarkerColor(2);
     hratio->Draw("E0 same");
  }

}