TRestRawSignalChannelActivityProcess.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 "TRestRawSignalChannelActivityProcess.h"
 
using namespace std;
 
ClassImp(TRestRawSignalChannelActivityProcess);
 
TRestRawSignalChannelActivityProcess::TRestRawSignalChannelActivityProcess() { Initialize(); }
 
TRestRawSignalChannelActivityProcess::~TRestRawSignalChannelActivityProcess() {}
 
void TRestRawSignalChannelActivityProcess::Initialize() {
    SetSectionName(this->ClassName());
    SetLibraryVersion(LIBRARY_VERSION);
 
    fInputEvent = nullptr;
}
 
void TRestRawSignalChannelActivityProcess::InitProcess() {
    if (!fReadOnly) {
        // Create a unique histogram name based on signal type
        string histogramName = "daqChannelActivityRaw_" + fChannelType;
 
        // Create the histogram using the unique name
        fDaqChannelsHisto = new TH1D(histogramName.c_str(), histogramName.c_str(), fDaqChannels,
                                     fDaqStartChannel, fDaqEndChannel);
    }
}
 
TRestEvent* TRestRawSignalChannelActivityProcess::ProcessEvent(TRestEvent* inputEvent) {
    fInputEvent = dynamic_cast<TRestRawSignalEvent*>(inputEvent);
 
    const auto run = GetRunInfo();
    if (run != nullptr) {
        fInputEvent->InitializeReferences(run);
    }
 
    if (fReadoutMetadata == nullptr) {
        fReadoutMetadata = fInputEvent->GetReadoutMetadata();
    }
 
    if (fReadoutMetadata == nullptr && !fChannelType.empty()) {
        cerr << "TRestRawSignalChannelActivityProcess::ProcessEvent: readout metadata is null, cannot filter "
                "the process by signal type"
             << endl;
        exit(1);
    }
 
    for (int s = 0; s < fInputEvent->GetNumberOfSignals(); s++) {
        const auto signal = fInputEvent->GetSignal(s);
        if (!fChannelType.empty()) {
            const auto channelType = fReadoutMetadata->GetTypeForChannelDaqId(signal->GetID());
            if (fChannelType != channelType) {
                continue;
            }
        }
        // Adding signal to the channel activity histogram
        if (!fReadOnly) {
            Int_t daqChannel = signal->GetID();
            fDaqChannelsHisto->Fill(daqChannel);
        }
    }
 
    if (GetVerboseLevel() >= TRestStringOutput::REST_Verbose_Level::REST_Debug)
        fAnalysisTree->PrintObservables();
 
    return fInputEvent;
}
 
void TRestRawSignalChannelActivityProcess::EndProcess() {
    if (!fReadOnly) {
        fDaqChannelsHisto->Write();
    }
}