CbmStsAlgoFindHitsOrtho.h

Go to the documentation of this file.

/* Copyright (C) 2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese [committer] */
 
#ifndef CBMSTSALGOFINDHITSORTHO_H
#define CBMSTSALGOFINDHITSORTHO_H 1
 
 
#include <Rtypes.h>
 
#include <cassert>
 
class TGeoHMatrix;
class CbmStsCluster;
class CbmStsHit;
class CbmStsParSensor;
 
 
class CbmStsAlgoFindHitsOrtho {
 
 public:
  CbmStsAlgoFindHitsOrtho();
 
 
  CbmStsAlgoFindHitsOrtho(const CbmStsAlgoFindHitsOrtho&) = delete;
 
 
  CbmStsAlgoFindHitsOrtho& operator=(const CbmStsAlgoFindHitsOrtho&) = delete;
 
 
  virtual ~CbmStsAlgoFindHitsOrtho(){};
 
 
  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, UInt_t nStripsF,
                UInt_t nStripsB, Double_t pitchF, Double_t pitchB, 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
  {
    assert(channel < fNofStripsF + fNofStripsB);
    return (channel < Double_t(fNofStripsF) ? 0 : 1);
  }
 
 
  std::pair<Int_t, Int_t> GetStrip(UInt_t channel) const;
 
 
  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 fNofStripsF   = 0;        
  UInt_t fNofStripsB   = 0;        
  Double_t fDx         = 0.;       
  Double_t fDy         = 0.;       
  Double_t fPitchF     = 0.;       
  Double_t fPitchB     = 0.;       
  Double_t fLorentzF   = 0.;       
  Double_t fLorentzB   = 0.;       
  TGeoHMatrix* fMatrix = nullptr;  
 
  std::vector<CbmStsHit>* fHits = nullptr;  
 
 
  ClassDef(CbmStsAlgoFindHitsOrtho, 1);
};
 
 
#endif