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

63 : 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:

Public Member Functions