LArSamples::Chi2Calc Class Reference

#include <Chi2Calc.h>

Inheritance diagram for LArSamples::Chi2Calc:

Collaboration diagram for LArSamples::Chi2Calc:

Public Member Functions
	Chi2Calc (int pars=0)
	Constructor (takes ownership of LArCellInfo object)
virtual	~Chi2Calc ()
TVectorD	deltas (const AbsShape &data, const AbsShape &reference, CovMatrix &errors, const ScaledErrorData *shapeError=0, int lwb=-1, int upb=-1, bool noDataError=false)
double	chi2 (const AbsShape &data, const AbsShape &reference, const ScaledErrorData *shapeError=0, int lwb=-1, int upb=-1)
double	scalarProduct (const TVectorD &values1, const TVectorD &values2, const CovMatrix &invCovMat) const
int	lwb () const
int	upb () const
unsigned int	nDof () const
bool	bestRescale (const AbsShape &data, const AbsShape &reference, double &k, double &chi2, double deltaT=0, const ScaledErrorData *sed=0, unsigned int minNDof=0) const
bool	isInRange (int i) const
bool	hasSameRange (int lw, int up) const
bool	hasSameRange (const IndexRange &other) const
bool	hasSameRange (const TVectorD &v) const
bool	providesRange (int lw, int up) const
bool	providesRange (const IndexRange &other) const
bool	providesRange (const TVectorD &v) const
TString	rangeStr () const
int	commonLwb (const IndexRange &other) const
int	commonUpb (const IndexRange &other) const
bool	checkRange (int &l, int &h) const

Private Member Functions
bool	useCorrs () const

Private Attributes
int	m_pars
int	m_lwb
int	m_upb

Detailed Description

Definition at line 26 of file Chi2Calc.h.

Constructor & Destructor Documentation

Chi2Calc()

LArSamples::Chi2Calc::Chi2Calc ( int pars = 0 )

inline

Constructor (takes ownership of LArCellInfo object)

Definition at line 31 of file Chi2Calc.h.

      : m_pars(pars), m_lwb(-1), m_upb(-1) { }

~Chi2Calc()

virtual LArSamples::Chi2Calc::~Chi2Calc ( )

inlinevirtual

Definition at line 34 of file Chi2Calc.h.

{ }

Member Function Documentation

bestRescale()

bool Chi2Calc::bestRescale	(	const AbsShape &	data,
		const AbsShape &	reference,
		double &	k,
		double &	chi2,
		double	deltaT = 0,
		const ScaledErrorData *	sed = 0,
		unsigned int	minNDof = 0 ) const

Definition at line 88 of file Chi2Calc.cxx.

{
  k = 1;
  chi2 = -1;
  ScaledShiftedShape shiftedData(data, 1, -deltaT);
  TVectorD r;
  CovMatrix errors;  
  bool result = reference.interpolate(shiftedData, r, errors, -1, -1); // lwb and upb will come from data
  if (!result || r.GetNrows() < (int)minNDof) {
    k = 0;
    chi2 = 1E99;
    return true;
  }
  
  if (sed) {
    ScaledErrorData sed2(*sed, sed->scale(), sed->time() + deltaT);
    if (sed2.providesRange(r)) r += sed2.offsets(r.GetLwb(), r.GetUpb());
  }
  
  TVectorD v = data.values(r.GetLwb(), r.GetUpb());
  CovMatrix invCovMat = data.invCovarianceMatrix(r.GetLwb(), r.GetUpb(), errors, useCorrs());
  double sumR2 = scalarProduct(r, r, invCovMat);
  double sumRV = scalarProduct(r, v, invCovMat);
  double sumV2 = scalarProduct(v, v, invCovMat);
  k = sumR2==0 ? 0 : sumRV/sumR2; // optimal K-factor
  chi2 = sumV2 - 2*k*sumRV + k*k*sumR2; // chi2
 
  return true;
}

checkRange()

bool IndexRange::checkRange ( int & l, int & h ) const

inherited

Definition at line 14 of file IndexRange.cxx.

{  
  if (l < 0) 
    l = lwb();
  else if (!isInRange(l)) {
    cout << "IndexRange::checkRange : lower bound " << l << " is out of bounds" << endl;
    return false;
  }
 
  if (h < 0) 
    h = upb();
  else if (!isInRange(h)) {
    cout << "IndexRange::checkRange : upper bound " << h << " is out of bounds" << endl;
    return false;
  }
 
  return true;
}

chi2()

double Chi2Calc::chi2	(	const AbsShape &	data,
		const AbsShape &	reference,
		const ScaledErrorData *	shapeError = 0,
		int	lwb = -1,
		int	upb = -1 )

Definition at line 60 of file Chi2Calc.cxx.

{ 
  double chi2Val = 0;
  CovMatrix errors;
  TVectorD dv = deltas(data, reference, errors, shapeError, lw, up, true);
  if (dv.GetNrows() == 0) return -1;
  
  // OFC version (gives correct error terms for individual samplings after OFC fit is performed
  if (m_pars & OFCChi2) {
    OFC ofc(reference, data, errors, lwb(), upb(), shapeError, useCorrs());
    return scalarProduct(dv, dv, ofc.invGamma());
  }
  // Basic version --- fast, simple
  if (m_pars & BasicChi2) {
    for (int k = lwb(); k < upb(); k++)
      chi2Val += dv(k)*dv(k)/(TMath::Power(data.error(k),2) + errors(k,k));
    return chi2Val;
  }
  // Default version : 
  CovMatrix invCovMat = data.invCovarianceMatrix(lwb(), upb(), errors, useCorrs());
/*  dv.Print();
  data.covarianceMatrix(lwb(), upb()).Print();
  data.covarianceMatrix(lwb(), upb(), errors, useCorrs()).Print();*/
  return scalarProduct(dv, dv, invCovMat);
}

commonLwb()

int LArSamples::IndexRange::commonLwb ( const IndexRange & other ) const

inlineinherited

Definition at line 39 of file IndexRange.h.

{ return (other.lwb() > lwb() ? other.lwb() : lwb()); }

commonUpb()

int LArSamples::IndexRange::commonUpb ( const IndexRange & other ) const

inlineinherited

Definition at line 40 of file IndexRange.h.

{ return (other.upb() < upb() ? other.upb() : upb()); }

deltas()

TVectorD Chi2Calc::deltas	(	const AbsShape &	data,
		const AbsShape &	reference,
		CovMatrix &	errors,
		const ScaledErrorData *	shapeError = 0,
		int	lwb = -1,
		int	upb = -1,
		bool	noDataError = false )

Definition at line 36 of file Chi2Calc.cxx.

{
  TVectorD refVals;
  if (!reference.interpolate(data, refVals, errors, lw, up)) {
    cout << "Chi2Calc::deltas : Could not find overlap region between data and reference in [" << lw << ", " << up << "]." << endl;
    return TVectorD();
  }
  m_lwb = refVals.GetLwb();
  m_upb = refVals.GetUpb();
 
  TVectorD deltas = data.values(lwb(), upb());
  deltas -= refVals;
  if (shapeError && shapeError->providesRange(*this)) {
    TVectorD offsets = shapeError->offsets(lwb(), upb());
    deltas -= offsets;
    errors += shapeError->errors(lwb(), upb());
  }
 
  if (!noDataError) errors = data.covarianceMatrix(lwb(), upb(), errors, useCorrs());
  return deltas;
}

hasSameRange() [1/3]

bool LArSamples::IndexRange::hasSameRange ( const IndexRange & other ) const

inlineinherited

Definition at line 30 of file IndexRange.h.

{ return hasSameRange(other.lwb(), other.upb()); }

hasSameRange() [2/3]

bool LArSamples::IndexRange::hasSameRange ( const TVectorD & v ) const

inlineinherited

Definition at line 31 of file IndexRange.h.

{ return hasSameRange(v.GetLwb(), v.GetUpb()); }

hasSameRange() [3/3]

bool LArSamples::IndexRange::hasSameRange ( int lw, int up ) const

inlineinherited

Definition at line 29 of file IndexRange.h.

{ return (lwb() == lw && upb() == up); }

isInRange()

bool LArSamples::IndexRange::isInRange ( int i ) const

inlineinherited

Definition at line 27 of file IndexRange.h.

{ return (i >= lwb() && i <= upb()); }

lwb()

int LArSamples::Chi2Calc::lwb ( ) const

inlinevirtual

Implements LArSamples::IndexRange.

Definition at line 45 of file Chi2Calc.h.

{ return m_lwb; }

nDof()

unsigned int LArSamples::Chi2Calc::nDof ( ) const

inline

Definition at line 47 of file Chi2Calc.h.

{ return upb() - lwb() + 1; }  

providesRange() [1/3]

bool LArSamples::IndexRange::providesRange ( const IndexRange & other ) const

inlineinherited

Definition at line 34 of file IndexRange.h.

{ return providesRange(other.lwb(), other.upb()); }

providesRange() [2/3]

bool LArSamples::IndexRange::providesRange ( const TVectorD & v ) const

inlineinherited

Definition at line 35 of file IndexRange.h.

{ return providesRange(v.GetLwb(), v.GetUpb()); }

providesRange() [3/3]

bool LArSamples::IndexRange::providesRange ( int lw, int up ) const

inlineinherited

Definition at line 33 of file IndexRange.h.

{ return (lwb() <= lw && upb() >= up); }

rangeStr()

TString LArSamples::IndexRange::rangeStr ( ) const

inlineinherited

Definition at line 37 of file IndexRange.h.

{ return Form("[%d, %d]", lwb(), upb()); }

scalarProduct()

double Chi2Calc::scalarProduct	(	const TVectorD &	values1,
		const TVectorD &	values2,
		const CovMatrix &	invCovMat ) const

Definition at line 19 of file Chi2Calc.cxx.

{
  if (values1.GetLwb() != invCovMat.GetRowLwb() || values1.GetUpb() != invCovMat.GetRowUpb()) return -1;
  if (values2.GetLwb() != invCovMat.GetColLwb() || values2.GetUpb() != invCovMat.GetColUpb()) return -1;
 
  double chi2Val = 0;
  
  for (int i = values1.GetLwb(); i <= values1.GetUpb(); i++) {
    for (int j = values2.GetLwb(); j <= values2.GetUpb(); j++) {
      chi2Val += invCovMat(i, j)*values1(i)*values2(j);
    }
  }
  return chi2Val;
}

upb()

int LArSamples::Chi2Calc::upb ( ) const

inlinevirtual

Implements LArSamples::IndexRange.

Definition at line 46 of file Chi2Calc.h.

{ return m_upb; }

useCorrs()

bool LArSamples::Chi2Calc::useCorrs ( ) const

inlineprivate

Definition at line 54 of file Chi2Calc.h.

{ return !(m_pars & NoCorrs); }

Member Data Documentation

m_lwb

int LArSamples::Chi2Calc::m_lwb

private

Definition at line 58 of file Chi2Calc.h.

m_pars

int LArSamples::Chi2Calc::m_pars

private

Definition at line 56 of file Chi2Calc.h.

m_upb

int LArSamples::Chi2Calc::m_upb

private

Definition at line 58 of file Chi2Calc.h.

---

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

• Chi2Calc.h
• Chi2Calc.cxx