LArSamples::Residuals Class Reference

#include <Residual.h>

Inheritance diagram for LArSamples::Residuals:

Collaboration diagram for LArSamples::Residuals:

Public Member Functions
	Residuals ()
	Residuals (const std::vector< Residual > &r)
virtual	~Residuals ()
unsigned int	size () const
const Residual *	residual (unsigned int i) const
int	lwb () const
int	upb () const
bool	add (const Residual &residual)
bool	medianVars (TVectorD &medians, TVectorD &widths) const
std::unique_ptr< Residuals >	truncate (double nWidthsRes, double nWidthsTime=-1, unsigned int nMax=0) const
ResidualCalculator *	calculator (bool weigh=false) const
TH1D *	histogram (short sample, const TString &name, int nBins, double xMin, double xMax) const
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 Attributes
std::vector< Residual >	m_residuals

Detailed Description

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

Constructor & Destructor Documentation

Residuals() [1/2]

LArSamples::Residuals::Residuals ( )

inline

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

{ }

Residuals() [2/2]

LArSamples::Residuals::Residuals ( const std::vector< Residual > & r )

inline

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

          : m_residuals(r) { }

~Residuals()

virtual LArSamples::Residuals::~Residuals ( )

inlinevirtual

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

{ }

Member Function Documentation

add()

bool LArSamples::Residuals::add ( const Residual & residual )

inline

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

{ m_residuals.push_back(residual); return true; }

calculator()

ResidualCalculator * Residuals::calculator

(

bool

weigh = false

)

const

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

{           
  if (size() == 0) return nullptr;
 
  ResidualCalculator* calc = new ResidualCalculator(lwb(), upb(), weigh);
  for (const Residual& residual : m_residuals)
    calc->fill(residual);
  
  return calc;
}

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

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

histogram()

TH1D * Residuals::histogram	(	short	sample,
		const TString &	name,
		int	nBins,
		double	xMin,
		double	xMax ) const

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

{
  TH1D* h = new TH1D(name, "", nBins, xMin, xMax);
  for (const Residual& residual : m_residuals)
    h->Fill(sample <= upb() ? residual.scaledDelta(sample) : residual.time());
  return h;
}

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

inlinevirtual

Implements LArSamples::IndexRange.

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

{ return (size() > 0 ? residual(0)->lwb() : 0); }

medianVars()

bool Residuals::medianVars	(	TVectorD &	medians,
		TVectorD &	widths ) const

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

{
  if (size() == 0) return false;
  medians.ResizeTo(lwb(), upb() + 1);
  widths.ResizeTo (lwb(), upb() + 1);
  
  double halfSigmaQuantile = TMath::Prob(1,1)/2;
  double medianQuantile = 0.5;
  
  std::vector<Residual> sortedResiduals = m_residuals;
    
  for (short i = lwb(); i <= upb() + 1; i++) {
//     cout << "test-sorting sample " << i << " (n = " << m_residuals.size() << ")" << endl;
//     std::vector<const Residual*> testSort(m_residuals.size());
//     ResidualCompare comparer(i);
//     for (unsigned int i1 = 0; i1 < m_residuals.size(); i1++) {
//       const Residual& res1 = m_residuals[i1];
//       unsigned int nBefore = 0;
//       for (unsigned int i2 = 0; i2 < m_residuals.size(); i2++) {
//         if (i1 == i2) continue;
//         const Residual& res2 = m_residuals[i2];
//         //cout << "comparing " << i1 << " " << i2 << endl;
//         if (comparer(res2, res1)) nBefore++;
//         if ((i1 == 3617 || i1 == 3886) && comparer(res2, res1)) 
//           cout << i2 << " is bfore " << i1 << " " << res2.scaledDelta(i) << " < " << res1.scaledDelta(i) 
//                << (res2.scaledDelta(i) < res1.scaledDelta(i)) << " " << (res2.scaledDelta(i) < res1.scaledDelta(i) - 1E-3) << endl;
//       }
//       cout << "index " << i1 << " : nBefore = " << nBefore << " " << res1.scaledDelta(i) << endl;
//       while (testSort[nBefore]) { nBefore++; if (nBefore == testSort.size()) cout << "ERROR!!! " << i1  << endl; }
//       testSort[nBefore] = &res1;
//     }
//     cout << "done sorting" << endl;
//     for (unsigned int i1 = 0; i1 < testSort.size(); i1++) {
//       const Residual* res1 = testSort[i1];
//       if (!res1) cout << "ERROR : null at index " << i1 << endl;
//       for (unsigned int i2 = 0; i2 < i1; i2++) {
//         const Residual* res2 = testSort[i2];
//         if (comparer(*res1, *res2)) cout << "Wrong order " << i1 << " " << i2 << endl;
//       }
//     }
//     cout << "done checking" << endl;
    std::sort(sortedResiduals.begin(), sortedResiduals.end(), ResidualCompare(i));
    Residual& medianRes      = sortedResiduals[(unsigned int)(sortedResiduals.size()*medianQuantile)];
    Residual& loHalfSigmaRes = sortedResiduals[(unsigned int)(sortedResiduals.size()*halfSigmaQuantile)];
    Residual& hiHalfSigmaRes = sortedResiduals[(unsigned int)(sortedResiduals.size()*(1 - halfSigmaQuantile))];
    double median      = (i < upb() + 1 ? medianRes.scaledDelta(i) : medianRes.time()) ;
    double loHalfSigma = (i < upb() + 1 ? loHalfSigmaRes.scaledDelta(i) : loHalfSigmaRes.time());
    double hiHalfSigma = (i < upb() + 1 ? hiHalfSigmaRes.scaledDelta(i) : hiHalfSigmaRes.time());
    //cout << loHalfSigma << " " << median << " " << hiHalfSigma << endl;
    medians[i] = median;
    widths[i]  = (hiHalfSigma - median > median - loHalfSigma ? median - loHalfSigma : hiHalfSigma - median);
  }
  return true;
}

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

residual()

const Residual * LArSamples::Residuals::residual ( unsigned int i ) const

inline

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

{ return (i < size() ? &m_residuals[i] : 0); } 

size()

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

inline

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

{ return m_residuals.size(); }

truncate()

std::unique_ptr< Residuals > Residuals::truncate	(	double	nWidthsRes,
		double	nWidthsTime = -1,
		unsigned int	nMax = 0 ) const

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

{
  if (size() == 0) return std::make_unique<Residuals>();
  TVectorD medians, widths;
 
  if (nMax > 0) {
    auto original = std::make_unique<Residuals>();
    for (unsigned int i = 0; i < nMax; i++) original->add(*residual(i));
    if (!original->medianVars(medians, widths)) return nullptr;
  }
  else {
    if (!medianVars(medians, widths)) return nullptr;
  }
  auto truncated = std::make_unique<Residuals>();
 
  for (const Residual& residual : m_residuals) {
    bool pass = true;
    if (nMax > 0 && truncated->size() == nMax) break;
    for (short i = lwb(); i <= upb(); i++) {
      if (nWidthsRes > 0 && TMath::Abs(residual.scaledDelta(i) - medians[i]) > nWidthsRes*widths[i]) {
        pass = false;
        break;
      }
    }
    if (!pass) continue;
    if (nWidthsTime > 0 && TMath::Abs(residual.time() - medians[upb() + 1]) > nWidthsTime*widths[upb() + 1]) continue;
    truncated->add(residual);
  }
  return truncated;
}

upb()

int LArSamples::Residuals::upb ( ) const

inlinevirtual

Implements LArSamples::IndexRange.

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

{ return (size() > 0 ? residual(0)->upb() : 0); }

Member Data Documentation

m_residuals

std::vector<Residual> LArSamples::Residuals::m_residuals

private

Definition at line 96 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