SpacePointAlg_TimeSort.h

Go to the documentation of this file.

#ifndef SPACEPOINTALG_TIMESORT_H
#define SPACEPOINTALG_TIMESORT_H

// LArSoft Includes
#include "lardataobj/RecoBase/Hit.h"
#include "lardataobj/RecoBase/SpacePoint.h"
#include "canvas/Persistency/Common/Ptr.h"
#include "fhiclcpp/ParameterSet.h"

//boost includes
#include "boost/multi_array.hpp"

namespace sppt{

  bool  HitTimeComparison(art::Ptr<recob::Hit> a, art::Ptr<recob::Hit> b) { return a->PeakTime() < b->PeakTime(); }

  class SpacePointAlg_TimeSort {

  public:
    SpacePointAlg_TimeSort(fhicl::ParameterSet const& pset);
    ~SpacePointAlg_TimeSort();

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

    void createSpacePoints(std::vector< art::Ptr<recob::Hit> > &hitVec_U,
                           std::vector< art::Ptr<recob::Hit> > &hitVec_V,
                           std::vector< art::Ptr<recob::Hit> > &hitVec_Y,
                           std::unique_ptr<std::vector<recob::SpacePoint> > &spptCollection,
                           std::unique_ptr<std::vector<std::vector<art::Ptr<recob::Hit> > > > &spptAssociatedHits);

  private:

    float          fTimeDiffMax;        
    float          fYDiffMax;           
    float          fZDiffMax;           

    bool TIME_OFFSET_SET;
    bool COORDINATES_FILLED;

    double TIME_OFFSET_U;
    double TIME_OFFSET_V;
    double TIME_OFFSET_Y;
    double TICKS_TO_X;

    boost::multi_array<double, 2> coordinates_UV_y;
    boost::multi_array<double, 2> coordinates_UV_z;
    boost::multi_array<double, 2> coordinates_UY_y;
    boost::multi_array<double, 2> coordinates_UY_z;

    void sortHitsByTime(std::vector< art::Ptr<recob::Hit> > &hits_handle);
    
  }; //class SpacePointAlg_TimeSort

} //end sppt namespace

#endif