CBAlgoCenterOfMassSmall.h

Go to the documentation of this file.

#ifndef RECOTOOL_CBALGOCENTEROFMASSSMALL_H
#define RECOTOOL_CBALGOCENTEROFMASSSMALL_H

#include <iostream>
#include "larreco/RecoAlg/CMTool/CMToolBase/CBoolAlgoBase.h"
#include <math.h>

namespace cmtool {
  class CBAlgoCenterOfMassSmall: public CBoolAlgoBase {
    
  public:
    
    CBAlgoCenterOfMassSmall();
    
    virtual ~CBAlgoCenterOfMassSmall(){};
 
    virtual bool Bool(const ::cluster::ClusterParamsAlg &cluster1,
                      const ::cluster::ClusterParamsAlg &cluster2);

    virtual void Reset(){}

    virtual void Report();

    void SetDebug(bool on) { _debug = on; }

    void SetMaxHitsSmallClus(size_t n) { _maxHits = n; }

    void SetMaxDistance(double d) { _MaxDist = d; }

    void SetMaxCOMDistance(double d) { _MaxCOMDistSquared = d; }

    void UseCOMInPoly(bool on) { _COMinPolyAlg = on; }

    void UseCOMClose(bool on) { _COMsClose = on; }

    void UseCOMNearClus(bool on) { _COMNearClus = on; }

    void SetLengthReach(double frac){ _lengthReach = frac; }

    double ShortestDistanceSquared(double point_x, double point_y, 
                                   double start_x, double start_y,
                                   double end_x,   double end_y  ) const;
    
    
  protected:

    bool _debug;
    size_t _maxHits;
    double _MaxDist; //Max distance between COM and start-end line of cluster
    double _MaxCOMDistSquared; //Max distance^2 between COMs in cm^2
    double _lengthReach; 
    bool _COMinPolyAlg;
    bool _COMsClose;
    bool _COMNearClus;
  };
}

#endif
 // end of doxygen group