CbmTaskTofClusterizerParWrite.h

Go to the documentation of this file.

/* Copyright (C) 2023-2025 Facility for Antiproton and Ion Research in Europe, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Dominik Smith [committer], Pierre-Alain Loizeau, Norbert Herrmann */
 
#ifndef CBMTASKTOFCLUSTERIZERPARWRITE_H
#define CBMTASKTOFCLUSTERIZERPARWRITE_H 1
 
// TOF Classes and includes
// Geometry
class CbmTofDetectorId;
class CbmTofDigiPar;
class CbmTofDigiBdfPar;
class CbmTofCell;
 
#include "algo/detectors/tof/Clusterizer.h"
 
// ROOT Classes and includes
#include "Math/Rotation3D.h"
#include "Math/Vector3Dfwd.h"
 
// FAIR classes and includes
#include "FairTask.h"
 
// ROOT Classes and includes
class TClonesArray;
 
// C++ Classes and includes
#include <vector>
 
class CbmTaskTofClusterizerParWrite : public FairTask {
 
 public:
  static constexpr int kBmonAssignedSmType = 5;  
  static constexpr int kBmonAssignedSide   = 0;  

  CbmTaskTofClusterizerParWrite();

  CbmTaskTofClusterizerParWrite(const char* name, int32_t verbose = 1, bool writeDataInOut = true);
  virtual ~CbmTaskTofClusterizerParWrite();

  virtual InitStatus Init();

  virtual void SetParContainers();

  virtual void Exec(Option_t* /*option*/){};

  virtual void Finish(){};
 
  void SetCalMode(int32_t iMode) { fCalMode = iMode; }
  void PosYMaxScal(double val) { fPosYMaxScal = val; }
  void SetTotMax(double val) { fTotMax = val; }
  void SetTotMin(double val) { fTotMin = val; }
  void SetTotMean(double val) { fTotMean = val; }
  void SetMaxTimeDist(double val) { fMaxTimeDist = val; }
  void SetChannelDeadtime(double val) { fdChannelDeadtime = val; }
  void SetCalParFileName(TString CalParFileName) { fCalParFileName = CalParFileName; }
  void SetRunId(uint32_t runId) { fRunId = runId; }
 
  double GetTotMean() { return fTotMean; }
 
  void SwapChannelSides(bool bSwap) { fbSwapChannelSides = bSwap; }
  void SetDeadStrips(int32_t iDet, uint32_t ival);

  void SetRecoParOutputDir(const std::string& dirname) { fsRecoParOutputDir = dirname; }
 
 
 protected:
 private:
  bool bAddBeamCounterSideDigi = true;
  const int32_t nbClWalkBinX   = 50;  // was 100 (11.10.2018)
  const int32_t nbClWalkBinY   = 41;  // choose odd number to have central bin symmetric around 0

  CbmTaskTofClusterizerParWrite(const CbmTaskTofClusterizerParWrite&);
  CbmTaskTofClusterizerParWrite& operator=(const CbmTaskTofClusterizerParWrite&);
 
  // Functions common for all clusters approximations
  bool InitParameters();
  bool InitCalibParameter();

  bool InitAlgosTof();

  bool InitAlgosBmon();
 
  // ToF geometry variables
  CbmTofDetectorId* fTofId;
  CbmTofDigiPar* fDigiPar;
  CbmTofDigiBdfPar* fDigiBdfPar;
 
  //Calibration parameters
  std::vector<std::vector<std::vector<std::vector<double>>>> fvCPTOff;     //[nSMT][nRpc][nCh][nbSide]
  std::vector<std::vector<std::vector<std::vector<double>>>> fvCPTotGain;  //[nSMT][nRpc][nCh][nbSide]
  std::vector<std::vector<std::vector<std::vector<double>>>> fvCPTotOff;   //[nSMT][nRpc][nCh][nbSide]
  std::vector<std::vector<std::vector<std::vector<std::vector<double>>>>>
    fvCPWalk;                                               //[nSMT][nRpc][nCh][nbSide][nbWalkBins]
  std::vector<std::vector<std::vector<double>>> fvCPTOffY;  //[nSMT][nRpc][nBin]
  std::vector<std::vector<double>> fvCPTOffYBinWidth;       //[nSMT][nRpc]
  std::vector<std::vector<double>> fvCPTOffYRange;          //[nSMT][nRpc]
 
  std::vector<uint32_t> fvDeadStrips;  //[nbDet]
 
  // Calib
  int32_t fCalMode;
  uint32_t fRunId{0};  
 
  double fPosYMaxScal;
  double fTotMax;
  double fTotMin;
  double fTotOff;
  double fTotMean;
  double fMaxTimeDist;
  double fdChannelDeadtime;
 
  TString fCalParFileName;  // name of the file name with Calibration Parameters
  std::string fsRecoParOutputDir = ".";  // name of the output parameter directory
 
  double fdTOTMax;
  double fdTOTMin;
  double fdTTotMean;
 
  double fdMaxTimeDist;   // Isn't this just a local variable? Why make it global and preset?!?
  double fdMaxSpaceDist;  // Isn't this just a local variable? Why make it global and preset?!?
  double fdModifySigvel;
 
  bool fbSwapChannelSides;
 
  ClassDef(CbmTaskTofClusterizerParWrite, 1);
};
 
#endif  // CBMTASKTOFCLUSTERIZER_H