/*
 * CbmClusteringGeometry.h
 *
 *  Created on: Feb 22, 2012
 *      Author: kozlov
 */

#ifndef CBMCLUSTERINGGEOMETRY_H_
#define CBMCLUSTERINGGEOMETRY_H_

#include "FairTask.h"
#include "CbmMuchGeoScheme.h"

class CbmClusteringGeometry
{
public:
	CbmClusteringGeometry();
	CbmClusteringGeometry(Int_t nStation, Int_t nLayer, Bool_t nSide, Int_t nModule, CbmMuchGeoScheme* scheme);
	void CbmClusteringSetPad(Int_t nPad, Float_t x, Float_t y, Float_t dx, Float_t dy, Double_t phi1, Double_t phi2, Float_t r1, Float_t r2,
			   Int_t digiNum, UInt_t charge, Long64_t chID);            //Addition of a single pad
   virtual ~CbmClusteringGeometry();

   //void SetMuchModuleGeometryRectangular(Int_t nStation, Int_t nLayer, Bool_t nSide, Int_t nModule, CbmMuchGeoScheme* scheme);  //Old geometry version
   void SetMuchModuleGeometryRadial(Int_t nStation, Int_t nLayer, Bool_t nSide, Int_t nModule, CbmMuchGeoScheme* scheme);
   void SetMuchModuleGeometryRadialFast(Int_t nStation, Int_t nLayer, Bool_t nSide, Int_t nModule, CbmMuchGeoScheme* scheme);

   Int_t GetDetId()			const { return fDetId; }
   Int_t GetNPads()			const { return fNofPads; }
   Int_t GetAPadsNom()		const { return fNofActivePads;}
   void SetAPadsNom(Int_t nPads);
   void SetAPadsPlusOne();

   Float_t GetDx(Int_t iPad);
   Float_t GetDy(Int_t iPad);
   Float_t GetX0(Int_t iPad);
   Float_t GetY0(Int_t iPad);

   Int_t GetDigiNum(Int_t iPad);
   void SetDigiNum(Int_t iPad, Int_t iDigi);
   Int_t GetNeighborsNum(Int_t iPad);
   Int_t GetGoodNeighborsNum(Int_t iPad);
   Int_t GetNeighbor(Int_t iPad, Int_t iNeighbor);
   Long64_t GetPadID(Int_t iPad);
   UInt_t GetPadCharge(Int_t iPad);
   void SetPadCharge(Int_t iPad, UInt_t iCharge);
   Int_t GetPadByChannelId(Long64_t chId);
   Double_t GetPhi1(Int_t iPad);
   Double_t GetPhi2(Int_t iPad);
   Float_t GetR1(Int_t iPad);
   Float_t GetR2(Int_t iPad);
   vector<Int_t> GetNeighbors(Int_t iPad);
   Long64_t GetChannelID(Int_t iPad);

private:

   template <typename T1>
   T1 GetMin(T1& a, T1& b);
   template <typename T2>
   T2 GetMax(T2& a, T2& b);
   Bool_t SubEqual(Double_t x1, Double_t x2, Double_t l);

   Long64_t fNofPads;                       //Number of pads in module
   Int_t fDetId;                            //Detector Id of MuchModule
   Int_t fNofActivePads;                    //Number of charged pads in module

   map <Long64_t, Int_t> fPadByChannelId;
   struct PadInformation{                   //Description of single pad
	   Float_t fX, fY;
	   Float_t fDx, fDy;
	   Double_t fPhi1, fPhi2;
	   Float_t fR1, fR2;
	   Int_t fDigiNum;
	   UInt_t fCharge;
	   Int_t fNofNeighbors;
	   Int_t fNofGoodNeighbors;
	   vector<Int_t> fNeighbors;
	   Long64_t channelID;
	   Int_t nSector;
   };
   PadInformation* fPadList;

   CbmClusteringGeometry(const CbmClusteringGeometry&);
   CbmClusteringGeometry& operator=(const CbmClusteringGeometry&);
};

#endif