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

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.

104 : 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.

107 : 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.

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

◆ ~ResidualCalculator()

virtual LArSamples::ResidualCalculator::~ResidualCalculator ( )

inlinevirtual

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

112 { }

Member Function Documentation

◆ add()

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

Definition at line 205 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 277 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.

39 { 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.

40 { return (other.upb() < upb() ? other.upb() : upb()); }

◆ description()

TString ResidualCalculator::description ( ) const

Definition at line 288 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.

132 { return m_events[i]; }

◆ fill()

bool ResidualCalculator::fill ( const Residual & residual )

Definition at line 213 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 270 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 197 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.

30 { 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.

31 { 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.

29 { return (lwb() == lw && upb() == up); }

◆ isInRange()

bool LArSamples::IndexRange::isInRange ( int i ) const

inlineinherited

Definition at line 27 of file IndexRange.h.

27 { 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.

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

◆ providesRange() [1/3]

bool LArSamples::IndexRange::providesRange ( const IndexRange & other ) const

inlineinherited

Definition at line 34 of file IndexRange.h.

34 { 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.

35 { 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.

33 { return (lwb() <= lw && upb() >= up); }

◆ rangeStr()

TString LArSamples::IndexRange::rangeStr ( ) const

inlineinherited

Definition at line 37 of file IndexRange.h.

37 { 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.

119 { return &m_regresser; }

◆ remove()

bool ResidualCalculator::remove ( const Residual & residual )

Definition at line 221 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.

131 { return m_runs[i]; }

◆ shapeErrorData()

ShapeErrorData * ResidualCalculator::shapeErrorData ( ) const

Definition at line 236 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 new 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 new 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.

117 { 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.

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

◆ weigh()

bool LArSamples::ResidualCalculator::weigh ( ) const

inline

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

120 { return m_weigh; }

◆ weight()

double ResidualCalculator::weight ( const Residual & residual ) const

Definition at line 263 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:

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ResidualCalculator() [1/3]

◆ ResidualCalculator() [2/3]

◆ ResidualCalculator() [3/3]

◆ ~ResidualCalculator()

Member Function Documentation

◆ add()

◆ append()

◆ checkRange()

◆ commonLwb()

◆ commonUpb()

◆ description()

◆ event()

◆ fill()

◆ fill_with_weight()

◆ find()

◆ hasSameRange() [1/3]

◆ hasSameRange() [2/3]

◆ hasSameRange() [3/3]

◆ isInRange()

◆ lwb()

◆ providesRange() [1/3]

◆ providesRange() [2/3]

◆ providesRange() [3/3]

◆ rangeStr()

◆ regresser()

◆ remove()

◆ run()

◆ shapeErrorData()

◆ size()

◆ upb()

◆ weigh()

◆ weight()

Member Data Documentation

◆ m_events

◆ m_regresser

◆ m_runs

◆ m_weigh