APAGeometryAlg.h

Go to the documentation of this file.

//  \fileAPAGeometryAlg.h
//  tylerdalion@gmail.com
#ifndef APAGeometryALG_H
#define APAGeometryALG_H
#include <vector>
#include <cmath>
#include <iostream>
#include <stdint.h>

#include "art/Framework/Services/Registry/ServiceHandle.h"
#include "fhiclcpp/ParameterSet.h"
#include "canvas/Persistency/Common/Ptr.h"
#include "canvas/Persistency/Common/PtrVector.h"

#include "larcore/Geometry/Geometry.h"
#include "larcoreobj/SimpleTypesAndConstants/geo_types.h"
#include "lardataobj/RecoBase/Wire.h"
#include "lardataobj/RecoBase/Hit.h"

#include "TMatrixD.h"
#include "TVectorD.h"
#include "TVector3.h"



namespace apa{

  // each APA has 4 separate views
  typedef enum _apa_plane_proj {
    kU,      
    kV,      
    kZ0,     
    kZ1,     
    kUnknown
  } APAView_t;


  //--------------------------------------------------------------- 
  class APAGeometryAlg {
  public:
    
    
    APAGeometryAlg(fhicl::ParameterSet const& pset);
    APAGeometryAlg();
    virtual ~APAGeometryAlg();
    
    void                 reconfigure(fhicl::ParameterSet const& p);

    void                 Init();                       

    bool                 APAChannelsIntersect( uint32_t chan1,
                                               uint32_t chan2,
                                               std::vector< geo::WireIDIntersection > & IntersectVector);

    bool                 LineSegChanIntersect(     TVector3 xyzStart,
                                                   TVector3 xyzEnd,
                                                   uint32_t chan,
                                                   std::vector< geo::WireID >& widsCrossed,
                                                   bool ExtendLine );

    std::vector<geo::WireID> ChanSegsPerSide(uint32_t chan, unsigned int side);
    std::vector<geo::WireID> ChanSegsPerSide(std::vector<geo::WireID> wids, unsigned int side);

    std::vector<double>  ThreeChanPos( uint32_t u, uint32_t v, uint32_t z );

    geo::WireID          NearestWireIDOnChan( const double WorldLoc[3],
                                              uint32_t chan,
                                              unsigned int const plane, 
                                              unsigned int const tpc=0, 
                                              unsigned int const cstat=0);

    unsigned int         ChannelToAPA(uint32_t chan);  
    void                 ChannelToAPA(         uint32_t chan, 
                                               unsigned int & apa, 
                                               unsigned int & cryo);
    APAView_t            APAView(uint32_t chan);       
    unsigned int         ChannelsInView( geo::View_t geoview );
    uint32_t             FirstChannelInView(   geo::View_t geoview, 
                                               unsigned int apa, 
                                               unsigned int cryo );
    uint32_t             FirstChannelInView(   geo::View_t geoview, 
                                               uint32_t chan );
    uint32_t             FirstChannelInView(   uint32_t chan );
    unsigned int         ChannelsInAPAView( APAView_t apaview );
    unsigned int         ChannelsPerAPA(){ return fChannelsPerAPA; };


  private:

    art::ServiceHandle<geo::Geometry> fGeom;           // handle to geometry service

    unsigned int fChannelsPerAPA;                      
    unsigned int fAPAsPerCryo;

    // channel boundaries to avoid calling ChannelToWire repetitively
    uint32_t fFirstU;
    uint32_t fLastU;
    uint32_t fFirstV;
    uint32_t fLastV;
    uint32_t fFirstZ0;
    uint32_t fLastZ0;
    uint32_t fFirstZ1;
    uint32_t fLastZ1;

    double fChannelRange[2]; // for each induction view: U=0, V=1

  }; // class APAGeometryAlg

} // namespace apa

#endif // ifndef APAGeometryALG_H