PXDClusterPositionEstimator.h

/**************************************************************************
* BASF2 (Belle Analysis Framework 2)                                     *
* Copyright(C) 2010 - Belle II Collaboration                             *
*                                                                        *
* Author: The Belle II Collaboration                                     *
* Contributors: Benjamin Schwenker                                       *
*                                                                        *
* This software is provided "as is" without any warranty.                *
**************************************************************************/
 
#pragma once
 
#include <pxd/dbobjects/PXDClusterPositionEstimatorPar.h>
#include <pxd/dbobjects/PXDClusterShapeIndexPar.h>
#include <pxd/dataobjects/PXDCluster.h>
#include <vxd/dataobjects/VxdID.h>
#include <framework/database/DBObjPtr.h>
#include <set>
#include <vector>
#include <pxd/reconstruction/Pixel.h>
#include <memory>
 
namespace Belle2 {
  namespace PXD {
    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();
 
      float computeEta(const std::set<Pixel>& pixels, int vStart, int vSize, double thetaU, double thetaV) const;
 
      int computeShapeIndex(const std::set<Pixel>& pixels, int uStart, int vStart, int vSize, double thetaU,
                            double thetaV) const;
 
      const std::string getShortName(const std::set<Pixel>& pixels, int uStart, int vStart, int vSize, double thetaU,
                                     double thetaV) const;
 
      const std::string getMirroredShortName(const std::set<Pixel>& pixels, int uStart, int vStart, int vSize, double thetaU,
                                             double thetaV) const;
 
      const std::string getFullName(const std::set<Pixel>& pixels, int uStart, int vStart) const;
 
      int getClusterkind(const PXDCluster& cluster) const;
 
      int getClusterkind(const std::vector<Pixel>& pixels, const VxdID& sensorID) const;
 
      int getSectorIndex(double thetaU, double thetaV) const;
 
    private:
 
      const Pixel& getHeadPixel(const std::set<Pixel>& pixels, int vStart, int vSize, double thetaU, double thetaV) const;
 
      const Pixel& getTailPixel(const std::set<Pixel>& pixels, int vStart, int vSize, double thetaU, double thetaV) const;
 
      const Pixel& getLastPixelWithVOffset(const std::set<Pixel>& pixels, int vStart, int vOffset) const;
 
      const Pixel& getFirstPixelWithVOffset(const std::set<Pixel>& pixels, int vStart, int vOffset) const;
 
      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