SimpleFacetFilter.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/trackFindingCDC/filters/facet/SimpleFacetFilter.h>
 
#include <tracking/trackingUtilities/eventdata/hits/CDCFacet.h>
 
#include <tracking/trackingUtilities/utilities/StringManipulation.h>
 
#include <framework/core/ModuleParamList.templateDetails.h>
#include <framework/geometry/VectorUtil.h>
 
#include <cmath>
 
using namespace Belle2;
using namespace TrackFindingCDC;
using namespace TrackingUtilities;
 
SimpleFacetFilter::SimpleFacetFilter()
  : m_param_deviationCosCut(std::cos(M_PI / 180.0 * 9))
{
}
 
SimpleFacetFilter::SimpleFacetFilter(double deviationCosCut)
  : m_param_deviationCosCut(deviationCosCut)
{
}
 
void SimpleFacetFilter::exposeParameters(ModuleParamList* moduleParamList,
                                         const std::string& prefix)
{
  Super::exposeParameters(moduleParamList, prefix);
  moduleParamList->addParameter(prefixed(prefix, "deviationCosCut"),
                                m_param_deviationCosCut,
                                "Acceptable deviation cosine in the angle of adjacent tangents "
                                "to the drift circles.",
                                m_param_deviationCosCut);
}
 
Weight SimpleFacetFilter::operator()(const CDCFacet& facet)
{
  facet.adjustFitLine();
 
  const ParameterLine2D& startToMiddle = facet.getStartToMiddleLine();
  const ParameterLine2D& startToEnd = facet.getStartToEndLine();
  const ParameterLine2D& middleToEnd = facet.getMiddleToEndLine();
 
  const double startCos  = VectorUtil::CosPhi(startToMiddle.tangential(), startToEnd.tangential());
  // const double middleCos = VectorUtil::CosPhi(startToMiddle.tangential(), middleToEnd.tangential());
  const double endCos = VectorUtil::CosPhi(startToEnd.tangential(), middleToEnd.tangential());
 
  // cut on the angle of
  if (startCos > m_param_deviationCosCut and endCos > m_param_deviationCosCut) {
    // Good facet contains three points of the track
    // the amount carried by this facet can the adjusted more realistically
    return 3;
  } else {
    return NAN;
  }
}