LArSamples::ShapeFitter Class Reference

#include <ShapeFitter.h>

Collaboration diagram for LArSamples::ShapeFitter:

Public Member Functions
	ShapeFitter (Chi2Params chi2Params=DefaultChi2)
bool	fit (const LArSamples::AbsShape &data, const AbsShape &reference, double &k, double &deltaT, double &chi2, const ScaledErrorData *sed=0) const
const Chi2Calc &	chi2Calc () const

Static Public Member Functions
static void	adjusted_reference (Int_t &nPar, Double_t *grad, Double_t &f, Double_t *par, Int_t iflag)

Private Attributes
Chi2Calc	m_c2c

Detailed Description

Definition at line 24 of file ShapeFitter.h.

Constructor & Destructor Documentation

ShapeFitter()

LArSamples::ShapeFitter::ShapeFitter ( Chi2Params  chi2Params = DefaultChi2 )

inline

Definition at line 28 of file ShapeFitter.h.

: m_c2c(chi2Params) { }

Member Function Documentation

adjusted_reference()

void ShapeFitter::adjusted_reference ( Int_t &  nPar, Double_t *  grad, Double_t &  f, Double_t *  par, Int_t  iflag  )

static

Definition at line 77 of file ShapeFitter.cxx.

{
  double deltaT = par[0], k = 1;
  //std::cout << "Eval at " << deltaT << std::endl;
  if (!FitterData::fitter->chi2Calc().bestRescale(*FitterData::data, *FitterData::reference, k, f, deltaT, 
                                                  FitterData::sed, FitterData::minNDof)) f = 1E99; 
  //std::cout << "-> " << k << " " << deltaT << " " << f << std::endl;
}

chi2Calc()

const Chi2Calc& LArSamples::ShapeFitter::chi2Calc ( ) const

inline

Definition at line 36 of file ShapeFitter.h.

{ return m_c2c; }

fit()

bool ShapeFitter::fit	(	const LArSamples::AbsShape &	data,
		const AbsShape &	reference,
		double &	k,
		double &	deltaT,
		double &	chi2,
		const ScaledErrorData *	sed = 0
	)		const

Definition at line 32 of file ShapeFitter.cxx.

{
  TMinuit* minuit = new TMinuit(1);  //initialize TMinuit with 1 param
  gMinuit->SetFCN(adjusted_reference);
  
  FitterData::fitter = this;
  FitterData::data = &data;
  FitterData::reference = &reference;
  FitterData::sed = sed;
 
  double nData = (data.time(data.nPoints() - 1) - data.time(0))/Definitions::samplingInterval;
  double nRef  = (reference.time(data.nPoints() - 1) - reference.time(0))/Definitions::samplingInterval;
 
  FitterData::minNDof =  (nData < nRef ? ((unsigned int)nData+1)/2 : ((unsigned int)nRef+1)/2); // we want at least n/2 samples in range
  //std::cout << minNDof << std::endl;
  Double_t arglist[10];
  Int_t ierflg = 0;
  
  arglist[0] = -1;
  gMinuit->mnexcm("SET PRINT",arglist,1,ierflg);
  arglist[0] = 0;
  gMinuit->mnexcm("SET NOW",  arglist,0,ierflg);      
 
  arglist[0] = 1;
  minuit->mnexcm("SET ERR", arglist, 1, ierflg);
  minuit->mnparm(0, "deltaT", 0, Definitions::samplingInterval/24, -5*(double)Definitions::samplingInterval, 5*(double)Definitions::samplingInterval, ierflg);
 
  arglist[0] = 500;
  arglist[1] = 1.;
  minuit->mnexcm("MIGRAD", arglist, 2, ierflg);
 
  Double_t edm, errdef;
  Int_t nvpar, nparx, icstat;
  minuit->mnstat(chi2, edm, errdef, nvpar, nparx, icstat);
 
  double err, low, high;
  int iuint;
  TString name;
  minuit->mnpout(0, name, deltaT, err, low, high, iuint);
  delete minuit;
  return m_c2c.bestRescale(data, reference, k, chi2, deltaT, sed);
}

Member Data Documentation

m_c2c

Chi2Calc LArSamples::ShapeFitter::m_c2c

private

Definition at line 40 of file ShapeFitter.h.

---

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

• ShapeFitter.h
• ShapeFitter.cxx