StBbcGeom.h

#ifndef _StBbcGeom
#define _StBbcGeom

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

/*************************************
 * \author Mike Lisa
 * \date 1 Mar 2018
 *
 * \description:
 *  class defining geometrical aspects of the inner BBC
 *  (position of tile center, tile corners, etc.)
 *
 * The BBC is a holy nightmare because one PMT connects
 *  to *either* one or two disjoint tiles.  Also, the
 *  PMT number is the same as the tile number sometimes,
 *  but not all the time.  Plus the numbering starts
 *  at zero sometimes and not other times.
 *
 * In this class, we are only talking about tiles
 *  and tile geometry.  Another class will give you the
 *  map from PMT number to tile(s).
 *
 *
 * Because some users will use the sign of the tileId
 *  to denote east/west, and some will use an explicit
 *  integer, we provide overloaded functions, to allow both.
 *
 ******************************************/



class StBbcGeom{
 private:

  TVector3 mEastTileCenter[18];
  TVector3 mWestTileCenter[18];
  double mX[6];   // x-coordinates of hexagonal tile corners, if tile center is at (0,0)
  double mY[6];   // y-coordinates of hexagonal tile corners, if tile center is at (0,0)
  unsigned short mPMT[18];  // phototube number corresponding to a given inner BBC tile.  Same for east and west BBCs

  unsigned short mNtilesOfPmt[16];          // how many tiles correspond to the phototube?  (either 1 or 2)
  unsigned short mTileNumbersOfPmt[16][2];  // WHICH tiles correspond to the phototube?


// mEastNumberOfOverlappingEpdTiles[BB-1]  gives the number of EPD tiles that have any overlap with BBC tile BB
  static short mEastNumberOfOverlappingEpdTiles[18];
// mEastEpdTilesWhichOverlap[BB-1][i]  gives the EPD tile id of the jth EPD tile that overlaps BBC inner tile BB
  static short mEastEpdTilesWhichOverlap[18][10];

// mWestNumberOfOverlappingEpdTiles[BB-1]  gives the number of EPD tiles that have any overlap with BBC tile BB
  static short mWestNumberOfOverlappingEpdTiles[18];
// mWestEpdTilesWhichOverlap[BB-1][i]  gives the EPD tile id of the jth EPD tile that overlaps BBC inner tile BB
  static short mWestEpdTilesWhichOverlap[18][10];





  void InitializeGeometry();

 public:
  StBbcGeom();
  ~StBbcGeom();

  //------------------------------------------------------

  TVector3 TileCenter(short tileId);

  TVector3 TileCenter(unsigned short absValueTileNumber, short eastwest);

  //------------------------------------------------------

  void GetCorners(short tileId, double* x, double* y);

  void GetCorners(unsigned short absValueTileNumber, short eastwest, double* x, double* y);

  //------------------------------------------------------

  void GetOverlappingEpdTiles(short tileId,
                              int* nOverlappingEpdTiles, short* EpdTileIDs);


  void GetOverlappingEpdTiles(unsigned short absValueTileNumber, short eastwest,
                              int* nOverlappingEpdTiles, short* EpdTileIDs);

  //------------------------------------------------------

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

  bool IsInTile(unsigned short absValueTileNumber, short eastwest, double x, double y);


  unsigned short PmtOfTile(unsigned short tileNumber);

  void GetTilesOfPmt(unsigned short pmtNumber, unsigned short* nTiles, unsigned short* tileNumbers);

  ClassDef(StBbcGeom,0)
};


inline unsigned short StBbcGeom::PmtOfTile(unsigned short tile){return mPMT[tile-1];}

#endif