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;
  }

}