BmnTrigWaveDigit.h

Go to the documentation of this file.

#ifndef BMNTRIGWAVEDIGIT_H
#define BMNTRIGWAVEDIGIT_H
 
#include <iostream>
#include <vector>
 
using std::vector;
 
#include <TObject.h>
 
/**
 * class for TQDC data including waveform
 * can be initialized both from an std::vector and raw pointer to an array
 */
class BmnTrigWaveDigit : public TObject
{
  public:
    BmnTrigWaveDigit();
    BmnTrigWaveDigit(uint16_t iMod,
                     vector<int16_t>&& iValue,
                     Double_t trigTimestamp,
                     Double_t adcTimestamp,
                     Double_t time = 0);
 
    BmnTrigWaveDigit(uint16_t iMod,
                     int16_t* iValue,
                     uint32_t nVals,
                     Double_t trigTimestamp,
                     Double_t adcTimestamp,
                     Double_t time = 0);
 
    uint16_t GetMod() const { return fMod; }
 
    size_t GetNSamples() const
    {
        if (fValueI)
            return fNsmpl;   // backward compatibility
        else
            return fValue.size();
    }
 
    int32_t GetIntegral(size_t start = 0, size_t stop = SIZE_MAX) const
    {
        stop = std::min(stop, fValue.size());
        start = std::min(start, stop);
        int32_t spectra = 0;
        for (size_t i = start; i < stop; ++i)
            spectra += fValue[i];
        return spectra;
    }
 
    int16_t GetPeak(size_t start = 0, size_t stop = SIZE_MAX) const
    {
        stop = std::min(stop, fValue.size());
        start = std::min(start, stop);
        int16_t peak = INT16_MIN;
        for (size_t i = start; i < stop; ++i)
            if (fValue[i] > peak)
                peak = fValue[i];
        return peak;
    }
 
    //    int16_t *GetShortValue() const {
    //        return fValueI;
    //    }
 
    vector<int16_t>& GetValue() { return fValue; }
 
    Double_t GetAdcTimestamp() const { return fAdcTimestamp; }
 
    Double_t GetTrigTimestamp() const { return fTrigTimestamp; }
 
    Double_t GetTime() const { return fTime; }
 
    vector<Double_t>& TdcVector() { return fTdcTimes; }
 
    virtual ~BmnTrigWaveDigit();
 
  protected:
    uint16_t fMod;
    /**
     * ROOT schema evolution does not work as stated in their manual
     * So we have to hold old fields in order to transition from the old version.
     * https://root.cern/manual/io/ even says it works automatically,
     * actually it does not work neither automatically nor even by linkdef
     * pragma read function if the fields are not present in the current version.
     * And the custom Streamer method is impossible to make being called during
     * object deserialization from a tree
     */
    uint32_t fNsmpl;
    int16_t* fValueI;   //[fNsmpl]
    vector<int16_t> fValue;
    Double_t fTrigTimestamp;
    Double_t fAdcTimestamp;
    Double_t fTime;
    vector<Double_t> fTdcTimes;
 
    ClassDef(BmnTrigWaveDigit, 6);
};
 
#endif /* BMNTRIGWAVEDIGIT_H */