FourHitFilters.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Jakob Lettenbichler                                      *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <TVector3.h>
#include "tracking/vxdCaTracking/ThreeHitFilters.h"
 
 
 
namespace Belle2 {
  class FourHitFilters {
 
  public:
 
    FourHitFilters():
      m_radiusABC(0.),
      m_radiusBCD(0.),
      m_circleCalculated(false) {}
 
    FourHitFilters(TVector3& outer, TVector3& outerCenter, TVector3& innerCenter, TVector3& inner, 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(TVector3& outer, TVector3& outerCenter, TVector3& innerCenter, TVector3& inner)
    {
      m_hitA = outer;
      m_hitB = outerCenter;
      m_hitC = innerCenter;
      m_hitD = inner;
      m_radiusABC = 0.;
      m_radiusBCD = 0.;
      m_circleCalculated = false;
    }
 
    void resetMagneticField(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;
    }
 
    TVector3 m_hitA; 
    TVector3 m_hitB; 
    TVector3 m_hitC; 
    TVector3 m_hitD;  
    TVector3 m_centerABC;  
    TVector3 m_centerBCD; 
    double m_radiusABC; 
    double m_radiusBCD; 
    ThreeHitFilters m_threeHitFilter;  
    bool m_circleCalculated; 
  }; //end class FourHitFilters
} //end namespace Belle2