CBAlgoAngleAlign.cxx

Go to the documentation of this file.

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

#include <cmath>

namespace cmtool {

  CBAlgoAngleAlign::CBAlgoAngleAlign() : CBoolAlgoBase()
  {

    //this just sets default values
    SetDebug(true);
    SetAngleCut(10.); // in degrees
    SetAllow180Ambig(false);
    SetMinNHits(30);

  } //end constructor

  bool CBAlgoAngleAlign::Bool(const ::cluster::ClusterParamsAlg& cluster1,
                              const ::cluster::ClusterParamsAlg& cluster2)
  {

    double angle1 = cluster1.GetParams().angle_2d;
    double angle2 = cluster2.GetParams().angle_2d;
    size_t hits1 = cluster1.GetHitVector().size();
    size_t hits2 = cluster2.GetHitVector().size();

    //if don't make hit cut return aflse
    if ((hits1 < _MinNHits) or (hits2 < _MinNHits)) return false;

    if (_debug) {
      std::cout << "Cluster1:" << std::endl;
      std::cout << "\tAngle: " << angle1 << std::endl;
      std::cout << "\t Start: (" << cluster1.GetParams().start_point.w << ", "
                << cluster1.GetParams().start_point.t << ")" << std::endl;
      std::cout << "Cluster2:" << std::endl;
      std::cout << "\tAngle: " << angle2 << std::endl;
      std::cout << "\t Start: (" << cluster2.GetParams().start_point.w << ", "
                << cluster2.GetParams().start_point.t << ")" << std::endl;
      std::cout << std::endl;
    }

    //for some reason angles are frequently -999.99.
    //if either angle is this, clearly the cluster 2d angle is not well defined
    //and this algorithm does not apply
    if (angle1 < -998 || angle2 < -998) return false;

    bool compatible = false;

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

    if (_verbose) {
      if (compatible) std::cout << "These two clusters are compatible in angle." << std::endl;
    }

    return compatible;

  } // end Merge function

} //end namespace cmtool