displaying bg histograms (Michele)

[u/mrichter/AliRoot.git] / PWG0 / eventStats / show.C

void show(const char* fileName = "event_stat.root")
{
  if (TString(fileName).BeginsWith("alien://"))
    TGrid::Connect("alien://");

  TFile::Open(fileName);
  
  c = new TCanvas;
  hist = (TH1*) gFile->Get("physics_selection/fHistStatistics");
  hist->SetStats(0);
  hist->Draw("TEXT");
  c->SetLeftMargin(0.25);
  c->SetBottomMargin(0.2);
  c->SaveAs("stat.png");
  
  c = new TCanvas;
  c->SetLeftMargin(0.25);
  hist = (TH1*) gFile->Get("physics_selection/fHistBunchCrossing");
  hist->SetStats(0);
  hist->Draw("TEXT");
  c->SaveAs("bc.png");

  // BG
  gStyle->SetStatX(0.87);
  gStyle->SetStatY(0.93);

  TH2F * hBB  = (TH2F*) gDirectory->Get("physics_selection/background_identification/h_CINT1B-ABCE-NOPF-ALL");
  TH2F * hBE  = (TH2F*) gDirectory->Get("physics_selection/background_identification/h_CINT1C-ABCE-NOPF-ALL");
  TH2F * hBEA = (TH2F*) gDirectory->Get("physics_selection/background_identification/h_CINT1A-ABCE-NOPF-ALL");

  if(!hBB || !hBE || !hBEA) {
    printf("WARNING: no BG histos\n");
  }
  else
  {
    hBE->Add(hBEA);
    hBE->SetTitle("CINT1A-ABCE-NOPF-ALL + CINT1C-ABCE-NOPF-ALL");
    c = new TCanvas("cbg","cbg",800,400);
    c->Divide(2,1);
    c->cd(1);
    hBB->Draw("colz");
    c->cd(2);
    hBE->Draw("colz");
    c->Print("bg.png");
  }
}