#include<iostream>
#include<iomanip>
#include<TFile.h>
#include<TTree.h>
#include<TChain.h>
#include<TCanvas.h>
#include<TLegend.h>
#include "TString.h"
#include "TROOT.h"

using namespace std;

void E_had() 
{
   const int nbins = 40;
   const int min   = 0;
   const int max   = 3;
    
   TFile* file_C12 = new TFile("default_ccqe_Carbon.root","OPEN");
   if (!file_C12 || file_C12->IsZombie() ) { delete file_C12; return; }
   
   TTree* t_C12    = (TTree*) file_C12->Get("gst");
   if (! t_C12) { delete file_C12; return;}
   
  
   auto c = new TCanvas();
   c->cd();
   
   TH1D* h_C12_default_KE_nucleon = new TH1D("h_C12_default_KE_nucleon "," ", nbins, min,max);
   TH1D* h_C12_default_pion_plus = new TH1D ("h_C12_default_pion_plus "," ", nbins, min,max);
   TH1D* h_C12_default_pion_0 = new TH1D ("h_C12_default_pion_0 "," ", nbins, min,max);
    
   h_C12_default_KE_nucleon->Sumw2(kTRUE);
 //  t_C12->Draw("(Ef-0.938)>>h_C12_default_KE_nucleon","cc&&pdgf==2212"," goff");
   t_C12->Draw("(Ef-0.13957)>>h_C12_default_pion_plus","cc&&pdgf==211","goff");
   t_C12->Draw("(Ef-0.13498)>>h_C12_default_pion_0","cc&&pdgf==111","goff");

    TString cut_nucleon; 
    cut_nucleon = "cc&&pdgf==2212";
    
    t_C12->Project("h_C12_default_KE_nucleon", "(Ef-0.938)", cut_nucleon);
  
    h_C12_default_KE_nucleon->Draw("hist");
 
   c->Update();
  
}