StEEmcPointMaker.h

#ifndef __StEEmcPointMaker_h__
#define __StEEmcPointMaker_h__

#include "StMaker.h"
#include "TString.h"
#include "StEEmcPool/StEEmcClusterMaker/StEEmcCluster.h"
#include "StEEmcPool/StEEmcClusterMaker/StEEmcSmdCluster.h"
#include "StEEmcPoint.h" 


#include <map>

class StEEmcA2EMaker;
class StEEmcClusterMaker;

class EEmcGeomSimple;
class EEmcSmdMap;
class EEmcSmdGeom; 

class StEEmcPointMaker : public StMaker {

 public:

  StEEmcPointMaker( const Char_t *name );
  ~StEEmcPointMaker(){ /* nada */ };

  Int_t Init();
  Int_t Make();
  void  Clear(Option_t *opts="");
  void analysis(const Char_t *name);
  void clusters(const Char_t *name);

  void towerThreshold( Float_t t) { mTowerThreshold=t; }

  void smdMatch( Float_t s ){ mSmdMatch = s; }

  StEEmcPointVec_t points();
  Int_t numberOfPoints(); 
  StEEmcPoint point(Int_t ipoint); 

  StEEmcPointVec_t smdPoints();
  Int_t numberOfSmdPoints();
  StEEmcPoint smdPoint(Int_t ip); 

  Float_t energySeen(){ return mEseen; }

  void setFillStEvent(){ mFillStEvent=true; }

  void setEnergyMode(Int_t mode){ mEnergyMode=mode; }
  void setLimit(Int_t l){ mLimit=l; }


  StEEmcPoint point( StEmcPoint *p ){ return mEtoEE[p]; }


 private:
 protected:

  TString mNameAnalysis; 
  TString mNameClusters; 
  StEEmcA2EMaker     *mEEanalysis; 
  StEEmcClusterMaker *mEEclusters; 

  EEmcGeomSimple *mEEtow; 
  EEmcSmdGeom    *mEEsmd; 
  EEmcSmdMap     *mEEmap; 

  void fillStEvent();
  void verifyStEvent();

  StEEmcPointVec_t mPoints;  

  StEEmcPointVec_t mSmdPoints;  

  std::map< StEEmcSmdCluster, StEEmcPointVec_t > mUclusters2points; 
  std::map< StEEmcSmdCluster, StEEmcPointVec_t > mVclusters2points; 
  
  StEEmcPointVec_t buildSmdPoints( Int_t sector, StEEmcSmdClusterVec_t &u, StEEmcSmdClusterVec_t &v );
  Bool_t findPoints ( Int_t sector, 
                      StEEmcSmdClusterVec_t u, 
                      StEEmcSmdClusterVec_t v, 
                      StEEmcPointVec_t &p );

  void shareEnergy();
  void shareEnergySimple();
  void shareEnergySmd();
  void countRelatives();
  void removeCluster( StEEmcSmdClusterVec_t &clusters, Int_t key );

  Float_t  fracp2t( StEEmcPoint &p, StEEmcTower &t );

  Float_t mTowerThreshold;
  Float_t mSmdMatch;
  Float_t mEseen;
  Bool_t mFillStEvent;
  Int_t mEnergyMode;
  Int_t mLimit;

  std::map<StEmcPoint *, StEEmcPoint> mEtoEE;    

  ClassDef(StEEmcPointMaker,1);

};

inline void StEEmcPointMaker::analysis(const Char_t *n){ mNameAnalysis=n; }
inline void StEEmcPointMaker::clusters(const Char_t *n){ mNameClusters=n; }

inline StEEmcPointVec_t StEEmcPointMaker::points(){ return mPoints ; }
inline Int_t StEEmcPointMaker::numberOfPoints(){ return mPoints.size(); } 
inline StEEmcPoint StEEmcPointMaker::point(Int_t ip){ return mPoints[ip]; } 

inline StEEmcPointVec_t StEEmcPointMaker::smdPoints(){ return mSmdPoints; }
inline Int_t StEEmcPointMaker::numberOfSmdPoints(){ return mSmdPoints.size(); } 
inline StEEmcPoint StEEmcPointMaker::smdPoint(Int_t ip){ return mSmdPoints[ip]; } 
#endif