void setOPT_canvas(TCanvas *canvas)
{
    gStyle->SetPalette(1);
    gStyle->SetOptTitle(0);
    gStyle->SetOptStat(0);
    gStyle->SetOptFit(1111);
    gStyle->SetCanvasColor(10);
    gStyle->SetPadColor(10);
    canvas->SetRightMargin(0.01);
    canvas->SetTopMargin(0.01);
    canvas->SetBottomMargin(0.01);
    canvas->Range(-194.483,-10.3682,1041.38,-2.08469);
    canvas->SetFillColor(0);
    canvas->SetBorderMode(0);
    canvas->SetBorderSize(0);
    canvas->SetFrameLineWidth(2);
    canvas->SetFrameBorderMode(0);
    canvas->SetFrameBorderSize(0);
}
void setOPT_hists(TH2F *hist, TString xAxisTitle=" ", TString yAxisTitle=" ", Int_t color=1)
{
    hist->SetStats(1);
    hist->SetDirectory(0);
    hist->GetXaxis()->SetTitle(xAxisTitle);
    hist->GetYaxis()->SetTitle(yAxisTitle);
    hist->GetXaxis()->SetTitleSize(0.06);
    hist->GetYaxis()->SetTitleSize(0.06);
    hist->GetXaxis()->SetTitleFont(42);
    hist->GetYaxis()->SetTitleFont(42);
    hist->GetYaxis()->SetTitleOffset(1.2);
    hist->GetZaxis()->SetTitleOffset(1.2);
    hist->GetYaxis()->SetTitleOffset(1.25);
    hist->GetXaxis()->CenterTitle();
    hist->GetYaxis()->CenterTitle();
    hist->GetXaxis()->SetLabelFont(42);
    hist->GetYaxis()->SetLabelFont(42);
    hist->GetZaxis()->SetLabelFont(42);
    hist->GetXaxis()->SetLabelSize(0.05);
    hist->GetYaxis()->SetLabelSize(0.05);
    hist->GetZaxis()->SetLabelSize(0.05);
}

void setOPT_text(TLatex *text, Int_t color=1)
{   text->SetTextFont(42);
    text->SetTextSize(0.05);
    text->SetTextColor(color);
}

void plots_efficiency()
{

    TFile *input = new TFile("apr06_protons_tot_2DIM.root","read");

    TH2F *hist[6];

    TCanvas *ala[6];

    char name[256];

    for(Int_t i=0; i<6; i++)
    {
	sprintf(name,"ala_%i",i);
	ala[i]= new TCanvas(name,name,158,28,900,918);
	setOPT_canvas(ala[i]);
	ala[i]->Divide(3,3);
    }

    input->cd();

    hist[0]  = (TH2F*)input->Get("ProtonsEff_1");
    hist[1]  = (TH2F*)input->Get("ProtonsEff_2");
    hist[2]  = (TH2F*)input->Get("ProtonsEff_3");
    hist[3]  = (TH2F*)input->Get("ProtonsEff_4");
    hist[4]  = (TH2F*)input->Get("ProtonsEff_5");
    hist[5]  = (TH2F*)input->Get("ProtonsEff_6");

    TLatex *tex_ele_lost[6];

    for(Int_t i=0;i<6;i++)
    {
	setOPT_hists(hist[i]   ,"#theta [rad]","momentum [MeV/c]");
	sprintf(name,"sector %i",i);
	tex_ele_lost[i] = new TLatex(0.6,2050,name);
	setOPT_text(tex_ele_lost[i],1);
    }

    for(Int_t i=0;i<6;i++)
    {
	ala[i]->cd();
	gPad->SetTopMargin(0.1);
	gPad->SetLeftMargin(0.15);
	gPad->SetRightMargin(0.15);
	gPad->SetBottomMargin(0.15);
	hist[i]->SetMaximum(1.1);
	hist[i]->Draw("colz");
        tex_ele_lost[i]->Draw();
    }

}