StEEmcA2EMaker.h

#ifndef __StEEmcA2EMaker_h__
#define __StEEmcA2EMaker_h__

#include <StMaker.h>
#include <TString.h>
#include <vector>

#include "StEEmcUtil/EEmcGeom/EEmcGeomDefs.h"
#include "StEEmcUtil/EEmcGeom/EEmcGeomSimple.h"

#include "StEEmcTower.h"
#include "StEEmcStrip.h"

class StEEmcDb;

class StMuEmcCollection;
class StEmcCollection;

class StEEmcA2EMaker : public StMaker {
public:
  StEEmcA2EMaker( const Char_t *name="mEEanalysis" );
  virtual ~StEEmcA2EMaker();

  virtual Int_t Init();
  virtual Int_t Make();
  virtual void  Clear(Option_t *opts="");

  void threshold(Float_t cut, Int_t layer) { mSigmaPed[layer]=cut; }

  void database(const Char_t * /*dbname*/) { /*mDbName=dbname;*/ } /* obsolete */

  void source(const Char_t * /*name*/, Int_t /*type*/ = 0) { /*mInputName=name; mInputType=type;*/ } /* obsolete */

  StEEmcTowerVec_t &towers(Int_t layer=0){ return mHitTowers[layer]; }
  const StEEmcTowerVec_t &towers(Int_t layer=0) const { return mHitTowers[layer]; }

  StEEmcStripVec_t &strips(Int_t sec, Int_t pln ){ return mHitStrips[sec][pln]; }
  const StEEmcStripVec_t &strips(Int_t sec, Int_t pln ) const { return mHitStrips[sec][pln]; }

  Int_t numberOfHitTowers(Int_t layer) const { return (Int_t)mHitTowers[layer].size(); }

  StEEmcTower &hittower(Int_t hit, Int_t layer){ return mHitTowers[layer][hit]; } 
  const StEEmcTower &hittower(Int_t hit, Int_t layer) const { return mHitTowers[layer][hit]; } 

  StEEmcTower &tower(Int_t index, Int_t layer=0) { return mTowers[index][layer]; }
  const StEEmcTower &tower(Int_t index, Int_t layer=0) const { return mTowers[index][layer]; }

  StEEmcTower &tower(Int_t sector, Int_t subsector, Int_t etabin, Int_t layer=0 ){return tower(index(sector,subsector,etabin),layer);}
  const StEEmcTower &tower(Int_t sector, Int_t subsector, Int_t etabin, Int_t layer=0) const {return tower(index(sector,subsector,etabin),layer);}

  StEEmcTower *tower(TVector3 &r, Int_t layer);
  const StEEmcTower *tower(TVector3 &r, Int_t layer) const;
  
  Int_t phibin( Int_t sector, Int_t subsector) const { return sector * 5 + subsector; } 
  Int_t index( Int_t sector, Int_t subsector, Int_t etabin) const { return 12 * phibin(sector, subsector) + etabin; }  
  

  StEEmcTower &hightower(Int_t layer=0) { return *mHighTower[layer]; }
  const StEEmcTower &hightower(Int_t layer=0) const { return *mHighTower[layer]; }
  
  Int_t numberOfHitStrips(Int_t sector, Int_t plane) const { return (Int_t)mHitStrips[sector][plane].size(); }

  StEEmcStrip &hitstrip(Int_t sec,Int_t pl, Int_t hit){ return mHitStrips[sec][pl][hit]; } 
  const StEEmcStrip &hitstrip(Int_t sec,Int_t pl, Int_t hit) const { return mHitStrips[sec][pl][hit]; } 
  
  StEEmcStrip &strip(Int_t sector, Int_t plane, Int_t strip) { return mStrips[sector][plane][strip]; }
  const StEEmcStrip &strip(Int_t sector, Int_t plane, Int_t strip) const { return mStrips[sector][plane][strip]; }
  
  Float_t energy(Int_t sector, Int_t layer) const { return mEnergy[sector][layer]; }
  Float_t energy(Int_t layer) const; 

  Int_t numberOfHits(Int_t sector, Int_t layer) const { return mHits[sector][layer]; } 

  void scale(Float_t s){ mScale=s; }

  virtual const char* GetCVS() const
    {static const char cvs[]="Tag $Name:  $ $Id: StEEmcA2EMaker.h,v 1.8 2014/08/06 11:42:59 jeromel Exp $ built " __DATE__ " " __TIME__; return cvs;}

protected:
  
  Float_t mScale;         
  Float_t mSigmaPed[6];   
  const StEEmcDb      *mDbMaker;      
  const EEmcGeomSimple *mEEgeom;      

  Bool_t readData();     
  Bool_t fillFromMuDst( const StMuEmcCollection *emc );
  Bool_t fillFromSt( const StEmcCollection *emc );

  Float_t addTowerHit(Int_t sec, Int_t sub, Int_t eta, Float_t adc, Int_t layer);
  Float_t addSmdHit(Int_t sec, Int_t plane, Int_t str, Float_t adc );

  StEEmcTower mTowers[kEEmcNumSectors*kEEmcNumSubSectors*kEEmcNumEtas][4];
  StEEmcStrip mStrips[kEEmcNumSectors][kEEmcNumSmdUVs][kEEmcNumStrips];

  std::vector< StEEmcTowerVec_t > mHitTowers;

  StEEmcTower *mHighTower[4]; 
  std::vector< std::vector< StEEmcStripVec_t > > mHitStrips;


  Float_t mEnergy[kEEmcNumSectors][6];
  Int_t   mHits[kEEmcNumSectors][6]; 

  ClassDef(StEEmcA2EMaker,1);

};

#endif