// prepare final histograms for human analysis
// of reconstructed info
//data from HIMster
// 12/05/2014 A.Karavdina

void DrawAna(TString path="/panda/myResults/hc_5pi_phsp_scrut2014/mom_5_61/samp0/",TString fileIN="TOT_plotRes_PHSPModel_OmegaEtaAna.root"){
  TString file_sig_name = path+fileIN;
  TFile *finSIG = new TFile(file_sig_name.Data(),"READ");
  TString fout_root = file_sig_name+"_FINALE.root";

  // TH1 *hhcTh_mc = (TH1D*)finSIG->GetObject("hhcTh_mc");
  // TH1 *hhcTh_raw =  (TH1D*)finSIG->GetObject("hhcTh_raw");
  // hhcTh_raw->Divide(hhcTh_mc);
  // TH1 *hhcPh_mc = (TH1D*)finSIG->GetObject("hhcPh_mc");
  // TH1 *hhcPh_raw =  (TH1D*)finSIG->GetObject("hhcPh_raw");
  // hhcPh_raw->Divide(hhcPh_mc);
  TH1 *hrhoTh_mc = (TH1D*)finSIG->Get("hrhoTh_mc");
  TH1 *hrhoTh_raw = (TH1D*)finSIG->Get("hrhoTh_raw");
  hrhoTh_raw->Divide(hrhoTh_mc);
 TH1 *hrhoPh_mc = (TH1D*)finSIG->Get("hrhoPh_mc");
  TH1 *hrhoPh_raw = (TH1D*)finSIG->Get("hrhoPh_raw");
  hrhoPh_raw->Divide(hrhoPh_mc);
 TH1 *hetaTh_mc = (TH1D*)finSIG->Get("hetaTh_mc");
  TH1 *hetaTh_raw = (TH1D*)finSIG->Get("hetaTh_raw");
  hetaTh_raw->Divide(hetaTh_mc);
 TH1 *hetaPh_mc = (TH1D*)finSIG->Get("hetaPh_mc");
  TH1 *hetaPh_raw = (TH1D*)finSIG->Get("hetaPh_raw");
  hetaPh_raw->Divide(hetaPh_mc);
  TH1 *hpi0Th_mc = (TH1D*)finSIG->Get("hpi0Th_mc");
  TH1 *hpi0Th_raw = (TH1D*)finSIG->Get("hpi0Th_raw");
  hpi0Th_raw->Divide(hpi0Th_mc);
 TH1 *hpi0Ph_mc = (TH1D*)finSIG->Get("hpi0Ph_mc");
  TH1 *hpi0Ph_raw = (TH1D*)finSIG->Get("hpi0Ph_raw");
  hpi0Ph_raw->Divide(hpi0Ph_mc);
  TH1 *hpiplTh_mc = (TH1D*)finSIG->Get("hpiplTh_mc");
  TH1 *hpiplTh_raw = (TH1D*)finSIG->Get("hpiplTh_raw");
  hpiplTh_raw->Divide(hpiplTh_mc);
 TH1 *hpiplPh_mc = (TH1D*)finSIG->Get("hpiplPh_mc");
  TH1 *hpiplPh_raw = (TH1D*)finSIG->Get("hpiplPh_raw");
  hpiplPh_raw->Divide(hpiplPh_mc);
 TH1 *hpimnTh_mc = (TH1D*)finSIG->Get("hpimnTh_mc");
  TH1 *hpimnTh_raw = (TH1D*)finSIG->Get("hpimnTh_raw");
  hpimnTh_raw->Divide(hpimnTh_mc);
  TH1 *hpimnPh_mc = (TH1D*)finSIG->Get("hpimnPh_mc");
  TH1 *hpimnPh_raw = (TH1D*)finSIG->Get("hpimnPh_raw");
  hpimnPh_raw->Divide(hpimnPh_mc);

  TH1 *hgammaPh_mc = (TH1D*)finSIG->Get("hgammaPh_mc");
  TH1 *hgammaPh_raw = (TH1D*)finSIG->Get("hgammaPh_raw");
  hgammaPh_raw->Divide(hgammaPh_mc);
  TH1 *hgammaTh_mc = (TH1D*)finSIG->Get("hgammaTh_mc");
  TH1 *hgammaTh_raw = (TH1D*)finSIG->Get("hgammaTh_raw");
  hgammaTh_raw->Divide(hgammaTh_mc);

  // hhcTh_raw->Write();
  // hhcPh_raw->Write();
  TFile *fout = new TFile(fout_root,"RECREATE");
  hrhoTh_raw->Write();
  hrhoPh_raw->Write();
  hetaTh_raw->Write();
  hetaPh_raw->Write();
  hpi0Th_raw->Write();
  hpi0Ph_raw->Write();
  hgammaTh_raw->Write();
  hgammaPh_raw->Write();
  hpiplTh_raw->Write();
  hpiplPh_raw->Write();
  hpimnTh_raw->Write();
  hpimnPh_raw->Write();

 //  TH1 *hhcmraw = new TH1D("hhcmraw","Signal ;m, GeV",1e3,3.4,3.7);
 //  tree_sig_raw->Project("hhcmraw","hcm");
 //  TH2 *hDalitz = new TH2D("hDalitz","; m^{2}_{#pi^{+}#pi^{-}#pi^{0}}, (GeV/c^{2})^{2};m^{2}_{#pi^{+}#pi^{-}}, (GeV/c^{2})^{2}",1e3,0.,1.,1e3,0.,1.0);
 //  tree_sig_raw->Project("hDalitz","hcd0m*hcd0m:hcd1m*hcd1m");
 //  TH1 *hhcmfit = new TH1D("hhcmfit","Signal ;m, GeV",2e3,3.3,3.7);
 //  tree_sig_4Cfit->Project("hhcmfit","hcfitm");

 //  TH1 *hrhom = new TH1D("hrhom","Signal ; m, GeV",1e3,0.,5.0);
 //  tree_sig_rho->Project("hrhom","rhom");
 //  TH2 *hPeyrou_rho = new TH2D("hPeyrou_rho","Peyrou Diagram: #rho (SIG); p_{z}, GeV/c ; p_{#perp} , GeV/c",1e3,-1,7,1e3,0,2);
 //  tree_sig_rho->Project("hPeyrou_rho","rhopt:rhopz");
 //  TH1 *hpi0m = new TH1D("hpi0m","Signal ; m, GeV",1e3,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",1e3,-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",1e3,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",1e3,5,7,1e3,0,4e-11);
 //  //  tree_sig_4Cfit->Project("hPeyrou_hc_fit","hcfitpt:hcfitpz");

 //  TH1 *hetam = new TH1D("hetam","Signal ; m, GeV",1e3,0,5.0);
 //  tree_sig_eta->Project("hetam","etam");
 //  TH2 *hPeyrou_eta = new TH2D("hPeyrou_eta","Peyrou Diagram: #eta (SIG); p_{z}, GeV/c; p_{#perp}, GeV/c",1e3,-1,7,1e3,0,2);
 //  tree_sig_eta->Project("hPeyrou_eta","etapt:etapz");
 
 //  TCanvas c1;
 //  c1.SetLogy(1);
 //  c1.Divide(3,2);
 //  c1.cd(1);
 //  hpi0m->Draw();
 //  c1.cd(2);
 //  hetam->Draw();
 //  c1.cd(3);
 //  hrhom->Draw();
 //  c1.cd(4);
 //  hhcmraw->Draw();
 //  c1.cd(5);
 //  hhcmfit->Draw();
 //  c1.cd(6);
 //  hchi2->Draw();
 //  c1.Print(fout_pdf_o);

 //  TCanvas c2;
 //  c2.Divide(2,2);
 //  c2.cd(1);
 //  //hPeyrou_pi0->Draw("colz");
 //  hPeyrou_rho->Draw("colz");
 //  c2.cd(2);
 //  hPeyrou_eta->Draw("colz");
 //  c2.cd(3);
 //  hPeyrou_hc->Draw("colz");
 //  c2.cd(4);
 //  hDalitz->Draw("colz");
 //  //  hPeyrou_hc_fit->Draw("colz");
 //  c2.Print(fout_pdf_o);

 //  TCanvas c3;
 //  c3.Divide(2,3);
 //  c3.cd(1);
 //  hpi0Th_raw->SetMinimum(0);
 //  hpi0Th_raw->SetMaximum(1.2);
 //  hpi0Th_raw->Draw();
 //  c3.cd(2);
 //  hpi0Ph_raw->SetMinimum(0);
 //  hpi0Ph_raw->SetMaximum(1.2);
 //  hpi0Ph_raw->Draw();
 //  c3.cd(3);
 //  hpiplTh_raw->SetMinimum(0);
 //  hpiplTh_raw->SetMaximum(1.2);
 //  hpiplTh_raw->Draw();
 //  c3.cd(4);
 //  hpiplPh_raw->SetMinimum(0);
 //  hpiplPh_raw->SetMaximum(1.2);
 //  hpiplPh_raw->Draw();
 //  c3.cd(5);
 //  hpimnTh_raw->SetMinimum(0);
 //  hpimnTh_raw->SetMaximum(1.2);
 //  hpimnTh_raw->Draw();
 //  c3.cd(6);
 //  hpimnPh_raw->SetMinimum(0);
 //  hpimnPh_raw->SetMaximum(1.2);
 //  hpimnPh_raw->Draw();
 //  c3.Print(fout_pdf_o);

 // TCanvas c4;
 // c4.Divide(2,2);
 // c4.cd(1);
 // hetaTh_raw->SetMinimum(0);
 // hetaTh_raw->SetMaximum(1.2);
 //  hetaTh_raw->Draw();
 //  c4.cd(2);
 //  hetaPh_raw->SetMinimum(0);
 //  hetaPh_raw->SetMaximum(1.2);
 //  hetaPh_raw->Draw();
 //  c4.cd(3);
 //  hrhoTh_raw->SetMinimum(0);
 //  hrhoTh_raw->SetMaximum(1.2);
 //  hrhoTh_raw->Draw();
 //  c4.cd(4);
 //  hrhoPh_raw->SetMinimum(0);
 //  hrhoPh_raw->SetMaximum(1.2);
 //  hrhoPh_raw->Draw();
 //  c4.Print(fout_pdf_c);


 //  // Compare reconstructed info with true one
 //  TH2 *hPeyrou_rho_tr = new TH2D("hPeyrou_rho_tr","Peyrou Diagram: #rho (SIG); p_{z}, GeV/c ; p_{#perp} , GeV/c",1e3,-1,7,1e3,0,2);
 //  tree_sig_rho->Project("hPeyrou_rho_tr","rhopt:rhopz","trrhom>0");
 //  TH2 *hPeyrou_eta_tr = new TH2D("hPeyrou_eta_tr","Peyrou Diagram: #eta (SIG); p_{z}, GeV/c; p_{#perp}, GeV/c",1e3,-1,7,1e3,0,2);
 //  tree_sig_eta->Project("hPeyrou_eta_tr","etapt:etapz","tretam>0");
 //  // TCanvas c4;
 //  // c4.Divide(2,2);
 //  // c4.cd(1);
 //  // //hPeyrou_pi0->Draw("colz");
 //  // hPeyrou_rho_tr->Draw("colz");
 //  // c4.cd(2);
 //  // hPeyrou_eta_tr->Draw("colz");
 //  // c4.Print(fout_pdf_c);


 //  //  TString out="plotRes.root";
 //  TFile *fout = new TFile(fout_root,"RECREATE");
 //   hpi0m->Write();
 //  hetam->Write();
 //  hrhom->Write();
 //  hhcmraw->Write();
 //  hhcmfit->Write();
 //  hPeyrou_pi0->Write();
 //  hPeyrou_rho->Write();
 //  hPeyrou_eta->Write();
 //  hPeyrou_hc->Write();
 //  //  hPeyrou_hc_fit->Write();
 //  hDalitz->Write();
 //  hhcTh_mc->Write();
 //  hhcTh_raw->Write();
 //  hhcPh_mc->Write();
 //  hhcPh_raw->Write();
 //  hrhoTh_mc->Write();
 //  hrhoPh_mc->Write();
 //  hrhoTh_raw->Write();
 //  hrhoPh_raw->Write();
 //  hetaTh_mc->Write();
 //  hetaPh_mc->Write();
 //  hetaTh_raw->Write();
 //  hetaPh_raw->Write();
 //  hpi0Th_mc->Write();
 //  hpi0Ph_mc->Write();
 //  hpi0Th_raw->Write();
 //  hpi0Ph_raw->Write();

 //  // hrhoTh_fit->Write();
 //  // hrhoPh_fit->Write(); 
 //  // hetaTh_fit->Write();
 //  // hetaPh_fit->Write();
 //  // hpi0Th_fit->Write();
 //  // hpi0Ph_fit->Write();
 //  hPeyrou_rho_tr->Write();
 //  hPeyrou_eta_tr->Write();
 //  hpiplPh_raw->Write();
 //  hpiplPh_mc->Write();
 //  hpimnPh_raw->Write();
 //  hpimnPh_mc->Write();
 //  hpiplTh_raw->Write();
 //  hpiplTh_mc->Write();
 //  hpimnTh_raw->Write();
 //  hpimnTh_mc->Write();
 //  c1.Write();
 //  c2.Write();
 //  c3.Write();
 //  c4.Write();
  fout->Close();
}