PXDClusterPositionEstimator.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <pxd/dbobjects/PXDClusterPositionEstimatorPar.h>
#include <pxd/dbobjects/PXDClusterShapeIndexPar.h>
#include <framework/database/DBObjPtr.h>
#include <set>
#include <vector>
#include <memory>
 
namespace Belle2 {
  class PXDCluster;
  class VxdID;
 
  namespace PXD {
    class Pixel;

    class PXDClusterPositionEstimator {
 
    public:

      void initialize();

      void setShapeIndexFromDB();

      void setPositionEstimatorFromDB();

      const PXDClusterPositionEstimatorPar& getPositionEstimatorParameters() const {return  m_positionEstimatorPar;}

      const PXDClusterShapeIndexPar& getShapeIndexParameters() const {return  m_shapeIndexPar;}

      const PXDClusterOffsetPar* getClusterOffset(const PXDCluster& cluster, double tu, double tv) const;

      float getShapeLikelyhood(const PXDCluster& cluster, double tu, double tv) const;

      static PXDClusterPositionEstimator& getInstance();

      static float computeEta(const std::set<Pixel>& pixels, int vStart, int vSize, double thetaU, double thetaV);

      int computeShapeIndex(const std::set<Pixel>& pixels, int uStart, int vStart, int vSize, double thetaU,
                            double thetaV) const;

      static const std::string getShortName(const std::set<Pixel>& pixels, int uStart, int vStart, int vSize, double thetaU,
                                            double thetaV);

      static const std::string getMirroredShortName(const std::set<Pixel>& pixels, int uStart, int vStart, int vSize, double thetaU,
                                                    double thetaV);

      static const std::string getFullName(const std::set<Pixel>& pixels, int uStart, int vStart);

      int getClusterkind(const PXDCluster& cluster) const;

      int getClusterkind(const std::vector<Belle2::PXD::Pixel>& pixels, const Belle2::VxdID& sensorID) const;

      static int getSectorIndex(double thetaU, double thetaV);
 
    private:

      static const Pixel& getHeadPixel(const std::set<Pixel>& pixels, int vStart, int vSize, double thetaU, double thetaV);

      static const Pixel& getTailPixel(const std::set<Pixel>& pixels, int vStart, int vSize, double thetaU, double thetaV);

      static const Pixel& getLastPixelWithVOffset(const std::set<Pixel>& pixels, int vStart, int vOffset);

      static const Pixel& getFirstPixelWithVOffset(const std::set<Pixel>& pixels, int vStart, int vOffset);

      PXDClusterPositionEstimator() {};
      PXDClusterPositionEstimator(const PXDClusterPositionEstimator&) = delete;
      PXDClusterPositionEstimator& operator=(const PXDClusterPositionEstimator&) = delete;

      std::unique_ptr<DBObjPtr<PXDClusterShapeIndexPar>> m_shapeIndexFromDB;
      std::unique_ptr<DBObjPtr<PXDClusterPositionEstimatorPar>> m_positionEstimatorFromDB;

      PXDClusterShapeIndexPar m_shapeIndexPar;
      PXDClusterPositionEstimatorPar m_positionEstimatorPar;
    };
  }
} //Belle2 namespace