TRestSummaryProcess.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 "TRestSummaryProcess.h"
using namespace std;
 
ClassImp(TRestSummaryProcess);
 
TRestSummaryProcess::TRestSummaryProcess() { Initialize(); }
 
TRestSummaryProcess::TRestSummaryProcess(const char* configFilename) {
    Initialize();
 
    LoadConfig(configFilename);
}
 
TRestSummaryProcess::~TRestSummaryProcess() {}
 
void TRestSummaryProcess::LoadDefaultConfig() {
    SetName(this->ClassName());
    SetTitle("Default config");
}
 
void TRestSummaryProcess::LoadConfig(const string& configFilename, const string& name) {
    if (LoadConfigFromFile(configFilename, name)) LoadDefaultConfig();
}
 
void TRestSummaryProcess::InitProcess() {
    // For example, try to initialize a pointer to existing metadata
    // accessible from TRestRun
}
 
void TRestSummaryProcess::Initialize() { SetSectionName(this->ClassName()); }
 
TRestEvent* TRestSummaryProcess::ProcessEvent(TRestEvent* inputEvent) {
    fEvent = inputEvent;
 
    // This process does nothing at the event-by-event level
 
    return fEvent;
}
 
void TRestSummaryProcess::EndProcess() {
    Int_t nEntries = GetFullAnalysisTree()->GetEntries();
    Double_t startTime = fRunInfo->GetStartTimestamp();
    Double_t endTime = fRunInfo->GetEndTimestamp();
 
    fMeanRate = nEntries / (endTime - startTime);
    fMeanRateSigma = TMath::Sqrt(nEntries) / (endTime - startTime);
 
    for (auto const& x : fAverage) {
        TVector2 range = fAverageRange[x.first];
        fAverage[x.first] = this->GetFullAnalysisTree()->GetObservableAverage(x.first, range.X(), range.Y());
    }
 
    for (auto const& x : fRMS) {
        TVector2 range = fRMSRange[x.first];
        fRMS[x.first] = this->GetFullAnalysisTree()->GetObservableRMS(x.first, range.X(), range.Y());
    }
 
    for (auto const& x : fMaximum) {
        TVector2 range = fMaximumRange[x.first];
        fMaximum[x.first] = this->GetFullAnalysisTree()->GetObservableMaximum(x.first, range.X(), range.Y());
    }
 
    for (auto const& x : fMinimum) {
        TVector2 range = fMinimumRange[x.first];
        fMinimum[x.first] = this->GetFullAnalysisTree()->GetObservableMinimum(x.first, range.X(), range.Y());
    }
 
    if (GetVerboseLevel() >= TRestStringOutput::REST_Verbose_Level::REST_Info) PrintMetadata();
}
 
void TRestSummaryProcess::InitFromConfigFile() {
    string definition;
    size_t pos = 0;
    while ((definition = GetKEYDefinition("average", pos)) != "") {
        TString obsName = GetFieldValue("obsName", definition);
 
        fAverage[obsName] = 0;
        fAverageRange[obsName] = StringTo2DVector(GetFieldValue("range", definition));
    }
 
    pos = 0;
    while ((definition = GetKEYDefinition("rms", pos)) != "") {
        TString obsName = GetFieldValue("obsName", definition);
 
        fRMS[obsName] = 0;
        fRMSRange[obsName] = StringTo2DVector(GetFieldValue("range", definition));
    }
 
    pos = 0;
    while ((definition = GetKEYDefinition("maximum", pos)) != "") {
        TString obsName = GetFieldValue("obsName", definition);
 
        fMaximum[obsName] = 0;
        fMaximumRange[obsName] = StringTo2DVector(GetFieldValue("range", definition));
    }
 
    pos = 0;
    while ((definition = GetKEYDefinition("minimum", pos)) != "") {
        TString obsName = GetFieldValue("obsName", definition);
 
        fMinimum[obsName] = 0;
        fMinimumRange[obsName] = StringTo2DVector(GetFieldValue("range", definition));
    }
}
 
void TRestSummaryProcess::PrintMetadata() {
    BeginPrintProcess();
 
    RESTMetadata << " - Mean rate : " << fMeanRate << " Hz" << RESTendl;
    RESTMetadata << " - Mean rate sigma : " << fMeanRateSigma << " Hz" << RESTendl;
    for (auto const& x : fAverage) {
        RESTMetadata << " " << RESTendl;
        RESTMetadata << x.first << " average:" << x.second << RESTendl;
        TVector2 a = fAverageRange[x.first];
        if (a.X() != -1 && a.Y() != -1)
            RESTMetadata << "    range : (" << a.X() << ", " << a.Y() << ")" << RESTendl;
    }
    for (auto const& x : fRMS) {
        RESTMetadata << " " << RESTendl;
        RESTMetadata << x.first << " RMS:" << x.second << RESTendl;
        TVector2 a = fRMSRange[x.first];
        if (a.X() != -1 && a.Y() != -1)
            RESTMetadata << "    range : (" << a.X() << ", " << a.Y() << ")" << RESTendl;
    }
    for (auto const& x : fMaximum) {
        RESTMetadata << " " << RESTendl;
        RESTMetadata << x.first << " Maximum:" << x.second << RESTendl;
        TVector2 a = fMaximumRange[x.first];
        if (a.X() != -1 && a.Y() != -1)
            RESTMetadata << "    range : (" << a.X() << ", " << a.Y() << ")" << RESTendl;
    }
    for (auto const& x : fMinimum) {
        RESTMetadata << " " << RESTendl;
        RESTMetadata << x.first << " Minimum:" << x.second << RESTendl;
        TVector2 a = fMinimumRange[x.first];
        if (a.X() != -1 && a.Y() != -1)
            RESTMetadata << "    range : (" << a.X() << ", " << a.Y() << ")" << RESTendl;
    }
    EndPrintProcess();
}