StEEmcIUCluster.h

#ifndef __StEEmcIUCluster_h__
#define __StEEmcIUCluster_h__

#include <TObject.h>
#include <TVector3.h>
#include "StEEmcPool/StEEmcA2EMaker/StEEmcTower.h"

class StEmcCluster;

class StEEmcIUCluster : public TObject {

 public:

  StEEmcIUCluster();
  ~StEEmcIUCluster();

  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 StEEmcIUCluster &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(StEEmcIUCluster,1);

};

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

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

typedef std::vector<StEEmcIUCluster> StEEmcIUClusterVec_t;

#endif