LArSamples::ScaledShiftedShape Class Reference

#include <ScaledShiftedShape.h>

Inheritance diagram for LArSamples::ScaledShiftedShape:

Collaboration diagram for LArSamples::ScaledShiftedShape:

Public Member Functions
	ScaledShiftedShape (const AbsShape &base, double scaling=1, double shift=0)
	Constructor.
	ScaledShiftedShape (const ScaledShiftedShape &other)
virtual	~ScaledShiftedShape ()
const AbsShape &	base () const
double	scaling () const
double	shift () const
unsigned int	nPoints () const
double	value (unsigned int i) const
double	error (unsigned int i) const
double	time (unsigned int i) const
double	covariance (unsigned int i, unsigned int j) const
TVectorD	values (int lwb, int upb) const
int	findTimeInterval (double time) const
int	interpolate (double time, double &value, double &error) const
bool	interpolate (const AbsShape &other, TVectorD &values, CovMatrix &errors, int lwb=-1, int upb=-1) const
int	interpolateDiff (double time, double &diff) 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	covarianceMatrix (unsigned int nPoints, bool withCorrs=true) const
CovMatrix	invCovarianceMatrix (int lwb=-1, int upb=-1, const CovMatrix &refErr=CovMatrix(), 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

Private Attributes
const AbsShape *	m_base
double	m_scaling
double	m_shift

Detailed Description

Definition at line 17 of file ScaledShiftedShape.h.

Constructor & Destructor Documentation

ScaledShiftedShape() [1/2]

LArSamples::ScaledShiftedShape::ScaledShiftedShape ( const AbsShape & base, double scaling = 1, double shift = 0 )

inline

Constructor.

Definition at line 22 of file ScaledShiftedShape.h.

      : m_base(&base), m_scaling(scaling), m_shift(shift) { }

ScaledShiftedShape() [2/2]

LArSamples::ScaledShiftedShape::ScaledShiftedShape ( const ScaledShiftedShape & other )

inline

Definition at line 25 of file ScaledShiftedShape.h.

      : AbsShape(), m_base(&other.base()), m_scaling(other.scaling()), m_shift(other.shift()) { }

~ScaledShiftedShape()

virtual LArSamples::ScaledShiftedShape::~ScaledShiftedShape ( )

inlinevirtual

Definition at line 28 of file ScaledShiftedShape.h.

{ }

Member Function Documentation

base()

const AbsShape & LArSamples::ScaledShiftedShape::base ( ) const

inline

Definition at line 30 of file ScaledShiftedShape.h.

{ return *m_base; }

covariance()

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

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 38 of file ScaledShiftedShape.h.

{ return base().covariance(i, j)*scaling()*scaling(); }

covarianceMatrix() [1/2]

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

inherited

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

inherited

Definition at line 223 of file AbsShape.cxx.

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

error()

double LArSamples::ScaledShiftedShape::error ( unsigned int i ) const

inlinevirtual

Reimplemented from LArSamples::AbsShape.

Definition at line 36 of file ScaledShiftedShape.h.

{ return scaling()*base().error(i); }

findTimeInterval()

int AbsShape::findTimeInterval ( double time ) const

inherited

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

inherited

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

inherited

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

inherited

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

inherited

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

inherited

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

inherited

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

inherited

Definition at line 230 of file AbsShape.cxx.

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

maxPosition()

int AbsShape::maxPosition ( ) const

inherited

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

inherited

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

inherited

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

inherited

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()

unsigned int LArSamples::ScaledShiftedShape::nPoints ( ) const

inlinevirtual

Returns

data points

Implements LArSamples::AbsShape.

Definition at line 34 of file ScaledShiftedShape.h.

{ return base().nPoints(); }  

resample()

SimpleShape * AbsShape::resample ( unsigned int nPts ) const

inherited

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);
}

scaling()

double LArSamples::ScaledShiftedShape::scaling ( ) const

inline

Definition at line 31 of file ScaledShiftedShape.h.

{ return m_scaling; } 

shift()

double LArSamples::ScaledShiftedShape::shift ( ) const

inline

Definition at line 32 of file ScaledShiftedShape.h.

{ return m_shift; }

time()

double LArSamples::ScaledShiftedShape::time ( unsigned int i ) const

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 37 of file ScaledShiftedShape.h.

{ return base().time(i) + shift(); }

value()

double LArSamples::ScaledShiftedShape::value ( unsigned int i ) const

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 35 of file ScaledShiftedShape.h.

{ return scaling()*base().value(i); }

values()

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

inherited

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;
}

Member Data Documentation

m_base

const AbsShape* LArSamples::ScaledShiftedShape::m_base

private

Definition at line 42 of file ScaledShiftedShape.h.

m_scaling

double LArSamples::ScaledShiftedShape::m_scaling

private

Definition at line 43 of file ScaledShiftedShape.h.

m_shift

double LArSamples::ScaledShiftedShape::m_shift

private

Definition at line 43 of file ScaledShiftedShape.h.

---

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

• ScaledShiftedShape.h