#ifndef rcfitter_rdf_h
#define rcfitter_rdf_h

#include <vector>
#include <TH1D.h>
#include <TH2D.h>
#include <TFitResultPtr.h>
#include <TFitResult.h>

// Fits a single peak on top of a flat background in one rocking-curve
// histogram. Exposed here in case calling code wants to fit/plot a single
// histogram directly (mirrors the old rcfitter::FitRC(TH1D*)).
TFitResultPtr FitRC(TH1D *rc);

// Runs the full rocking-curve fit analysis over the tree "rctree" found
// in treeUrl (a local path or an xrootd URL, exactly as TChain::Add()
// would accept). Returns the 8 final result histograms in the order
// base, amp, mu, sigma, mean, rms, max, peak.
//
// If start/count restrict to a sub-range of entries, RDataFrame::Range()
// is used internally, which forces single-threaded execution -- pass
// count=-1 (the default) for the normal full-range, parallel path.
std::vector<TH2D> RunRcFitter(const char *treeUrl,
                              Long64_t start = 0,
                              Long64_t count = -1);

#endif // rcfitter_rdf_h