Xsec_final.C

Go to the documentation of this file.

{
 
  //Gargamelle Experiment
  double x[6] = {1.46513,2.4752,3.44589,4.44633,5.42126,7.93998};
  double y[6] = {13.53,17.8427,24.5361,33.1036,37.9249,36.5359};
  double exl[6] = {0.917728,0.98571,0.91773,0.95133,1.46119,2.10699};
  double exh[6] = {1.08768,0.98571,0.98571,0.91773,1.52993,3.94362};
  double eyl[6] = {2.2105,3.9125,4.5932,5.9587,8.5096,6.9779};
  double eyh[6] = {3.576,3.4038,4.5971,5.4423,8.3394,6.4654};

  TGraphAsymmErrors* gargamelle = new TGraphAsymmErrors(6,x,y,exl,exh,eyl,eyh);
  gargamelle->SetTitle("Gargamelle #nu_{e} CC data");
  gargamelle->SetMarkerStyle(21);
  gargamelle->SetMarkerColor(kViolet-1);
  gargamelle->SetLineColor(kViolet-1);
  gargamelle->SetLineWidth(2);
  //gargamelle->Draw("AP");
  //gargamelle->GetXaxis()->SetTitle("E_{#nu} (GeV)");
  //gargamelle->GetYaxis()->SetTitle("#nu_{e} CC #it{#sigma} (#times10^{-39} cm^{2}/nucleon)");
  //gargamelle->Draw("AP");

  //T2K results
  double x_t2k[1]={1.343};
  double y_t2k[1]={11.1};
  double exl_t2k[1]={0.918};
  double exh_t2k[1]={1.782};
  double eyl_t2k[1]={2.1};
  double eyh_t2k[1]={2.1};

  TGraphAsymmErrors *t2k=new TGraphAsymmErrors(1,x_t2k,y_t2k,exl_t2k,exh_t2k,eyl_t2k,eyh_t2k);
  t2k->SetMarkerStyle(24);
  t2k->SetLineColor(4);
  t2k->SetMarkerColor(4);
  t2k->SetLineWidth(2);

  //NOvA results
  double x_nova[4]={1.25,1.75,2.25,2.75};
  double y_nova[4]={22.0921,28.8235,28.16,37.6431};
  double ex_nova[4]={0.25,0.25,0.25,0.25};
  double ey_nova[4]={6.58498,5.99768,5.53289,6.66383};
 
  TGraphErrors *nova=new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova); 
  //nova->SetMarkerStyle(8);
  //nova->SetMarkerSize(1.2);
  nova->SetLineColor(2);
  //nova->SetMarkerColor(2);
  nova->SetLineWidth(3);

  double ey_nova_stat[4]={4.88965,3.12006,2.64657,2.54683};
  TGraphErrors *nova_stat=new TGraphErrors(4,x_nova,y_nova,ex_nova,ey_nova_stat);
  nova_stat->SetMarkerStyle(8);
  nova_stat->SetMarkerSize(1.2);
  nova_stat->SetLineWidth(3);

  //GENIE prediction
  TFile *fgenie=new TFile("xsec_genie.root");
  TGraph *gr_genie=(TGraph*)fgenie->Get("nu_e_C12/tot_cc");
  TSpline3* spline_nue=new TSpline3("mysplinenue",gr_genie);
  TH1F *Xnue=new TH1F("Xnue","",100,0,10);
  double x_genie[100]={0.};
  double y_genie[100]={0.};

  for( int ibin=1; ibin<Xnue->GetNbinsX()+1; ++ibin ){
    double binctr=Xnue->GetBinCenter(ibin);
    double xnue=(spline_nue->Eval(binctr)/12.)*10.0;
    Xnue->SetBinContent(ibin,xnue);

    x_genie[ibin-1]=binctr;
    y_genie[ibin-1]=xnue;
    
   }

  //get genie error band
  TFile *fband=new TFile("genie_band.root");
  TH1F *hplus=(TH1F*)fband->Get("ccnue_1p_total");
  TH1F *hminus=(TH1F*)fband->Get("ccnue_1m_total");
  double x_band[49]={0.};
  double y_band[49]={0.};
  double exl_band[49]={0.};
  double exh_band[49]={0.};
  double eyl_band[49]={0.};
  double eyh_band[49]={0.};
  TH1F *Xnue2=new TH1F("Xnue2","",50,0,10);
  for( int ibin=1; ibin<Xnue2->GetNbinsX()+1; ++ibin ){
    double binctr=Xnue2->GetBinCenter(ibin);
    //if( binctr>8. ) continue;
    double xnue2=(spline_nue->Eval(binctr)/12.)*10.0;
    Xnue2->SetBinContent(ibin,xnue2);

   if( ibin>1 ){
      x_band[ibin-2]=binctr;
      exl_band[ibin-2]=0.1;
      exh_band[ibin-2]=0.1;
      y_band[ibin-2]=xnue2;
      eyl_band[ibin-2]=xnue2*(hminus->GetBinContent(ibin));
      eyh_band[ibin-2]=xnue2*(hplus->GetBinContent(ibin));
    }
  }//loop bins

  TGraphAsymmErrors *band=new TGraphAsymmErrors(49,x_band,y_band,exl_band,exh_band,eyl_band,eyh_band);
  band->SetFillColor(3);
  band->SetLineColor(3);

 
/*
  TBox *errband[50]; 
  for( int ibin=1; ibin<Xnue2->GetNbinsX()+1; ++ibin ){
    float xlow=(ibin-1)*0.2;
    float xhigh=ibin*0.2;
    float ylow=Xnue2->GetBinContent(ibin)*(1.0-hminus->GetBinContent(ibin));
    float yhigh=Xnue2->GetBinContent(ibin)*(1.0+hplus->GetBinContent(ibin));
    errband[ibin-1]=new TBox(xlow,ylow,xhigh,yhigh);
  }
*/
  
  TGraph *genie=new TGraph(100,x_genie,y_genie);
  genie->SetLineStyle(2);
  genie->SetLineWidth(3);
  genie->SetLineColor(1);
  genie->SetFillColor(3);

  gStyle->SetOptStat(0);
  //TH2F template
  TCanvas *c1=new TCanvas("Xsec_Enue_withNOvA","",700,600);
  TH2F *h2f=new TH2F("h2f","",10,0,10,10,0,59);
  //h2f->SetXTitle("E_{#nu} (GeV)");
  h2f->SetXTitle("Electron Neutrino Energy (GeV)");
  h2f->SetYTitle("#nu_{e} CC #it{#sigma} (#times10^{-39} cm^{2}/nucleon)");
  h2f->GetXaxis()->CenterTitle();
  h2f->GetYaxis()->CenterTitle();
  h2f->Draw();

  TLatex* prelim = new TLatex(.9, .95, "NO#nuA Preliminary");
     prelim->SetTextColor(kBlue+1);
     prelim->SetNDC();
     prelim->SetTextSize(2/30.);
     prelim->SetTextAlign(32);
     //prelim->Draw();

/*
  for( int i=1; i<40; ++i ){
    errband[i]->SetFillColor(3);
    errband[i]->SetLineColor(3);
    errband[i]->Draw("c same");
  }
*/

  band->Draw("E3 same");

  gargamelle->Draw("P same");
  t2k->Draw("P same");
  nova->Draw("P z same");
  
  gStyle->SetEndErrorSize(4);

  nova_stat->Draw("P same");

  genie->Draw("c same");
  //Xnue->Draw("HISt same");
  //Xnue2->SetLineColor(2);
  //Xnue2->Draw("HIST same");

  TLegend *leg=new TLegend(0.6,0.6,0.89,0.89);
  leg->SetLineColor(0);
  leg->SetFillColor(0);
  leg->AddEntry(genie,"GENIE 2.8.0 #pm 1#sigma on C12","fl");
  leg->AddEntry(gargamelle,"Gargamelle #nu_{e} data","pl");
  leg->AddEntry(t2k,"T2K #nu_{e} data","pl");
  leg->AddEntry(nova_stat,"NOvA #nu_{e} data","pl");
  leg->AddEntry(nova_stat,"stat. uncert.","E");
  leg->AddEntry(nova,"stat. + syst. uncert.","E");
  leg->Draw("same");

/*
   TLatex* prelim = new TLatex(.9, .95, "NO#nuA Preliminary");
   prelim->SetTextColor(kBlue+1);
   prelim->SetNDC();
   prelim->SetTextSize(2/30.);
   prelim->SetTextAlign(32);
   prelim->Draw();
*/

}