/** CbmMuchModuleGem.h
 *
 *@author  M.Ryzhinskiy <m.ryzhinskiy@gsi.de>
 *@version 1.0
 *@since   11.02.08
 *
 * This class holds geometry parameters of much modules
 *
 * The channel ID consists of:
 *   sector number  (0-2047),       bits 0-10
 *   channel number (0-1024),       bits 11-20
 *
 */

#ifndef CBMMUCHMODULEGEM_H
#define CBMMUCHMODULEGEM_H 1

#include "CbmMuchGeoScheme.h"
#include "CbmMuchModule.h"

#include "TObject.h"
#include "TObjArray.h"
#include "TVector3.h"
#include "TPolyLine.h"
#include "TClonesArray.h"
#include "TPave.h"

#include <vector>
#include <map>

using std::vector;
using std::multimap;
using std::pair;

// Length of the index of the corresponding volume
#define WL_SECTOR 2147483647LL //2047
#define WL_CHANNEL 1023LL

// Number of a start bit for each volume
#define SB_SECTOR 0LL
#define SB_CHANNEL 32LL //11


class CbmMuchSector;

class CbmMuchModuleGem : public CbmMuchModule
{

  public:

    /** Default constructor **/
    CbmMuchModuleGem();

    /** Standard constructor
     *@param detId     Detector ID
     *@param position  Location of the module center in global c.s. (all dimensions in [cm])
     *@param size      Size of the module (all dimensions in [cm])
     *@param cutRadius Radius of the cut (if any, otherwise = -1.) [cm]
     **/
    CbmMuchModuleGem(Int_t detId, TVector3 position, TVector3 size, Double_t cutRadius);

    /** Standard constructor
     *@param iStation  Station index
     *@param iLayer    Layer index
     *@param iSide     Defines side of the layer (0 - Front, 1 - Back)
     *@param iModule   Module index
     *@param position  Location of the module center in global c.s. (all dimensions in [cm])
     *@param size      Size of the module (all dimensions in [cm])
     *@param cutRadius Radius of the cut (if any, otherwise = -1.) [cm]
     **/
    CbmMuchModuleGem(Int_t iStation, Int_t iLayer, Bool_t iSide, Int_t iModule, TVector3 position, TVector3 size, Double_t cutRadius);

    /** Destructor **/
    virtual ~CbmMuchModuleGem();

    /**
     * Gets sector index for the given channel Id.
     * @param channelId   Channel Id.
     * @return            Sector index within the module.
     */
    static Int_t GetSectorIndex(Long64_t channelId) {
      return (channelId & (WL_SECTOR << SB_SECTOR)) >> SB_SECTOR;
    }
    /**
     * Gets channel index for the given channel Id.
     * @param channelId   Channel Id.
     * @return            Channel index within the sector.
     */
    static Int_t GetChannelIndex(Long64_t channelId) {
      return (channelId & (WL_CHANNEL << SB_CHANNEL)) >> SB_CHANNEL;
    }
    /**
     * Gets channel Id by sector index and channel index.
     * @param iSector    Sector index within the module.
     * @param iChannel   Channel index within the sector.
     * @return           Detector Id.
     */
    static Long64_t GetChannelId(Int_t iSector, Int_t iChannel) {
      return ((Long64_t)iSector << SB_SECTOR) | ((Long64_t)iChannel << SB_CHANNEL);
    }

    /** Gets sector by the given channel Id. **/
    CbmMuchSector* GetSector(Long64_t channelId);
    /** Gets sector by the given numbers of column and row in the grid. **/
    CbmMuchSector* GetSector(Int_t iGridColumn, Int_t iGridRow);
    /** Gets sector by the given coordinates in global c.s. */
    CbmMuchSector* GetSector(Double_t x, Double_t y);
    /** Gets pad by the given Id. */
    CbmMuchPad* GetPad(Long64_t channelId);
    /** Gets array of pads for this module. */
    vector<CbmMuchPad*> GetPads();
    Int_t GetNPads();

    /** Adds a given sector to the array.
     *@param  sector   CbmMuchSector which should be added to the array.**/
    void AddSector(CbmMuchSector* sector);

    virtual Bool_t InitModule();
    virtual void DrawModule(Color_t color);

    Int_t GetNSectorChannels() const { return fNSectorChannels; }
    void SetNSectorChannels(Int_t nSectorChannels) { fNSectorChannels = nSectorChannels; }

    void ClearDigis() { fDigis.clear(); }
    void AddDigi(Double_t time,Int_t id) { fDigis.insert(pair<Double_t,Int_t>(time,id)); }
    multimap<Double_t,Int_t> GetDigis() { return fDigis; } 
    
  private:
    Bool_t                 fUseModuleDesign;       // Whether to use module or monolithic design
    Double_t               fGridDx, fGridDy;       // Size of the grid cell (for fast search purpose)
    Int_t                  fGridCols;              // Number of columns in the grid
    Int_t                  fGridRows;              // Number of rows in the grid
    vector<vector<Int_t> > fGridIndices;           // 2D-vector of sector numbers
    multimap<Double_t,Int_t>  fDigis;              //!
    Int_t                  fNSectorChannels;       // Number of channels per "complete" sector

    Double_t GetGridCellY(Int_t iRow);
    Double_t GetGridCellX(Int_t iCol);
    Double_t GetInitX(CbmMuchSector* sector);
    Double_t GetInitY(CbmMuchSector* sector);
    /** Grid initialization. **/
    Bool_t InitGrid(Bool_t useModuleDesign);
    /** Adds to each sector info about neighbour sectors **/
    void InitNeighbourSectors();
    /** Adds to each pad info about neighbour pads **/
    void InitNeighbourPads();

    ClassDef(CbmMuchModuleGem,2);
};
#endif