LArSamples::GraphShape Class Reference

#include <GraphShape.h>

Inheritance diagram for LArSamples::GraphShape:

Collaboration diagram for LArSamples::GraphShape:

Public Member Functions
	GraphShape (const std::vector< double > &times, const std::vector< double > &values, const std::vector< double > &errors)
	Constructor.
	GraphShape (const GraphShape &other)
	GraphShape (const AbsShape &other, double scale=1, double shift=0)
virtual	~GraphShape ()
bool	add (const AbsShape &other)
unsigned int	nPoints () const
double	value (unsigned int i) const
double	time (unsigned int i) const
double	covariance (unsigned int i, unsigned int j) const
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
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
std::vector< GraphPoint >	m_points

Detailed Description

Definition at line 32 of file GraphShape.h.

Constructor & Destructor Documentation

GraphShape() [1/3]

GraphShape::GraphShape	(	const std::vector< double > &	times,
		const std::vector< double > &	values,
		const std::vector< double > &	errors )

Constructor.

Definition at line 11 of file GraphShape.cxx.

{
  for (unsigned int k = 0; k < times.size(); k++)
    m_points.emplace_back(times[k], values[k], errors[k]);
}

GraphShape() [2/3]

LArSamples::GraphShape::GraphShape ( const GraphShape & other )

inline

Definition at line 39 of file GraphShape.h.

                                        :
        AbsShape(),
        m_points(other.m_points) { }

GraphShape() [3/3]

GraphShape::GraphShape	(	const AbsShape &	other,
		double	scale = 1,
		double	shift = 0 )

Definition at line 19 of file GraphShape.cxx.

{
  for (unsigned int k = 0; k < shape.nPoints(); k++)
    m_points.emplace_back(shape.time(k) + shift, shape.value(k)*scale, shape.error(k)*scale);
}

~GraphShape()

virtual LArSamples::GraphShape::~GraphShape ( )

inlinevirtual

Definition at line 45 of file GraphShape.h.

{ }

Member Function Documentation

add()

bool GraphShape::add

(

const AbsShape &

other

)

Definition at line 26 of file GraphShape.cxx.

{
  for (unsigned int k = 0; k < other.nPoints(); k++)
    m_points.emplace_back(other.time(k), other.value(k), other.error(k));
 
  std::sort(m_points.begin(), m_points.end());
  return true;
}

covariance()

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

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 53 of file GraphShape.h.

{ return (i == j ? m_points[i].error*m_points[i].error : 0); }

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 AbsShape::error ( unsigned int i ) const

virtualinherited

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

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::GraphShape::nPoints ( ) const

inlinevirtual

Returns

size

Implements LArSamples::AbsShape.

Definition at line 50 of file GraphShape.h.

{ return m_points.size(); }  

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

time()

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

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 52 of file GraphShape.h.

{ return m_points[i].time; }

value()

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

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 51 of file GraphShape.h.

{ return m_points[i].value; }

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_points

std::vector<GraphPoint> LArSamples::GraphShape::m_points

private

Definition at line 57 of file GraphShape.h.

---

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

• GraphShape.h
• GraphShape.cxx