void effSUMplot(){
  double scen[5]={1,2,3,4,5};
  double NsimSIG=1e4;
  double NsimBKG=1e8;
  double NsimBKG2=1e8;
  //FwdSpec EmcBar Drc Dsc MvdGem
  //FwdSpec Drc Dsc MvdGem
  //EmcBar Drc Dsc MvdGem
  //FwdSpec EmcBar Drc MvdGem
  //FwdSpec EmcBar Drc Dsc
  double eff_sig_ph_ph[5]={55.96, 0.09, 34.05, 56.08, 31.07};
  double eff_sig_ph_oe[5]={0.06,0.002, 0.04, 0.06, 0.03};
  double eff_sig_oe_ph[5]={39.92,1.93,  31.17,39.87,  20.52 };
  double eff_sig_oe_oe[5]={39.86, 1.93, 31.11, 39.79, 20.48};
  double eff_bkg_ph_all[5]={0.496, 0.004,  0.13, 0.497, 0.25};
  double eff_bkg_oe_all[5]={0.005, 0.0004, 0.002, 0.005, 0.002};
  // double eff_bkg_ph_filt[5]={};
  // double eff_bkg_oe_filt[5]={};
  // double err_eff_bkg_ph_filt[5];  double err_eff_bkg_oe_filt[5]; 
  double err_eff_sig_ph_ph[5];
  double err_eff_sig_ph_oe[5];
  double err_eff_sig_oe_ph[5];
  double err_eff_sig_oe_oe[5];
  double err_eff_bkg_ph_all[5];   double err_eff_bkg_oe_all[5];
  for(int i=0;i<5;i++){
    eff_sig_ph_ph[i] *=NsimSIG/100.;     eff_sig_oe_ph[i] *=NsimSIG/100.;     eff_sig_ph_oe[i] *=NsimSIG/100.;     eff_sig_oe_oe[i] *=NsimSIG/100.;
    eff_bkg_ph_all[i] *=double(NsimBKG)/100.;     eff_bkg_oe_all[i] *= double(NsimBKG)/100.;

    //Binomial errors
    err_eff_sig_ph_ph[i] = 2e3*(1./NsimSIG)*sqrt(eff_sig_ph_ph[i]*(1-eff_sig_ph_ph[i]/NsimSIG));
    err_eff_sig_ph_oe[i] = 2e3*(1./NsimSIG)*sqrt(eff_sig_ph_oe[i]*(1-eff_sig_ph_oe[i]/NsimSIG));
    err_eff_sig_oe_ph[i] = 2e3*(1./NsimSIG)*sqrt(eff_sig_oe_ph[i]*(1-eff_sig_oe_ph[i]/NsimSIG));
    err_eff_sig_oe_oe[i] = 2e3*(1./NsimSIG)*sqrt(eff_sig_oe_oe[i]*(1-eff_sig_oe_oe[i]/NsimSIG));

    err_eff_bkg_ph_all[i] = 2e3*(1./NsimBKG)*sqrt(eff_bkg_ph_all[i]*(1-eff_bkg_ph_all[i]/NsimBKG));
    err_eff_bkg_oe_all[i] = 2e3*(1./NsimBKG)*sqrt(eff_bkg_oe_all[i]*(1-eff_bkg_oe_all[i]/NsimBKG));

    eff_sig_ph_ph[i] *=100./NsimSIG ;     eff_sig_oe_ph[i] *=100./NsimSIG;     eff_sig_ph_oe[i] *=100./NsimSIG;     eff_sig_oe_oe[i] *=100./NsimSIG;
    eff_bkg_ph_all[i] *=100./double(NsimBKG);     eff_bkg_oe_all[i] *= 10./double(NsimBKG);
  }
  TGraphErrors *gr_sig_ph_ph = new TGraphErrors(5,scen,eff_sig_ph_ph,0,err_eff_sig_ph_ph);
  gr_sig_ph_ph->SetMarkerStyle(24);
  gr_sig_ph_ph->SetMarkerColor(8);
  gr_sig_ph_ph->SetMarkerSize(2.5);
  TGraphErrors *gr_sig_ph_oe = new TGraphErrors(5,scen,eff_sig_ph_oe,0,err_eff_sig_ph_oe);
  gr_sig_ph_oe->SetMarkerStyle(20);
  gr_sig_ph_oe->SetMarkerColor(8);
  gr_sig_ph_oe->SetMarkerSize(2.5);
 TGraphErrors *gr_sig_oe_ph = new TGraphErrors(5,scen,eff_sig_oe_ph,0,err_eff_sig_oe_ph);
  gr_sig_oe_ph->SetMarkerStyle(24);
  gr_sig_oe_ph->SetMarkerColor(2);
  gr_sig_oe_ph->SetMarkerSize(2.5);
  TGraphErrors *gr_sig_oe_oe = new TGraphErrors(5,scen,eff_sig_oe_oe,0,err_eff_sig_oe_oe);
  gr_sig_oe_oe->SetMarkerStyle(20);
  gr_sig_oe_oe->SetMarkerColor(2);
  gr_sig_oe_oe->SetMarkerSize(2.5);
 TLegend *legsh = new TLegend(0.7,0.75,0.98,0.97);
 legsh->SetFillColor(0);
 legsh->AddEntry(gr_sig_ph_ph,"PHSP sim, PHSP analysis","lep");
 legsh->AddEntry(gr_sig_ph_oe,"PHSP sim, #eta analysis","lep");
 legsh->AddEntry(gr_sig_oe_ph,"#omega#eta sim, PHSP analysis","lep");
 legsh->AddEntry(gr_sig_oe_oe,"#omega#eta sim, #eta analysis","lep");
 TMultiGraph *mgr_sig = new  TMultiGraph();
  mgr_sig->Add(gr_sig_ph_ph);
  mgr_sig->Add(gr_sig_ph_oe);
  mgr_sig->Add(gr_sig_oe_ph);
  mgr_sig->Add(gr_sig_oe_oe);
  // mgr_sig->Draw("AP");
  // legsh->Draw();

  // TPaveText *pave1= new TPaveText(0.45,-7,1.55,-1);
  // pave1->AddText("Full set-up");
  // pave1->SetTextSize(0.025);
  // pave1->SetFillColor(kWhite);
  // pave1->Draw();

  // TPaveText *pave2= new TPaveText(1.6,-7,2.5,-1);
  // pave2->AddText("w/o EmcBar");
  // pave2->SetTextSize(0.025);
  // pave2->SetFillColor(kWhite);
  // pave2->Draw();

  // TPaveText *pave3= new TPaveText(2.6,-7,3.5,-1);
  // pave3->AddText("w/o FwdSpec");
  // pave3->SetTextSize(0.025);
  // pave3->SetFillColor(kWhite);
  // pave3->Draw();

  // TPaveText *pave4= new TPaveText(3.6,-7,4.5,-1);
  // pave4->AddText("w/o Disc DIRC");
  // pave4->SetTextSize(0.025);
  // pave4->SetFillColor(kWhite);
  // pave4->Draw();

  // TPaveText *pave5= new TPaveText(4.6,-7,5.5,-1);
  // pave5->AddText("STT only");
  // pave5->SetTextSize(0.025);
  // pave5->SetFillColor(kWhite);
  // pave5->Draw();


  TGraphErrors *gr_bkg_ph_all = new TGraphErrors(5,scen,eff_bkg_ph_all,0,err_eff_bkg_ph_all);
  gr_bkg_ph_all->SetMarkerStyle(24);
  gr_bkg_ph_all->SetMarkerColor(8);
  gr_bkg_ph_all->SetMarkerSize(2.5);
  TGraphErrors *gr_bkg_oe_all = new TGraphErrors(5,scen,eff_bkg_oe_all,0,err_eff_bkg_oe_all);
  gr_bkg_oe_all->SetMarkerStyle(20);
  gr_bkg_oe_all->SetMarkerColor(8);
  gr_bkg_oe_all->SetMarkerSize(2.5);
  TMultiGraph *mgr_bkg = new  TMultiGraph();
  mgr_bkg->Add(gr_bkg_ph_all);
  mgr_bkg->Add(gr_bkg_oe_all);
  mgr_bkg->Draw("AP");
  TLegend *legsh2 = new TLegend(0.7,0.75,0.98,0.97);
  legsh2->SetFillColor(0);
  legsh2->AddEntry(gr_bkg_ph_all,"DPM all, PHSP analysis","lep");
  legsh2->AddEntry(gr_bkg_oe_all,"DPM all, #eta analysis","lep");
  legsh2->Draw();

TPaveText *pave1= new TPaveText(0.45,0,1.55,-0.05);
  pave1->AddText("Full set-up");
  pave1->SetTextSize(0.025);
  pave1->SetFillColor(kWhite);
  pave1->Draw();

  TPaveText *pave2= new TPaveText(1.6,0,2.5,-0.05);
  pave2->AddText("w/o EmcBar");
  pave2->SetTextSize(0.025);
  pave2->SetFillColor(kWhite);
  pave2->Draw();

  TPaveText *pave3= new TPaveText(2.6,0,3.5,-0.05);
  pave3->AddText("w/o FwdSpec");
  pave3->SetTextSize(0.025);
  pave3->SetFillColor(kWhite);
  pave3->Draw();

  TPaveText *pave4= new TPaveText(3.6,0,4.5,-0.05);
  pave4->AddText("w/o Disc DIRC");
  pave4->SetTextSize(0.025);
  pave4->SetFillColor(kWhite);
  pave4->Draw();

  TPaveText *pave5= new TPaveText(4.6,0,5.5,-0.05);
  pave5->AddText("STT only");
  pave5->SetTextSize(0.025);
  pave5->SetFillColor(kWhite);
  pave5->Draw();
}