combinePlotsEtac.C

Go to the documentation of this file.

// --------------------------------------------------------------------

void confgraph(TGraph *g, TString tit, int col=1, int marker=20, double shift=0.)
{
        if (tit!="") g->GetHistogram()->SetTitle(tit);
        g->GetHistogram()->SetMinimum(0);
        g->SetLineColor(col);
        g->SetMarkerColor(col);
        g->SetMarkerStyle(marker);
        g->SetMarkerSize(1.5);
        g->SetLineWidth(2);
        
        TF1 *f1 = g->GetFunction("f1");
        
        if (f1)
        {
                f1->SetLineColor(col);
                f1->SetLineStyle(2);
        }
        
        for (int i=0;i<g->GetN();++i)
        {
                double x,y;
                g->GetPoint(i,x,y);
                g->SetPoint(i,x+shift,y);
                
        //      cout <<x<<" "<<y<<endl;
        }
}

// --------------------------------------------------------------------


TH1F* createHistoGraph(TGraph *g, TString tit="", double yymax=0, double xmin=0, double xmax=0)
{
        static int cnt=0;
        if (tit=="") tit=g->GetTitle();
        
        double ymax  = TMath::MaxElement(g->GetN(), g->GetY());
        double dymax = g->GetErrorY(TMath::LocMax(g->GetN(), g->GetEY()));
        
        if (xmin>=xmax)
        {
                xmin = TMath::MinElement(g->GetN(), g->GetX());
                xmax = TMath::MaxElement(g->GetN(), g->GetX());
        }       
        TH1F *h=new TH1F(Form("h%03d",cnt++),tit,g->GetN(),xmin-0.5*(xmax-xmin)/g->GetN(), xmax+0.5*(xmax-xmin)/g->GetN());
        
        h->GetXaxis()->SetLabelSize(0.05);
        h->GetXaxis()->SetTitleSize(0.05);
        h->GetYaxis()->SetLabelSize(0.05);
        h->GetYaxis()->SetTitleSize(0.05);
        h->GetYaxis()->SetTitleOffset(1.4);
        
        if (yymax==0) 
                h->SetMaximum((ymax+dymax)*1.05);
        else 
                h->SetMaximum(yymax);
        h->SetMinimum(0);
        
        h->SetStats(0);
        
        return h;
}

// --------------------------------------------------------------------

void combinePlotsEtac(TString fname="ana_etac.root")
{
        gStyle->SetOptStat(0);
        gStyle->SetOptFit(0);
        
        TFile *f=new TFile(fname);
        
        TString plotnam[2]={"etac1_A","etac1_B"};
        TString lab[4] = {  "signal efficiency;EMC missing (polar #theta) [%];signal efficiency [%]", "background efficiency;EMC missing (polar #theta) [%];background efficiency [%]", "signal to noise;EMC missing (polar #theta) [%];S/B", "significance;EMC missing (polar #theta) [%];significance [#sigma]"};
        
        TGraphErrors *geffs[2];
        TGraphErrors *geffb[2];
        TGraphErrors *gston[2];
        TGraphErrors *gsign[2];

        TCanvas *c1=new TCanvas("c1","c1",10,10,1300,1000);
        c1->Divide(2,2,0.0001,0.0001);

        int colors[2] = {kBlue,kCyan-2};
        //int colors[2] = {2,kRed-7};

        TH1F *h[4]={0};
        double hmaxy[4] = {20., 1e-5, 0.25, 7.};

        for (int i=0;i<2;++i)
        {
                c1->cd(1);
                geffs[i]=(TGraphErrors*)f->Get(Form("effs_%s",plotnam[i].Data()));
                confgraph(geffs[i],"",colors[i],20+(i%2)*3,i);
                if (h[0]==0) h[0]=createHistoGraph(geffs[i],lab[0], hmaxy[0]);
                if (i==0) h[0]->Draw();
                geffs[i]->Draw("P same");
                
                c1->cd(2);      gPad->SetLogy();
                geffb[i]=(TGraphErrors*)f->Get(Form("effb_%s",plotnam[i].Data()));
                confgraph(geffb[i],"",colors[i],20+(i%2)*3,i);
                if (h[1]==0) h[1]=createHistoGraph(geffb[i],lab[1], hmaxy[1]);
                h[1]->SetMinimum(1e-8);
                if (i==0) h[1]->Draw();
                geffb[i]->Draw("P same");
                
                c1->cd(3);      
                gston[i]=(TGraphErrors*)f->Get(Form("sn_%s",plotnam[i].Data()));
                confgraph(gston[i],"",colors[i],20+(i%2)*3,i);
                if (h[2]==0) h[2]=createHistoGraph(gston[i],lab[2], hmaxy[2]);
                if (i==0) h[2]->Draw();
                gston[i]->Draw("P same");
                
                c1->cd(4);      
                gsign[i]=(TGraphErrors*)f->Get(Form("sign_%s",plotnam[i].Data()));      
                confgraph(gsign[i],"",colors[i],20+(i%2)*3,i);
                if (h[3]==0) h[3]=createHistoGraph(gsign[i],lab[3], hmaxy[3]);
                if (i==0) h[3]->Draw();
                gsign[i]->Draw("P same");
        }
        
        c1->cd(1);
        TLegend *leg1=new TLegend(0.6,0.8,0.94,0.92);
        leg1->AddEntry(geffs[0],"#eta_{c1}#eta - Setup A","lep");
        leg1->AddEntry(geffs[1],"#eta_{c1}#eta - Setup B","lep");
        leg1->Draw();

        c1->cd(2);
        TLegend *leg2=new TLegend(0.6,0.8,0.94,0.92);
        leg2->AddEntry(geffb[0],"#eta_{c1}#eta - Setup A","lep");
        leg2->AddEntry(geffb[1],"#eta_{c1}#eta - Setup B","lep");
        leg2->Draw();
        
        c1->cd(3);
        TLegend *leg3=new TLegend(0.6,0.8,0.94,0.92);
        leg3->AddEntry(gston[0],"#eta_{c1}#eta - Setup A","lep");
        leg3->AddEntry(gston[1],"#eta_{c1}#eta - Setup B","lep");
        leg3->Draw();
        
        c1->cd(4);
        TLegend *leg4=new TLegend(0.6,0.8,0.94,0.92);
        leg4->AddEntry(gsign[0],"#eta_{c1}#eta - Setup A","lep");
        leg4->AddEntry(gsign[1],"#eta_{c1}#eta - Setup B","lep");
        leg4->Draw();
        
        c1->SaveAs("fig/comb_etac.gif");
}