void effSUMplot2(){
  double scen[5]={1,2,3,4,5};
  double NsimSIG=1e4;
  double NsimBKG=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]={42.44,8.38,33.35,43.14,11.96};
  double eff_sig_ph_oe[5]={0.21,0.01,0.19,0.14,0.01};
  double eff_sig_F0F0RhoOmegaEta_ph[5]={43.83,7.69,36.38,43.83,14.32};
  double eff_sig_F0F0RhoOmegaEta_oe[5]={44.18,7.79,36.95,44.39,14.48};
  double eff_sig_OmegaEta_ph[5]={44.11,7.94,36.21,42.9,13.58};
  double eff_sig_OmegaEta_oe[5]={44.11,7.96,36.2,42.89,13.59};
  double eff_sig_RhoEta_ph[5]={43.67,7.85,35.73,42.81,13.45};
  double eff_sig_RhoEta_oe[5]={43.74,7.86,35.82,42.91,13.51};
  double eff_sig_F0Eta_ph[5]={44.4,8.13,37.64,45.63,13.97};
  double eff_sig_F0Eta_oe[5]={45.01,8.21,38.2,46.21,14.24};
  double eff_sig_F01500Eta_ph[5]={46.04,7.64,37.11,44.47,13.91};
  double eff_sig_F01500Eta_oe[5]={47.43,7.85,38.14,45.71,14.45};


  double eff_bkg_ph_all[5]={1.30442,0.410251,0.56972,1.30575,0.325867};
  double eff_bkg_oe_all[5]={0.017256,0.005956,0.007461,0.017344,0.004289};

  double err_eff_sig_ph_ph[5];
  double err_eff_sig_ph_oe[5];
 double err_eff_sig_F0F0RhoOmegaEta_ph[5];
  double err_eff_sig_F0F0RhoOmegaEta_oe[5];
  double err_eff_sig_OmegaEta_ph[5];
  double err_eff_sig_OmegaEta_oe[5];
  double err_eff_sig_RhoEta_ph[5];
  double err_eff_sig_RhoEta_oe[5];
  double err_eff_sig_F0Eta_ph[5];
  double err_eff_sig_F0Eta_oe[5];
  double err_eff_sig_F01500Eta_ph[5];
  double err_eff_sig_F01500Eta_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_ph_oe[i] *=NsimSIG/100.;     
    eff_sig_F0F0RhoOmegaEta_ph[i] *=NsimSIG/100.; 
    eff_sig_F0F0RhoOmegaEta_oe[i] *=NsimSIG/100.; 
    eff_sig_OmegaEta_ph[i] *=NsimSIG/100.; 
    eff_sig_OmegaEta_oe[i] *=NsimSIG/100.; 
    eff_sig_RhoEta_ph[i] *=NsimSIG/100.; 
    eff_sig_RhoEta_oe[i] *=NsimSIG/100.; 
    eff_sig_F0Eta_ph[i] *=NsimSIG/100.; 
    eff_sig_F0Eta_oe[i] *=NsimSIG/100.; 
    eff_sig_F01500Eta_ph[i] *=NsimSIG/100.; 
    eff_sig_F01500Eta_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] = 3e2*(1./NsimSIG)*sqrt(eff_sig_ph_ph[i]*(1-eff_sig_ph_ph[i]/NsimSIG));
    err_eff_sig_ph_oe[i] = 3e2*(1./NsimSIG)*sqrt(eff_sig_ph_oe[i]*(1-eff_sig_ph_oe[i]/NsimSIG));
    err_eff_sig_F0F0RhoOmegaEta_ph[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_F0F0RhoOmegaEta_ph[i]*(1-eff_sig_F0F0RhoOmegaEta_ph[i]/NsimSIG)); 
    err_eff_sig_F0F0RhoOmegaEta_oe[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_F0F0RhoOmegaEta_oe[i]*(1-eff_sig_F0F0RhoOmegaEta_oe[i]/NsimSIG)); 
    err_eff_sig_OmegaEta_ph[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_OmegaEta_ph[i]*(1-eff_sig_OmegaEta_ph[i]/NsimSIG)); 
    err_eff_sig_OmegaEta_oe[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_OmegaEta_oe[i]*(1-eff_sig_OmegaEta_oe[i]/NsimSIG)); 
    err_eff_sig_RhoEta_ph[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_RhoEta_ph[i]*(1-eff_sig_RhoEta_ph[i]/NsimSIG)); 
    err_eff_sig_RhoEta_oe[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_RhoEta_oe[i]*(1-eff_sig_RhoEta_oe[i]/NsimSIG)); 
    err_eff_sig_F0Eta_ph[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_F0Eta_ph[i]*(1-eff_sig_F0Eta_ph[i]/NsimSIG)); 
    err_eff_sig_F0Eta_oe[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_F0Eta_oe[i]*(1-eff_sig_F0Eta_oe[i]/NsimSIG)); 
    err_eff_sig_F01500Eta_ph[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_F01500Eta_ph[i]*(1-eff_sig_F01500Eta_ph[i]/NsimSIG)); 
    err_eff_sig_F01500Eta_oe[i]  = 3e2*(1./NsimSIG)*sqrt(eff_sig_F01500Eta_oe[i]*(1-eff_sig_F01500Eta_oe[i]/NsimSIG)); 

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

    eff_sig_ph_ph[i] *=100./NsimSIG ;     eff_sig_ph_oe[i] *=100./NsimSIG;   
    eff_sig_F0F0RhoOmegaEta_ph[i] *=100./NsimSIG; 
    eff_sig_F0F0RhoOmegaEta_oe[i] *=100./NsimSIG; 
    eff_sig_OmegaEta_ph[i] *=100./NsimSIG; 
    eff_sig_OmegaEta_oe[i] *=100./NsimSIG; 
    eff_sig_RhoEta_ph[i] *=100./NsimSIG; 
    eff_sig_RhoEta_oe[i] *=100./NsimSIG; 
    eff_sig_F0Eta_ph[i] *=100./NsimSIG; 
    eff_sig_F0Eta_oe[i] *=100./NsimSIG; 
    eff_sig_F01500Eta_ph[i] *=100./NsimSIG; 
    eff_sig_F01500Eta_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_F0F0RhoOmegaEta_ph = new TGraphErrors(5,scen,eff_sig_F0F0RhoOmegaEta_ph,0,err_eff_sig_F0F0RhoOmegaEta_ph);
  gr_sig_F0F0RhoOmegaEta_ph->SetMarkerStyle(24);
  gr_sig_F0F0RhoOmegaEta_ph->SetMarkerColor(1);
  gr_sig_F0F0RhoOmegaEta_ph->SetMarkerSize(2.5);
  TGraphErrors *gr_sig_F0F0RhoOmegaEta_oe = new TGraphErrors(5,scen,eff_sig_F0F0RhoOmegaEta_oe,0,err_eff_sig_F0F0RhoOmegaEta_oe);
  gr_sig_F0F0RhoOmegaEta_oe->SetMarkerStyle(20);
  gr_sig_F0F0RhoOmegaEta_oe->SetMarkerColor(1);
  gr_sig_F0F0RhoOmegaEta_oe->SetMarkerSize(2.5);

TGraphErrors *gr_sig_OmegaEta_ph = new TGraphErrors(5,scen,eff_sig_OmegaEta_ph,0,err_eff_sig_OmegaEta_ph);
  gr_sig_OmegaEta_ph->SetMarkerStyle(24);
  gr_sig_OmegaEta_ph->SetMarkerColor(2);
  gr_sig_OmegaEta_ph->SetMarkerSize(2.5);
  TGraphErrors *gr_sig_OmegaEta_oe = new TGraphErrors(5,scen,eff_sig_OmegaEta_oe,0,err_eff_sig_OmegaEta_oe);
  gr_sig_OmegaEta_oe->SetMarkerStyle(20);
  gr_sig_OmegaEta_oe->SetMarkerColor(2);
  gr_sig_OmegaEta_oe->SetMarkerSize(2.5);


TGraphErrors *gr_sig_RhoEta_ph = new TGraphErrors(5,scen,eff_sig_RhoEta_ph,0,err_eff_sig_RhoEta_ph);
  gr_sig_RhoEta_ph->SetMarkerStyle(24);
  gr_sig_RhoEta_ph->SetMarkerColor(4);
  gr_sig_RhoEta_ph->SetMarkerSize(2.5);
  TGraphErrors *gr_sig_RhoEta_oe = new TGraphErrors(5,scen,eff_sig_RhoEta_oe,0,err_eff_sig_RhoEta_oe);
  gr_sig_RhoEta_oe->SetMarkerStyle(20);
  gr_sig_RhoEta_oe->SetMarkerColor(4);
  gr_sig_RhoEta_oe->SetMarkerSize(2.5);

TGraphErrors *gr_sig_F0Eta_ph = new TGraphErrors(5,scen,eff_sig_F0Eta_ph,0,err_eff_sig_F0Eta_ph);
  gr_sig_F0Eta_ph->SetMarkerStyle(24);
  gr_sig_F0Eta_ph->SetMarkerColor(6);
  gr_sig_F0Eta_ph->SetMarkerSize(2.5);
  TGraphErrors *gr_sig_F0Eta_oe = new TGraphErrors(5,scen,eff_sig_F0Eta_oe,0,err_eff_sig_F0Eta_oe);
  gr_sig_F0Eta_oe->SetMarkerStyle(20);
  gr_sig_F0Eta_oe->SetMarkerColor(6);
  gr_sig_F0Eta_oe->SetMarkerSize(2.5);

  TGraphErrors *gr_sig_F01500Eta_ph = new TGraphErrors(5,scen,eff_sig_F01500Eta_ph,0,err_eff_sig_F01500Eta_ph);
  gr_sig_F01500Eta_ph->SetMarkerStyle(24);
  gr_sig_F01500Eta_ph->SetMarkerColor(12);
  gr_sig_F01500Eta_ph->SetMarkerSize(2.5);
  TGraphErrors *gr_sig_F01500Eta_oe = new TGraphErrors(5,scen,eff_sig_F01500Eta_oe,0,err_eff_sig_F01500Eta_oe);
  gr_sig_F01500Eta_oe->SetMarkerStyle(20);
  gr_sig_F01500Eta_oe->SetMarkerColor(12);
  gr_sig_F01500Eta_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_F0F0RhoOmegaEta_ph,"(f_{0}(980)#eta+f_{0}(1500)#eta+#rho#eta+#omega#eta) sim, PHSP analysis","lep");
 legsh->AddEntry(gr_sig_F0F0RhoOmegaEta_oe,"(f_{0}(980)#eta+f_{0}(1500)#eta+#rho#eta+#omega#eta) sim, #eta analysis","lep");
 legsh->AddEntry(gr_sig_OmegaEta_ph,"#omega#eta sim, PHSP analysis","lep");
 legsh->AddEntry(gr_sig_OmegaEta_oe,"#omega#eta sim, #eta analysis","lep");
 legsh->AddEntry(gr_sig_RhoEta_ph,"#rho#eta sim, PHSP analysis","lep");
 legsh->AddEntry(gr_sig_RhoEta_oe,"#rho#eta sim, #eta analysis","lep");
 legsh->AddEntry(gr_sig_F0Eta_ph,"f_{0}(980)#eta sim, PHSP analysis","lep");
 legsh->AddEntry(gr_sig_F0Eta_oe,"f_{0}(980)#eta sim, #eta analysis","lep");
 legsh->AddEntry(gr_sig_F01500Eta_ph,"f_{0}(1500)#eta sim, PHSP analysis","lep");
 legsh->AddEntry(gr_sig_F01500Eta_oe,"f_{0}(1500)#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_F0F0RhoOmegaEta_ph);
  mgr_sig->Add(gr_sig_F0F0RhoOmegaEta_oe);
  mgr_sig->Add(gr_sig_RhoEta_ph);
  mgr_sig->Add(gr_sig_RhoEta_oe);
  mgr_sig->Add(gr_sig_F0Eta_ph);
  mgr_sig->Add(gr_sig_F0Eta_oe);
  mgr_sig->Add(gr_sig_F01500Eta_ph);
  mgr_sig->Add(gr_sig_F01500Eta_oe);
  mgr_sig->Add(gr_sig_OmegaEta_ph);
  mgr_sig->Add(gr_sig_OmegaEta_oe);
  mgr_sig->Draw("AP");
  legsh->Draw();

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

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

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

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

  TPaveText *pave5= new TPaveText(4.6,-5,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.05,1.55,-0.15);
  // pave1->AddText("Full set-up");
  // pave1->SetTextSize(0.025);
  // pave1->SetFillColor(kWhite);
  // pave1->Draw();

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

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

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

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