AdaptiveResidualSmoothing.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// 29.07.2008, AUTHOR: OLIVER KORTNER
// Modified: 20.03.2009 by O. Kortner, additional method ,,performSmoothing''
//                      with better performance than the old method added;
//                      softer chi^2 cuts in addResidualsFromSegment.
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
#ifndef MuonCalib_AdaptiveResidualSmoothingH
#define MuonCalib_AdaptiveResidualSmoothingH
 
//:::::::::::::::::::::::::::::::::::::
//:: CLASS AdaptiveResidualSmoothing ::
//:::::::::::::::::::::::::::::::::::::

 
//::::::::::::::::::
//:: HEADER FILES ::
//::::::::::::::::::
 
// STL //
#include <vector>
 
// MuonCalib //
#include "MdtCalibData/RtRelationLookUp.h"
#include "MdtCalibFitters/CurvedPatRec.h"
#include "MdtCalibFitters/StraightPatRec.h"
#include "MuonCalibMath/DataPoint.h"
 
namespace MuonCalib {
 
    class IRtRelation;
 
    class AdaptiveResidualSmoothing {
    public:
        // Constructor //
        AdaptiveResidualSmoothing();
 
        // Methods //
        void clear();
        void addResidual(const double  radius, const double  residual);
        bool addResidualsFromSegment(MuonCalibSegment& seg, bool curved, double road_width);
        RtRelationLookUp performSmoothing(const IRtRelation& rt_rel, unsigned int nb_entries_per_bin, bool fix_t_min, bool fix_t_max);
        RtRelationLookUp performSmoothing(const IRtRelation& rt_rel, const bool& fix_t_min, const bool& fix_t_max);
 
    private:
        std::vector<DataPoint> m_residual_point;  // vector of residual points
        StraightPatRec m_sfitter;                 // straight-line fitter
        CurvedPatRec m_cfitter;                   // curved-line fitter
        double t_from_r(const IRtRelation& rt_rel, const double  r);
        // get t(r) for the given r-t relationship,
        // the method is auxiliary and not optimized;
        // it will disappear when the t(r) will be
        // available in the MuonCalib framework
    };
 
}  // namespace MuonCalib
 
#endif