void effResDPM(TString path="/home/karavdin/datastorage/DPMbkg_SCRUT_FwdSpec_EmcBar_Drc_Dsc_MvdGem/mom_5_61/samp0/", int nsim=2e7){
  //  const int nsim = 2e7;
  double i_nsim = 100./nsim;//%
  TString file_PH = path+"/TOT_plotRes_PHSP.root";
  TFile *finAna_PH = new TFile(file_PH.Data(),"READ");
  TH1D *hhcmfit_ph = (TH1D*)finAna_PH->Get("hhcmfit");
  int n_ph_raw = hhcmfit_ph->GetEntries();
  double eff_ph_raw = i_nsim*n_ph_raw;
  cout<<"Eff for PHSP  ana = "<<eff_ph_raw<<" %"<<endl;
  
  TString file_OE = path+"/TOT_plotRes_OmegaEta.root";
  TFile *finAna_OE = new TFile(file_OE.Data(),"READ");
  TH1D *hhcmfit_oe = (TH1D*)finAna_OE->Get("hhcmfit");
  int n_oe_raw = hhcmfit_oe->GetEntries();
  double eff_oe_raw = i_nsim*n_oe_raw;
  cout<<"Eff for OE  ana = "<<eff_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");
}