void effLOCRes(TString path="/home/karavdin/datastorage/hc_5pi_scrut2014_FwdSpec_EmcBar_Drc_Dsc_MvdGem/mom_5_61/samp0/", int nsim=1e6, TString pbeam="5_61"){
  //  const int nsim = 1e6;
  double i_nsim = 100./nsim;//%hhcmfit
  // double i_nsim = 1./nsim;//%hhcmfit

  // TString file_PHmod_PHana = path+"/plotRes_sig_phsp_5_61_PHSP.root";
  // TFile *fin_PHmod_PHana = new TFile(file_PHmod_PHana.Data(),"READ");
  // TH1D *hhcmfit_ph_ph = (TH1D*)fin_PHmod_PHana->Get("hhcmfit");
  // int n_ph_ph_raw = hhcmfit_ph_ph->GetEntries();
  // double eff_ph_ph_raw = i_nsim*n_ph_ph_raw;
  // cout<<"Eff for PHSP model and PHSP ana = "<<eff_ph_ph_raw<<" %"<<endl;
  // TString file_PHmod_OEana = path+"/plotRes_sig_phsp_5_61_OmegaEta.root";
  // TFile *fin_PHmod_OEana = new TFile(file_PHmod_OEana.Data(),"READ");
  // TH1D *hhcmfit_ph_oe = (TH1D*)fin_PHmod_OEana->Get("hhcmfit");
  // int n_ph_oe_raw = hhcmfit_ph_oe->GetEntries();
  // double eff_ph_oe_raw = i_nsim*n_ph_oe_raw;
  // cout<<"Eff for PHSP model and OE ana = "<<eff_ph_oe_raw<<" %"<<endl;
  // cout<<" "<<endl;

  // TString file_PHmod_PHana = path+"/plotRes_sig_F0F0RhoOmegaEta_5_61_PHSP.root";
  // TFile *fin_PHmod_PHana = new TFile(file_PHmod_PHana.Data(),"READ");
  // TH1D *hhcmfit_ph_ph = (TH1D*)fin_PHmod_PHana->Get("hhcmfit");
  // int n_ph_ph_raw = hhcmfit_ph_ph->GetEntries();
  // double eff_ph_ph_raw = i_nsim*n_ph_ph_raw;
  // cout<<"Eff for F0F0RhoOmegaEta model and PHSP ana = "<<eff_ph_ph_raw<<" %"<<endl;
  // TString file_PHmod_OEana = path+"/plotRes_sig_F0F0RhoOmegaEta_5_61_OmegaEta.root";
  // TFile *fin_PHmod_OEana = new TFile(file_PHmod_OEana.Data(),"READ");
  // TH1D *hhcmfit_ph_oe = (TH1D*)fin_PHmod_OEana->Get("hhcmfit");
  // int n_ph_oe_raw = hhcmfit_ph_oe->GetEntries();
  // double eff_ph_oe_raw = i_nsim*n_ph_oe_raw;
  // cout<<"Eff for F0F0RhoOmegaEta model and OE ana = "<<eff_ph_oe_raw<<" %"<<endl;
  // cout<<" "<<endl;

  TString file_PHmod_PHana = path+"/plotRes_sig_OmegaEta_";
  file_PHmod_PHana +=pbeam;
  file_PHmod_PHana += "_PHSP.root";
  TFile *fin_PHmod_PHana = new TFile(file_PHmod_PHana.Data(),"READ");
  TH1D *hhcmfit_ph_ph = (TH1D*)fin_PHmod_PHana->Get("hhcmfit");
  int n_ph_ph_raw = hhcmfit_ph_ph->GetEntries();
  double eff_ph_ph_raw = i_nsim*n_ph_ph_raw;
  cout<<"Eff for OmegaEta model and PHSP ana = "<<eff_ph_ph_raw<<" %"<<endl;
  TString file_PHmod_OEana = path+"/plotRes_sig_OmegaEta_";
  file_PHmod_OEana+=pbeam;
  file_PHmod_OEana +="_OmegaEta.root";
  TFile *fin_PHmod_OEana = new TFile(file_PHmod_OEana.Data(),"READ");
  TH1D *hhcmfit_ph_oe = (TH1D*)fin_PHmod_OEana->Get("hhcmfit");
  int n_ph_oe_raw = hhcmfit_ph_oe->GetEntries();
  double eff_ph_oe_raw = i_nsim*n_ph_oe_raw;
  cout<<"Eff for OmegaEta model and OE ana = "<<eff_ph_oe_raw<<" %"<<endl;
  cout<<" "<<endl;

//   TString file_PHmod_PHana = path+"/plotRes_sig_RhoEta_5_61_PHSP.root";
//   TFile *fin_PHmod_PHana = new TFile(file_PHmod_PHana.Data(),"READ");
//   TH1D *hhcmfit_ph_ph = (TH1D*)fin_PHmod_PHana->Get("hhcmfit");
//   int n_ph_ph_raw = hhcmfit_ph_ph->GetEntries();
//   double eff_ph_ph_raw = i_nsim*n_ph_ph_raw;
//   cout<<"Eff for RhoEta model and PHSP ana = "<<eff_ph_ph_raw<<" %"<<endl;
//   TString file_PHmod_OEana = path+"/plotRes_sig_RhoEta_5_61_OmegaEta.root";
//   TFile *fin_PHmod_OEana = new TFile(file_PHmod_OEana.Data(),"READ");
//   TH1D *hhcmfit_ph_oe = (TH1D*)fin_PHmod_OEana->Get("hhcmfit");
//   int n_ph_oe_raw = hhcmfit_ph_oe->GetEntries();
//   double eff_ph_oe_raw = i_nsim*n_ph_oe_raw;
//   cout<<"Eff for RhoEta model and OE ana = "<<eff_ph_oe_raw<<" %"<<endl;
//   cout<<" "<<endl;


//  TString file_PHmod_PHana = path+"/plotRes_sig_F0Eta_5_61_PHSP.root";
//   TFile *fin_PHmod_PHana = new TFile(file_PHmod_PHana.Data(),"READ");
//   TH1D *hhcmfit_ph_ph = (TH1D*)fin_PHmod_PHana->Get("hhcmfit");
//   int n_ph_ph_raw = hhcmfit_ph_ph->GetEntries();
//   double eff_ph_ph_raw = i_nsim*n_ph_ph_raw;
//   cout<<"Eff for F0Eta model and PHSP ana = "<<eff_ph_ph_raw<<" %"<<endl;
//   TString file_PHmod_OEana = path+"/plotRes_sig_F0Eta_5_61_OmegaEta.root";
//   TFile *fin_PHmod_OEana = new TFile(file_PHmod_OEana.Data(),"READ");
//   TH1D *hhcmfit_ph_oe = (TH1D*)fin_PHmod_OEana->Get("hhcmfit");
//   int n_ph_oe_raw = hhcmfit_ph_oe->GetEntries();
//   double eff_ph_oe_raw = i_nsim*n_ph_oe_raw;
//   cout<<"Eff for F0Eta model and OE ana = "<<eff_ph_oe_raw<<" %"<<endl;
//   cout<<" "<<endl;


// TString file_PHmod_PHana = path+"/plotRes_sig_F01500Eta_5_61_PHSP.root";
//   TFile *fin_PHmod_PHana = new TFile(file_PHmod_PHana.Data(),"READ");
//   TH1D *hhcmfit_ph_ph = (TH1D*)fin_PHmod_PHana->Get("hhcmfit");
//   int n_ph_ph_raw = hhcmfit_ph_ph->GetEntries();
//   double eff_ph_ph_raw = i_nsim*n_ph_ph_raw;
//   cout<<"Eff for F01500Eta model and PHSP ana = "<<eff_ph_ph_raw<<" %"<<endl;
//   TString file_PHmod_OEana = path+"/plotRes_sig_F01500Eta_5_61_OmegaEta.root";
//   TFile *fin_PHmod_OEana = new TFile(file_PHmod_OEana.Data(),"READ");
//   TH1D *hhcmfit_ph_oe = (TH1D*)fin_PHmod_OEana->Get("hhcmfit");
//   int n_ph_oe_raw = hhcmfit_ph_oe->GetEntries();
//   double eff_ph_oe_raw = i_nsim*n_ph_oe_raw;
//   cout<<"Eff for F01500Eta model and OE ana = "<<eff_ph_oe_raw<<" %"<<endl;
//   cout<<" "<<endl;

}