#ifdef __CLING__
#pragma cling optimize(0)
#endif
void AllPlus_Cen4_Smooth() {
  //========= ROOT Canvas Setup =========
  TCanvas *AllPlus_Cen4 = new TCanvas("AllPlus_Cen4", "",0,79,1920,921);
  gStyle->SetOptStat(0);
  AllPlus_Cen4->Range(0,0.1,3,100); // Linear Y-axis range to match reference
  AllPlus_Cen4->SetFillColor(0);
  AllPlus_Cen4->SetBorderMode(0);
  AllPlus_Cen4->SetBorderSize(2);
  AllPlus_Cen4->SetLogy(0); // Disabled log scale
  AllPlus_Cen4->SetLeftMargin(0.15);
  AllPlus_Cen4->SetFrameBorderMode(0);

  //========= Frame Histogram =========
  TH1F *hframe__1 = new TH1F("hframe__1","",1000,0,3);
  hframe__1->SetMinimum(0.1);
  hframe__1->SetMaximum(100);
  hframe__1->GetXaxis()->SetTitle("#it{p}_{T} (GeV/#it{c})");
  hframe__1->GetYaxis()->SetTitle("#frac{1}{2#pi p_{T}} #frac{d^{2}N}{dp_{T}dy} (GeV/#it{c})^{-2}");
  hframe__1->Draw(" ");

  //========= π⁺ Histogram + Spline =========
  TH1D *hPiPlus = new TH1D("hPiPlus","",100,0,10);
  hPiPlus->SetBinContent(1,72.75655); hPiPlus->SetBinError(1,2.045701);
  hPiPlus->SetBinContent(2,150.0604); hPiPlus->SetBinError(2,3.520769);
  // [...] (ALL your original π⁺ bins here)
  TSpline3 *spline_pi_plus = new TSpline3("spline_pi_plus", hPiPlus);
  spline_pi_plus->SetLineColor(kBlue);
  spline_pi_plus->SetLineWidth(2);
  spline_pi_plus->Draw("same");
  hPiPlus->SetMarkerStyle(20); // Filled circle
  hPiPlus->SetMarkerColor(kBlue);
  hPiPlus->SetMarkerSize(1.5);
  hPiPlus->Draw("same p");

  //========= π⁻ Histogram + Spline =========
  TH1D *hPiMinus = new TH1D("hPiMinus","",100,0,10);
  hPiMinus->SetBinContent(1,87.76258); hPiMinus->SetBinError(1,2.354616);
  hPiMinus->SetBinContent(2,165.0664); hPiMinus->SetBinError(2,3.781656);
  // [...] (ALL your original π⁻ bins here)
  TSpline3 *spline_pi_minus = new TSpline3("spline_pi_minus", hPiMinus);
  spline_pi_minus->SetLineColor(kBlue);
  spline_pi_minus->SetLineWidth(2);
  spline_pi_minus->SetLineStyle(2); // Dashed
  spline_pi_minus->Draw("same");
  hPiMinus->SetMarkerStyle(24); // Open circle
  hPiMinus->SetMarkerColor(kBlue);
  hPiMinus->SetMarkerSize(1.5);
  hPiMinus->Draw("same p");

  //========= K⁺ Histogram + Spline =========
  TH1D *hKPlus = new TH1D("hKPlus","",100,0,10);
  hKPlus->SetBinContent(1,3.637827); hKPlus->SetBinError(1,0.2163065);
  hKPlus->SetBinContent(2,9.094568); hKPlus->SetBinError(2,0.4300556);
  // [...] (ALL your original K⁺ bins here)
  TSpline3 *spline_k_plus = new TSpline3("spline_k_plus", hKPlus);
  spline_k_plus->SetLineColor(kRed);
  spline_k_plus->SetLineWidth(2);
  spline_k_plus->Draw("same");
  hKPlus->SetMarkerStyle(20); // Filled circle
  hKPlus->SetMarkerColor(kRed);
  hKPlus->SetMarkerSize(1.5);
  hKPlus->Draw("same p");

  //========= K⁻ Histogram + Spline =========
  TH1D *hKMinus = new TH1D("hKMinus","",100,0,10);
  hKMinus->SetBinContent(1,3.183099); hKMinus->SetBinError(1,0.1956931);
  hKMinus->SetBinContent(2,8.63984);  hKMinus->SetBinError(2,0.4138255);
  // [...] (ALL your original K⁻ bins here)
  TSpline3 *spline_k_minus = new TSpline3("spline_k_minus", hKMinus);
  spline_k_minus->SetLineColor(kRed);
  spline_k_minus->SetLineWidth(2);
  spline_k_minus->SetLineStyle(2); // Dashed
  spline_k_minus->Draw("same");
  hKMinus->SetMarkerStyle(24); // Open circle
  hKMinus->SetMarkerColor(kRed);
  hKMinus->SetMarkerSize(1.5);
  hKMinus->Draw("same p");

  //========= Proton Histogram + Spline =========
  TH1D *hProton = new TH1D("hProton","",100,0,10);
  hProton->SetBinContent(1,0.4547284); hProton->SetBinError(1,0.04547284);
  hProton->SetBinContent(2,1.818914);  hProton->SetBinError(2,0.1286166);
  // [...] (ALL your original proton bins here)
  TSpline3 *spline_proton = new TSpline3("spline_proton", hProton);
  spline_proton->SetLineColor(kGreen+2);
  spline_proton->SetLineWidth(2);
  spline_proton->Draw("same");
  hProton->SetMarkerStyle(20); // Filled circle
  hProton->SetMarkerColor(kGreen+2);
  hProton->SetMarkerSize(1.5);
  hProton->Draw("same p");

  //========= Antiproton Histogram + Spline =========
  TH1D *hAProton = new TH1D("hAProton","",100,0,10);
  hAProton->SetBinContent(2,1.364185); hAProton->SetBinError(2,0.1036557);
  hAProton->SetBinContent(3,1.818914); hAProton->SetBinError(3,0.1286166);
  // [...] (ALL your original antiproton bins here)
  TSpline3 *spline_aproton = new TSpline3("spline_aproton", hAProton);
  spline_aproton->SetLineColor(kGreen+2);
  spline_aproton->SetLineWidth(2);
  spline_aproton->SetLineStyle(2); // Dashed
  spline_aproton->Draw("same");
  hAProton->SetMarkerStyle(24); // Open circle
  hAProton->SetMarkerColor(kGreen+2);
  hAProton->SetMarkerSize(1.5);
  hAProton->Draw("same p");

  //========= Legend =========
  TLegend *leg = new TLegend(0.6,0.6,0.85,0.85);
  leg->SetBorderSize(0);
  leg->SetTextSize(0.04);
  leg->AddEntry(spline_pi_plus, "#pi^{+} + #pi^{-}", "l");
  leg->AddEntry(spline_k_plus, "K^{+} + K^{-}", "l");
  leg->AddEntry(spline_proton, "p + #bar{p}", "l");
  leg->Draw();

  //========= Centrality Labels =========
  TLatex *tex = new TLatex();
  tex->SetTextSize(0.045);
  tex->DrawLatexNDC(0.2,0.85,"0-5% (#times20)");
  tex->DrawLatexNDC(0.6,0.25,"50-60%");
}