//draw results of fast simulation analysis
//in nice and complete form
//16/04/2014 A.Karavdina

void DrawAnaResults(){
 
  TString path = "/panda/scrut14/macro/scrut/fastsim_hc/";
  TString file_sig = path+"/hc_rho0pi0_fast.root_taskana.root";
  TString file_bkg = path+"/dpm_5_6119_fast.root_taskana.root";

  TFile *finSIG = new TFile(file_sig.Data(),"READ");
  TTree *tree_sig_raw = (TTree*)finSIG->Get("ntp3");
  TTree *tree_sig_4Cfit = (TTree*)finSIG->Get("ntp4");
  TTree *tree_sig_rho = (TTree*)finSIG->Get("ntp2");
  TTree *tree_sig_pi0 = (TTree*)finSIG->Get("ntp1");

  TH1 *hhcmraw = new TH1D("hhcmraw","Signal ;m, GeV",2e2,3.4,3.7);
  tree_sig_raw->Project("hhcmraw","hcm");
  TH1 *hhcmfit = new TH1D("hhcmfit","Signal ;m, GeV",2e2,3.5255,3.527);
  tree_sig_4Cfit->Project("hhcmfit","hcfitm");

  TH1 *hrhom = new TH1D("hrhom","Signal ; m, GeV",1e2,0,2.0);
  tree_sig_rho->Project("hrhom","pichpairm");
TH2 *hPeyrou_rho = new TH2D("hPeyrou_rho","Peyrou Diagram: #rho (SIG); p_{z}, GeV/c ; p_{#perp} , GeV/c",1e2,-1,7,1e3,0,2);
  tree_sig_rho->Project("hPeyrou_rho","pichpairpt:pichpairpz");
  TH1 *hpi0m = new TH1D("hpi0m","Signal ; m, GeV",1e2,0,0.3);
  tree_sig_pi0->Project("hpi0m","pi0m");
TH2 *hPeyrou_pi0 = new TH2D("hPeyrou_pi0","Peyrou Diagram: #pi^{0} (SIG); p_{z}, GeV/c; p_{#perp} , GeV/c",1e2,-1,7,1e3,0,2);
  tree_sig_pi0->Project("hPeyrou_pi0","pi0pt:pi0pz");
  TH2 *hPeyrou_hc = new TH2D("hPeyrou_hc","Peyrou Diagram: h_{c} (SIG); p_{z}, GeV/c; p_{#perp} , GeV/c",1e2,5,7,1e3,0,0.2);
  tree_sig_raw->Project("hPeyrou_hc","hcpt:hcpz");
  TH2 *hPeyrou_hc_fit = new TH2D("hPeyrou_hc_fit","Peyrou Diagram: h_{c} (SIG, 4C fit); p_{z}, GeV/c; p_{#perp} , GeV/c",1e2,5,7,1e3,0,4e-11);
  tree_sig_4Cfit->Project("hPeyrou_hc_fit","hcfitpt:hcfitpz");

  TFile *finBKG = new TFile(file_bkg.Data(),"READ");
  TTree *tree_bkg_raw = (TTree*)finBKG->Get("ntp3");
  TTree *tree_bkg_4Cfit = (TTree*)finBKG->Get("ntp4");
  TTree *tree_bkg_rho = (TTree*)finBKG->Get("ntp2");
  TTree *tree_bkg_pi0 = (TTree*)finBKG->Get("ntp1");
  TH1 *hhcmraw_bkg = new TH1D("hhcmraw_bkg","Background ;m, GeV",2e2,3.4,3.7);
  tree_bkg_raw->Project("hhcmraw_bkg","hcm");
  TH1 *hhcmfit_bkg = new TH1D("hhcmfit_bkg","Background ;m, GeV",2e2,3.5255,3.527);
  tree_bkg_4Cfit->Project("hhcmfit_bkg","hcfitm");
  TH1 *hrhom_bkg = new TH1D("hrhom_bkg","Background ; m, GeV",1e2,0,2.0);
  tree_bkg_rho->Project("hrhom_bkg","pichpairm");
  TH1 *hpi0m_bkg = new TH1D("hpi0m_bkg","Background ; m, GeV",1e2,0,0.3);
  tree_bkg_pi0->Project("hpi0m_bkg","pi0m");
  TH2 *hPeyrou_pi0_bkg = new TH2D("hPeyrou_pi0_bkg","Peyrou Diagram: #pi^{0} (BKG);  p_{z}, GeV/c; p_{#perp} , GeV/c",1e2,-1,7,1e3,0,2);
  tree_bkg_pi0->Project("hPeyrou_pi0_bkg","pi0pt:pi0pz");
TH2 *hPeyrou_rho_bkg = new TH2D("hPeyrou_rho_bkg","Peyrou Diagram: #rho (BKG);  p_{z}, GeV/c; p_{#perp} , GeV/c ",1e2,-1,7,1e3,0,2);
  tree_bkg_rho->Project("hPeyrou_rho_bkg","pichpairpt:pichpairpz");
TH2 *hPeyrou_hc_bkg = new TH2D("hPeyrou_hc_bkg","Peyrou Diagram: h_{c} (SIG); p_{z}, GeV/c ; p_{#perp} , GeV/c",1e2,5,7,1e3,0,0.2);
  tree_bkg_raw->Project("hPeyrou_hc_bkg","hcpt:hcpz");
  TH2 *hPeyrou_hc_fit_bkg = new TH2D("hPeyrou_hc_fit_bkg","Peyrou Diagram: h_{c} (SIG, 4C fit); p_{z}, GeV/c; p_{#perp} , GeV/c",1e2,5,7,1e3,0,4e-11);
  tree_bkg_4Cfit->Project("hPeyrou_hc_fit_bkg","hcfitpt:hcfitpz");
  TCanvas c1;
  c1.SetLogy(1);
  c1.Divide(4,2);
  c1.cd(1);
  hpi0m->Draw();
  c1.cd(2);
  hrhom->Draw();
  c1.cd(3);
  hhcmraw->Draw();
  c1.cd(4);
  hhcmfit->Draw();
  c1.cd(5);
  hpi0m_bkg->Draw();
  c1.cd(6);
  hrhom_bkg->Draw();
  c1.cd(7);
  hhcmraw_bkg->Draw();
  c1.cd(8);
  hhcmfit_bkg->Draw();

  c1.Print("plotRes.pdf(");
  TCanvas c2;
  c2.Divide(4,2);
  c2.cd(1);
  hPeyrou_pi0->Draw("colz");
  c2.cd(2);
  hPeyrou_rho->Draw("colz");
  c2.cd(3);
  hPeyrou_hc->Draw("colz");
  c2.cd(4);
  hPeyrou_hc_fit->Draw("colz");
  c2.cd(5);
  hPeyrou_pi0_bkg->Draw("colz");
  c2.cd(6);
  hPeyrou_rho_bkg->Draw("colz");
  c2.cd(7);
  hPeyrou_hc_bkg->Draw("colz");
  c2.cd(8);
  hPeyrou_hc_fit_bkg->Draw("colz");
  c2.Print("plotRes.pdf)");

  TString out="plotRes.root";
  TFile *fout = new TFile(out,"RECREATE");
  hpi0m->Write();
  hrhom->Write();
  hhcmraw->Write();
  hhcmfit->Write();
  hpi0m_bkg->Write();
  hrhom_bkg->Write();
  hhcmraw_bkg->Write();
  hhcmfit_bkg->Write();
  hPeyrou_pi0->Write();
  hPeyrou_rho->Write();
  hPeyrou_hc->Write();
  hPeyrou_hc_fit->Write();
  hPeyrou_pi0_bkg->Write();
  hPeyrou_rho_bkg->Write();
  hPeyrou_hc_bkg->Write();
  hPeyrou_hc_fit_bkg->Write();
  c1.Write();
  c2.Write();
  fout->Close();
}