CBAlgoStartNearEnd.cxx

Go to the documentation of this file.

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

namespace cmtool {

  //----------------------------------------
  CBAlgoStartNearEnd::CBAlgoStartNearEnd() : CBoolAlgoBase()
  //----------------------------------------
  {
    SetMaxStartEndSeparation(10);   //cm^2
    SetMinHits(40);                 //cm^2
    SetMaxAngle(20 * (3.14 / 180)); //rad: max angle of end-point cluster
    // Nothing to be done in the base class
  }

  //--------------------------------------------------------
  bool CBAlgoStartNearEnd::Bool(const ::cluster::ClusterParamsAlg& cluster1,
                                const ::cluster::ClusterParamsAlg& cluster2)
  //--------------------------------------------------------
  {

    double start_w1 = cluster1.GetParams().start_point.w;
    double start_t1 = cluster1.GetParams().start_point.t;
    double end_w1 = cluster1.GetParams().end_point.w;
    double end_t1 = cluster1.GetParams().end_point.t;

    double start_w2 = cluster2.GetParams().start_point.w;
    double start_t2 = cluster2.GetParams().start_point.t;
    double end_w2 = cluster2.GetParams().end_point.w;
    double end_t2 = cluster2.GetParams().end_point.t;

    double angle_1 = cluster1.GetParams().opening_angle;
    double angle_2 = cluster2.GetParams().opening_angle;

    size_t hits_1 = cluster1.GetHitVector().size();
    size_t hits_2 = cluster2.GetHitVector().size();

    if ((angle_1 < _maxopeningangle) and (hits_1 > _MinHits) and
        (((start_w2 - end_w1) * (start_w2 - end_w1) + (start_t2 - end_t1) * (start_t2 - end_t1)) <
         _separation)) {
      if (_verbose) std::cout << "Start in End!" << std::endl;
      return true;
    }

    if ((angle_2 < _maxopeningangle) and (hits_2 > _MinHits) and
        (((start_w1 - end_w2) * (start_w1 - end_w2) + (start_t1 - end_t2) * (start_t1 - end_t2)) <
         _separation)) {
      if (_verbose) std::cout << "Start in End!" << std::endl;
      return true;
    }

    return false;
  }

  //-----------------------
  void CBAlgoStartNearEnd::Report()
  //-----------------------
  {}

}