HighResTimingData.hpp

Go to the documentation of this file.

#ifndef __HIGHRESTIMINGDATA_HPP__
#define __HIGHRESTIMINGDATA_HPP__

#include "ChanEvent.hpp"
#include "Globals.hpp"
#include "Trace.hpp"

class HighResTimingData {
public:
    HighResTimingData() {};
    virtual ~HighResTimingData() {};

    HighResTimingData(ChanEvent *chan) {
        trace_ = &chan->GetTrace();
        highResTime_ = chan->GetHighResTime()*1e9;
        pixieEnergy_ = chan->GetEnergy();

        snr_ =
            20*log10(trace_->GetValue("maxval") /
                     trace_->GetValue("sigmaBaseline"));

        if(!std::isnan(trace_->GetValue("maxval")) &&
           !std::isnan(trace_->GetValue("phase")) &&
           !std::isnan(trace_->GetValue("tqdc")) &&
           !std::isnan(trace_->GetValue("sigmaBaseline")) ) {
            isValidData_ = true;
        }else
            isValidData_ = false;
    };

    double CalcEnergy(const double &tof, const double &z0) {
        return((0.5*Globals::get()->neutronMass()*
        pow((z0/tof)/Globals::get()->speedOfLight(), 2)));
    }

    bool GetIsValidData() const { return isValidData_; };
    double GetAveBaseline() const { return trace_->GetValue("baseline"); };
    double GetDiscrimination() const { return trace_->GetValue("discrim"); };
    double GetHighResTime() const { return highResTime_; };
    double GetMaximumPosition() const { return trace_->GetValue("maxpos"); }
    double GetMaximumValue() const { return trace_->GetValue("maxval"); };
    int GetNumAboveThresh() const { return trace_->GetValue("numAboveThresh");; };
    double GetPhase() const { return trace_->GetValue("phase")*
                        (Globals::get()->clockInSeconds()*1e+9); };
    double GetPixieEnergy() const { return(pixieEnergy_); };
    double GetSignalToNoiseRatio() const { return snr_; };
    double GetStdDevBaseline() const { return trace_->GetValue("sigmaBaseline"); };
    Trace* GetTrace() const { return trace_; };
    double GetTraceQdc() const { return trace_->GetValue("tqdc") /
        Globals::get()->qdcCompression(); };
    double GetWalk() const { return trace_->GetValue("walk"); };
    double GetWalkCorrectedTime() const
        { return highResTime_ - trace_->GetValue("walk"); };

    void SetDataValid(const bool &val) { isValidData_ = val; };

private:
    Trace *trace_; 

    bool isValidData_; 
    double highResTime_; 
    double snr_; 
    double pixieEnergy_; 
};

typedef std::map<TimingDefs::TimingIdentifier, HighResTimingData> TimingMap;
#endif // __HIGHRESTIMINGDATA_HPP__