CbmStsAlgoFindHits.h

Go to the documentation of this file.

 
#ifndef CBMSTSALGOFINDHITS_H
#define CBMSTSALGOFINDHITS_H 1
 
 
#include <Rtypes.h>
 
class TGeoHMatrix;
class CbmStsCluster;
class CbmStsHit;
class CbmStsParSensor;
 
 
class CbmStsAlgoFindHits {
 
public:
  CbmStsAlgoFindHits();
 
 
  CbmStsAlgoFindHits(const CbmStsAlgoFindHits&) = delete;
 
 
  CbmStsAlgoFindHits& operator=(const CbmStsAlgoFindHits&) = delete;
 
 
  virtual ~CbmStsAlgoFindHits() {};
 
 
  Long64_t Exec(const std::vector<CbmStsCluster>& clustersF,
                const std::vector<CbmStsCluster>& clustersB,
                std::vector<CbmStsHit>& hits,
                UInt_t address,
                Double_t timeCutSig,
                Double_t timeCutAbs,
                Double_t dY,
                UInt_t nStrips,
                Double_t pitch,
                Double_t stereoF,
                Double_t stereoB,
                Double_t lorentzF,
                Double_t lorentzB,
                TGeoHMatrix* matrix);
 
 
private:
  void CreateHit(Double_t xLocal,
                 Double_t yLocal,
                 Double_t varX,
                 Double_t varY,
                 Double_t varXY,
                 const CbmStsCluster& clusterF,
                 const CbmStsCluster& clusterB,
                 UInt_t indexF,
                 UInt_t indexB,
                 Double_t du = 0.,
                 Double_t dv = 0.);
 
 
  void
  GetClusterPosition(Double_t ClusterCentre, Double_t& xCluster, Int_t& side);
 
 
  Int_t GetSide(Double_t channel) const {
    return (channel < Double_t(fNofStrips) ? 0 : 1);
  }
 
 
  std::pair<Int_t, Int_t> GetStrip(UInt_t channel) const;
 
 
  Bool_t Intersect(Double_t xF,
                   Double_t exF,
                   Double_t xB,
                   Double_t exB,
                   Double_t& x,
                   Double_t& y,
                   Double_t& varX,
                   Double_t& varY,
                   Double_t& varXY);
 
 
  Int_t IntersectClusters(const CbmStsCluster& clusterF,
                          const CbmStsCluster& clusterB,
                          UInt_t indexF,
                          UInt_t indexB);
 
 
  Bool_t IsInside(Double_t x, Double_t y);
 
 
private:
  UInt_t fAddress      = 0;        
  Double_t fTimeCutSig = 0.;       
  Double_t fTimeCutAbs = 0.;       
  UInt_t fNofStrips    = 0;        
  Double_t fDx         = 0.;       
  Double_t fDy         = 0.;       
  Double_t fPitch      = 0.;       
  Double_t fStereoF    = 0.;       
  Double_t fStereoB    = 0.;       
  Double_t fLorentzF   = 0.;       
  Double_t fLorentzB   = 0.;       
  TGeoHMatrix* fMatrix = nullptr;  
  Double_t fTanStereoF = 0.;       
  Double_t fTanStereoB = 0.;       
  Double_t fErrorFac   = 0.;       
 
  std::vector<CbmStsHit>* fHits = nullptr;  
 
 
  ClassDef(CbmStsAlgoFindHits, 1);
};
 
 
#endif