// $Id$
//
//    File: DBCALGeometry.h
// Created: Thu Nov 17 15:10:51 CST 2005
// Creator: gluexuser (on Linux hydra.phys.uregina.ca 2.4.20-8smp i686)
//

#ifndef _DBCALGeometry_
#define _DBCALGeometry_

#include <JANA/JObject.h>
#include <JANA/JFactory.h>
using namespace jana;

// create a single number channel id which is useful in algorithms
// if M L S are module layer sector the bit map looks like:
// MMMM MMMM LLLL SSSS

#define MODULE_SHIFT 8
#define MODULE_MASK 0xFF00
#define LAYER_SHIFT 4
#define LAYER_MASK 0x00F0
#define SECTOR_SHIFT 0
#define SECTOR_MASK 0X000F

#define PI 3.14159265358

// with this set one will utilize the default summing -- see DBCALGeometry.cc
// for a description of summed geometry.
//
// Author's note:  this is seen as a temporary feature to study effects of
// two different summing schemes.  A preprocessor macro is not the best way
// to change the functionality of code, but mcsmear doesn't use factories to
// provide objects so its parameters aren't helpful.

#define BCAL_SUM_CELL

class DBCALGeometry : public JObject {
  
public:
  
  JOBJECT_PUBLIC( DBCALGeometry );
  
  DBCALGeometry();
  
  enum End { kUpstream, kDownstream };
  
  static const int NBCALMODS=48;         ///< number of modules

  //the distinction between inner layers and outer layers is important, since only the inner layers have TDC readout
#ifdef BCAL_SUM_CELL
  static const int NBCALLAYSIN=3;        ///< number of readout layers in inner BCAL (first 6 SiPM layers)
  static const int NBCALLAYSOUT=1;       ///< number of readout layers in outer BCAL (outer 4 SiPM layers)
#else
  static const int NBCALLAYSIN=6;
  static const int NBCALLAYSOUT=4;
#endif
  static int NSUMLAYSIN[NBCALLAYSIN];          ///< number of radial SiPM layers summed for digitization in each inner readout layer
  static int NSUMLAYSOUT[NBCALLAYSOUT];        ///< number of radial SiPM layers summed for digitization in each outer readout layer
  static int NSUMSECSIN;         ///< for the inner layers, the number of SiPM that will be summed in the azimuthal direction
  static int NSUMSECSOUT;        ///< for the outer layer(s), the number of SiPM that will be summed in the azimuthal direction
  static int NBCALSECSIN;        ///<number of sectors in inner region
  static int NBCALSECSOUT;       ///<number of sectors in outer region
  static float BCALINNERRAD;     ///< innner radius of BCAL in cm
  static float BCAL_PHI_SHIFT;   ///< overall phi roation of BCAL in radians

  // Enter the index of the SiPM that designates the first
  // (counting radially outward) of the outer cells (default 7)
  static const int BCALMID=7;         ///< first outer layer (default 7)

  static float m_radius[11];
  static float fADC_radius[5];
  static float BCALMIDRAD;       ///< mid radius of BCAL in cm (boundary between inner and outer layers)
  static float BCALOUTERRAD;     ///< outer radius of BCAL in cm
  static float BCALFIBERLENGTH;  ///< BCAL Scintilator fiber lenth in cm
  static float GLOBAL_CENTER;    ///< center of BCAL in gloobal coordinate system
  
  static float ATTEN_LENGTH;     ///< attenuation length
  static float C_EFFECTIVE;      ///< speed of light in fibers 
  
  static bool summingOn() {
  
#ifdef BCAL_SUM_CELL
    return true;
#else
    return false;
#endif
  }

  ///these functions are about encoding/decoding module/layer/sector info in a cellId
  static int cellId( int module, int layer, int sector );  ///< This object can be used for the SiPM ID or for the fADC ID since they are defined in the same way (4 bits for sector then 4 bits for layer then 8 bits for module.)
  static int module( int cellId );  ///< This method can be used for the SiPM ID or for the fADC ID since they are defined in the same way
  static int layer( int cellId );   ///< This method can be used for the SiPM ID or for the fADC ID since they are defined in the same way
  static int sector( int cellId );  ///< This method can be used for the SiPM ID or for the fADC ID since they are defined in the same way

  ///these functions are about finding which readout cell contains a specific SiPM cell
  static int fADC_layer( int SiPM_cellId );
  static int fADC_sector( int SiPM_cellId );
  static int fADCId( int module, int SiPM_layer, int SiPM_sector );
  static int NSiPMs(int fADCId);

  ///these functions are about the physical location and dimensions of a readout cell
  static float phi( int fADC_cellId );
  static float phiSize( int fADC_cellId );  
  static float r( int fADC_cellId );
  static float rSize( int fADC_cellId );

  ///these are missing functions that fill in some previous gaps.
  static int fADCcellId_rphi( float r, float phi );  ///< Method to get the fADC cell ID from an (R, phi) combination.\n  R in cm and phi in radians.
  static int getglobalsector(int module, int sector);
  static int getsector(int globalsector);
  static int getmodule(int globalsector);

};

#endif // _DBCALGeometry_