CbmStsUtils.h

Go to the documentation of this file.

/* Copyright (C) 2023 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Dario Ramirez [committer] */
 
#ifndef CBMSTSUTILS_H
#define CBMSTSUTILS_H
 
#include "CbmStsHit.h"
#include "TClonesArray.h"
 
#include <set>
 
class CbmStsParModule;

namespace cbm_sts_utils
{
  static const double kStsClock = 3.125;
  static const double kStsDx    = 6. / 1024;
  static const double kStsDy    = kStsDx / std::tan(7.5 * TMath::Pi() / 180);
  static const double kStsErrX  = kStsDx / sqrt(12);
  static const double kStsErrY  = kStsDy / sqrt(12);

  int32_t GetHitCluSizeF(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  int32_t GetHitCluSizeB(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  double GetHitCharge(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  double GetHitChargeB(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  double GetHitChargeF(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  double GetHitTimeB(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  double GetHitTimeF(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  double GetHitSizeAsy(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  double GetHitChargeAsy(CbmStsHit* hit = nullptr, TClonesArray* clusters = nullptr);

  struct HBinning {
    uint32_t n_of_bins{0};
    double x_min{0};
    double x_max{0};
  };

  std::pair<HBinning, HBinning> ChargeBinning(const CbmStsParModule& par_module, const uint32_t max_clu_size = 1);

  std::set<int> GetUnits(const std::vector<int32_t> addresses);

  std::vector<int32_t> GetUnitModules(const std::vector<int32_t> addresses, const uint32_t unit);
};  // namespace cbm_sts_utils
#endif