Cluster.h

Go to the documentation of this file.

// $Id: Cluster.h,v 1.7 2009/08/25 00:52:55 messier Exp $
//
// Definition of reco Cluster of CellHits
//
// \author $Author: messier $
// \date   $Date: 2009/08/25 00:52:55 $
//

#ifndef CLUSTER_H
#define CLUSTER_H

#include "TObject.h"
#include "TRefArray.h"
#include "RecoBase/CellHit.h"
#include <vector>

using namespace std;

namespace recobase {
  class Cluster : public TObject
  {
  public:
    Cluster();
    Cluster(vector<const recobase::CellHit*>& cellvec);
    ~Cluster();
    
    // Note: the redundant 'recobase::' is not required by C++ but by
    // CINT. Best to just keep it until CINT catches up.
    bool Add(const recobase::CellHit* cell);
    int  NPlane()  const {return fXPlane.size() + fYPlane.size(); }
    int  NXPlane() const {return fXPlane.size(); }
    int  NYPlane() const {return fYPlane.size(); }
    
    int  Plane(int i) const;
    int  XPlane(int i) const;
    int  YPlane(int i) const;
    
    int  NCell()   const;
    int  NXCell()  const;
    int  NYCell()  const;
    int  NCell(geo::View_t view)   const;
    
    bool Is3D() const {return (this->NXPlane()>0 && this->NYPlane()>0); }
    
    const recobase::CellHit* XCell(int i) const;
    const recobase::CellHit* YCell(int i) const;
    const recobase::CellHit* Cell(geo::View_t view, int i) const;
    
    const recobase::CellHit* XCell(int icell, int iplane) const;
    const recobase::CellHit* YCell(int icell, int iplane) const;
    const recobase::CellHit* Cell(geo::View_t view, int icell, int iplane) const;
    
    const double* XYZT() const { return f4Vec; }
    double Z()       const {return f4Vec[2];}
    double dZ()      const {return fd4Vec[2];}
    double T()       const {return f4Vec[3];}
    double dT()      const {return fd4Vec[3];}
    double X()       const {return f4Vec[0];}
    double dX()      const {return fd4Vec[0];}
    double Y()       const {return f4Vec[1];}
    double dY()      const {return fd4Vec[1];}
    double Q()       const {return fQ;}
    double dQ()      const {return fdQ;}
    
    void SetQ(double q)        { fQ = q; }
    void SetdQ(double dq)      { fdQ = dq; }
    void SetXYZT(double x[])   { for (int i=0; i<4; ++i) f4Vec[i] = x[i]; }
    void SetdXYZT(double dx[]) { for (int i=0; i<4; ++i) fd4Vec[i] = dx[i]; }
    
    void Clear();
    
  private:
    vector<const recobase::CellHit*> fXCell;
    vector<const recobase::CellHit*> fYCell;
    
    vector<int> fXPlane; 
    vector<int> fYPlane; 
    vector<int> fPlane;  
    
    short  fStatus;
    double f4Vec[4];
    double fd4Vec[4];
    double fQ, fdQ;
    
    ClassDef(Cluster, 2);
    
  };
}

#endif // CLUSTER_H

---