Linear_fit_of_Hits_t Struct Reference

Public Member Functions
template<typename Container_T >
auto	operator() (const Container_T &container) const

Detailed Description

Definition at line 102 of file klm_trig_linear_fit.cc.

Member Function Documentation

operator()()

auto operator() ( const Container_T &  container ) const

inline

Definition at line 104 of file klm_trig_linear_fit.cc.

    {
 
 
 
      int64_t sumX = 0, sumY = 0, sumXX = 0, sumXY = 0, sumYY = 0;
      for (const auto& e : container) {
        sumX += e.x_pos;
        sumXX += e.x_pos * e.x_pos;
 
        sumY += e.y_pos;
        sumYY += e.y_pos * e.y_pos;
 
        sumXY += e.x_pos * e.y_pos;
      }
      int64_t nHits = container.size();
      int denom = sumXX * nHits - sumX * sumX;
      if (denom == 0) {
        return nt::ntuple(slopeXY(1e9), interceptXY(1e9), ipXY(1e9), chisqXY(1e9), Nhits(nHits));
      }
 
      auto slopeXY_ = slopeXY((double)(sumXY * nHits - sumX * sumY) / (double)denom);
      auto  interceptXY_ = interceptXY((double)(sumXX * sumY - sumX * sumXY) / (double)denom);
 
      auto  ipXY_ = ipXY(interceptXY_ * interceptXY_ * (1.0 - slopeXY_ * slopeXY_));
 
      auto  chisqXY_ = chisqXY(slopeXY_ * slopeXY_ * sumXX
                               + interceptXY_ * interceptXY_ * nHits
                               + sumYY
                               + 2.0 * slopeXY_ * interceptXY_ * sumX
                               - 2.0 * slopeXY_ * sumXY
                               - 2.0 * interceptXY_ * sumY);
 
      return  nt::ntuple(slopeXY_, interceptXY_, ipXY_, chisqXY_, Nhits(nHits));
    }

---

The documentation for this struct was generated from the following file:

• trg/klm/modules/klmtrigger/src/klm_trig_linear_fit.cc