spacepoint.hh

Go to the documentation of this file.

#ifndef SPACEPOINT_HH
#define SPACEPOINT_HH

#include "data_base.hh"
#include "TVector3.h"
#include "TMatrixD.h"

namespace larlight {
  class spacepoint : public data_base {
    
  public:
    
    spacepoint(DATA::DATA_TYPE type=DATA::SpacePoint);
    
    spacepoint(Int_t id, 
           Double_t x,  Double_t y,  Double_t z,
           Double_t ex, Double_t ey, Double_t ez,
           Double_t chi2) : data_base(DATA::SpacePoint)
    {
      clear_data();
      fID=id;
      fXYZ[0]=x;
      fXYZ[1]=y;
      fXYZ[2]=z;
      fErrXYZ[0]=ex;
      fErrXYZ[1]=ey;
      fErrXYZ[2]=ez;
      fChisq=chi2;
    }
    
    spacepoint(const spacepoint& original) : data_base(original),
                         fID(original.fID),
                         fChisq(original.fChisq)
    {
      fXYZ[0]=original.fXYZ[0];
      fXYZ[1]=original.fXYZ[1];
      fXYZ[2]=original.fXYZ[2];
      fErrXYZ[0]=original.fErrXYZ[0];
      fErrXYZ[1]=original.fErrXYZ[1];
      fErrXYZ[2]=original.fErrXYZ[2];
    }
    
    virtual ~spacepoint(){}
    
    Int_t           ID()     const { return fID;     }
    const Double_t* XYZ()    const { return fXYZ;    }
    const Double_t* ErrXYZ() const { return fErrXYZ; }
    Double_t        Chisq()  const { return fChisq;  }
    
    void clear_data() {
      data_base::clear_data();
      fID = -1;
      fXYZ[0]=-99999;
      fXYZ[1]=-99999;
      fXYZ[2]=-99999;
      fErrXYZ[0]=-99999;
      fErrXYZ[1]=-99999;
      fErrXYZ[2]=-99999;
      fChisq=-1;
    }
    
  protected:
    
    Int_t                        fID;        
    Double_t                     fXYZ[3];    
    Double_t                     fErrXYZ[6]; 
    Double_t                     fChisq;     
    
  private:
    
    ClassDef(spacepoint,1)
      
  };
  
  class event_sps : public std::vector<larlight::spacepoint>,
            public event_base {
    
  public:
    
    event_sps(DATA::DATA_TYPE type=DATA::SpacePoint);
    
    event_sps(const event_sps& original) : std::vector<larlight::spacepoint>(original),
                       event_base(original)
    {}
    
    ~event_sps(){}
    
    virtual void clear_data(){event_base::clear_data(); clear();}
    
    void get_axis_range (Double_t &max, Double_t &min, const int axis) const;
    
    void get_axis_range (Double_t &xmax, Double_t &xmin,
             Double_t &ymax, Double_t &ymin,
             Double_t &zmax, Double_t &zmin) const;
    
  private:
    
    ClassDef(event_sps,3)
  };
}
#endif
 // end of doxygen group 

---