CbmStsRecoModule.h

Go to the documentation of this file.

/* Copyright (C) 2020-2021 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese [committer] */
 
#ifndef CBMSTSRECOMODULE_H
#define CBMSTSRECOMODULE_H 1
 
 
#include "CbmStsCluster.h"
#include "CbmStsHit.h"
 
#include <TNamed.h>
 
#include <mutex>
 
class TGeoHMatrix;
class CbmStsAlgoAnaCluster;
class CbmStsAlgoFindClusters;
class CbmStsAlgoFindHits;
class CbmStsAlgoFindHitsOrtho;
class CbmStsHit;
class CbmStsDigi;
class CbmStsModule;
class CbmStsParModule;
class CbmStsParSensor;
class CbmStsParSensorCond;
 
 
class CbmStsRecoModule : public TNamed {
 
 public:
  struct Timings {
    double timeSortDigi    = 0;
    double timeCluster     = 0;
    double timeSortCluster = 0;
    double timeHits        = 0;
  };
 
  CbmStsRecoModule();
 
 
  CbmStsRecoModule(CbmStsModule* setupModule, const CbmStsParModule& parModule, const CbmStsParSensor& parSensor,
                   Double_t lorentzShiftF, Double_t lorentzShiftB);
 
 
  CbmStsRecoModule(const CbmStsRecoModule&) = delete;
 
 
  CbmStsRecoModule& operator=(const CbmStsRecoModule&) = delete;
 
 
  virtual ~CbmStsRecoModule();
 
 
  void AddDigiToQueue(const CbmStsDigi* digi, Int_t digiIndex);
 
 
  const std::vector<CbmStsCluster>& GetClustersF() const { return fClustersF; }
 
 
  const std::vector<CbmStsCluster>& GetClustersB() const { return fClustersB; }
 
 
  const std::vector<CbmStsHit>& GetHits() const { return fHits; }
 
  Timings GetTimings() const { return fTimings; }
 
 
  void Reconstruct();
 
  void SortDigis();
 
  void FindClusters();
 
  void SortClusters();
 
  void FindHits();
 
 
  void Reset();
 
 
  TGeoHMatrix* getMatrix() { return fMatrix; }
 
 
  std::string ToString() const;
 
  void SetTimeCutClustersAbs(Double_t value) { fTimeCutClustersAbs = value; }
 
 
  void SetTimeCutClustersSig(Double_t value) { fTimeCutClustersSig = value; }
 
 
  void SetTimeCutDigisAbs(Double_t value) { fTimeCutDigisAbs = value; }
 
 
  void SetTimeCutDigisSig(Double_t value) { fTimeCutDigisSig = value; }
 
 
 private:
  void Init();
 
 
 private:
  // --- Algorithms
  CbmStsAlgoAnaCluster* fClusterAna        = nullptr;  
  CbmStsAlgoFindClusters* fClusterFinder   = nullptr;  
  CbmStsAlgoFindHits* fHitFinder           = nullptr;  
  CbmStsAlgoFindHitsOrtho* fHitFinderOrtho = nullptr;  
  std::mutex fLock{};
 
  // --- Parameters
  CbmStsModule* fSetupModule        = nullptr;  
  const CbmStsParModule* fParModule = nullptr;  
  const CbmStsParSensor* fParSensor = nullptr;  
  Double_t fDyActive                = 0.;       
  UInt_t fNofStripsF                = 0;        
  UInt_t fNofStripsB                = 0;        
  Double_t fStripPitchF             = 0.;       
  Double_t fStripPitchB             = 0.;       
  Double_t fStereoFront             = 0.;       
  Double_t fStereoBack              = 0.;       
  TGeoHMatrix* fMatrix              = nullptr;  
  Double_t fLorentzShiftF           = 0.;       
  Double_t fLorentzShiftB           = 0.;       
 
  // --- Data
  std::vector<std::pair<const CbmStsDigi*, Long64_t>> fDigisF{};  
  std::vector<std::pair<const CbmStsDigi*, Long64_t>> fDigisB{};  
  std::vector<CbmStsCluster> fClustersF{};                        
  std::vector<CbmStsCluster> fClustersB{};                        
  std::vector<CbmStsHit> fHits{};                                 
 
  // --- Settings
  Double_t fTimeCutDigisSig    = 3.;      
  Double_t fTimeCutDigisAbs    = -1.;     
  Double_t fTimeCutClustersSig = 4.;      
  Double_t fTimeCutClustersAbs = -1.;     
  Bool_t fConnectEdgeFront     = kFALSE;  
  Bool_t fConnectEdgeBack      = kFALSE;  
 
  // --- Time measurement
  Timings fTimings;
 
 
  ClassDef(CbmStsRecoModule, 1);
};
 
#endif /* CBMSTSRECOMODULE_H */