GeometryHelper.h

#ifndef GEOMETRYHELPER_H
#define GEOMETRYHELPER_H

#include "HelperBase.h"

#include "larcore/Geometry/Geometry.h"
#include "lardataobj/RecoBase/PFParticle.h"
#include "lardataobj/RecoBase/SpacePoint.h"
#include "lardataobj/RecoBase/Vertex.h"
#include "lardataobj/RecoBase/Shower.h"

#include "canvas/Persistency/Common/FindManyP.h"
#include "canvas/Persistency/Common/FindOneP.h"

#include "TVector3.h"

namespace lee {

class GeometryHelper : public HelperBase {
public:
  GeometryHelper() = default;
  ~GeometryHelper() = default;

  bool isInside( std::vector<double> P, std::vector< std::vector<double> > V);

  bool isFiducial(const std::vector<double> &x) const;

  bool isFiducial(const TVector3 &x) const;

  bool isFiducial(const double x[3]) const;

  bool isActive(const std::vector<double> &x) const;

  bool isActive(const double x[3]) const;

  double distance(const std::vector<double> &a,
                  const std::vector<double> &b) const;

  double distance(const TVector3 &a, const TVector3 &b) const;

  double getPitch(const TVector3 &direction, const int &pl) const;

  void setFiducialVolumeCuts(float m_fidvolXstart, float m_fidvolXend,
                             float m_fidvolYstart, float m_fidvolYend,
                             float m_fidvolZstart, float m_fidvolZend);


  TVector3
  getAveragePosition(std::vector<art::Ptr<recob::SpacePoint>> &spcpnts);

  void buildRectangle(double length, double width, std::vector<double> &start,
                      std::vector<double> &axis,
                      std::vector<std::vector<double>> &points);


private:
  float m_fidvolXstart;
  float m_fidvolXend;
  float m_fidvolYstart;
  float m_fidvolYend;
  float m_fidvolZstart;
  float m_fidvolZend;
};
}

#endif