TRestRawFindResponseSignalProcess.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 "TRestRawFindResponseSignalProcess.h"
 
using namespace std;
 
#include <TFile.h>
 
ClassImp(TRestRawFindResponseSignalProcess);
 
TRestRawFindResponseSignalProcess::TRestRawFindResponseSignalProcess() { Initialize(); }
 
TRestRawFindResponseSignalProcess::TRestRawFindResponseSignalProcess(const char* configFilename) {
    Initialize();
 
    if (LoadConfigFromFile(configFilename) == -1) LoadDefaultConfig();
 
    PrintMetadata();
    // TRestRawFindResponseSignalProcess default constructor
}
 
TRestRawFindResponseSignalProcess::~TRestRawFindResponseSignalProcess() { delete fOutputSignalEvent; }
 
void TRestRawFindResponseSignalProcess::LoadDefaultConfig() {
    SetName("findResponseSignal-Default");
    SetTitle("Default config");
}
 
void TRestRawFindResponseSignalProcess::Initialize() {
    SetSectionName(this->ClassName());
    SetLibraryVersion(LIBRARY_VERSION);
 
    fInputSignalEvent = nullptr;
    fOutputSignalEvent = new TRestRawSignalEvent();
}
 
void TRestRawFindResponseSignalProcess::LoadConfig(const string& configFilename, const string& name) {
    if (LoadConfigFromFile(configFilename, name) == -1) LoadDefaultConfig();
}
 
void TRestRawFindResponseSignalProcess::InitProcess() {}
 
TRestEvent* TRestRawFindResponseSignalProcess::ProcessEvent(TRestEvent* inputEvent) {
    fInputSignalEvent = (TRestRawSignalEvent*)inputEvent;
 
    // We accept signals that are inside a given condition.
    // TODO: Now it is also possible to use ApplyCut and <cut definitions?
    if (fInputSignalEvent->GetNumberOfSignals() <= 0) return nullptr;
    if (fInputSignalEvent->GetNumberOfSignals() > 8) return nullptr;
 
    Int_t dominantSignal = -1;
    Double_t maxPeak = 0;
    //   Double_t maxTime = 0;
    for (int n = 0; n < fInputSignalEvent->GetNumberOfSignals(); n++) {
        if (fInputSignalEvent->GetSignal(n)->GetMaxPeakValue() > maxPeak) {
            maxPeak = fInputSignalEvent->GetSignal(n)->GetMaxPeakValue();
            // maxTime = fInputSignalEvent->GetSignal(n)->GetMaxPeakBin();
            dominantSignal = n;
        }
    }
 
    if (maxPeak < 400 || maxPeak > 600) return nullptr;
 
    TRestRawSignal* sgnl = fInputSignalEvent->GetSignal(dominantSignal);
    sgnl->Scale(1000. / maxPeak);
 
    fOutputSignalEvent->AddSignal(*sgnl);
 
    return fOutputSignalEvent;
}
 
void TRestRawFindResponseSignalProcess::EndProcess() {}