void effRes(TString path="/home/karavdin/datastorage/hc_5pi_scrut2014_FwdSpec_EmcBar_Drc_Dsc_MvdGem/mom_5_61/samp0/", int nsim=1e6){
  //  const int nsim = 1e6;
  double i_nsim = 100./nsim;//%hhcmfit
  TString file_PHmod_PHana = path+"/TOT_plotRes_PHSPModel_PHSPAna.root";
  TFile *finAna_PHmod_PHana = new TFile(file_PHmod_PHana.Data(),"READ");
  TH1D *hhcmfit_ph_ph = (TH1D*)finAna_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+"/TOT_plotRes_PHSPModel_OmegaEtaAna.root";
  TFile *finAna_PHmod_OEana = new TFile(file_PHmod_OEana.Data(),"READ");
  TH1D *hhcmfit_ph_oe = (TH1D*)finAna_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;
  TString file_OEmod_PHana = path+"/TOT_plotRes_OmegaEtaModel_PHSPAna.root";
  TFile *finAna_OEmod_PHana = new TFile(file_OEmod_PHana.Data(),"READ");
  TH1D *hhcmfit_oe_ph = (TH1D*)finAna_OEmod_PHana->Get("hhcmfit");
  int n_oe_ph_raw = hhcmfit_oe_ph->GetEntries();
  double eff_oe_ph_raw = i_nsim*n_oe_ph_raw;
  cout<<"Eff for OE model and PHSP ana = "<<eff_oe_ph_raw<<" %"<<endl;
  TString file_OEmod_OEana = path+"/TOT_plotRes_OmegaEtaModel_OmegaEtaAna.root";
  TFile *finAna_OEmod_OEana = new TFile(file_OEmod_OEana.Data(),"READ");
  TH1D *hhcmfit_oe_oe = (TH1D*)finAna_OEmod_OEana->Get("hhcmfit");
  int n_oe_oe_raw = hhcmfit_oe_oe->GetEntries();
  double eff_oe_oe_raw = i_nsim*n_oe_oe_raw;
  cout<<"Eff for OE model and OE ana = "<<eff_oe_oe_raw<<" %"<<endl;

  // TString filens = path+"/plotRes_";
  // const int numPoints=13;
  // TString fext[numPoints]={"4_6","4_9","5_2","5_3","5_4","5_5","5_612", "5_65", "5_7", "5_8","6_0", "6_3", "6_6"};
  // double pbeam[numPoints]={4.6, 4.9, 5.2, 5.3, 5.4, 5.5, 5.612, 5.65, 5.7, 5.8, 6.0, 6.3, 6.6};
  // double nev_hcraw_sig[numPoints]; double nev_hcfit_sig[numPoints];
  // double nev_hcraw_bkg[numPoints]; double nev_hcfit_bkg[numPoints];
  // for(int i=0;i<numPoints;i++){
  //   TString fname = filens+ fext[i];
  //   fname +=".root";
  //   TFile *finAna = new TFile(fname.Data(),"READ");
  //   TH1D *hhcmfit = (TH1D*)finAna->Get("hhcmfit");
  //   TH1D *hhcmfit = (TH1D*)finAna->Get("hhcmfit");
  //   TH1D *hhcmfit_bkg = (TH1D*)finAna->Get("hhcmfit_bkg");
  //   TH1D *hhcmfit_bkg =  (TH1D*)finAna->Get("hhcmfit_bkg");
  //   double norm_sig = 1./10000;
  //   double norm_bkg = 1./100000;
  //   nev_hcraw_sig[i] = hhcmfit->GetEntries();
  //   nev_hcfit_sig[i] = hhcmfit->GetEntries();
  //   nev_hcraw_bkg[i] = hhcmfit_bkg->GetEntries();
  //   nev_hcfit_bkg[i] = hhcmfit_bkg->GetEntries();
  // }
  // TGraph *grhcraw = new TGraph(numPoints,pbeam,nev_hcraw_sig);
  // grhcraw->SetMarkerStyle(20);
  // grhcraw->SetMarkerColor(2);
  // TGraph *grhcfit = new TGraph(numPoints,pbeam,nev_hcfit_sig);
  // grhcfit->SetMarkerStyle(20);
  // grhcfit->SetMarkerColor(3);
  // TGraph *grhcraw_bkg = new TGraph(numPoints,pbeam,nev_hcraw_bkg);
  // grhcraw_bkg->SetMarkerStyle(20);
  // grhcraw_bkg->SetMarkerColor(4);
  // TGraph *grhcfit_bkg = new TGraph(numPoints,pbeam,nev_hcfit_bkg);
  // grhcfit_bkg->SetMarkerStyle(20);
  // grhcfit_bkg->SetMarkerColor(5);
  // TMultiGraph *mgr = new TMultiGraph();
  // mgr->Add(grhcraw);
  // mgr->Add(grhcfit);
  // mgr->Add(grhcraw_bkg);
  // mgr->Add(grhcfit_bkg);
  // mgr->Draw("AP");
}