TRestRawSignalIdTaggingProcess.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 "TRestRawSignalIdTaggingProcess.h"
 
using namespace std;
 
ClassImp(TRestRawSignalIdTaggingProcess);
 
TRestRawSignalIdTaggingProcess::TRestRawSignalIdTaggingProcess() { Initialize(); }
 
TRestRawSignalIdTaggingProcess::~TRestRawSignalIdTaggingProcess() {}
 
void TRestRawSignalIdTaggingProcess::Initialize() {
    SetSectionName(this->ClassName());
    SetLibraryVersion(LIBRARY_VERSION);
 
    fSignalEvent = nullptr;
}
 
void TRestRawSignalIdTaggingProcess::InitProcess() {}
 
void TRestRawSignalIdTaggingProcess::InitFromConfigFile() {
    // This line is to exploit the retrieval of parameter as it is done at any process
    TRestEventProcess::InitFromConfigFile();
 
    RESTDebug << "fPointThreshold: " << fPointThreshold << RESTendl;
    RESTDebug << "fSignalThreshold: " << fSignalThreshold << RESTendl;
    RESTDebug << "fPointsOverThreshold: " << fPointsOverThreshold << RESTendl;
    RESTDebug << "fBaseLineRange.X(): " << fBaseLineRange.X() << RESTendl;
    RESTDebug << "fBaseLineRange.Y(): " << fBaseLineRange.Y() << RESTendl;
 
    if (fPointThreshold != -1 && fSignalThreshold != -1 && fPointsOverThreshold != -1 &&
        fBaseLineRange.X() != -1 && fBaseLineRange.Y() != -1) {
        fGoodSignalsOnly = true;
    }
 
    // This is the additional code required by the process to read tags
    TiXmlElement* tagDefinition = GetElement("tag");
    while (tagDefinition) {
        fTagNames.push_back(GetFieldValue("name", tagDefinition));
        fIdRanges.push_back(Get2DVectorParameterWithUnits("ids", tagDefinition));
 
        tagDefinition = GetNextElement(tagDefinition);
    }
}
 
TRestEvent* TRestRawSignalIdTaggingProcess::ProcessEvent(TRestEvent* evInput) {
    fSignalEvent = (TRestRawSignalEvent*)evInput;
 
    std::vector<int> tagCode;
 
    for (int j = 0; j < fSignalEvent->GetNumberOfSignals(); j++) {
        TRestRawSignal* singleSignal = fSignalEvent->GetSignal(j);
 
        if (fGoodSignalsOnly == true) {
            singleSignal->CalculateBaseLine(fBaseLineRange.X(), fBaseLineRange.Y());
            singleSignal->InitializePointsOverThreshold(TVector2(fPointThreshold, fSignalThreshold),
                                                        fPointsOverThreshold);
        }
 
        if (fGoodSignalsOnly == false ||
            singleSignal->GetPointsOverThreshold().size() >= (unsigned int)fPointsOverThreshold) {
            for (unsigned int n = 0; n < fIdRanges.size(); n++) {
                if (singleSignal->GetID() >= fIdRanges[n].X() && singleSignal->GetID() <= fIdRanges[n].Y()) {
                    // If it is not already in the vector, adds it. n+1 to avoid 0.
                    if (std::find(tagCode.begin(), tagCode.end(), n + 1) == tagCode.end()) {
                        tagCode.push_back(n + 1);
                    }
                }
            }
        }
    }
    std::sort(tagCode.begin(), tagCode.end());
 
    int result = 0;
    for (auto d : tagCode) {
        result = result * 10 + d;
    }
    SetObservableValue("tagId", result);
 
    // If cut condition matches the event will be not registered.
    if (ApplyCut()) return nullptr;
 
    return fSignalEvent;
}