TRestRawToDetectorSignalProcess.h

/*************************************************************************
 * 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.                  *
 *************************************************************************/
 
#ifndef RestCore_TRestRawToDetectorSignalProcess
#define RestCore_TRestRawToDetectorSignalProcess
 
#include <TRestDetectorSignalEvent.h>
#include <TRestRawSignalEvent.h>
 
#include "TRestEventProcess.h"
 
class TRestRawToDetectorSignalProcess : public TRestEventProcess {
   private:
    TRestRawSignalEvent* fInputSignalEvent;  
 
    TRestDetectorSignalEvent* fOutputSignalEvent;  
 
    void Initialize() override;
 
   protected:
    Double_t fSampling = 0.1;
 
    Double_t fTriggerStarts = 0;
 
    Double_t fGain = 1;
 
    Double_t fThreshold = 0.1;
 
    // Perform Zero suppression to the data
    Bool_t fZeroSuppression = false;
 
    TVector2 fBaseLineRange = TVector2(5, 55);
 
    TVector2 fIntegralRange = TVector2(10, 500);
 
    Double_t fPointThreshold = 3;
 
    Double_t fSignalThreshold = 5;
 
    Int_t fNPointsOverThreshold = 5;
 
    Bool_t fBaseLineCorrection = false;
 
   public:
    RESTValue GetInputEvent() const override { return fInputSignalEvent; }
    RESTValue GetOutputEvent() const override { return fOutputSignalEvent; }
 
    TRestEvent* ProcessEvent(TRestEvent* inputEvent) override;
 
    void ZeroSuppresion(TRestRawSignal* rawSignal, TRestDetectorSignal& signal);
 
    void PrintMetadata() override {
        BeginPrintProcess();
 
        RESTMetadata << "Sampling time : " << fSampling << " us" << RESTendl;
        RESTMetadata << "Trigger starts : " << fTriggerStarts << " us" << RESTendl;
        RESTMetadata << "Gain : " << fGain << RESTendl;
 
        if (fZeroSuppression) {
            RESTMetadata << "Base line range definition : ( " << fBaseLineRange.X() << " , "
                         << fBaseLineRange.Y() << " ) " << RESTendl;
            RESTMetadata << "Integral range : ( " << fIntegralRange.X() << " , " << fIntegralRange.Y()
                         << " ) " << RESTendl;
            RESTMetadata << "Point Threshold : " << fPointThreshold << " sigmas" << RESTendl;
            RESTMetadata << "Signal threshold : " << fSignalThreshold << " sigmas" << RESTendl;
            RESTMetadata << "Number of points over threshold : " << fNPointsOverThreshold << RESTendl;
        }
 
        if (fBaseLineCorrection)
            RESTMetadata << "BaseLine correction is enabled for TRestRawSignalAnalysisProcess" << RESTendl;
 
        EndPrintProcess();
    }
 
    TRestEventProcess* Maker() { return new TRestRawToDetectorSignalProcess; }
 
    const char* GetProcessName() const override { return "rawSignalToSignal"; }
 
    // Constructor
    TRestRawToDetectorSignalProcess();
 
    // Destructor
    ~TRestRawToDetectorSignalProcess();
 
    ClassDefOverride(TRestRawToDetectorSignalProcess, 2);
};
#endif