TRestEventSelectionProcess.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 "TRestEventSelectionProcess.h"
 
using namespace std;
 
ClassImp(TRestEventSelectionProcess);
 
TRestEventSelectionProcess::TRestEventSelectionProcess() { Initialize(); }
 
void TRestEventSelectionProcess::Initialize() {
    SetSectionName(this->ClassName());
 
    fEvent = nullptr;
}
 
void TRestEventSelectionProcess::InitProcess() {
    if (TRestTools::GetFileNameExtension(fFileWithIDs) == "txt") {
        string line;
        ifstream File(fFileWithIDs);
 
        if (File.is_open()) {
            while (getline(File, line)) {
                fList.push_back(stoi(line));
            }
            File.close();
        }
    } else if (TRestTools::GetFileNameExtension(fFileWithIDs) == "root") {
        TRestRun run(fFileWithIDs);
        fList = run.GetEventIdsWithConditions(fConditions);
    } else {
        RESTDebug << "TRestEventSelectionProcess: using the processing file itself." << RESTendl;
    }
}
 
TRestEvent* TRestEventSelectionProcess::ProcessEvent(TRestEvent* inputEvent) {
    fEvent = inputEvent;
 
    if (fFileWithIDs.empty()) {
        if (this->GetAnalysisTree()->EvaluateCuts(fConditions)) {
            return fEvent;
        }
    }
 
    for (auto id : fList) {
        if (id == fEvent->GetID()) {
            return fEvent;
        }
    }
 
    return nullptr;
}
 
void TRestEventSelectionProcess::PrintMetadata() {
    BeginPrintProcess();
 
    RESTMetadata << "File with IDs: " << fFileWithIDs << RESTendl;
    RESTMetadata << "Conditions: " << fConditions << RESTendl;
 
    EndPrintProcess();
}