TRestEventRateAnalysisProcess.cxx

/*************************************************************************
 * This file is part of the REST software framework.                     *
 *                                                                       *
 * Copyright (C) 2016 GIFNA/TREX (University of Zaragoza)                *
 * For more information see http://gifna.unizar.es/trex                  *
 *                                                                       *
 * REST is free software: you can redistribute it and/or modify          *
 * it under the terms of the GNU General Public License as published by  *
 * the Free Software Foundation, either version 3 of the License, or     *
 * (at your option) any later version.                                   *
 *                                                                       *
 * REST is distributed in the hope that it will be useful,               *
 * but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the          *
 * GNU General Public License for more details.                          *
 *                                                                       *
 * You should have a copy of the GNU General Public License along with   *
 * REST in $REST_PATH/LICENSE.                                           *
 * If not, see http://www.gnu.org/licenses/.                             *
 * For the list of contributors see $REST_PATH/CREDITS.                  *
 *************************************************************************/
 
#include "TRestEventRateAnalysisProcess.h"
 
#include "TRestDataBase.h"
using namespace std;
 
ClassImp(TRestEventRateAnalysisProcess);
 
TRestEventRateAnalysisProcess::TRestEventRateAnalysisProcess() { Initialize(); }
 
TRestEventRateAnalysisProcess::~TRestEventRateAnalysisProcess() {}
 
void TRestEventRateAnalysisProcess::Initialize() {
    SetSectionName(this->ClassName());
 
    fEvent = NULL;
 
    fSingleThreadOnly = true;
 
    fFirstEventTime = -1;
    fPreviousEventTime.clear();
}
 
void TRestEventRateAnalysisProcess::InitProcess() {
    if (fRunInfo != NULL && fRunInfo->GetStartTimestamp() != 0) {
        fFirstEventTime = fRunInfo->GetStartTimestamp();
    } else {
        fFirstEventTime = -1;
    }
 
    // if is run under single thread mode, we add rate observables
    fRateAnalysis = GetNumberOfParallelProcesses() <= 1;
}
 
TRestEvent* TRestEventRateAnalysisProcess::ProcessEvent(TRestEvent* inputEvent) {
    fEvent = inputEvent;
 
    if (fFirstEventTime == -1) fFirstEventTime = fEvent->GetTime();
 
    Double_t secondsFromStart = fEvent->GetTime() - fFirstEventTime;
    SetObservableValue("SecondsFromStart", secondsFromStart);
    SetObservableValue("HoursFromStart", secondsFromStart / 3600.);
 
    if (fRateAnalysis) {
        Double_t evTimeDelay = 0;
        if (fPreviousEventTime.size() > 0) evTimeDelay = fEvent->GetTime() - fPreviousEventTime.back();
        SetObservableValue("EventTimeDelay", evTimeDelay);
 
        Double_t meanRate = 0;
        if (fPreviousEventTime.size() == 10)
            meanRate = 10. / (fEvent->GetTime() - fPreviousEventTime.front());
        SetObservableValue("MeanRate_InHz", meanRate);
 
        if (GetVerboseLevel() >= TRestStringOutput::REST_Verbose_Level::REST_Debug) {
            for (auto i : fObservablesDefined) {
                fAnalysisTree->PrintObservable(i.second);
            }
        }
 
        fPreviousEventTime.push_back(fEvent->GetTimeStamp());
        if (fPreviousEventTime.size() > 10) fPreviousEventTime.erase(fPreviousEventTime.begin());
    }
    // If cut condition matches the event will be not registered.
    if (ApplyCut()) return NULL;
 
    return fEvent;
}