ana::FitToFourier Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Core/Utilities.h"

Public Member Functions
	FitToFourier (TH1 *h, double xlo, double xhi, int NOsc)
	~FitToFourier ()
TF1 *	Fit () const
double	operator() (double *x, double *par) const

Private Attributes
const TH1 *	fHist
const double	fxlo
const double	fxhi
const int	fNOsc

Detailed Description

Definition at line 210 of file Utilities.h.

Constructor & Destructor Documentation

ana::FitToFourier::FitToFourier	(	TH1 *	h,
		double	xlo,
		double	xhi,
		int	NOsc
	)

Definition at line 482 of file Utilities.cxx.

    : fHist(h), fxlo(xlo), fxhi(xhi), fNOsc(NOsc)
  {
  }

ana::FitToFourier::~FitToFourier

(

)

Definition at line 488 of file Utilities.cxx.

  {
  }

Member Function Documentation

TF1 * ana::FitToFourier::Fit

(

)

const

Definition at line 503 of file Utilities.cxx.

References plot_validation_datamc::c, std::cos(), MakeMiniprodValidationCuts::f, fHist, fNOsc, fxhi, fxlo, MECModelEnuComparisons::i, M_PI, cet::sqlite::max(), min(), getGoodRuns4SAM::n, plotROC::nBins, std::sin(), ana::UniqueName(), submit_syst::x, and submit_syst::y.

Referenced by SmoothWithFourierFit().

  {
    std::vector<double> s(fNOsc, 0.);
    std::vector<double> c(fNOsc, 0.);
    int nBins = 0;
    for(int i = 1; i <= fHist->GetNbinsX(); ++i){
      const double x = M_PI * fHist->GetXaxis()->GetBinCenter(i);
      const double y = fHist->GetBinContent(i);

      if(y == 0) continue;
      ++nBins;

      for(int n = 0; n <= fNOsc; ++n){
        s[n] += y * sin(n*x);
        c[n] += y * cos(n*x);
      }
    }

    for(int n = 0; n <= fNOsc; ++n){
      s[n] *= 2./nBins;
      c[n] *= 2./nBins;
    }

    TF1* f = new TF1(UniqueName().c_str(), this, fxlo, fxhi, 2*fNOsc+1);

    f->SetParameter(0, c[0]/2);
    for(int n = 1; n <= fNOsc; ++n){
      f->SetParameter(n*2-1, s[n]);
      f->SetParameter(n*2,   c[n]);
    }

    // Because ROOT is having problems drawing f if I don't
    double min = fHist->GetMinimum();
    double max = fHist->GetMaximum();
    f->GetYaxis()->SetRangeUser(0.8*min, 1.2*max);
    return f;
  }

double ana::FitToFourier::operator()	(	double *	x,
		double *	par
	)		const

Definition at line 493 of file Utilities.cxx.

References std::cos(), fNOsc, MECModelEnuComparisons::i, M_PI, std::sin(), and febshutoff_auto::val.

  {
    double x0 = x[0];
    double val = par[0];
    for (int i = 1; i <= fNOsc; i++)
      val += par[2*i-1]*sin(i*M_PI*x0) + par[2*i]*cos(i*M_PI*x0);
    return val;
  }

Member Data Documentation

const TH1* ana::FitToFourier::fHist

private

Definition at line 219 of file Utilities.h.

Referenced by Fit().

const int ana::FitToFourier::fNOsc

private

Definition at line 222 of file Utilities.h.

Referenced by Fit(), and operator()().

const double ana::FitToFourier::fxhi

private

Definition at line 221 of file Utilities.h.

Referenced by Fit().

const double ana::FitToFourier::fxlo

private

Definition at line 220 of file Utilities.h.

Referenced by Fit().

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Core/Utilities.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-02-25/CAFAna/Core/Utilities.cxx