FourHitFilters.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 <framework/geometry/B2Vector3.h>
#include "tracking/vxdCaTracking/ThreeHitFilters.h"
 
 
 
namespace Belle2 {
  class FourHitFilters {
 
  public:
 
    FourHitFilters():
      m_radiusABC(0.),
      m_radiusBCD(0.),
      m_circleCalculated(false) {}
 
    FourHitFilters(const B2Vector3D& outer, const B2Vector3D& outerCenter, const B2Vector3D& innerCenter, const B2Vector3D& inner,
                   const double magneticFieldStrength = 1.5):
      m_hitA(outer),
      m_hitB(outerCenter),
      m_hitC(innerCenter),
      m_hitD(inner),
      m_radiusABC(0.),
      m_radiusBCD(0.),
      m_circleCalculated(false) { m_threeHitFilter.resetMagneticField(magneticFieldStrength); }
 
 
    ~FourHitFilters() {}
 
    void resetValues(const B2Vector3D& outer, const B2Vector3D& outerCenter, const B2Vector3D& innerCenter, const B2Vector3D& inner)
    {
      m_hitA = outer;
      m_hitB = outerCenter;
      m_hitC = innerCenter;
      m_hitD = inner;
      m_radiusABC = 0.;
      m_radiusBCD = 0.;
      m_circleCalculated = false;
    }
 
    void resetMagneticField(const double magneticFieldStrength = 1.5) { m_threeHitFilter.resetMagneticField(magneticFieldStrength); }
 
 
    double calcDeltapT();
 
 
    double deltapT() { return calcDeltapT(); }
 
 
    double calcDeltaDistCircleCenter();
 
 
    double deltaDistCircleCenter() { return calcDeltaDistCircleCenter(); }
  protected:
 
    void calcCircle()
    {
      m_threeHitFilter.calcCircleCenter(m_hitA, m_hitB, m_hitC, m_centerABC);
      m_radiusABC = m_threeHitFilter.calcRadius(m_hitA, m_hitB, m_hitC, m_centerABC);
      m_threeHitFilter.calcCircleCenter(m_hitB, m_hitC, m_hitD, m_centerBCD);
      m_radiusBCD = m_threeHitFilter.calcRadius(m_hitB, m_hitC, m_hitD, m_centerBCD);
 
      m_circleCalculated = true;
    }
 
    B2Vector3D m_hitA; 
    B2Vector3D m_hitB; 
    B2Vector3D m_hitC; 
    B2Vector3D m_hitD;  
    B2Vector3D m_centerABC;  
    B2Vector3D m_centerBCD; 
    double m_radiusABC; 
    double m_radiusBCD; 
    ThreeHitFilters m_threeHitFilter;  
    bool m_circleCalculated; 
  }; //end class FourHitFilters
} //end namespace Belle2