EndPointAlg.h

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \file  EndPointAlg.h
/// \brief algorithm to find 2D endpoints
///
/// \author  joshua.spitz@yale.edu
////////////////////////////////////////////////////////////////////////

#ifndef ENDPOINTALG_H
#define ENDPOINTALG_H

namespace art {
  class Event;
}
#include "canvas/Persistency/Common/PtrVector.h"
namespace fhicl { class ParameterSet; }

#include <string>
#include <vector>

#include "lardataobj/RecoBase/Cluster.h"
#include "lardataobj/RecoBase/EndPoint2D.h"
#include "lardataobj/RecoBase/Hit.h"

namespace cluster {

  ///Algorithm to find 2D end points
 class EndPointAlg {

  public:

    explicit EndPointAlg(fhicl::ParameterSet const& pset);

    void   reconfigure(fhicl::ParameterSet const& pset);

    size_t EndPoint(const art::PtrVector<recob::Cluster>           & clusIn,
                    std::vector<recob::EndPoint2D>                 & vtxcol,
                    std::vector< art::PtrVector<recob::Hit> >      & vtxHitsOut,
                    art::Event                                const& evt,
                    std::string                               const& label) const;

  private:

    double Gaussian(int x, int y, double sigma) const;
    double GaussianDerivativeX(int x, int y) const;
    double GaussianDerivativeY(int x, int y) const;
    void VSSaveBMPFile(const char *fileName, unsigned char *pix, int dx, int dy) const;


    int          fTimeBins;
    int          fMaxCorners;
    double       fGsigma;
    int          fWindow;
    double       fThreshold;
    int          fSaveVertexMap;
  };

}



#endif // ENDPOINTALG_H