TRestRawToDetectorSignalProcess.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 "TRestRawToDetectorSignalProcess.h"
 
using namespace std;
 
ClassImp(TRestRawToDetectorSignalProcess);
 
TRestRawToDetectorSignalProcess::TRestRawToDetectorSignalProcess() { Initialize(); }
 
TRestRawToDetectorSignalProcess::~TRestRawToDetectorSignalProcess() { delete fOutputSignalEvent; }
 
void TRestRawToDetectorSignalProcess::Initialize() {
    SetSectionName(this->ClassName());
    SetLibraryVersion(LIBRARY_VERSION);
 
    fInputSignalEvent = nullptr;
    fOutputSignalEvent = new TRestDetectorSignalEvent();
}
 
TRestEvent* TRestRawToDetectorSignalProcess::ProcessEvent(TRestEvent* inputEvent) {
    fInputSignalEvent = (TRestRawSignalEvent*)inputEvent;
 
    Int_t rejectedSignal = 0;
 
    if (fZeroSuppression) {
        fInputSignalEvent->SetBaseLineRange(fBaseLineRange);
        fInputSignalEvent->SetRange(fIntegralRange);
    }
 
    for (int n = 0; n < fInputSignalEvent->GetNumberOfSignals(); n++) {
        TRestDetectorSignal signal;
        TRestRawSignal* rawSignal = fInputSignalEvent->GetSignal(n);
        signal.SetID(rawSignal->GetID());
 
        if (fZeroSuppression) {
            ZeroSuppresion(rawSignal, signal);
        } else {
            for (int p = 0; p < int(rawSignal->GetNumberOfPoints()); p++) {
                if (rawSignal->GetData(p) > fThreshold) {
                    signal.NewPoint(fTriggerStarts + fSampling * p, fGain * rawSignal->GetData(p));
                }
            }
        }
 
        if (signal.GetNumberOfPoints() > 0) {
            fOutputSignalEvent->AddSignal(signal);
        } else {
            rejectedSignal++;
        }
    }
 
    SetObservableValue("NSignalsRejected", rejectedSignal);
 
    if (fOutputSignalEvent->GetNumberOfSignals() <= 0) {
        return nullptr;
    }
 
    return fOutputSignalEvent;
}
 
void TRestRawToDetectorSignalProcess::ZeroSuppresion(TRestRawSignal* rawSignal, TRestDetectorSignal& signal) {
    rawSignal->InitializePointsOverThreshold(TVector2(fPointThreshold, fSignalThreshold),
                                             fNPointsOverThreshold, 512);
 
    std::vector<Int_t> pOver = rawSignal->GetPointsOverThreshold();
    for (unsigned int n = 0; n < pOver.size(); n++) {
        int j = pOver[n];
        signal.NewPoint(fTriggerStarts + fSampling * j, fGain * rawSignal->GetData(j));
    }
}