CbmTrdHitMC.h

Go to the documentation of this file.

/* Copyright (C) 2018-2022 Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Alexandru Bercuci [committer] */
 
#ifndef CBMTRDHITMC_H
#define CBMTRDHITMC_H
 
#include "CbmTrdCluster.h"
#include "CbmTrdHit.h"
#include "CbmTrdPoint.h"
 
#include <string>  // for ToString
#include <vector>  // for fTrdPoints
 
class CbmTrdDigi;
class CbmTrdHitMC : public CbmTrdHit {
 public:
  enum eCbmTrdHitMCshape
  {
    kRT = 0,  // open left - open right shape
    kRR,      // open left - closed right shape
    kTT,      // closed left - open right shape
    kTR       // closed left - closed right shape
  };
 
  CbmTrdHitMC();
 
  CbmTrdHitMC(const CbmTrdHit& hit);
 
  virtual ~CbmTrdHitMC();
 
  void AddCluster(const CbmTrdCluster* c);
 
  size_t AddPoint(const CbmTrdPoint* p, double t, int id);
  size_t AddSignal(const CbmTrdDigi* d, uint64_t t0);
 
  std::string GetErrorMsg() const { return fErrMsg; }
 
  const CbmTrdPoint* GetPoint(uint idx = 0) const;
 
  double GetSignal(uint idx = 0) const;
 
  size_t GetNPoints() const { return fTrdPoints.size(); }
 
  size_t GetNSignals() const { return fTrdSignals.size(); }
 
  eCbmTrdHitMCshape GetClShape() const;
 
  double GetDx() const;
 
  double GetSx() const;
 
  double GetDy() const;
 
  double GetSy() const;
 
  double GetDt() const;
 
  void SetErrorMsg(std::string msg) { fErrMsg = msg; }
 
  size_t PurgeSignals();
 
  virtual std::string ToString() const;
 
 private:
  CbmTrdHitMC(const CbmTrdHitMC&) = default;
  CbmTrdHitMC& operator=(const CbmTrdHitMC&) = default;
 
  std::string fErrMsg                             = "";  //< error message from the QA task
  std::vector<std::pair<double, int>> fTrdSignals = {};  //< list of signal/time in cluster
  std::vector<std::tuple<CbmTrdPoint, double, int>> fTrdPoints =
    {};                     //< list of MC points together with the event time and particle PDG code producing them
  CbmTrdCluster fCluster;   //< data from the cluster
  static int fSx[5][2];     //< x error parametrization as function of cluster size and incident direction
  ClassDef(CbmTrdHitMC, 1)  // Hit to MC point data correlation
};
 
#endif