CBAlgoAngleIncompat.cxx

Go to the documentation of this file.

#include "larreco/RecoAlg/CMTool/CMTAlgMerge/CBAlgoAngleIncompat.h"

#include <cmath>

namespace cmtool {

  //-------------------------------------------------------
  CBAlgoAngleIncompat::CBAlgoAngleIncompat() : CBoolAlgoBase()
  //-------------------------------------------------------
  {
    SetAngleCut(30.);

    SetMinHits(0);

    SetAllow180Ambig(false);

    SetUseOpeningAngle(false);

    SetMinLength(20);

    SetDebug(false);
  }

  //----------------------------------------------------------------
  bool CBAlgoAngleIncompat::Bool(const ::cluster::ClusterParamsAlg& cluster1,
                                 const ::cluster::ClusterParamsAlg& cluster2)
  //----------------------------------------------------------------
  {
    //if either cluster is too small, do not prohibit them
    if (cluster1.GetNHits() < _minHits || cluster2.GetNHits() < _minHits) return false;

    //if either cluster is too short, do not prohibit them
    if (cluster1.GetParams().length < _min_length || cluster2.GetParams().length < _min_length)
      return false;

    double angle1 = cluster1.GetParams().angle_2d;
    double angle2 = cluster2.GetParams().angle_2d;

    if (angle1 < -998 || angle2 < -998) return false;

    bool compatible = false;

    double my_cut_value = _max_allowed_2D_angle_diff;

    //if using opening angle, have angle cutoff be the smaller of the two opening angles
    if (_use_opening_angle)
      my_cut_value =
        std::min(cluster1.GetParams().opening_angle, cluster2.GetParams().opening_angle);

    //if you don't care if clusters have been reconstructed backwards
    if (_allow_180_ambig)
      compatible =
        (abs(angle1 - angle2) < my_cut_value || abs(angle1 - angle2 - 180) < my_cut_value ||
         abs(angle1 - angle2 + 180) < my_cut_value);
    else
      compatible = (abs(angle1 - angle2) < my_cut_value);

    if (_debug) {
      if (!compatible)

        std::cout << "CBAlgoAngleIncompat is PROHIBITING clusters with angles " << angle1 << " and "
                  << angle2 << std::endl;
    }

    return !compatible;
  }

}