TModuleDecoderTimestamp.cc

Go to the documentation of this file.

/**
 * @file    TModuleDecoderTimestamp.cc
 * @brief
 * @author  Kodai Okawa <okawa@cns.s.u-tokyo.ac.jp>
 * @date    2022?
 * @note    last modified: 2024-08-23 20:59:52
 * @details
 */
 
#include "TModuleDecoderTimestamp.h"
 
#include <TRawDataSimple.h>
 
using art::crib::TModuleDecoderTimestamp;
 
typedef art::TRawDataSimple<ULong64_t> TimestampRaw_t;
 
ClassImp(TModuleDecoderTimestamp);
 
TModuleDecoderTimestamp::TModuleDecoderTimestamp()
    : TModuleDecoder(kID, TimestampRaw_t::Class()) {
    fHitData = new TObjArray;
}
 
TModuleDecoderTimestamp::~TModuleDecoderTimestamp() {
    delete fHitData;
    fHitData = nullptr;
}
 
Int_t TModuleDecoderTimestamp::Decode(char *buf, const Int_t &size, TObjArray *seg) {
    UShort_t *evtdata = reinterpret_cast<UShort_t *>(buf);
    UInt_t evtsize = size / sizeof(UShort_t);
    Int_t igeo, ich;
    ULong64_t ts;
    TimestampRaw_t *data;
 
    // clear old hits
    fHitData->Clear();
 
    for (UInt_t i = 0; i < evtsize; i++) {
        igeo = 0;
        ich = 0;
 
        ts = (ULong64_t)evtdata[3] << 48;
        ts += (ULong64_t)evtdata[2] << 32;
        ts += (ULong64_t)evtdata[1] << 16;
        ts += (ULong64_t)evtdata[0];
 
        // check if the data object exists.
        if (fHitData->GetEntriesFast() <= ich || !fHitData->At(ich)) {
            // if no data object is available, create one
            TimestampRaw_t *obj = static_cast<TimestampRaw_t *>(this->New());
            obj->SetSegInfo(seg->GetUniqueID(), igeo, ich);
            fHitData->AddAtAndExpand(obj, ich);
            seg->Add(obj);
        }
 
        data = static_cast<TimestampRaw_t *>(fHitData->At(ich));
        data->Set(ts);
 
        fHitData->AddAt(nullptr, ich);
    }
 
    return 0;
}