StEEmcPoint.h

#ifndef __StEEmcPoint_h__
#define __StEEmcPoint_h__

#include "TObject.h"
#include "TVector3.h"

#include "StEEmcTower.h"
#include "StEEmcSmdCluster.h"

class StEmcPoint;

class StEEmcPoint : public TObject {

 public:

  StEEmcPoint();
  ~StEEmcPoint(){ /* nada */ };

  StEEmcPoint( const StEEmcPoint &p );

  void position( TVector3 p ) { mPosition=p; }
  void energy( Float_t e ) { mEnergy=e; }
  void fraction ( Float_t f ) { mFraction=f; }
  void tower( StEEmcTower t, Float_t w=1. ) { mTowers.push_back(t); mWeights.push_back(w); }
  void cluster( StEEmcSmdCluster c, Int_t plane ){mSmdClusters[plane]=c;} 

  void numberOfRelatives( Int_t r ){ mRelatives=r; }

  TVector3 position(){ return mPosition; }
  Float_t  energy(){ return mEnergy; }
  Float_t  fraction(){ return mFraction; }
  Float_t  energy() const { return mEnergy; }
  Int_t numberOfTowers(){ return (Int_t)mTowers.size(); }
  StEEmcTower tower(Int_t t){ return mTowers[t]; }
  Float_t weight(Int_t t){ return mWeights[t]; }
  StEEmcSmdCluster cluster(Int_t c){ return mSmdClusters[c]; } 

  Int_t numberOfRelatives(){ return mRelatives; }

  StEmcPoint *stemc();

  Int_t sector();

  void print();
  
  Bool_t operator<( const StEEmcPoint &other ) const; 
  Bool_t chiSquare( const StEEmcPoint &other ) const;

 private:
 protected:

  TVector3 mPosition;
  Float_t mEnergy;
  Float_t mFraction;
 
  StEEmcSmdCluster mSmdClusters[2];

  StEEmcTowerVec_t mTowers;       
  std::vector<Float_t> mWeights;  

  Int_t mRelatives;

  StEmcPoint *mEmcPoint;

  ClassDef(StEEmcPoint,1);

};

typedef std::vector<StEEmcPoint> StEEmcPointVec_t;

inline Bool_t StEEmcPoint::operator<( const StEEmcPoint &other ) const { return this->energy() < other.energy(); }  
inline Bool_t chiSquare( const StEEmcPoint &me, const StEEmcPoint &other ) { return me.chiSquare( other ); }

inline Int_t StEEmcPoint::sector(){ if ( mSmdClusters[0].sector()==mSmdClusters[1].sector() ) return mSmdClusters[0].sector(); else return -1; }

#endif