PndGemSensor.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                     PndGemSensor header file                  -----
// -----                  Created 12/02/2009  by R. Karabowicz         -----
// -------------------------------------------------------------------------

#ifndef PNDGEMSENSOR_H
#define PNDGEMSENSOR_H 1

#include "TNamed.h"
#include "PndDetectorList.h"

#include <map>
#include <list>
#include <set>
#include <vector>


class PndGemSensor : public TNamed
{

 public:

  PndGemSensor();


  PndGemSensor(TString tempName, Int_t detId, Int_t iType, 
               Double_t x0, Double_t y0, Double_t z0,
               Double_t rotation, 
               Double_t innerRad, Double_t outerRad,
               Double_t d, 
               Double_t stripAngle0, Double_t stripAngle1,
               Double_t pitch0, Double_t pitch1);
  PndGemSensor(TString tempName, Int_t stationNr, Int_t sectorNr, Int_t iType, 
               Double_t x0, Double_t y0, Double_t z0,
               Double_t rotation, 
               Double_t innerRad, Double_t outerRad,
               Double_t d, 
               Double_t stripAngle0, Double_t stripAngle1,
               Double_t pitch0, Double_t pitch1);
  PndGemSensor(const PndGemSensor &tempSensor);

  virtual ~PndGemSensor();


  TString  GetDetectorName()          const { return fName.Data(); }
  void     SetDetectorId(Int_t stationNr,  Int_t sensorNr) {
    fDetectorId = kGEM << 27 | 0 << 21 | stationNr << 8 | sensorNr << 6; }
  Int_t    GetDetectorId()            const { 
    return fDetectorId; }
  Int_t    GetSystemId()   const { 
    return ( ( fDetectorId & (  31<<27) ) >> 27); }
  Int_t    GetStationNr()  const { 
    return ( ( fDetectorId & (8191<< 8) ) >>  8 ); }
  Int_t    GetSensorNr()   const {  // sensor number within station
    return ( ( fDetectorId & (   3<< 6) ) >>  6 ); }
  Int_t    GetType()                  const { return fType; }
  Double_t GetX0()                    const { return fPosition[0]; }
  Double_t GetY0()                    const { return fPosition[1]; }
  Double_t GetZ0()                    const { return fPosition[2]; }
  Double_t GetRotation()              const { return fRotation; }
  Double_t GetInnerRadius()           const { return fInnerRadius; }
  Double_t GetOuterRadius()           const { return fOuterRadius; }
  Double_t GetD()                     const { return fD;  }
  Double_t GetStripAngle(Int_t index) const { return fStripAngle[index]; }
  Double_t GetPitch(Int_t index)      const { return fPitch[index]; }
  Double_t GetSigmaX()                const { return fSigmaX; }
  Double_t GetSigmaY()                const { return fSigmaY; }
  Double_t GetSigmaXY()               const { return fSigmaXY; }
  Int_t    GetNChannels()             const { return fNChannelsFront+fNChannelsBack; }
  Int_t    GetNChannelsFront()        const { return fNChannelsFront; }
  Int_t    GetNChannelsBack ()        const { return fNChannelsBack;  }

  Int_t    GetSideChannels (Int_t si) const { if ( si==0 ) return fNChannelsFront; return fNChannelsBack;} // will return NChFront for si==0, and NChBack for anything else

  Int_t GetChannel(Double_t x, Double_t y, Int_t iSide);

  Int_t GetChannel2(Double_t x, Double_t y, Int_t iSide, Double_t& feeDist); 
  //Double_t GetChannel(Double_t x, Double_t y, Int_t iSide, Double_t& stripWidth);
 
  Bool_t Inside(Double_t x, Double_t y);
  Bool_t Inside(Double_t radius);


  //Bool_t ActivateChannels(Int_t ipt, Double_t x, Double_t y);  // not implemented

   Double_t GetStripOrientation(Double_t x, Double_t y, Int_t iSide);
   Double_t GetDistance(Int_t iSide, Double_t chan1, Double_t chan2);
   Int_t    GetDistance(Int_t iSide, Int_t chanMin, Int_t chanMax, Int_t chanTest);
   Double_t GetDistance2(Int_t iSide, Double_t chan1, Double_t chan2);
   
   // returns the part of sensor in which chan is located.
   // for most sensors it should return -1
   // only for the sensors where the channel order is broken
   Int_t    GetSensorPart(Int_t iSide, Int_t chan);

   Double_t GetMeanChannel(Int_t iSide, Double_t chan1, Double_t weight1, Double_t chan2, Double_t weight2);

   Int_t    GetNeighbours(Int_t iSide, Int_t iChan, Int_t& nChan1, Int_t& nChan2, Int_t& nChan3);

  //Int_t Intersect(Int_t iFStrip, Int_t iBStrip, std::vector<Double_t>& xCross, std::vector<Double_t>& yCross);  // not implemented

  Int_t Intersect(Double_t iFStrip, Double_t iBStrip, Double_t& xCross, Double_t& yCross, Double_t& zCross);

  Int_t Intersect(Double_t iFStrip, Double_t iBStrip, Double_t& xCross, Double_t& yCross, Double_t& zCross,
                  Double_t& dr, Double_t& dp);

  Int_t Intersect(Double_t iFStrip, Double_t iBStrip, Double_t& xCross, Double_t& yCross, Double_t& zCross,
                  Double_t& dx, Double_t& dy, Double_t& dr, Double_t& dp);

  //Int_t IntersectClusters(Double_t fChan, Double_t bChan, Double_t& xCross, Double_t& yCross, Double_t& zCross);  // not implemented

  //Int_t PointIndex(Int_t iFStrip, Int_t IBStrip);  // not implemented


  void Reset();


  void Print();
  


 private:

  Int_t    fDetectorId;   // Unique detector ID
  Int_t    fType;         // Sensor type 
  Double_t fPosition[3];  // Coordinates of the sensor centre [cm]
  Double_t fRotation;     // Rotation angle in global c.m. [rad]
  Double_t fInnerRadius;  // Inner radius of the sensor [cm]
  Double_t fOuterRadius;  // Outer radius of the sensor [cm]
  Double_t fD;            // thickness of the sensor [cm]
  Double_t fStripAngle[2];     // Strips angle
  Double_t fPitch[2];     // Strip readout pitch or pixel size in x/y

  Int_t fNChannelsFront;
  Int_t fNChannelsBack;

  // Logically, this belongs to the HitFinder, but it is here
  // for performance reasons: to be executed once per sensor, not for
  // each MCPoint
  Double_t fSigmaX;   // RMS in x, global c.s. [cm]
  Double_t fSigmaY;   // RMS in y, global c.s. [cm]
  Double_t fSigmaXY;  // Covariance in global c.s. [cm**2]

  Int_t FrontStripNumber(Double_t x, Double_t y) const;


  Int_t BackStripNumber(Double_t x, Double_t y) const;


  Bool_t IntCoord(Double_t x, Double_t y, Double_t z, 
                  Double_t& xint, Double_t& yint, Double_t& zint) const;

  
  Bool_t IsInside(Double_t xint, Double_t yint) const;




  ClassDef(PndGemSensor,1);

};




#endif