void appearance_nu_mode_0() {
  
   TCanvas *c1 = new TCanvas("c1","",200,10,700,500);

//nu_e appearance
TH1F *h1 = new TH1F("h1", "", 26,0,1.30);
h1->Fill(0.125,1.96212);
h1->Fill(0.175,3.96698);
h1->Fill(0.225,7.96632);
h1->Fill(0.275,15.8017);
h1->Fill(0.325,29.728);
h1->Fill(0.375,51.2277);
h1->Fill(0.425,79.6251);
h1->Fill(0.475,111.388);
h1->Fill(0.525,140.652);
h1->Fill(0.575,160.968);
h1->Fill(0.625,167.718);
h1->Fill(0.675,160.003);
h1->Fill(0.725,140.845);
h1->Fill(0.775,115.542);
h1->Fill(0.825,89.3378);
h1->Fill(0.875,65.8507);
h1->Fill(0.925,46.7573);
h1->Fill(0.975,32.2991);
h1->Fill(1.025,21.9334);
h1->Fill(1.075,14.7951);
h1->Fill(1.125,9.97224);
h1->Fill(1.175,6.68135);
h1->Fill(1.225,4.36624);
h1->SetLineColor(2);
h1->SetLineWidth(2);
h1->SetTitle("Appearance #nu mode (#delta_{CP} =0), T2HK");
h1->GetXaxis()->SetTitle("Reconstructed Energy E^{rec}_{#nu} (GeV)");
h1->GetYaxis()->SetTitle("Number of events/50MeV");
h1->GetYaxis()->SetRangeUser(0,250);
h1->GetXaxis()->SetRangeUser(0,1.30);
h1->Draw("");

//nu_e bar appearance
TH1F *h2 = new TH1F("h2", "", 26,0,1.30);
h2->Fill(0.125,0.0636863);
h2->Fill(0.175,0.086231);
h2->Fill(0.225,0.105766);
h2->Fill(0.275,0.129829);
h2->Fill(0.325,0.173571);
h2->Fill(0.375,0.249346);
h2->Fill(0.425,0.356767);
h2->Fill(0.475,0.481846);
h2->Fill(0.525,0.604162);
h2->Fill(0.575,0.705091);
h2->Fill(0.625,0.772653);
h2->Fill(0.675,0.802991);
h2->Fill(0.725,0.79962);
h2->Fill(0.775,0.770952);
h2->Fill(0.825,0.726768);
h2->Fill(0.875,0.675123);
h2->Fill(0.925,0.620853);
h2->Fill(0.975,0.565642);
h2->Fill(1.025,0.508763);
h2->Fill(1.075,0.44797);
h2->Fill(1.125,0.380878);
h2->Fill(1.175,0.307055);
h2->Fill(1.225,0.229948);
h2->SetLineColor(3);
h2->SetLineWidth(2);
h2->Draw("same");

//nu_mu+numu_bar beam
TH1F *h3 = new TH1F("h3", "", 26,0,1.30);
h3->Fill(0.125,0.00117504+0);
h3->Fill(0.175,0.00352266+0);
h3->Fill(0.225,0.00851674+0);
h3->Fill(0.275,0.0168964+0);
h3->Fill(0.325,0.0278187+0);
h3->Fill(0.375,0.0382866+0);
h3->Fill(0.425,0.0451024+0);
h3->Fill(0.475,0.0501771+0);
h3->Fill(0.525,0.0651834+0);
h3->Fill(0.575,0.108156+0);
h3->Fill(0.625,0.190819+0);
h3->Fill(0.675,0.306894+0);
h3->Fill(0.725,0.432797+0);
h3->Fill(0.775,0.539502+0);
h3->Fill(0.825,0.604905+0);
h3->Fill(0.875,0.620392+0);
h3->Fill(0.925,0.592213+0);
h3->Fill(0.975,0.537992+0);
h3->Fill(1.025,0.47784+0);
h3->Fill(1.075,0.423396+0);
h3->Fill(1.125,0.371542+0);
h3->Fill(1.175,0.308949+0);
h3->Fill(1.225,0.227929+0);
h3->SetLineColor(4);
h3->SetLineWidth(2);
h3->Draw("same");


//nue+nue_bar CC beam
TH1F *h4 = new TH1F("h4", "",26,0,1.30);
h4->Fill(0.125,0.551575+0.0225402);
h4->Fill(0.175,1.20627+0.0446868);
h4->Fill(0.225,2.24115+0.076913);
h4->Fill(0.275,3.61859+0.117816);
h4->Fill(0.325,5.21137+0.164572);
h4->Fill(0.375,6.87829+0.2145);
h4->Fill(0.425,8.52428+0.265999);
h4->Fill(0.475,10.1+0.318316);
h4->Fill(0.525,11.5688+0.370827);
h4->Fill(0.575,12.8847+0.422745);
h4->Fill(0.625,13.9963+0.473376);
h4->Fill(0.675,14.865+0.522322);
h4->Fill(0.725,15.4785+0.569425);
h4->Fill(0.775,15.8472+0.614848);
h4->Fill(0.825,15.991+0.659302);
h4->Fill(0.875,15.9331+0.703774);
h4->Fill(0.925,15.7009+0.748691);
h4->Fill(0.975,15.3245+0.793128);
h4->Fill(1.025,14.8163+0.833667);
h4->Fill(1.075,14.1059+0.860069);
h4->Fill(1.125,12.9536+0.847983);
h4->Fill(1.175,11.0103+0.762916);
h4->Fill(1.225,8.19187+0.591599);
h4->SetLineColor(5);
h4->SetLineWidth(2);
h4->Draw("same");

//NC bckg
TH1F *h5 = new TH1F("h5", "",26,0,1.30);
h5->Fill(0.125,0.371822);
h5->Fill(0.175,0.740551);
h5->Fill(0.225,1.32521);
h5->Fill(0.275,2.1789);
h5->Fill(0.325,3.34877);
h5->Fill(0.375,4.87485);
h5->Fill(0.425,6.7747);
h5->Fill(0.475,8.9818);
h5->Fill(0.525,11.2419);
h5->Fill(0.575,13.0772);
h5->Fill(0.625,13.9588);
h5->Fill(0.675,13.6257);
h5->Fill(0.725,12.258);
h5->Fill(0.775,10.3255);
h5->Fill(0.825,8.28329);
h5->Fill(0.875,6.40341);
h5->Fill(0.925,4.80385);
h5->Fill(0.975,3.53141);
h5->Fill(1.025,2.59161);
h5->Fill(1.075,1.93953);
h5->Fill(1.125,1.48231);
h5->Fill(1.175,1.11345);
h5->Fill(1.225,0.767099);
h5->SetLineColor(6);
h5->SetLineWidth(2);
h5->Draw("same");


//total
TH1F *h6 = new TH1F("h6", "",26,0,1.30);
h6->Fill(0.125,2.63907724);
h6->Fill(0.175,5.16912746);
h6->Fill(0.225,9.54414474);
h6->Fill(0.275,17.8844084);
h6->Fill(0.325,33.6963787);
h6->Fill(0.375,60.1854166);
h6->Fill(0.425,97.8405554);
h6->Fill(0.475,142.9426691);
h6->Fill(0.525,187.4372044);
h6->Fill(0.575,220.553345);
h6->Fill(0.625,233.020515);
h6->Fill(0.675,222.235296);
h6->Fill(0.725,193.857062);
h6->Fill(0.775,157.84289);
h6->Fill(0.825,122.719621);
h6->Fill(0.875,92.888811);
h6->Fill(0.925,69.489378);
h6->Fill(0.975,52.177423);
h6->Fill(1.025,40.041551);
h6->Fill(1.075,31.760773);
h6->Fill(1.125,25.669744);
h6->Fill(1.175,20.128826);
h6->Fill(1.225,14.273191);
h6->SetLineColor(1);
h6->SetLineWidth(2);
h6->Draw("same");

gStyle->SetOptStat(0);
TLegend*leg = new TLegend(0.7,0.6,0.9,0.85);
leg->SetTextSize(0.02);
leg->SetTextColor(1);
leg->SetFillColor(0);
leg->SetBorderSize(0);
leg ->SetTextSize(12);
leg ->SetTextFont(43);
leg->AddEntry(h6,"Total","l");
leg->AddEntry(h1,"SIG #nu_{#mu} #rightarrow #nu_{e}","l");
leg->AddEntry(h2,"SIG Wrong-Sign #bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e}","l");
leg->AddEntry(h3,"BG #nu_{#mu}+#bar{#nu}_{#mu}","l");
leg->AddEntry(h4,"BG #nu_{e} + #bar{#nu}_{e}","l");
leg->AddEntry(h5,"BG NC","l");
leg->Draw();

   // TCanvas::Update() draws the frame, after which one can change it
   c1->Update();
   c1->GetFrame()->SetFillColor(0);
   c1->GetFrame()->SetBorderSize(12);
   c1->Modified();
    gPad->Print("appearance_nu_mode_0.pdf");//pdf format
 

}