void scanRes(TString path="/panda/scrut14/macro/scrut/results"){
  //  TString filens = path+"/plotRes_sig_FwdSpec_EmcBar_Drc_Dsc_MvdGem_";
  TString filens = path+"/plotRes_sig_OE_FwdSpec_EmcBar_Drc_Dsc_MvdGem_";
  // const int numPoints=9;
  // TString fext[numPoints]={"5_4","5_5","5_55","5_61","5_65","5_7","5_75","5_8","5_9"};
  const int numPoints=10;
       
	TString fext[numPoints]={"5_58219","5_6035","5_60617","5_6075","5_60883", "5_61017",
				 "5_6115","5_61417", "5_63553", "5_87802"};
double Pbeam[numPoints]={5.58219,5.6035,5.60617,5.6075,5.60883, 5.61017,
				 5.6115,5.61417, 5.63553, 5.87802};
  double m_invraw[numPoints];
  double m_invfit[numPoints];
  double err_m_invraw[numPoints];
  double err_m_invfit[numPoints];
  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 +="_OmegaEta.root";
    fname +="_PHSP.root";
    //   fname +="_OmegaEta_Ang.root";
    TFile *finAna = new TFile(fname.Data(),"READ");
    TH1D *hhcmraw = (TH1D*)finAna->Get("hhcmraw");
    TH1D *hhcmfit = (TH1D*)finAna->Get("hhcmfit");
    double norm_sig = 1./50000;
    nev_hcraw_sig[i] = norm_sig*(hhcmraw->GetEntries());
    m_invraw[i] = hhcmraw->GetMean();
    err_m_invraw[i] = hhcmraw->GetRMS();
    //   cout<<" err_m_invraw[i] = "<<err_m_invraw[i]<<endl;
    nev_hcfit_sig[i] = norm_sig*(hhcmfit->GetEntries());
    m_invfit[i] = hhcmfit->GetMean();
    err_m_invfit[i] = hhcmfit->GetRMS();
    //    cout<<" err_m_invfit[i] = "<<err_m_invfit[i]<<endl;
  }
  TGraphErrors *grhcraw = new TGraphErrors(numPoints,Pbeam,nev_hcraw_sig,err_m_invraw,0);
  grhcraw->SetMarkerStyle(20);
  grhcraw->SetMarkerColor(2);
  TGraphErrors *grhcfit = new TGraphErrors(numPoints,Pbeam,nev_hcfit_sig,err_m_invfit);
  grhcfit->SetMarkerStyle(21);
  grhcfit->SetMarkerColor(3);
  TMultiGraph *mgr = new TMultiGraph();
  mgr->Add(grhcraw);
  mgr->Add(grhcfit);
  mgr->Draw("AP");
  TLegend *legsh = new TLegend(0.7,0.75,0.98,0.97);
  legsh->SetFillColor(0);
  legsh->AddEntry(grhcraw,"raw inv. mass","lep");
  legsh->AddEntry(grhcfit,"inv. mass after 4C fit","lep");
  legsh->Draw();

}