pl_over_clu.C

void pl_over_clu(Int_t SmT = 0, Int_t iSm = 0, Int_t iRpc = 0) {
  //  TCanvas *can = new TCanvas("can22","can22");
  //  can->Divide(2,2);
  TCanvas* can = new TCanvas("can", "can", 50, 0, 800, 800);
  can->Divide(3, 4);

  gPad->SetFillColor(0);
  gStyle->SetPalette(1);
  gStyle->SetOptStat(kTRUE);

  gROOT->cd();
  gROOT->SetDirLevel(1);

  TH1* h;
  TH1* h1;
  TH2* h2;
  TH1* h2y;
  void fit_ybox(const char* hname);
  void fit_ybox(const char* hname, Double_t ysize);
  // if (h!=NULL) h->Delete();

  can->cd(1);
  gROOT->cd();
  TString hname1 = Form("cl_SmT%d_sm%03d_rpc%03d_Size", SmT, iSm, iRpc);
  h2             = (TH2*) gROOT->FindObjectAny(hname1);
  if (h2 != NULL) {
    h2->Draw("colz");
    h2->ProfileX()->Draw("same");
    gPad->SetLogz();
  } else {
    cout << hname1 << " not found -> return" << endl;
    return;
  }

  can->cd(2);
  gROOT->cd();

  TString hname2 = Form("cl_SmT%d_sm%03d_rpc%03d_Pos", SmT, iSm, iRpc);
  h2             = (TH2*) gROOT->FindObjectAny(hname2);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");

    can->cd(3);
    h2y = h2->ProjectionY();
    //gROOT->LoadMacro("fit_ybox.h");
    cout << " Fit with ybox " << h2y->GetName() << endl;
    fit_ybox((const char*) (h2y->GetName()));

  } else {
    cout << hname2 << " not found" << endl;
  }

  /*
can->cd(3);
 gROOT->cd();
 TString hname=Form("cl_SmT%d_sm%03d_rpc%03d_TOff",SmT,iSm,iRpc);
 h2=(TH2 *)gROOT->FindObjectAny(hname);
 if (h2!=NULL) {
  h2->Draw("colz");
  gPad->SetLogz();
  h2->ProfileX()->Draw("same");
 }else 
   {
     cout << hname << " not found" << endl;
   }
 */

  can->cd(4);
  gROOT->cd();
  TString hname4 = Form("cl_SmT%d_sm%03d_rpc%03d_Tot", SmT, iSm, iRpc);
  h2             = (TH2*) gROOT->FindObjectAny(hname4);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname4 << " not found" << endl;
  }

  can->cd(5);
  gROOT->cd();
  TString hname5 = Form("cl_SmT%d_sm%03d_rpc%03d_AvWalk", SmT, iSm, iRpc);
  h2             = (TH2*) gROOT->FindObjectAny(hname5);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname5 << " not found" << endl;
  }

  can->cd(6);
  gROOT->cd();
  TString hname6 = Form("cl_SmT%d_sm%03d_rpc%03d_AvLnWalk", SmT, iSm, iRpc);
  h2             = (TH2*) gROOT->FindObjectAny(hname6);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname6 << " not found" << endl;
  }

  can->cd(7);
  gROOT->cd();
  TString hname7 = Form("cl_SmT%d_sm%03d_rpc%03d_Mul", SmT, iSm, iRpc);
  h1             = (TH1*) gROOT->FindObjectAny(hname7);
  if (h1 != NULL) {
    h1->Draw();
    gPad->SetLogy();
  } else {
    cout << hname7 << " not found" << endl;
  }

  can->cd(8);
  gROOT->cd();
  TString hname8 = Form("cl_SmT%d_sm%03d_rpc%03d_Trms", SmT, iSm, iRpc);
  h2             = (TH2*) gROOT->FindObjectAny(hname8);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname8 << " not found" << endl;
  }

  can->cd(9);
  gROOT->cd();
  TString hname9 = Form("cl_SmT%d_sm%03d_rpc%03d_DelPos", SmT, iSm, iRpc);
  h2             = (TH2*) gROOT->FindObjectAny(hname9);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname9 << " not found" << endl;
  }

  can->cd(10);
  gROOT->cd();
  TString hname10 = Form("cl_SmT%d_sm%03d_rpc%03d_DelTOff", SmT, iSm, iRpc);
  h2              = (TH2*) gROOT->FindObjectAny(hname10);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname10 << " not found" << endl;
  }

  can->cd(11);
  gROOT->cd();
  TString hname11 = Form("cl_SmT%d_sm%03d_rpc%03d_DelMatPos", SmT, iSm, iRpc);
  h2              = (TH2*) gROOT->FindObjectAny(hname11);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname11 << " not found" << endl;
  }

  can->cd(12);
  gROOT->cd();
  TString hname12 = Form("cl_SmT%d_sm%03d_rpc%03d_DelMatTOff", SmT, iSm, iRpc);
  h2              = (TH2*) gROOT->FindObjectAny(hname12);
  if (h2 != NULL) {
    h2->Draw("colz");
    gPad->SetLogz();
    h2->ProfileX()->Draw("same");
  } else {
    cout << hname12 << " not found" << endl;
  }

  can->SaveAs(Form("pl_over_clu%01d_%01d_%01d.pdf", SmT, iSm, iRpc));
}