LArSamples::ResidualCalculator Class Reference

#include <Residual.h>

Inheritance diagram for LArSamples::ResidualCalculator:

Collaboration diagram for LArSamples::ResidualCalculator:

Public Member Functions
	ResidualCalculator ()
	ResidualCalculator (unsigned int l, unsigned int u, bool weigh=false)
	ResidualCalculator (const Averager &regresser, bool weigh=false)
virtual	~ResidualCalculator ()
int	lwb () const
int	upb () const
unsigned int	size () const
const Averager *	regresser () const
bool	weigh () const
double	weight (const Residual &residual) const
bool	fill (const Residual &residual)
bool	remove (const Residual &residual)
bool	append (const ResidualCalculator &other)
std::unique_ptr< ShapeErrorData >	shapeErrorData () const
TString	description () const
int	run (unsigned int i) const
int	event (unsigned int i) const
int	find (int run, int event) const
bool	add (int run, int event)
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	fill_with_weight (const Residual &residual, double w)

Private Attributes
Averager	m_regresser
std::vector< int >	m_runs
std::vector< int >	m_events
bool	m_weigh

Detailed Description

Definition at line 100 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

Constructor & Destructor Documentation

ResidualCalculator() [1/3]

LArSamples::ResidualCalculator::ResidualCalculator ( )

inline

Definition at line 104 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

: m_weigh(false) { }

ResidualCalculator() [2/3]

LArSamples::ResidualCalculator::ResidualCalculator ( unsigned int l, unsigned int u, bool weigh = false )

inline

Definition at line 106 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

        : m_regresser(l, u + 1), m_weigh(weigh) { }

ResidualCalculator() [3/3]

LArSamples::ResidualCalculator::ResidualCalculator ( const Averager & regresser, bool weigh = false )

inline

Definition at line 109 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

        : m_regresser(regresser), m_weigh(weigh) { }

~ResidualCalculator()

virtual LArSamples::ResidualCalculator::~ResidualCalculator ( )

inlinevirtual

Definition at line 112 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ }

Member Function Documentation

add()

bool ResidualCalculator::add	(	int	run,
		int	event )

Definition at line 204 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{
  m_runs.push_back(r);
  m_events.push_back(e);
  return true;
}

append()

bool ResidualCalculator::append

(

const ResidualCalculator &

other

)

Definition at line 276 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{
  if (!m_regresser.append(*other.regresser())) return false;
  //for (std::map< std::pair<int, int >, bool >::const_iterator event = other.events().begin();
  //     event != other.events().end(); event++)
  //  m_events[event->first] = event->second;
  for (unsigned int i = 0; i < other.size(); i++) add(other.run(i), other.event(i));
  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;
}

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

description()

TString ResidualCalculator::description

(

)

const

Definition at line 287 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{
  TString str = "\n";
  /*for (std::map< std::pair<int, int>, bool>::const_iterator ev = m_events.begin();
       ev != m_events.end(); ev++)
   str += Form("run %6d event %10d\n", ev->first.first, ev->first.second);
  */
  for (unsigned int i = 0; i < size(); i++)
    str += Form("run %6d event %10d\n", run(i), event(i));
 
  return str;
}

event()

int LArSamples::ResidualCalculator::event ( unsigned int i ) const

inline

Definition at line 132 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ return m_events[i]; }

fill()

bool ResidualCalculator::fill

(

const Residual &

residual

)

Definition at line 212 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{ 
  if (find(residual.run(), residual.event()) >= 0) return true;
  add(residual.run(), residual.event());
  return fill_with_weight(residual, +1);
} 

fill_with_weight()

bool ResidualCalculator::fill_with_weight ( const Residual & residual, double w )

private

Definition at line 269 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{
  w *= weight(residual);
  return m_regresser.fill(residual.scaledDeltasAndTime(), w);
}

find()

int ResidualCalculator::find	(	int	run,
		int	event ) const

Definition at line 196 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{
  for (unsigned int i = 0; i < size(); i++)
    if (event(i) == e && run(i) == r) return i;
  return -1;
}

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::ResidualCalculator::lwb ( ) const

inlinevirtual

Implements LArSamples::IndexRange.

Definition at line 114 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ return regresser()->lwb(); }

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

regresser()

const Averager * LArSamples::ResidualCalculator::regresser ( ) const

inline

Definition at line 119 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ return &m_regresser; }

remove()

bool ResidualCalculator::remove

(

const Residual &

residual

)

Definition at line 220 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{ 
  //std::pair<int, int> ev(residual.run(), residual.event());
  //std::map< std::pair<int, int>, bool>::iterator findEv = m_events.find(ev);
  //if (findEv == m_events.end()) return true;
  //cout << "Will remove run = " << ev.first << ", event = " << ev.second << endl;  
  //m_events.erase(findEv);
  int index = find(residual.run(), residual.event());
  if (index == -1) return true;
  m_events[index] = -1;
  m_runs[index] = -1;
  return fill_with_weight(residual, -1); 
}

run()

int LArSamples::ResidualCalculator::run ( unsigned int i ) const

inline

Definition at line 131 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ return m_runs[i]; }

shapeErrorData()

std::unique_ptr< ShapeErrorData > ResidualCalculator::shapeErrorData

(

)

const

Definition at line 235 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{
  TVectorD xi = regresser()->means().GetSub(lwb(), upb(), "I");
  CovMatrix xiErr = regresser()->meanErrorMatrix().GetSub(lwb(), upb(), lwb(), upb(), "I");
  double tbar = regresser()->mean(upb() + 1);
  
  double denom = regresser()->covarianceMatrix()(upb() + 1, upb() + 1);   
  if (denom < 1E-6) {
    TVectorD xip(lwb(), upb());
    CovMatrix xipErr(lwb(), upb());
    // happens for size==2 if we are removing one of the residuals. 
    if (size() > 2) cout << "WARNING: variance of t < 1E-6, returning correction without derivative term. (V = " << denom << ", N = " << size() << ")" << endl;
    return std::make_unique<ShapeErrorData>(xi, xip, xiErr, xipErr, tbar, regresser()->nEntries());
  }  
   
  TVectorD xip = TVectorD(regresser()->covarianceMatrix()[upb() + 1]).GetSub(lwb(), upb(), "I");
  xip *= -1/denom;  
  
  TVectorD xipErrVect  = TVectorD(regresser()->covarianceMatrixErrors()[upb() + 1]).GetSub(lwb(), upb(), "I");
  xipErrVect *= 1/denom;
  CovMatrix xipErr(lwb(), upb());
  for (int k1 = lwb(); k1 <= upb(); k1++) xipErr(k1, k1) = TMath::Power(xipErrVect(k1), 2);
 
  return std::make_unique<ShapeErrorData>(xi, xip, xiErr, xipErr, tbar, regresser()->nEntries());
}

size()

unsigned int LArSamples::ResidualCalculator::size ( ) const

inline

Definition at line 117 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ return m_events.size(); }

upb()

int LArSamples::ResidualCalculator::upb ( ) const

inlinevirtual

Implements LArSamples::IndexRange.

Definition at line 115 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ return regresser()->upb() - 1; }

weigh()

bool LArSamples::ResidualCalculator::weigh ( ) const

inline

Definition at line 120 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

{ return m_weigh; }

weight()

double ResidualCalculator::weight

(

const Residual &

residual

)

const

Definition at line 262 of file LArCalorimeter/LArSamplesMon/src/Residual.cxx.

{
  if (!weigh()) return 1;
  return TMath::Power(residual.adcMax(), 2);
}

Member Data Documentation

m_events

std::vector<int> LArSamples::ResidualCalculator::m_events

private

Definition at line 143 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

m_regresser

Averager LArSamples::ResidualCalculator::m_regresser

private

Definition at line 141 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

m_runs

std::vector<int> LArSamples::ResidualCalculator::m_runs

private

Definition at line 143 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

m_weigh

bool LArSamples::ResidualCalculator::m_weigh

private

Definition at line 144 of file LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h.

---

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

• LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h
• LArCalorimeter/LArSamplesMon/src/Residual.cxx