CbmMcbm2020TrdTshiftPar.h

Go to the documentation of this file.

/* Copyright (C) 2021 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Pascal Raisig [committer] */
 
#ifndef CbmMcbm2020TrdTshiftPar_H
#define CbmMcbm2020TrdTshiftPar_H
 
#include "FairParGenericSet.h"  // mother class
 
// STD
#include <map>
 
// ROOT
#include <TArrayD.h>
#include <TArrayI.h>
#include <TFile.h>
#include <TH2.h>
#include <TH3.h>
 
class CbmMcbm2020TrdTshiftPar : public FairParGenericSet {
public:
  CbmMcbm2020TrdTshiftPar(const char* name    = "CbmMcbm2020TrdTshiftPar",
                          const char* title   = "TRD timeshift unpacker parameters mCbm 2020",
                          const char* context = "Default");
 
  virtual ~CbmMcbm2020TrdTshiftPar();
 
  virtual void printparams();
 
  virtual void clear();
 
  void putParams(FairParamList*);
  Bool_t getParams(FairParamList*);
 
  bool GetTimeshifts(std::shared_ptr<TH2> timeshiftHisto);
 
  double GetNEvents(std::shared_ptr<TFile> mcbmanafile);
 
 
  std::vector<Int_t> GetTimeshiftsVec(size_t tsidx);
 
  std::map<size_t, std::vector<Int_t>>* GetTimeshiftsMap() { return &fmapTimeshifts; }
 
 
  std::shared_ptr<TH3> GetCalibHisto(std::shared_ptr<TFile> mcbmanafile);
 
  std::shared_ptr<TH2> GetCalibHisto(std::shared_ptr<TH3> calibbasehisto);
 
  Int_t GetDominantShift(std::shared_ptr<TH3> calibbasehisto, size_t itsidx, size_t ichannel);
 
  bool FillTimeshiftArray(std::shared_ptr<TFile> mcbmanafile);
 
private:
  const std::vector<std::string> pararraynames = {"nTimeslices", "TrdTimeshifts"};
 
  const size_t fgNchannels = 768;
  // In principle this could be retrieved from the standard parameters. But since we hopefully need these ugly parameters only for mCbm2020 it is hardcoded, because a handling via the parameters would require full linking of the standard Trd parameter classes
 
  const UInt_t fgMicroSliceLength = 1.280e6;
 
  TArrayI fTimeshifts = {};
 
  TArrayI fNtimeslices;
 
  std::map<size_t, std::vector<Int_t>> fmapTimeshifts = {};  
 
  std::vector<Int_t> fvecCurrentTimeshifts = {};  
 
public:
  ClassDef(CbmMcbm2020TrdTshiftPar, 2);
};
 
#endif  // CbmMcbm2020TrdTshiftPar_H