// prepare final histograms for human analysis
// of reconstructed info
// rho-eta and omega-eta
//local data
// 19/05/2014 A.Karavdina

void DrawEff(TString path="/panda/scrut14/macro/scrut/results_2pi+2pi-pi0_ANG/FwdSpec_BarrelSpec/",TString fileIN1="plotRes_sig_PHSP_5_61_PHSP.root",TString fileIN2="plotRes_sig_OmegaEta_5_61_OmegaEta.root"){
  TString file_sig_name1 = path+fileIN1;
  TString file_sig_name2 = path+fileIN2;
  TFile *finSIG1 = new TFile(file_sig_name1.Data(),"READ");
  TFile *finSIG2 = new TFile(file_sig_name2.Data(),"READ");

  TString fout_root = path+"_OE_RE_FINALE.root";
  TH1 *hrhoTh_raw_1 = (TH1D*)finSIG1->Get("hrhoTh_raw2");
  hrhoTh_raw_1->SetLineColor(kRed+1);
  hrhoTh_raw_1->SetLineWidth(2);
  TH1 *hrhoPh_raw_1 = (TH1D*)finSIG1->Get("hrhoPh_raw2");
  hrhoPh_raw_1->SetLineColor(kRed+1);
  hrhoPh_raw_1->SetLineWidth(2);
  TH1 *hetaTh_raw_1 = (TH1D*)finSIG1->Get("hetaTh_raw2");
  hetaTh_raw_1->SetLineColor(kRed+1);
  hetaTh_raw_1->SetLineWidth(2);
  TH1 *hetaPh_raw_1 = (TH1D*)finSIG1->Get("hetaPh_raw2");
  hetaPh_raw_1->SetLineColor(kRed+1);
  hetaPh_raw_1->SetLineWidth(2);

  TH1 *hrhoTh_mc_1 = (TH1D*)finSIG1->Get("hrhoTh_mc");
  hrhoTh_mc_1->SetLineColor(kRed+1);
  hrhoTh_mc_1->SetLineWidth(2);
  TH1 *hrhoPh_mc_1 = (TH1D*)finSIG1->Get("hrhoPh_mc");
  hrhoPh_mc_1->SetLineColor(kRed+1);
  hrhoPh_mc_1->SetLineWidth(2);
  TH1 *hetaTh_mc_1 = (TH1D*)finSIG1->Get("hetaTh_mc");
  hetaTh_mc_1->SetLineColor(kRed+1);
  hetaTh_mc_1->SetLineWidth(2);
  TH1 *hetaPh_mc_1 = (TH1D*)finSIG1->Get("hetaPh_mc");
  hetaPh_mc_1->SetLineColor(kRed+1);
  hetaPh_mc_1->SetLineWidth(2);

  TH1 *hpi0Th_raw_1 = (TH1D*)finSIG1->Get("hpi0Th_raw2");
  hpi0Th_raw_1->SetLineColor(kRed+1);
  hpi0Th_raw_1->SetLineWidth(2);
  TH1 *hpi0Ph_raw_1 = (TH1D*)finSIG1->Get("hpi0Ph_raw2");
  hpi0Ph_raw_1->SetLineColor(kRed+1);
  hpi0Ph_raw_1->SetLineWidth(2);
   TH1 *hpiplTh_raw_1 = (TH1D*)finSIG1->Get("hpiplTh_raw2");
   hpiplTh_raw_1->SetLineColor(kRed+1);
   hpiplTh_raw_1->SetLineWidth(2);
  TH1 *hpiplPh_raw_1 = (TH1D*)finSIG1->Get("hpiplPh_raw2");
  hpiplPh_raw_1->SetLineColor(kRed+1);
  hpiplPh_raw_1->SetLineWidth(2);

TH1 *hrhoTh_raw_2 = (TH1D*)finSIG2->Get("hrhoTh_raw2");
  hrhoTh_raw_2->SetLineColor(kBlue+1);
  hrhoTh_raw_2->SetLineWidth(2);
  TH1 *hrhoPh_raw_2 = (TH1D*)finSIG2->Get("hrhoPh_raw2");
  hrhoPh_raw_2->SetLineColor(kBlue+1);
  hrhoPh_raw_2->SetLineWidth(2);
  TH1 *hetaTh_raw_2 = (TH1D*)finSIG2->Get("hetaTh_raw2");
  hetaTh_raw_2->SetLineColor(kBlue+1);
  hetaTh_raw_2->SetLineWidth(2);
  TH1 *hetaPh_raw_2 = (TH1D*)finSIG2->Get("hetaPh_raw2");
  hetaPh_raw_2->SetLineColor(kBlue+1);
  hetaPh_raw_2->SetLineWidth(2);
  
  TH1 *hrhoTh_mc_2 = (TH1D*)finSIG2->Get("hrhoTh_mc");
  hrhoTh_mc_2->SetLineColor(kBlue+1);
  hrhoTh_mc_2->SetLineWidth(2);
  TH1 *hrhoPh_mc_2 = (TH1D*)finSIG2->Get("hrhoPh_mc");
  hrhoPh_mc_2->SetLineColor(kBlue+1);
  hrhoPh_mc_2->SetLineWidth(2);
  TH1 *hetaTh_mc_2 = (TH1D*)finSIG2->Get("hetaTh_mc");
  hetaTh_mc_2->SetLineColor(kBlue+1);
  hetaTh_mc_2->SetLineWidth(2);
  TH1 *hetaPh_mc_2 = (TH1D*)finSIG2->Get("hetaPh_mc");
  hetaPh_mc_2->SetLineColor(kBlue+1);
  hetaPh_mc_2->SetLineWidth(2);


  TH1 *hpi0Th_raw_2 = (TH1D*)finSIG2->Get("hpi0Th_raw2");
  hpi0Th_raw_2->SetLineColor(kBlue+1);
  hpi0Th_raw_2->SetLineWidth(2);
  TH1 *hpi0Ph_raw_2 = (TH1D*)finSIG2->Get("hpi0Ph_raw2");
  hpi0Ph_raw_2->SetLineColor(kBlue+1);
  hpi0Ph_raw_2->SetLineWidth(2);
   TH1 *hpiplTh_raw_2 = (TH1D*)finSIG2->Get("hpiplTh_raw2");
   hpiplTh_raw_2->SetLineColor(kBlue+1);
   hpiplTh_raw_2->SetLineWidth(2);
  TH1 *hpiplPh_raw_2 = (TH1D*)finSIG2->Get("hpiplPh_raw2");
  hpiplPh_raw_2->SetLineColor(kBlue+1);
  hpiplPh_raw_2->SetLineWidth(2);
  


  TLegend *legsh = new TLegend(0.7,0.75,0.98,0.97);
  legsh->SetFillColor(0);
  legsh->AddEntry(hrhoTh_raw_1,"PHSP","l");
  legsh->AddEntry(hrhoTh_raw_2,"#eta#omega","l");

  TCanvas c1;
  hrhoTh_raw_1->SetMaximum(1.1);
  hrhoTh_raw_1->Draw();
  hrhoTh_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pipi_th_5_61.root");
  hrhoPh_raw_1->SetMaximum(1.1);
  hrhoPh_raw_1->Draw();
  hrhoPh_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pipi_ph_5_61.root");
  hetaTh_raw_1->SetMaximum(1.1);
  hetaTh_raw_1->Draw();
  hetaTh_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pipipi_th_5_61.root");
  hetaPh_raw_1->SetMaximum(1.1);
  hetaPh_raw_1->Draw();
  hetaPh_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pipipi_ph_5_61.root");

  //  hrhoTh_mc_1->SetMaximum(1.1);
  hrhoTh_mc_1->Draw();
  hrhoTh_mc_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/mc_hc_pipi_th_5_61.root");
  //  hrhoPh_mc_1->SetMaximum(1.1);
  hrhoPh_mc_1->Draw();
  hrhoPh_mc_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/mc_hc_pipi_ph_5_61.root");
  //hetaTh_mc_1->SetMaximum(1.1);
  hetaTh_mc_1->Draw();
  hetaTh_mc_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/mc_hc_pipipi_th_5_61.root");
  //hetaPh_mc_1->SetMaximum(1.1);
  hetaPh_mc_1->Draw();
  hetaPh_mc_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/mc_hc_pipipi_ph_5_61.root");

  hpi0Th_raw_1->SetMaximum(1.1);
  hpi0Th_raw_1->Draw();
  hpi0Th_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pi0_th_5_61.root");

  hpi0Ph_raw_1->SetMaximum(1.1);
  hpi0Ph_raw_1->Draw();
  hpi0Ph_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pi0_ph_5_61.root");



  hpiplTh_raw_1->SetMaximum(1.1);
  hpiplTh_raw_1->Draw();
  hpiplTh_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pipl_th_5_61.root");

  hpiplPh_raw_1->SetMaximum(1.1);
  hpiplPh_raw_1->Draw();
  hpiplPh_raw_2->Draw("same");
  legsh->Draw();
  c1.SaveAs("pics/hc_pipl_ph_5_61.root");
  // TFile *fout = new TFile(fout_root,"RECREATE");
  // fout->Close();
}