SensorPlane.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Peter Kvasnicka, Martin Ritter                           *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef VXD_SENSORPLANE_H
#define VXD_SENSORPLANE_H
 
#include <vxd/dataobjects/VxdID.h>
#include <vxd/geometry/SensorInfoBase.h>
#include <framework/logging/Logger.h>
//ROOT CINT has problems with the boost classes used by the GeoCache but it
//does not need to see them anyway
#ifndef __CINT__
#include <vxd/geometry/GeoCache.h>
#endif
 
#include <genfit/AbsFinitePlane.h>
#include <cmath>
 
namespace Belle2 {
  namespace VXD {
    class SensorPlane: public genfit::AbsFinitePlane {
    public:
      SensorPlane(VxdID sensorID = 0, double uTolerance = 0, double vTolerance = 0):
        m_sensorID(sensorID), m_uTolerance(uTolerance), m_vTolerance(vTolerance), m_cosPhi(1.), m_sinPhi(0.), m_sensorInfo(0) {}
 
      void setRotation(double phi)
      {
        m_cosPhi = cos(phi);
        m_sinPhi = sin(phi);
      }
      virtual ~SensorPlane() {}
 
      bool isInActive(double u, double v) const override
      {
#ifndef __CINT__
        //If running in ROOT CINT we do not know about GeoCache so we cannot get
        //the SensorInfo
        if (!m_sensorInfo) {
          m_sensorInfo = &VXD::GeoCache::get(m_sensorID);
        }
#endif
        //No sensorInfo set so we have to bail
        if (!m_sensorInfo) {
          B2FATAL("Could not find sensorInfo for VXD Sensor " << VxdID(m_sensorID));
        }
        double uRot = m_cosPhi * u - m_sinPhi * v;
        double vRot = m_sinPhi * u + m_cosPhi * v;
        return m_sensorInfo->inside(uRot, vRot, m_uTolerance, m_vTolerance);
      }
 
      void Print(const Option_t* option = "") const override;
 
      virtual genfit::AbsFinitePlane* clone() const override
      {
        return new SensorPlane(*this);
      }
 
    private:
      unsigned short m_sensorID;
      double m_uTolerance;
      double m_vTolerance;
      double m_cosPhi;
      double m_sinPhi;
      mutable const SensorInfoBase* m_sensorInfo; 
 
      ClassDefOverride(SensorPlane, 2)
    };
  } // vxd namespace
} // Belle2 namespace
 
#endif /* VXD_SENSORPLANE_H */