StEpdGeom.h

#ifndef _StEpdGeom
#define _StEpdGeom

//#include "Rtypes.h"
#include "TVector3.h"

/*************************************
 * \author Mike Lisa
 * \date 4 Jan 2018
 *
 * \description:
 *  class defining geometrical aspects of an EPD Tile
 *  (position of center, RandomPointOnTile(), etc.)
 *
 * The user may pass the PP/TT/EW _or_ the uniqueID to
 *   most functions.  No option to pass StEpdHit object,
 *   because we want to avoid StObject-dependence.
 *   Instead, simply use StEpdHit::id(),
 *   e.g. RandomPointOnTile(hit->id())
 *************************************/

class TRandom3;

class StEpdGeom{
 private:
  short   mPP;  
  short   mTT;  
  short   mEW;  
  TRandom3* pRan;

  double mPhiCenter[12][31][2];  // PP,TT,EW
  double mRmin[16];   // row
  double mRmax[16];   // row
  double mRave[16];   // row

  // EastNumberOfOverlappingBbcTiles[PP-1][TT-1] =  number of BBC tiles that overlap
  static short mEastNumberOfOverlappingBbcTiles[12][9];  
  // EastBbcTilesWhichOverlap[PP-1][TT-1][j] gives the BBC tile ID of the jth overlapping BBC tile
  static short mEastBbcTilesWhichOverlap[12][9][3];  // 

  // WestNumberOfOverlappingBbcTiles[PP-1][TT-1] =  number of BBC tiles that overlap
  static short mWestNumberOfOverlappingBbcTiles[12][9];  
  // WestBbcTilesWhichOverlap[PP-1][TT-1][j] gives the BBC tile ID of the jth overlapping BBC tile
  static short mWestBbcTilesWhichOverlap[12][9][3];  // 


  void InitializeGeometry();

  /* these methods are used internally */
  double GetZwheel();

  //--------- obsolete (always the plan; that's why it was private)  double GetPhiCenter();
  //--------- obsolete (always the plan; that's why it was private)  void     GetRminRmax(double *Rmin, double *Rmax);
  short  Row();

  void     SetPpTtEw(short uniqueID);

  TVector3 TileCenter();

  TVector3 RandomPointOnTile();

  void     GetCorners(int* nCorners, double* x, double* y);

  bool   IsInTile(double x, double y);

  void GetOverlappingBbcTiles(int* nOverlappingBbcTiles, short* BbcTileIDs);





 public:

  StEpdGeom();
  ~StEpdGeom();

  unsigned short position(short uniqueID);
  unsigned short tile(short uniqueID);


  TVector3 TileCenter(short uniqueID);

  TVector3 TileCenter(short position, short tilenumber, short eastwest);

  TVector3 RandomPointOnTile(short uniqueID);

  TVector3 RandomPointOnTile(short position, short tilenumber, short eastwest);

  void GetCorners(short uniqueID,
                  int* nCorners, double* x, double* y);
  void GetCorners(short position, short tilenumber, short eastwest,
                  int* nCorners, double* x, double* y);

  void GetOverlappingBbcTiles(short uniqueID,
                              int* nOverlappingBbcTiles, short* BbcTileIDs);

  void GetOverlappingBbcTiles(short position, short tilenumber, short eastwest,
                              int* nOverlappingBbcTiles, short* BbcTileIDs);



  bool IsInTile(short uniqueID,
                  double x, double y);

  bool IsInTile(short position, short tilenumber, short eastwest,
                  double x, double y);

  bool IsWest(short uniqueID);

  bool IsWest(short position, short tilenumber, short eastwest);

  bool IsEast(short uniqueID);
  bool IsEast(short position, short tilenumber, short eastwest);

  short Row(short uniqueID);

  short Row(short position, short tilenumber, short eastwest);
};

inline bool StEpdGeom::IsWest(short uniqueID){return uniqueID>0;}
inline bool StEpdGeom::IsWest(short position, short tilenumber, short eastwest){return eastwest>0;}
inline bool StEpdGeom::IsEast(short uniqueID){return uniqueID<0;}
inline bool StEpdGeom::IsEast(short position, short tilenumber, short eastwest){return eastwest<0;}
inline unsigned short StEpdGeom::position(short uniqueID){return abs(uniqueID/100);}
inline unsigned short StEpdGeom::tile(short uniqueID){return abs(uniqueID%100);}




#endif