LArSamples::AbsShape Class Reference

#include <AbsShape.h>

Inheritance diagram for LArSamples::AbsShape:

Collaboration diagram for LArSamples::AbsShape:

Public Member Functions
virtual	~AbsShape ()
virtual unsigned int	nPoints () const =0
virtual double	value (unsigned int i) const =0
virtual double	covariance (unsigned int i, unsigned int j) const =0
virtual double	time (unsigned int i) const =0
virtual double	error (unsigned int i) const
TVectorD	values (int lwb, int upb) const
int	findTimeInterval (double time) const
int	interpolate (double time, double &value, double &error) const
int	interpolateDiff (double time, double &diff) const
bool	interpolate (const AbsShape &other, TVectorD &values, CovMatrix &errors, int lwb=-1, int upb=-1) const
bool	interpolateDiff (const AbsShape &other, TVectorD &diffs, int lwb=-1, int upb=-1) const
CovMatrix	covarianceMatrix (int lwb=-1, int upb=-1, const CovMatrix &refErr=CovMatrix(), bool withCorrs=true) const
CovMatrix	invCovarianceMatrix (int lwb=-1, int upb=-1, const CovMatrix &refErr=CovMatrix(), bool withCorrs=true) const
CovMatrix	covarianceMatrix (unsigned int nPoints, bool withCorrs=true) const
CovMatrix	invCovarianceMatrix (unsigned int nPoints, bool withCorrs=true) const
double	maxValue (bool withErrors=false) const
double	minValue (bool withErrors=false) const
int	maxPosition () const
int	minPosition () const
TGraphErrors *	graph (bool timeInUnitOfSamples=false) const
SimpleShape *	resample (unsigned int nPts) const

Protected Member Functions
	AbsShape ()

Detailed Description

Definition at line 28 of file AbsShape.h.

Constructor & Destructor Documentation

~AbsShape()

virtual LArSamples::AbsShape::~AbsShape ( )

inlinevirtual

Definition at line 32 of file AbsShape.h.

{ }

AbsShape()

LArSamples::AbsShape::AbsShape ( )

inlineprotected

Definition at line 72 of file AbsShape.h.

{ }

Member Function Documentation

covariance()

virtual double LArSamples::AbsShape::covariance ( unsigned int  i, unsigned int  j  ) const

pure virtual

Implemented in LArSamples::Data, LArSamples::GraphShape, LArSamples::SimpleShape, and LArSamples::ScaledShiftedShape.

covarianceMatrix() [1/2]

CovMatrix AbsShape::covarianceMatrix	(	int	lwb = -1,
		int	upb = -1,
		const CovMatrix &	refErr = CovMatrix(),
		bool	withCorrs = true
	)		const

Definition at line 198 of file AbsShape.cxx.

{
  if (lwb < 0 || upb < 0) { lwb = 0; upb = nPoints() - 1; }
  CovMatrix c(lwb, upb);
  for (int i = lwb; i <= upb; i++) {
    for (int j = lwb; j <= upb; j++) {
      double cov = (withCorrs || i == j ? covariance(i, j) : 0);
      if (refErr.GetNrows() > 0) cov += refErr(i, j);
      c(i, j) = cov;
    }
  }
  return c;
}

covarianceMatrix() [2/2]

CovMatrix AbsShape::covarianceMatrix	(	unsigned int	nPoints,
		bool	withCorrs = true
	)		const

Definition at line 223 of file AbsShape.cxx.

{
  return covarianceMatrix(0, n - 1, CovMatrix(), withCorrs);
}

error()

double AbsShape::error ( unsigned int  i ) const

virtual

Reimplemented in LArSamples::ScaledShiftedShape.

Definition at line 24 of file AbsShape.cxx.

{ 
  return TMath::Sqrt(covariance(i, i)); 
}

findTimeInterval()

int AbsShape::findTimeInterval

(

double

time

)

const

Definition at line 71 of file AbsShape.cxx.

{
  if (t < time(0)) return -1;
  if (t > time(nPoints() - 1)) return int(nPoints()) - 1;
  
  for (unsigned int i = 0; i < nPoints() - 1; i++)
    if (time(i) <= t && t < time(i+1)) return i;
 
  return nPoints() - 1;
}

graph()

TGraphErrors * AbsShape::graph

(

bool

timeInUnitOfSamples = false

)

const

Definition at line 186 of file AbsShape.cxx.

{
  TGraphErrors* graph = new TGraphErrors(nPoints());
  double unit = (timeInUnitOfSamples ? Definitions::samplingInterval : 1);
  for (unsigned int i = 0; i < nPoints(); i++) {
    graph->SetPoint(i, time(i)/unit, value(i));    
    graph->SetPointError(i, 0, error(i));
  }
  return graph;
}

interpolate() [1/2]

bool AbsShape::interpolate	(	const AbsShape &	other,
		TVectorD &	values,
		CovMatrix &	errors,
		int	lwb = -1,
		int	upb = -1
	)		const

Definition at line 147 of file AbsShape.cxx.

{
  if (lwb < 0) lwb = 0;
  if (upb < 0) upb = other.nPoints() - 1;
  values.ResizeTo(lwb, upb);
  errors.ResizeTo(lwb, upb, lwb, upb);  
  int actualLwb = -1, actualUpb = upb;
  for (int i = lwb; i <= upb; i++) {
    double val, err;
    int stat = interpolate(other.time(i), val, err);
    if (stat ==  0 && actualLwb < 0) actualLwb = i;
    if (stat == +1) { actualUpb = i - 1; break; }
    values(i) = val;
    errors(i, i) = err*err;
  }
  values.ResizeTo(actualLwb, actualUpb);
  errors.ResizeTo(actualLwb, actualUpb, actualLwb, actualUpb);  
  return (actualLwb >= 0);
}

interpolate() [2/2]

int AbsShape::interpolate	(	double	time,
		double &	value,
		double &	error
	)		const

Definition at line 83 of file AbsShape.cxx.

{
  if (nPoints() < 2) return -1;
  
  int status = 0;
  int i1 = findTimeInterval(t);
  
  if (i1 < 0) {
    status = -1;
    i1 = 0;
  }
 
  if (i1 > int(nPoints()) - 2) {
    status = 1;
    i1 = nPoints() - 2;
  }
 
  double delta = (t - time(i1))/(time(i1+1) - time(i1));
  val = value(i1) + delta*(value(i1 + 1) - value(i1));
  err = error(i1) + delta*(error(i1 + 1) - error(i1));
  //cout << t << " " << i1 << " " << time(i1) << " " << delta << " " << value(i1) << " " << value(i1+1) << " " << val << " " << status << endl;
  return status;
}

interpolateDiff() [1/2]

bool AbsShape::interpolateDiff	(	const AbsShape &	other,
		TVectorD &	diffs,
		int	lwb = -1,
		int	upb = -1
	)		const

Definition at line 168 of file AbsShape.cxx.

{
  if (lwb < 0) lwb = 0;
  if (upb < 0) upb = other.nPoints() - 1;
  diffs.ResizeTo(lwb, upb);
  int actualLwb = -1, actualUpb = upb;  
  for (int i = lwb; i <= upb; i++) {
    double diff;
    int stat = interpolateDiff(other.time(i), diff);
    if (stat ==  0 && actualLwb < 0) actualLwb = i;
    if (stat == +1) { actualUpb = i - 1; break; }
    diffs(i) = diff;
  }
  diffs.ResizeTo(actualLwb, actualUpb);
  return (actualLwb >= 0);
}

interpolateDiff() [2/2]

int AbsShape::interpolateDiff	(	double	time,
		double &	diff
	)		const

Definition at line 108 of file AbsShape.cxx.

{
  int status = 0;
  int i1 = findTimeInterval(t);  
  if (i1 < 0) {
    i1 = 0;
    status = -1;
  }  
  if (i1 > int(nPoints()) - 1) {
    i1 = nPoints() - 1;
    status = 1;
  }  
  if (i1 > 1 && i1 < int(nPoints()) - 2)
    diff = (value(i1 - 2) - 8*value(i1 - 1) + 8*value(i1 + 1) - value(i1 + 2))/(time(i1 + 2) - time(i1 - 2))/3;
  else
    if (i1 > 0) 
      if (i1 <  int(nPoints()) - 1)
        diff = (value(i1 + 1) - value(i1 - 1))/(time(i1 + 1) - time(i1 - 1));
      else
        diff = (value(i1) - value(i1 - 1))/(time(i1) - time(i1 - 1));
    else 
      diff = (value(i1 + 1) - value(i1))/(time(i1 + 1) - time(i1));
 
  return status;
}

invCovarianceMatrix() [1/2]

CovMatrix AbsShape::invCovarianceMatrix	(	int	lwb = -1,
		int	upb = -1,
		const CovMatrix &	refErr = CovMatrix(),
		bool	withCorrs = true
	)		const

Definition at line 215 of file AbsShape.cxx.

{ 
  return covarianceMatrix(lwb, upb, refErr, withCorrs).Invert(); 
}

invCovarianceMatrix() [2/2]

CovMatrix AbsShape::invCovarianceMatrix	(	unsigned int	nPoints,
		bool	withCorrs = true
	)		const

Definition at line 230 of file AbsShape.cxx.

{
  return covarianceMatrix(n, withCorrs).Invert();
}

maxPosition()

int AbsShape::maxPosition

(

)

const

Definition at line 52 of file AbsShape.cxx.

{
  int maxPosition = -1;
  double maxValue = -1;
  for (unsigned int i = 0; i < nPoints(); i++)
    if (value(i) > maxValue) { maxPosition = i; maxValue = value(i); }
  return maxPosition;
}

maxValue()

double AbsShape::maxValue

(

bool

withErrors = false

)

const

Returns

sample max parameters

Definition at line 30 of file AbsShape.cxx.

{
  double maxValue = -DBL_MAX;
  for (unsigned int i = 0; i < nPoints(); i++) {
    double val = value(i) + (withErrors ? error(i) : 0);
    if (val > maxValue) maxValue = val;
  }
  return maxValue;
}

minPosition()

int AbsShape::minPosition

(

)

const

Definition at line 61 of file AbsShape.cxx.

{
  int minPosition = -1;
  double minValue = Definitions::none;
  for (unsigned int i = 0; i < nPoints(); i++)
    if (value(i) < minValue) { minPosition = i; minValue = value(i); }
  return minPosition;
}

minValue()

double AbsShape::minValue

(

bool

withErrors = false

)

const

Definition at line 41 of file AbsShape.cxx.

{
  double minValue = DBL_MAX;
  for (unsigned int i = 0; i < nPoints(); i++) {
    double val = value(i) - (withErrors ? error(i) : 0);
    if (val < minValue) minValue = val;
  }
  return minValue;
}

nPoints()

virtual unsigned int LArSamples::AbsShape::nPoints ( ) const

pure virtual

Returns

data points

Implemented in LArSamples::Data, LArSamples::GraphShape, LArSamples::SimpleShape, and LArSamples::ScaledShiftedShape.

resample()

SimpleShape * AbsShape::resample

(

unsigned int

nPts

)

const

Definition at line 237 of file AbsShape.cxx.

{
  std::vector<double> values, errors;
  double t0 = time(0);
  double t1 = time(nPoints() - 1);
  double dT = (t1 - t0)/nPts;
  double t = t0;
  for (unsigned int i = 0; i < nPts; i++, t += dT) {
    double val, err;
    int inRange = interpolate(t, val, err);
    if (inRange != 0) return nullptr;
    values.push_back(val);
    errors.push_back(err);
  }
  
  return new SimpleShape(values, errors, dT, t0);
}

time()

virtual double LArSamples::AbsShape::time ( unsigned int  i ) const

pure virtual

Implemented in LArSamples::Data, LArSamples::GraphShape, LArSamples::SimpleShape, and LArSamples::ScaledShiftedShape.

value()

virtual double LArSamples::AbsShape::value ( unsigned int  i ) const

pure virtual

Implemented in LArSamples::Data, LArSamples::GraphShape, LArSamples::SimpleShape, and LArSamples::ScaledShiftedShape.

values()

TVectorD AbsShape::values	(	int	lwb,
		int	upb
	)		const

Definition at line 135 of file AbsShape.cxx.

{
  if (lwb < 0) lwb = 0;
  if (upb < 0) upb = (int)nPoints() - 1;
  if (upb >= (int)nPoints()) upb = (int)nPoints() - 1;
  TVectorD vals(lwb, upb);
  for (int i = lwb; i <= upb; i++) 
    vals(i) = value(i);
  return vals;
}

---

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

• AbsShape.h
• AbsShape.cxx