LogChiSquareComparator.h

Go to the documentation of this file.

 
#ifndef MODELFITTING_LOGCHISQUARECOMPARATOR_H
#define MODELFITTING_LOGCHISQUARECOMPARATOR_H
 
#include <cmath>
 
namespace ModelFitting {
 
class LogChiSquareComparator {
 
public:
 
  LogChiSquareComparator(double u0=10.) : m_u0{u0} { }
 
  double operator()(double real, double model, double weight) const {
    double val =  weight * (real - model) / m_u0;
    return val>0. ? m_u0 * std::log1p(val) : -1. * m_u0 * std::log1p(-val);
  }
 
private:
 
  double m_u0;
 
}; // end of class LogChiSquareComparator
 
} // end of namespace ModelFitting
 
 
 
#endif  /* MODELFITTING_LOGCHISQUARECOMPARATOR_H */