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 = static_cast<int>(DetectorId::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