BmnTrigRaw2Digit.h

Go to the documentation of this file.

#ifndef BMNTRIGRAWTODIGIT_H
#define BMNTRIGRAWTODIGIT_H
// C++
#include <fstream>
#include <iostream>
#include <map>
#include <regex>
#include <sstream>
#include <unordered_map>
#include <vector>
// ROOT
#include <TCanvas.h>
#include <TClonesArray.h>
#include <TH1I.h>
#include <TMath.h>
#include <TPRegexp.h>
#include <TString.h>
#include <TTree.h>
// BmnRoot
#include "BmnAliases.h"
#include "BmnEnums.h"
#include "BmnTDCDigit.h"
#include "BmnTQDCADCDigit.h"
#include "BmnTrigDigit.h"
#include "BmnTrigWaveDigit.h"
 
#define CHANNEL_COUNT_MAX 72u
#define TDC72_CHANNEL_COUNT 72u
#define TQDC16_CHANNEL_COUNT 16u
#define TDC_BIN_COUNT 1024
#define HPTIMEBIN 0.02344
#define INVHPTIMEBIN 42.6666
#define ADC_CLOCK_TQDC16VS 8.0   // nano seconds
#define ADC_CLOCK 12.0
#define ADC_CLOCK_OLD 12.5
#define TDC_CLOCK 24.0
#define TDC_CLOCK_OLD 25.0
const UShort_t kNCHANNELS = 8;   // number of channels in one HPTDC
 
using std::map;
using std::vector;
 
struct BmnTrigChannelData
{
    TString name;
    UInt_t module;
    UInt_t serial;
    UShort_t slot;
    UShort_t channel;
    Bool_t isNegative;
    TClonesArray* branchArrayPtr;
};
 
struct BmnTrigParameters
{
    TString name;
    UInt_t BoardSerial;
    UInt_t CrateSerial;
    UShort_t slot;
    UShort_t ChannelCount;
    Double_t INL[CHANNEL_COUNT_MAX][TDC_BIN_COUNT];
    UShort_t ChannelMap[CHANNEL_COUNT_MAX];
    Bool_t NegativeMap[CHANNEL_COUNT_MAX];
    Bool_t IsT0;
    TClonesArray* branchArrayPtr[CHANNEL_COUNT_MAX];
    Double_t t[CHANNEL_COUNT_MAX];
    BmnTrigParameters();
};
 
class BmnTrigRaw2Digit
{
  public:
    BmnTrigRaw2Digit() {}
    BmnTrigRaw2Digit(TString mappingFile, TString INLFile, TTree* digiTree = NULL);
    // BmnTrigRaw2Digit(TString mappingFile, TString INLFile, FairRootManager* frm = NULL, bool saveOutput = true);
    BmnTrigRaw2Digit(TString mappingFile, TString INLFile, std::function<void(TString, TClonesArray*)> BranchRegFun);
 
    virtual ~BmnTrigRaw2Digit()
    {
        for (TClonesArray* ar : trigArrays)
            delete ar;
        //        for (auto el : fMap)
        //            delete el;
        fMap.clear();
        for (auto& el : fPlacementMap) {
            delete el.second;
        }
        fPlacementMap.clear();
    }
 
    vector<BmnTrigChannelData>* GetMap() { return &fMap; }
    T0ChMap& T0Map() { return fT0Map; }
    NameChMap& GetTrcMap() { return fTrcBitMap; }
    Name2SerCh& GetScalerMap() { return name2scaler_map; }
    SerCh2Name& GetScaler2NameMap() { return scalers2name_map; }
    map<std::pair<uint32_t, uint16_t>, uint16_t>& GetScaler2IdxMap() { return scalers2index_map; }
    map<uint16_t, uint16_t>& GetTrcIdx2ScalerIdxMap() { return trcIdx2scalerIdx; }
    NameChMap& GetScalerName2IdxMap() { return scalers_name2index_map; }
    vector<string>& GetScalerIdx2NameMap() { return scalers_index2name_map; }
 
    BmnStatus FillEvent(TClonesArray* tdc, unordered_map<UInt_t, Long64_t>& tsMap);
    BmnStatus FillEvent(TClonesArray* tdc, TClonesArray* adc, unordered_map<UInt_t, Long64_t>& tsMap);
    BmnStatus ReadINLFromFile(BmnTrigParameters* par);
    BmnStatus ReadChannelMap(TString mappingFile);
    BmnStatus ReadPlacementMap(TString mappingFile);
    BmnStatus ClearArrays();
 
    vector<TClonesArray*>* GetTrigArrays() { return &trigArrays; }
 
    UInt_t GetT0Serial()
    {
        //        if (fTrcSerial > 0)
        //            return fTrcSerial;
        for (auto& itPl : fPlacementMap)
            if (itPl.second->IsT0)
                return itPl.second->CrateSerial;
        return 0;
    }
 
    uint32_t GetTrcSerial() { return fTrcSerial; }
 
    void SetSetup(BmnSetup stp) { fSetup = stp; }
 
  private:
    TString Serial2FileName(TString boardName, UInt_t serial)
    {
        UShort_t h = serial >> 16;
        UShort_t l = serial & 0x0000FFFF;
        TString inlFileName = Form("%s-%04X-%04X.ini", boardName.Data(), h, l);
        return inlFileName;
    }
 
    UInt_t ChanCntByName(TString channelCountStr)
    {
        //    regex reBoardName("(\\D+)(\\d+)(.*)");
        TPRegexp reBoardName("(\\D+)(\\d+)(.*)");
        //        string channelCountStr = name;
        UInt_t channelCount = CHANNEL_COUNT_MAX;
        if (reBoardName.MatchB(channelCountStr)) {
            //        if (regex_match(name, reBoardName)){
            //            channelCountStr = regex_replace(name, reBoardName, "$2");
            reBoardName.Substitute(channelCountStr, "$2");
            channelCount = strtoul(channelCountStr.Data(), nullptr, 10);
            //            channelCount = strtoul(channelCountStr.c_str(), nullptr, 10);
        }
        return channelCount;
    }
 
    void ProcessWave(int16_t* iValue, const uint16_t& nVals, Bool_t& isNeg);
 
    Double_t GetBaseline(int16_t* iValue, const uint16_t& nVals);
 
    map<PlMapKey, BmnTrigParameters*> fPlacementMap;
    NameChMap fTrcBitMap;
    Name2SerCh name2scaler_map;
    T0ChMap fT0Map;
    vector<BmnTrigChannelData> fMap;
    BmnSetup fSetup;
    TString fMapFileName;
    vector<TClonesArray*> trigArrays;
 
    unique_ptr<TH1I> fWorkHist;
    uint32_t fTrcSerial;
 
    SerCh2Name scalers2name_map;
    // map: <serial, channel> -> index in the list (inside spillHeader->ScalersVec(), and MSC matrix)
    map<std::pair<uint32_t, uint16_t>, uint16_t> scalers2index_map;
    map<uint16_t, uint16_t> trcIdx2scalerIdx;
    // map: name -> index in the list
    NameChMap scalers_name2index_map;
    vector<string> scalers_index2name_map;
 
    ClassDef(BmnTrigRaw2Digit, 1);
};
 
#endif /* BMNTRIGRAWTODIGIT_H */