StEEmcCluster.h

#ifndef __StEEmcCluster_h__
#define __StEEmcCluster_h__

#include <TObject.h>
#include <TVector3.h>
#include "StEEmcTower.h"

class StEmcCluster;

class StEEmcCluster : public TObject {

 public:

  StEEmcCluster();
  ~StEEmcCluster();

  void add( StEEmcTower, Float_t weight=1.0 );
  
  Float_t energy();
  Float_t seedEnergy();
  TVector3 momentum();
  Int_t numberOfTowers(); 
  StEEmcTower tower(Int_t t); 
  StEEmcTower tower(Int_t t) const; 

  Float_t weight(Int_t t); 
  StEEmcTowerVec_t towers();

  void momentum( TVector3 p);

  StEmcCluster *stemc();

  void stemc(StEmcCluster *c){ mEmcCluster = c; }

  Int_t key(){ return mKey; }
  void  key(Int_t k){ mKey=k; }

  Bool_t operator==( const StEEmcCluster &other ) const { return this->tower(0).index() == other.tower(0).index(); }

  void print();
  
 private:
 protected:

  Int_t mKey;

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

  Float_t mEnergy;
  TVector3 mMomentum;

  StEmcCluster *mEmcCluster;

  ClassDef(StEEmcCluster,1);

};

inline Float_t StEEmcCluster::energy(){ return mEnergy; }
inline Float_t StEEmcCluster::seedEnergy(){ return (mTowers.size())?mTowers[0].energy():0.; }

inline void StEEmcCluster::momentum(TVector3 p){mMomentum=p;}
inline TVector3 StEEmcCluster::momentum(){ return mMomentum;}
inline StEEmcTower StEEmcCluster::tower(Int_t t){ return mTowers[t]; } 
inline StEEmcTower StEEmcCluster::tower(Int_t t)const{ return mTowers[t]; } 
inline Int_t StEEmcCluster::numberOfTowers(){ return (Int_t)mTowers.size(); } 
inline StEEmcTowerVec_t StEEmcCluster::towers(){ return mTowers; } 
inline Float_t StEEmcCluster::weight( Int_t t ) { return mWeights[t]; } 

typedef std::vector<StEEmcCluster> StEEmcClusterVec_t;

#endif