CbmStsSensor.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                     CbmStsSensor header file                  -----
// -----                  Created 01/07/2008  by R. Karabowicz         -----
// -------------------------------------------------------------------------
 
#ifndef CBMSTSSENSOR_H
#define CBMSTSSENSOR_H 1
 
#include "TNamed.h"
 
#include <map>
#include <list>
#include <set>
#include <vector>
 
#include "TGeoManager.h" //GP
#include "TGeoNode.h"//GP
#include "TGeoVolume.h"//GP
 
class CbmStsSensorDigiPar;
 
 
class CbmStsSensor : public TNamed
{
 
 public:
 
  CbmStsSensor();
 
 
  CbmStsSensor(Int_t detId, Int_t iType, Double_t x0, Double_t y0,
         Double_t rotation, Double_t lx, Double_t ly, 
         Double_t dx, Double_t dy, Double_t stereoF, Double_t stereoB);
  CbmStsSensor(TString tempName, Int_t detId, Int_t iType, Double_t x0, Double_t y0, Double_t z0,
         Double_t rotation, Double_t lx, Double_t ly, 
         Double_t d, Double_t dx, Double_t dy, Double_t stereoF, Double_t stereoB);
  
 
  virtual ~CbmStsSensor();
 
  TString GetDetectorName() const { return fName.Data(); }
  Int_t    GetDetectorId() const { 
    return fDetectorId; }
  Int_t    GetSystemId()   const {  
    return ( fDetectorId & (15<<24) ) >> 24; }
  Int_t    GetStationNr()  const { 
    return ( fDetectorId & (255<<16) ) >> 16; }
  Int_t    GetSectorNr()   const {  // sector number within station
    return ( fDetectorId & (4095<<4) ) >> 4; }
  /*   Int_t    GetSensorNr()   const {  // sensor number within sector 
     return ( fDetectorId & (32767<<1) ) >> 1; } */
  Int_t    GetSensorNr()   const {  // sensor number within sector
    return ( fDetectorId & (7<<1) ) >> 1; }
  Int_t    GetType()       const { return fType; }
  Double_t GetX0()         const { return fX0; }
  Double_t GetY0()         const { return fY0; }
  Double_t GetZ0()         const { return fZ0; }
  Double_t GetRotation()   const { return fRotation; }
  Double_t GetLx()         const { return fLx; }
  Double_t GetLy()         const { return fLy; }
  Double_t GetD()          const { return fD;  }
  Double_t GetDx()         const { return fDx; }
  Double_t GetDy()         const { return fDy; }
  Double_t GetStereoF()    const { return fStereoF; }
  Double_t GetStereoB()    const { return fStereoB; }
  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;  }
  Double_t GetBackStripShift()  const { return fBackStripShift; }
  Double_t GetFrontStripShift() const { return fFrontStripShift; }
 
  Int_t GetChannel(Double_t x, Double_t y, Int_t iSide);
 
  Float_t GetChannelPlus(Double_t x, Double_t y, Int_t iSide);
 
 // Int_t GetFrontChannel(Double_t x, Double_t y, Double_t z);
 // Int_t GetBackChannel (Double_t x, Double_t y, Double_t z);
  Int_t GetFrontChannel(Double_t x, Double_t y, Double_t z, Double_t &dPitch); //AZ
  Int_t GetBackChannel (Double_t x, Double_t y, Double_t z, Double_t &dPitch); //AZ
 
  Bool_t Inside(Double_t x, Double_t y);
 
 
  Bool_t ActivateChannels(Int_t ipt, Double_t x, Double_t y);
 
 
  Int_t Intersect(Int_t iFStrip, Int_t iBStrip,
      std::vector<Double_t>& xCross, std::vector<Double_t>& yCross);
 
  Int_t Intersect(Int_t iFStrip, Int_t iBStrip, Double_t& xCross, Double_t& yCross, Double_t& zCross);
 
  Int_t IntersectClusters(Double_t fChan, Double_t bChan, Double_t& xCross, Double_t& yCross, Double_t& zCross);
 
  Int_t PointIndex(Int_t iFStrip, Int_t IBStrip);
 
 
  void Reset();
 
 
  void Print();
  
 
 
 private:
  TString GeoMaskSi[3][8];
  TString GeoMaskGem[6][4];
  TGeoVolume* TopVolume;
  Double_t PIR = 3.1415926536 / 180.; // Ruf
 
  //TString    fName;             // Station name
  Int_t    fDetectorId;   // Unique detector ID
  Int_t    fType;         // Sensor type 
  Double_t fX0, fY0, fZ0; // Coordinates of the sensor centre [cm]
  Double_t fRotation;     // Rotation angle in global c.m. [rad]
  Double_t fLx, fLy;      // Dimensions of the sensor [cm]
  Double_t fDx;           // Strip readout pitch or pixel size in x [cm]
  Double_t fDy;           // Pixel size in y [cm] for pixel sensor
  Double_t fStereoF;      // StereoF angle [rad] for strip sensor.
  Double_t fStereoB;      // StereoB angle [rad] for strip sensor.
  Double_t fTanF,fTanB,fSinRot,fCosRot;
  
  Double_t fD;            // thickness of the sensor [cm]
 
  Int_t fNChannelsFront;
  Int_t fNChannelsBack;
 
  Double_t fBackStripShift;
  Double_t fFrontStripShift;
 
  // Logically, this belongs to the HitFinder, but it is here
  // for performance reasons: to be executed once per sensor, not for
  // each MCPoint
  Float_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]
 
  Double_t fXSmearWidth;
  Double_t fZSmearSlope;
 
  Double_t fFrontLorentzShift;
  Double_t fBackLorentzShift;
 
  std::set<Int_t> fFrontActive;       
  std::set<Int_t> fBackActive;        
 
 
  std::map<std::pair<Int_t,Int_t>, Int_t > fTrueHits;      
 
 
 
  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& xint, Double_t& yint) 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(CbmStsSensor,1);
 
};
 
 
 
 
#endif