CbmStsAlgoFindClusters.h

Go to the documentation of this file.

/* Copyright (C) 2017-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese [committer] */
 
#ifndef CBMSTSALGOFINDCLUSTERS_H
#define CBMSTSALGOFINDCLUSTERS_H 1
 
 
#include <Rtypes.h>
 
#include <vector>
 
class CbmStsCluster;
class CbmStsDigi;
class CbmStsParModule;
 
 
class CbmStsAlgoFindClusters {
 
 public:
  typedef std::pair<const CbmStsDigi*, Long64_t> InputData;
 
 
  CbmStsAlgoFindClusters(){};
 
 
  virtual ~CbmStsAlgoFindClusters(){};
 
 
  Long64_t Exec(const std::vector<InputData>& input, std::vector<CbmStsCluster>& output, UInt_t address,
                UShort_t nChannels, UShort_t channelOffset, Double_t timeCutSigma, Double_t timeCutAbs,
                Bool_t connectEdge, const CbmStsParModule* modPar);
 
 
 private:
  Short_t ChanLeft(UShort_t channel)
  {
    if (fConnectEdge) return (channel == 0 ? fNofChannels - 1 : channel - 1);
    return channel - 1;
  }
 
 
  Short_t ChanRight(UShort_t channel)
  {
    if (fConnectEdge) return (channel == fNofChannels - 1 ? 0 : channel + 1);
    return channel + 1;
  }
 
 
  Bool_t CheckChannel(Short_t channel, Double_t time);
 
 
  void CreateCluster(UShort_t channel);
 
 
  Bool_t IsActive(Short_t channel)
  {
    if (channel < 0 || channel >= fNofChannels) return kFALSE;
    return fStatus[channel].first > -1;
  }
 
 
  Bool_t ProcessDigi(UShort_t channel, Double_t time, Int_t index);
 
 
 private:
  std::vector<std::pair<Long64_t, Double_t>> fStatus{1024, {-1, 0.}};  
 
  std::vector<CbmStsCluster>* fOutput = nullptr;  
 
  // Required parameters
  UInt_t fAddress                = 0;        
  UShort_t fNofChannels          = 0;        
  UShort_t fChannelOffset        = 0;        
  Double_t fTimeCutSig           = 0.;       
  Double_t fTimeCutAbs           = 0.;       
  Bool_t fConnectEdge            = kFALSE;   
  const CbmStsParModule* fModPar = nullptr;  
};
 
#endif /* CBMSTSALGOFINDCLUSTERS_H */