CbmMcbm2018TofPar.h

Go to the documentation of this file.

/* Copyright (C) 2021 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Florian Uhlig [committer] */
 
// -------------------------------------------------------------------------
// -----                 CbmMcbm2018TofPar header file                 -----
// -----                 Created 09/09/18  by P.-A. Loizeau            -----
// -------------------------------------------------------------------------
 
#ifndef CBMMCBM2018TOFPAR_H
#define CBMMCBM2018TOFPAR_H
 
#include "FairParGenericSet.h"
#include "TArrayD.h"
#include "TArrayI.h"
 
class FairParIo;
class FairParamList;
 
 
class CbmMcbm2018TofPar : public FairParGenericSet {
 
 public:
  CbmMcbm2018TofPar(const char* name = "CbmMcbm2018TofPar", const char* title = "Tof unpacker parameters",
                    const char* context = "Default");
 

  virtual ~CbmMcbm2018TofPar();

  virtual void clear();
 
  void putParams(FairParamList*);
  Bool_t getParams(FairParamList*);
 
  static constexpr UInt_t GetNbByteMessage() { return kuBytesPerMessage; }
 
  Int_t Get4ChanToPadiChan(UInt_t uChannelInFee);
  Int_t PadiChanToGet4Chan(UInt_t uChannelInFee);
 
  Int_t ElinkIdxToGet4Idx(UInt_t uElink);
  Int_t Get4IdxToElinkIdx(UInt_t uGet4);
  Int_t ElinkIdxToGet4IdxA(UInt_t uElink);
  Int_t Get4IdxToElinkIdxA(UInt_t uGet4);
 
  static constexpr UInt_t GetNrOfPadiThrCodes() { return kuNbPadiThrCodes; }
  Double_t GetPadiThresholdVal(UInt_t uCode);
 
  //   static constexpr UInt_t GetNrOfChannelsPerGet4() { return kuNbChannelsPerGet4; }
  //   static constexpr UInt_t GetNrOfGet4PerFee()      { return kuNbGet4PerFee; }
  static constexpr UInt_t GetNrOfFeePerGbtx() { return kuNbFeePerGbtx; }
  static constexpr UInt_t GetNrOfGbtxPerGdpb() { return kuNbGbtxPerGdpb; }
  static constexpr UInt_t GetNrOfChannelsPerFee() { return kuNbChannelsPerFee; }
  static constexpr UInt_t GetNrOfChannelsPerGbtx() { return kuNbChannelsPerGbtx; }
  static constexpr UInt_t GetNrOfChannelsPerGdpb() { return kuNbChannelsPerGdpb; }
  static constexpr UInt_t GetNrOfGet4PerGbtx() { return kuNbGet4PerGbtx; }
  static constexpr UInt_t GetNrOfGet4PerGdpb() { return kuNbGet4PerGdpb; }
  static constexpr UInt_t GetNrOfFeePerGdpb() { return kuNbFeePerGdpb; }
  inline UInt_t GetNumberOfChannels() { return kuNbChannelsPerGdpb * fiNrOfGdpb; }

  inline Int_t GetNrOfGdpbs() { return fiNrOfGdpb; }
  inline Int_t GetGdpbId(Int_t i) { return fiGdpbIdArray[i]; }
  inline Int_t GetNrOfFeesPerGdpb() { return fiNrOfFeesPerGdpb; }
  inline Int_t GetNrOfGet4PerFee() { return fiNrOfGet4PerFee; }
  inline Int_t GetNrOfChannelsPerGet4() { return fiNrOfChannelsPerGet4; }
 
  inline Int_t GetNrOfGbtx() { return fiNrOfGbtx; }
  inline Int_t GetNrOfModules() { return fiNrOfModule; }
  inline Int_t GetNrOfRpc(Int_t i) { return fiNrOfRpc[i]; }
  inline Int_t GetRpcType(Int_t i) { return fiRpcType[i]; }
  inline Int_t GetRpcSide(Int_t i) { return fiRpcSide[i]; }
  inline Int_t GetModuleId(Int_t i) { return fiModuleId[i]; }
  inline std::vector<Int_t> GetRpcChUidMap() { return fviRpcChUId; }
 
  inline Int_t GetNbMsTot() { return fiNbMsTot; }
  inline Int_t GetNbMsOverlap() { return fiNbMsOverlap; }
  inline Double_t GetSizeMsInNs() { return fdSizeMsInNs; }
 
  inline Double_t GetStarTriggDeadtime(Int_t gdpb) { return fdStarTriggerDeadtime[gdpb]; }
  inline Double_t GetStarTriggDelay(Int_t gdpb) { return fdStarTriggerDelay[gdpb]; }
  inline Double_t GetStarTriggWinSize(Int_t gdpb) { return fdStarTriggerWinSize[gdpb]; }
  inline Double_t GetTsDeadtimePeriod() { return fdTsDeadtimePeriod; }
 
  inline bool CheckBmonComp(uint32_t uCompId) { return ((uCompId & 0xFFF0) == 0xABF0); }
  //inline bool CheckInnerComp(uint32_t uCompId) { return ((uCompId & 0xFFF0) == 0xBBC0); }
  inline bool CheckInnerComp(uint32_t uCompId) { return ((uCompId & 0xFF00) == 0xBB00); }
 
 private:
  void BuildChannelsUidMap();
  void BuildChannelsUidMapCbm(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapStar(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapBmon(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapBmon_2022(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapCern(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapCera(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapStar2(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapStar2A(UInt_t& uCh, UInt_t uGbtx);
  void BuildChannelsUidMapBuc(UInt_t& uCh, UInt_t uGbtx);

  static const uint32_t kuBytesPerMessage = 8;
  static const uint32_t kuNbChannelsPerGet4 = 4;
  static const uint32_t kuNbGet4PerFee      = 8;
  static const uint32_t kuNbFeePerGbtx      = 5;
  static const uint32_t kuNbGbtxPerGdpb     = 6;
  static const uint32_t kuNbChannelsPerFee  = kuNbChannelsPerGet4 * kuNbGet4PerFee;
  static const uint32_t kuNbChannelsPerGbtx = kuNbChannelsPerFee * kuNbFeePerGbtx;
  static const uint32_t kuNbChannelsPerGdpb = kuNbChannelsPerGbtx * kuNbGbtxPerGdpb;
  static const uint32_t kuNbGet4PerGbtx     = kuNbGet4PerFee * kuNbFeePerGbtx;
  static const uint32_t kuNbGet4PerGdpb     = kuNbGet4PerGbtx * kuNbGbtxPerGdpb;
  static const uint32_t kuNbFeePerGdpb      = kuNbFeePerGbtx * kuNbGbtxPerGdpb;
  const UInt_t kuGet4topadi[kuNbChannelsPerFee] = {
    3,  2,  1,  0,  7,  6,  5,  4,  11, 10, 9,  8,  15, 14, 13, 12,
    19, 18, 17, 16, 23, 22, 21, 20, 27, 26, 25, 24, 31, 30, 29, 28};  
 
  const UInt_t kuPaditoget4[kuNbChannelsPerFee] = {
    3,  2,  1,  0,  7,  6,  5,  4,  11, 10, 9,  8,  15, 14, 13, 12,
    19, 18, 17, 16, 23, 22, 21, 20, 27, 26, 25, 24, 31, 30, 29, 28};  
 
  const UInt_t kuElinkToGet4[kuNbGet4PerGbtx] = {27, 2,  7,  3,  31, 26, 30, 1,  33, 37, 32, 13, 9,  14,
                                                 10, 15, 17, 21, 16, 35, 34, 38, 25, 24, 0,  6,  20, 23,
                                                 18, 22, 28, 4,  29, 5,  19, 36, 39, 8,  12, 11};
  const UInt_t kuGet4ToElink[kuNbGet4PerGbtx] = {24, 7,  1,  3,  31, 33, 25, 2,  37, 12, 14, 39, 38, 11,
                                                 13, 15, 18, 16, 28, 34, 26, 17, 29, 27, 23, 22, 5,  0,
                                                 30, 32, 6,  4,  10, 8,  20, 19, 35, 9,  21, 36};
 
  const UInt_t kuElinkToGet4A[kuNbGet4PerGbtx] = {0,  16, 8,  17, 1,  18, 9,  19, 2,  20, 10, 21, 3,  22,
                                                  11, 23, 4,  24, 12, 25, 5,  26, 13, 27, 6,  28, 14, 29,
                                                  7,  30, 15, 31, 32, 33, 34, 35, 36, 37, 38, 39};
  const UInt_t kuGet4ToElinkA[kuNbGet4PerGbtx] = {0,  4,  8,  12, 16, 20, 24, 28, 2,  6,  10, 14, 18, 22,
                                                  26, 30, 1,  3,  5,  7,  9,  11, 13, 15, 17, 19, 21, 23,
                                                  25, 27, 29, 31, 32, 33, 34, 35, 36, 37, 38, 39};

  static const uint32_t kuNbPadiThrCodes        = 1024;  // 0x3FF + 1
  static const uint32_t kuNbThrMeasPoints       = 65;
  const UInt_t kuThrMeasCode[kuNbThrMeasPoints] = {
    0x000, 0x010, 0x020, 0x030, 0x040, 0x050, 0x060, 0x070, 0x080, 0x090, 0x0A0, 0x0B0, 0x0C0,
    0x0D0, 0x0E0, 0x0F0, 0x100, 0x110, 0x120, 0x130, 0x140, 0x150, 0x160, 0x170, 0x180, 0x190,
    0x1A0, 0x1B0, 0x1C0, 0x1D0, 0x1E0, 0x1F0, 0x200, 0x210, 0x220, 0x230, 0x240, 0x250, 0x260,
    0x270, 0x280, 0x290, 0x2A0, 0x2B0, 0x2C0, 0x2D0, 0x2E0, 0x2F0, 0x300, 0x310, 0x320, 0x330,
    0x340, 0x350, 0x360, 0x370, 0x380, 0x390, 0x3A0, 0x3B0, 0x3C0, 0x3D0, 0x3E0, 0x3F0, 0x3FF};
  const Double_t kdThrMeasVal[kuNbThrMeasPoints] = {
    -652.6, -631.2, -611.4, -590.6, -570.9, -550.0, -529.9, -509.4, -490.6, -469.5, -449.3, -428.5, -408.5,
    -388.2, -367.8, -347.2, -329.2, -308.2, -287.5, -266.8, -246.9, -226.0, -205.6, -185.0, -165.7, -144.9,
    -124.4, -103.8, -83.4,  -62.9,  -42.4,  -21.2,  -5.3,   15.5,   36.2,   56.8,   77.3,   97.8,   118.4,
    139.1,  158.7,  179.2,  199.7,  220.2,  240.8,  261.1,  281.7,  302.2,  321.3,  341.4,  362.0,  382.2,
    402.9,  422.8,  443.4,  463.7,  483.7,  503.7,  524.1,  544.3,  565.0,  585.0,  605.5,  626.0,  646.1};
  std::vector<Double_t> fvdPadiThrCodeToValue;
 
  Int_t fbMcbmTof2024 = 0;
 
  Int_t fiNrOfGdpb;       // Total number of GDPBs
  TArrayI fiGdpbIdArray;  // Array to hold the unique IDs for all Tof GDPBs
 
  Int_t fiNrOfFeesPerGdpb;      // Number of FEEs which are connected to one GDPB
  Int_t fiNrOfGet4PerFee;       // Number of GET4 chips which are connected to one FEB
  Int_t fiNrOfChannelsPerGet4;  // Number of channels per GET4
 
  Int_t fiNrOfGbtx;                     // Total number of Gbtx links
  Int_t fiNrOfModule;                   // Total number of Modules
  TArrayI fiNrOfRpc;                    // number of Rpcs connected to Gbtx link, i.e. 3 or 5
  TArrayI fiRpcType;                    // type of Rpcs connected to Gbtx link
  TArrayI fiRpcSide;                    // side of Rpcs connected to Gbtx link, i.e. 0 or 1
  TArrayI fiModuleId;                   // Module Identifier connected to Gbtx link, has to match geometry
  std::vector<Int_t> fviRpcChUId = {};  // UID/address for each channel, build from type, side and module
 
  Int_t fiNbMsTot;        // Total number of MS per link in TS
  Int_t fiNbMsOverlap;    // Number of overlap MS per TS
  Double_t fdSizeMsInNs;  // Size of the MS in ns, needed for MS border detection
 
  TArrayD fdStarTriggerDeadtime;  // STAR: Array to hold for each gDPB the deadtime between triggers in ns
  TArrayD
    fdStarTriggerDelay;  // STAR: Array to hold for each gDPB the Delay in ns to subtract when looking for beginning of coincidence of data with trigger window
  TArrayD fdStarTriggerWinSize;  // STAR: Array to hold for each gDPB the Size of the trigger window in ns
  Double_t
    fdTsDeadtimePeriod;  // Period (ns) in the first MS of each TS where events with missing triggers should be built using the overlap MS of previous TS (overlap events)
 
  ClassDef(CbmMcbm2018TofPar, 1);
};
 
class CbmMcbm2018BmonPar : public CbmMcbm2018TofPar {
 public:
  CbmMcbm2018BmonPar(const char* name = "CbmMcbm2018BmonPar", const char* title = "Bmon unpacker parameters",
                     const char* context = "Default")
    : CbmMcbm2018TofPar(name, title, context){};
 
 private:
  // just an alias for the Parameter container to allow two instances
  ClassDef(CbmMcbm2018BmonPar, 1);
};
#endif  // CBMMCBM2018TOFPAR_H