LArSamples::ResidualCompare Class Reference

#include <Residual.h>

Collaboration diagram for LArSamples::ResidualCompare:

Public Member Functions
	ResidualCompare (int sampling)
bool	operator() (const Residual &r1, const Residual &r2) const

Private Attributes
int	m_sampling

Detailed Description

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

Constructor & Destructor Documentation

ResidualCompare()

LArSamples::ResidualCompare::ResidualCompare ( int sampling )

inline

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

: m_sampling(sampling) { }

Member Function Documentation

operator()()

bool ResidualCompare::operator()	(	const Residual &	r1,
		const Residual &	r2 ) const

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

{
  if (!r1.hasSameRange(r2)) {
    cout << "WARNING : residuals have different ranges, this may crash the sorting..." << endl;
    return false;
  }
 
  // Very, very important! Since samples are integers, it *will* happen that some residuals will have *identical* deltas. 
  // It turns out that (at least in the way used here), operator<(double, double) is not a strict weak ordering: we will
  // have cases where residual1[4] == residual2[4], and then we will have residual1[4] < residual2[4] and residual2[4] < residual1[4]
  // both true, which will confuse std::sort and cause it to crash... so add this safety (lesson of 2 days of debugging)
  // Note: if we write a<<b iff a < b - epsilon, then << is indeed a strict weak ordering so all is fine. It does mean there
  // will be "sorting errors", i.e. residuals that differ by less than epsilon may be sorted the wrong way. But that's fine for our usage.
  const double epsilon = 1E-5;
//   cout << "Comparing sample " << m_sampling << ", = " 
//        << (m_sampling == r1.upb() + 1 ? r1.time() : r1.scaledDelta(m_sampling)) << " >/< "
//        << (m_sampling == r2.upb() + 1 ? r2.time() : r2.scaledDelta(m_sampling)) << " " 
//        << ((m_sampling == r1.upb() + 1 ? r1.time() : r1.scaledDelta(m_sampling)) < (m_sampling == r2.upb() + 1 ? r2.time() : r2.scaledDelta(m_sampling))) << endl;
 
  if (r1.isInRange(m_sampling)) return r1.scaledDelta(m_sampling) < r2.scaledDelta(m_sampling) - epsilon;
  if (m_sampling == r1.upb() + 1) return r1.time() < r2.time() - epsilon;
  cout << "WARNING : comparison on an unknown sample, this may crash the sorting..." << endl;
  return false;
}

Member Data Documentation

m_sampling

int LArSamples::ResidualCompare::m_sampling

private

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