LArSamples::SimpleShape Class Reference

#include <SimpleShape.h>

Inheritance diagram for LArSamples::SimpleShape:

Collaboration diagram for LArSamples::SimpleShape:

Public Member Functions
	SimpleShape (const std::vector< double > &values, const std::vector< double > &errors, double timeInterval=25, double startTime=0)
	Constructor.
	SimpleShape (unsigned int nPoints, double timeInterval=25, double startTime=0)
	SimpleShape (const SimpleShape &other)
	SimpleShape (const ShapeInfo &shapeInfo, double scale=1, double shift=0, bool samplingTimeOnly=false)
	SimpleShape (const AbsShape &other, double scale=1, double shift=0)
virtual	~SimpleShape ()
double	timeInterval () const
double	startTime () const
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
void	set (unsigned int i, double value, double error=-1)
void	setError (unsigned int i, double error)
SimpleShape *	diff () const
bool	add (unsigned int k, double value, double error)
SimpleShape *	add (const AbsShape &other, double scale=1, double shift=0)
SimpleShape *	createEmpty () const
TH1D *	histogram (const char *name="shape", const char *title="", bool timeInUnitOfSamples=false) 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

Static Public Member Functions
static bool	add (std::unique_ptr< SimpleShape > &s1, const AbsShape &s2)
static bool	scaleAndShift (std::unique_ptr< SimpleShape > &s1, double scale, double shift=0)

Private Member Functions
SimpleShape &	operator= (const SimpleShape &)

Private Attributes
std::vector< double >	m_values
std::vector< double >	m_errors
double	m_timeInterval {}
double	m_startTime {}

Detailed Description

Definition at line 25 of file SimpleShape.h.

Constructor & Destructor Documentation

SimpleShape() [1/5]

SimpleShape::SimpleShape	(	const std::vector< double > &	values,
		const std::vector< double > &	errors,
		double	timeInterval = 25,
		double	startTime = 0 )

Constructor.

Definition at line 15 of file SimpleShape.cxx.

  : m_values(values), m_errors(errors), 
    m_timeInterval(timeInterval), m_startTime(startTime)
{
}

SimpleShape() [2/5]

SimpleShape::SimpleShape	(	unsigned int	nPoints,
		double	timeInterval = 25,
		double	startTime = 0 )

Definition at line 23 of file SimpleShape.cxx.

  : m_values(nPoints), m_errors(nPoints), 
    m_timeInterval(timeInterval), m_startTime(startTime)
{
}

SimpleShape() [3/5]

LArSamples::SimpleShape::SimpleShape ( const SimpleShape & other )

inline

Definition at line 35 of file SimpleShape.h.

                                          :
      AbsShape(),
      m_values(other.m_values), m_errors(other.m_errors), 
      m_timeInterval(other.m_timeInterval), m_startTime(other.m_startTime) { }

SimpleShape() [4/5]

SimpleShape::SimpleShape	(	const ShapeInfo &	shapeInfo,
		double	scale = 1,
		double	shift = 0,
		bool	samplingTimeOnly = false )

Definition at line 42 of file SimpleShape.cxx.

{
  for (unsigned int i = 0; i < shapeInfo.nPoints(); i++) {
    if (samplingTimeOnly && shapeInfo.phase(i) != 0) continue;
    m_values.push_back(shapeInfo.value(i)*scale);
    m_errors.push_back(0);
  }
 
  double timeInterval = 1.0*Definitions::samplingInterval/(samplingTimeOnly ? 1 : shapeInfo.nIntervals());
  m_timeInterval = timeInterval;
  m_startTime = shift - shapeInfo.shift();
}

SimpleShape() [5/5]

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

Definition at line 30 of file SimpleShape.cxx.

{
  for (unsigned int i = 0; i < shape.nPoints(); i++) {
    m_values.push_back(shape.value(i)*scale);
    m_errors.push_back(shape.error(i)*scale);
  }
  
  m_timeInterval = (shape.time(shape.nPoints() - 1) - shape.time(0))/(shape.nPoints() - 1);
  m_startTime = shape.time(0) + shift;
}

~SimpleShape()

virtual LArSamples::SimpleShape::~SimpleShape ( )

inlinevirtual

Definition at line 43 of file SimpleShape.h.

{ }

Member Function Documentation

add() [1/3]

SimpleShape * SimpleShape::add	(	const AbsShape &	other,
		double	scale = 1,
		double	shift = 0 )

Definition at line 74 of file SimpleShape.cxx.

{
  std::vector<double> values, errors;
  for (unsigned int k = 0; k < nPoints(); ++k) {
    double val, err;
    int pb = other.interpolate(time(k) - shift, val, err);
    if (pb != 0) val = err = 0;
    //if (pb == 1) break; //overflow: stop. underflow is OK : shape is assumed to be 0
    values.push_back(this->value(k) + scale*val);
    errors.push_back(TMath::Sqrt(TMath::Power(this->error(k),2) + TMath::Power(scale*err, 2)));
  }
 
  return new SimpleShape(values, errors, timeInterval(), startTime());
}

add() [2/3]

bool SimpleShape::add ( std::unique_ptr< SimpleShape > & s1, const AbsShape & s2 )

static

Definition at line 124 of file SimpleShape.cxx.

{
  if (!s1) {
    s1.reset(new SimpleShape(s2));
    return true;
  }
  SimpleShape* sum = s1->add(s2);
  if (!sum) return false;
  s1.reset(sum);
  return true;
}

add() [3/3]

bool SimpleShape::add	(	unsigned int	k,
		double	value,
		double	error )

Definition at line 90 of file SimpleShape.cxx.

{
  m_values[k] = m_values[k] + val;
  m_errors [k] = TMath::Sqrt(m_errors[k]*m_errors[k] + err*err);
  return true;
}

covariance()

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

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 51 of file SimpleShape.h.

{ return (i == j ? m_errors[i]*m_errors[i] : 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);
}

createEmpty()

SimpleShape * SimpleShape::createEmpty

(

)

const

Definition at line 98 of file SimpleShape.cxx.

{
  return new SimpleShape(nPoints(), timeInterval(), startTime());
}

diff()

SimpleShape * SimpleShape::diff

(

)

const

Definition at line 62 of file SimpleShape.cxx.

{
  SimpleShape* diff = new SimpleShape(nPoints(), timeInterval(), startTime());
  for (unsigned int i = 0; i < nPoints(); i++) {
    double val;
    interpolateDiff(time(i), val);
    diff->set(i, val, 0);
  }
  return diff;
}

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

histogram()

TH1D * SimpleShape::histogram	(	const char *	name = "shape",
		const char *	title = "",
		bool	timeInUnitOfSamples = false ) const

Definition at line 104 of file SimpleShape.cxx.

{
  double xMin = (timeInUnitOfSamples ? -1.5 : time(0) - 1.5*timeInterval());
  double xMax = (timeInUnitOfSamples ? nPoints() + 0.5 : time(nPoints()) + 0.5*timeInterval());
  TH1D* h = new TH1D(name, title, nPoints() + 2, xMin, xMax);
  h->GetXaxis()->SetTitle(timeInUnitOfSamples ? "Sample Index" : "Time (ns)");
  h->GetYaxis()->SetTitle("ADC counts");
  h->GetYaxis()->SetTitleOffset(1.1);
  h->SetMarkerStyle(20);
  h->SetMarkerSize(1);
  
  for (unsigned int i = 0; i < nPoints(); i++) {
    h->SetBinContent(i + 2, value(i));
    if (error(i) > 0) h->SetBinError(i + 2, error(i));
  }
 
  return h;
}

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

inlinevirtual

Returns

data points

Implements LArSamples::AbsShape.

Definition at line 48 of file SimpleShape.h.

{ return m_values.size(); }  

operator=()

SimpleShape & LArSamples::SimpleShape::operator= ( const SimpleShape & )

private

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

scaleAndShift()

bool SimpleShape::scaleAndShift ( std::unique_ptr< SimpleShape > & s1, double scale, double shift = 0 )

static

Definition at line 137 of file SimpleShape.cxx.

{
  if (!s1) return false;
  s1.reset(new SimpleShape(*s1, scale, shift));
  return true;
}

set()

void LArSamples::SimpleShape::set ( unsigned int i, double value, double error = -1 )

inline

Definition at line 53 of file SimpleShape.h.

{ m_values[i] = value; if (error > 0) m_errors[i] = error; }

setError()

void LArSamples::SimpleShape::setError ( unsigned int i, double error )

inline

Definition at line 54 of file SimpleShape.h.

{ m_errors[i] = error; }

startTime()

double LArSamples::SimpleShape::startTime ( ) const

inline

Definition at line 46 of file SimpleShape.h.

{ return m_startTime; }

time()

double SimpleShape::time ( unsigned int i ) const

virtual

Implements LArSamples::AbsShape.

Definition at line 56 of file SimpleShape.cxx.

{ 
  return startTime() + i*timeInterval();
}

timeInterval()

double LArSamples::SimpleShape::timeInterval ( ) const

inline

Definition at line 45 of file SimpleShape.h.

{ return m_timeInterval; } 

value()

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

inlinevirtual

Implements LArSamples::AbsShape.

Definition at line 49 of file SimpleShape.h.

{ return (i < m_values.size() ? m_values[i] : -1E99); }

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_errors

std::vector<double> LArSamples::SimpleShape::m_errors

private

Definition at line 68 of file SimpleShape.h.

m_startTime

double LArSamples::SimpleShape::m_startTime {}

private

Definition at line 69 of file SimpleShape.h.

{}, m_startTime{};

m_timeInterval

double LArSamples::SimpleShape::m_timeInterval {}

private

Definition at line 69 of file SimpleShape.h.

{}, m_startTime{};

m_values

std::vector<double> LArSamples::SimpleShape::m_values

private

Definition at line 68 of file SimpleShape.h.

---

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

• SimpleShape.h
• SimpleShape.cxx