CDCRecoHit2D.cc

/**************************************************************************
 * 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.                  *
 **************************************************************************/
#include <tracking/trackingUtilities/eventdata/hits/CDCRecoHit2D.h>
 
#include <tracking/trackingUtilities/eventdata/hits/CDCRLWireHit.h>
 
#include <tracking/trackingUtilities/geometry/VectorUtil.h>
 
#include <tracking/trackingUtilities/numerics/ESign.h>
 
#include <cdc/dataobjects/CDCSimHit.h>
 
#include <ostream>
 
using namespace Belle2;
using namespace TrackingUtilities;
 
CDCRecoHit2D::CDCRecoHit2D(const CDCRLWireHit& rlWireHit) :
  m_rlWireHit(rlWireHit),
  m_recoDisp2D(0.0, 0.0)
{}
 
CDCRecoHit2D::CDCRecoHit2D(const CDCRLWireHit& rlWireHit,
                           const ROOT::Math::XYVector& recoDisp2D) :
  m_rlWireHit(rlWireHit),
  m_recoDisp2D(recoDisp2D)
{}
 
CDCRecoHit2D CDCRecoHit2D::fromSimHit(const CDCWireHit* wireHit,
                                      const CDCSimHit& simHit)
{
  // find out if the wire is right or left of the track ( view in flight direction )
  ROOT::Math::XYZVector trackPosToWire{simHit.getPosWire() - simHit.getPosTrack()};
  CDCRecoHit2D recoHit2D(CDCRLWireHit::fromSimHit(wireHit, simHit),
                         ROOT::Math::XYVector(-trackPosToWire.x(), -trackPosToWire.y()));
 
  recoHit2D.snapToDriftCircle();
  return recoHit2D;
}
 
CDCRecoHit2D CDCRecoHit2D::average(const CDCRecoHit2D& recoHit1,
                                   const CDCRecoHit2D& recoHit2)
{
  CDCRLWireHit rlWireHit =
    CDCRLWireHit::average(recoHit1.getRLWireHit(),
                          recoHit2.getRLWireHit());
 
  ROOT::Math::XYVector displacement =
    VectorUtil::average(recoHit1.getRecoDisp2D(),
                        recoHit2.getRecoDisp2D());
 
  CDCRecoHit2D result(rlWireHit, displacement);
  result.snapToDriftCircle();
 
  return result;
}
 
CDCRecoHit2D CDCRecoHit2D::average(const CDCRecoHit2D& recoHit1,
                                   const CDCRecoHit2D& recoHit2,
                                   const CDCRecoHit2D& recoHit3)
{
  CDCRLWireHit rlWireHit =
    CDCRLWireHit::average(recoHit1.getRLWireHit(),
                          recoHit2.getRLWireHit(),
                          recoHit3.getRLWireHit());
 
  ROOT::Math::XYVector displacement =
    VectorUtil::average(recoHit1.getRecoDisp2D(),
                        recoHit2.getRecoDisp2D(),
                        recoHit3.getRecoDisp2D());
 
  CDCRecoHit2D result(rlWireHit, displacement);
  result.snapToDriftCircle();
 
  return result;
}
 
CDCRecoHit2D CDCRecoHit2D::fromRecoPos2D(const CDCRLWireHit& rlWireHit,
                                         const ROOT::Math::XYVector& recoPos2D,
                                         bool snap)
{
  CDCRecoHit2D result(rlWireHit, recoPos2D - rlWireHit.getRefPos2D());
  if (snap) result.snapToDriftCircle();
  return result;
}
 
void CDCRecoHit2D::reverse()
{
  m_rlWireHit.reverse();
}
 
CDCRecoHit2D CDCRecoHit2D::reversed() const
{
  CDCRecoHit2D reversedRecoHit(*this);
  reversedRecoHit.reverse();
  return reversedRecoHit;
}
 
CDCRecoHit2D CDCRecoHit2D::getAlias() const
{
  return CDCRecoHit2D(getRLWireHit().getAlias(), -getRecoDisp2D());
}
 
 
void CDCRecoHit2D::setRefDriftLength(double driftLength, bool snapRecoPos)
{
  double oldDriftLength = m_rlWireHit.getRefDriftLength();
  m_rlWireHit.setRefDriftLength(driftLength);
  if (snapRecoPos) {
    bool switchSide = sign(oldDriftLength) != sign(driftLength);
    snapToDriftCircle(switchSide);
  }
}
 
void CDCRecoHit2D::snapToDriftCircle(bool switchSide)
{
  if (m_recoDisp2D.R() != 0.0) {
    m_recoDisp2D *= (getRLWireHit().getRefDriftLength() / m_recoDisp2D.R());
  }
  if (switchSide) {
    m_recoDisp2D = -m_recoDisp2D;
  }
}
 
ROOT::Math::XYZVector CDCRecoHit2D::reconstruct3D(const CDCTrajectory2D& trajectory2D, double z) const
{
  return getRLWireHit().reconstruct3D(trajectory2D, z);
}
 
std::ostream& TrackingUtilities::operator<<(std::ostream& output, const CDCRecoHit2D& recohit)
{
  output << "CDCRecoHit2D(" << recohit.getRLWireHit() << ","
         << recohit.getRecoDisp2D() << ")" ;
  return output;
}