18 vector<double>
extractBinning(
const unique_ptr<TH1> &h3D,
const TString &axis)
20 const TAxis *naxis =
nullptr;
21 if (axis ==
"X") naxis = h3D->GetXaxis();
22 else if (axis ==
"Y") naxis = h3D->GetYaxis();
23 else if (axis ==
"Z") naxis = h3D->GetZaxis();
25 int ndiv = naxis->GetNbins();
26 vector<double> binning(ndiv+1);
27 naxis->GetLowEdge(&binning[0]);
28 binning[ndiv] = naxis->GetBinUpEdge(ndiv);