CBAlgoAngleCompat.cxx

Go to the documentation of this file.

#ifndef RECOTOOL_CBALGOANGLECOMPAT_CXX
#define RECOTOOL_CBALGOANGLECOMPAT_CXX

#include "CBAlgoAngleCompat.h"

namespace cmtool {

  CBAlgoAngleCompat::CBAlgoAngleCompat() : CBoolAlgoBase() {

    //this just sets default values

    SetDebug(true);

    SetAngleCut(30.); // in degrees

    SetMinHits(0);
    
    SetAllow180Ambig(false);

    SetUseOpeningAngle(false);

    //    angle_dist_histo = 0;

    //    angle_dist_histo = new TH1F("angle_dist_histo","Cluster Angle Differences",100,-360,360);

  } //end constructor

  bool CBAlgoAngleCompat::Bool(const ::cluster::ClusterParamsAlg &cluster1,
                               const ::cluster::ClusterParamsAlg &cluster2)
  {
    
    //if number of hits not large enough skip
    if ( _minHits and ((cluster1.GetHitVector().size() < _minHits) or (cluster2.GetHitVector().size() < _minHits)) ) {
      return false;
    }

    //pretty sure we don't need conversion factors here. 
    //already in cm/cm units, degrees? need to check that
    double angle1 = cluster1.GetParams().angle_2d;// * _time_2_cm / _wire_2_cm;
    double angle2 = cluster2.GetParams().angle_2d;// * _time_2_cm / _wire_2_cm;

    double w_start1 = cluster1.GetParams().start_point.w;// * _wire_2_cm;
    double t_start1 = cluster1.GetParams().start_point.t;// * _time_2_cm;
    double w_start2 = cluster2.GetParams().start_point.w;// * _wire_2_cm;
    double t_start2 = cluster2.GetParams().start_point.t;// * _time_2_cm;

    if (_debug){

      std::cout << "Cluster 1:" << std::endl;
      std::cout << "\tAngle: " << angle1 << std::endl;
      std::cout << "\tStart: ( " << w_start1 << ", " << t_start1 << " )" << std::endl;
      std::cout << "Cluster 2:" << std::endl;
      std::cout << "\tAngle: " << angle2 << std::endl;
      std::cout << "\tStart: ( " << w_start2 << ", " << t_start2 << " )" << 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;

    //    if(angle_dist_histo){
    //      angle_dist_histo->Fill(angle1-angle2);
    //    }
    //    else
    //      std::cout<<"\n\n\nSOMETHING WENT HORRIBLY WRONG IN CBALGOANGLECOMPAT\n\n\n\n\n\n\n"<<std::endl;

    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(_verbose) {
        if(compatible) std::cout<<"These two clusters are compatible in angle."<<std::endl;
        else std::cout<<"These two clusters are NOT compatible in angle."<<std::endl;
    }
    
    return compatible;
    
  } // end Merge function 
  

}//end namespace cmtool
#endif