TileMuonTrackDistance.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
/********************************************************************
 *
 * NAME:     TileMuonTrackDistance
 * PACKAGE:  offline/TileCalorimeter/TileCosmicAlgs
 *
 * AUTHOR :  J. Maneira
 * CREATED:  10-Jul-2006 
 *
 * PURPOSE:  Minuit fit function for TileMuonFitter
 *           Supplies summed distance from cells to track
 *
 *  Input: Cells
 *  Output: Track parameters
 *   
 ********************************************************************/
 
// C++ STL includes
#include <vector>
 
// external packages declarations
#include "Minuit2/FCNBase.h"
#include "RVersion.h"
#include "CLHEP/Vector/ThreeVector.h"
 
namespace ROOT {
  namespace Minuit2 {

 
    class TileMuonTrackDistance: public FCNBase {
 
      public:

        TileMuonTrackDistance(const std::vector<double>& X, const std::vector<double>& Y, const std::vector<double>& Z,
            const std::vector<double>& E)
            : m_theX(X), m_theY(Y), m_theZ(Z), m_theE(E), m_doWeighted(false), m_meanX(0.0), m_meanY(0.0), m_meanZ(0.0) {
        }
 
        virtual ~TileMuonTrackDistance() {
        }

        virtual double operator()(const std::vector<double>&) const override;
        virtual double Up() const override {
          return 1.;
        }

        void Means();
        double Distance2SinglePoint(double x, double y, double z, const std::vector<double>& par) const;
        CLHEP::Hep3Vector ClosestPoint(CLHEP::Hep3Vector *dataPoint, const std::vector<double>& par) const;

        double GetMeanX() {
          return m_meanX;
        }

        double GetMeanY() {
          return m_meanY;
        }

        double GetMeanZ() {
          return m_meanZ;
        }

        void SetWeighted(bool weight) {
          m_doWeighted = weight;
        }

        bool GetWeighted() {
          return m_doWeighted;
        }
 
      private:
        std::vector<double> m_theX;
        std::vector<double> m_theY;
        std::vector<double> m_theZ;
        std::vector<double> m_theE;

        bool m_doWeighted;
 
        double m_meanX;
        double m_meanY;
        double m_meanZ;
 
    };
 
  }
}