SimpleFacetRelationFilter.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/facetRelation/SimpleFacetRelationFilter.h>
 
#include <tracking/trackingUtilities/eventdata/hits/CDCFacet.h>
#include <tracking/trackingUtilities/eventdata/hits/CDCWireHit.h>
 
#include <tracking/trackingUtilities/utilities/StringManipulation.h>
 
#include <framework/core/ModuleParamList.templateDetails.h>
#include <framework/geometry/VectorUtil.h>
 
using namespace Belle2;
using namespace TrackFindingCDC;
using namespace TrackingUtilities;
 
SimpleFacetRelationFilter::SimpleFacetRelationFilter()
  : m_param_deviationCosCut(cos(0.4))
{
}
 
SimpleFacetRelationFilter::SimpleFacetRelationFilter(double deviationCosCut)
  : m_param_deviationCosCut(deviationCosCut)
{
}
 
void SimpleFacetRelationFilter::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 SimpleFacetRelationFilter::operator()(const CDCFacet& fromFacet,
                                             const CDCFacet& toFacet)
{
  if (fromFacet.getStartWireHit().isOnWire(toFacet.getEndWire())) return NAN;
 
  // the compatibility of the short legs or all?
  // start end to continuation middle end
  // start middle to continuation start end
 
  const ParameterLine2D& fromStartToMiddle = fromFacet.getStartToMiddleLine();
  const ParameterLine2D& fromStartToEnd = fromFacet.getStartToEndLine();
 
  const ParameterLine2D& toStartToEnd = toFacet.getStartToEndLine();
  const ParameterLine2D& toMiddleToEnd = toFacet.getMiddleToEndLine();
 
  const double fromMiddleCos = VectorUtil::CosPhi(fromStartToMiddle.tangential(), toStartToEnd.tangential());
  const double toMiddleCos = VectorUtil::CosPhi(fromStartToEnd.tangential(), toMiddleToEnd.tangential());
 
  // check both
  if (fromMiddleCos > m_param_deviationCosCut and toMiddleCos > m_param_deviationCosCut) {
    // the weight must be -2 because the overlap of the facets is two points
    // so the amount of two facets is 4 points hence the cellular automat
    // must calculate 3 + (-2) + 3 = 4 as cellstate
    // this can of course be adjusted for a more realistic information measure
    // ( together with the facet creator filter)
    return -2;
  } else {
    return NAN;
  }
}