TRestDetectorSignalToRawSignalProcess.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_TRestDetectorSignalToRawSignalProcess
#define RestCore_TRestDetectorSignalToRawSignalProcess
 
#include <TRestDetectorSignalEvent.h>
#include <TRestEventProcess.h>
#include <TRestRawSignalEvent.h>
 
class TRestDetectorSignalToRawSignalProcess : public TRestEventProcess {
   private:
    TRestDetectorSignalEvent* fInputSignalEvent;  
 
    TRestRawSignalEvent* fOutputRawSignalEvent;  
 
    void Initialize() override;
 
   protected:
    Double_t fSampling = 1.0;  // ns
 
    Int_t fNPoints = 512;
 
    TString fTriggerMode = "firstDeposit";
 
    Int_t fTriggerDelay = 100;
 
    Int_t fTriggerFixedStartTime = 0;
 
    Double_t fCalibrationGain = 100.0;
    Double_t fCalibrationOffset = 0.0;  // adc units
 
    Double_t fIntegralThreshold = 1229.0;
 
    TVector2 fCalibrationEnergy = TVector2(0.0, 0.0);
    TVector2 fCalibrationRange = TVector2(0.0, 0.0);
 
    Double_t fShapingTime = 0.0;  // us
 
    Double_t fTimeStart;  
    Double_t fTimeEnd;    
 
   public:
    inline Double_t GetSampling() const { return fSampling; }
    inline void SetSampling(Double_t sampling) { fSampling = sampling; }
 
    inline Int_t GetNPoints() const { return fNPoints; }
    inline void SetNPoints(Int_t nPoints) { fNPoints = nPoints; }
 
    inline TString GetTriggerMode() const { return fTriggerMode; }
    inline void SetTriggerMode(const TString& triggerMode) { fTriggerMode = triggerMode; }
 
    inline Int_t GetTriggerDelay() const { return fTriggerDelay; }
    inline void SetTriggerDelay(Int_t triggerDelay) { fTriggerDelay = triggerDelay; }
 
    inline Double_t GetGain() const { return fCalibrationGain; }
    inline void SetGain(Double_t gain) { fCalibrationGain = gain; }
 
    inline Double_t GetCalibrationOffset() const { return fCalibrationOffset; }
    inline void SetCalibrationOffset(Double_t offset) { fCalibrationOffset = offset; }
 
    inline Double_t GetIntegralThreshold() const { return fIntegralThreshold; }
    inline void SetIntegralThreshold(Double_t integralThreshold) { fIntegralThreshold = integralThreshold; }
 
    inline Double_t GetShapingTime() const { return fShapingTime; }
    inline void SetShapingTime(Double_t shapingTime) { fShapingTime = shapingTime; }
 
    inline bool IsShapingEnabled() const { return fShapingTime > 0; }
 
    inline bool IsLinearCalibration() const {
        // Will return true if two points have been given for calibration
        return (fCalibrationEnergy.Mod() != 0 && fCalibrationRange.Mod() != 0);
    }
 
    inline TVector2 GetCalibrationEnergy() const { return fCalibrationEnergy; }
    inline void SetCalibrationEnergy(TVector2 calibrationEnergy) { fCalibrationEnergy = calibrationEnergy; }
 
    inline TVector2 GetCalibrationRange() const { return fCalibrationRange; }
    inline void SetCalibrationRange(TVector2 calibrationRange) { fCalibrationRange = calibrationRange; }
 
    any GetInputEvent() const override { return fInputSignalEvent; }
    any GetOutputEvent() const override { return fOutputRawSignalEvent; }
 
    void InitProcess() override;
 
    TRestEvent* ProcessEvent(TRestEvent* inputEvent) override;
 
    void LoadConfig(const std::string& configFilename, const std::string& name = "");
 
    void PrintMetadata() override {
        BeginPrintProcess();
 
        RESTMetadata << "Sampling time : " << fSampling << " us" << RESTendl;
        RESTMetadata << "Points per channel : " << fNPoints << RESTendl;
        RESTMetadata << "Trigger mode : " << fTriggerMode << RESTendl;
        RESTMetadata << "Trigger delay : " << fTriggerDelay << " time units" << RESTendl;
 
        if (IsLinearCalibration()) {
            RESTMetadata << "Calibration energy : (" << fCalibrationEnergy.X() << ", "
                         << fCalibrationEnergy.Y() << ") keV" << RESTendl;
            RESTMetadata << "Calibration range : (" << fCalibrationRange.X() << ", " << fCalibrationRange.Y()
                         << ")" << RESTendl;
        }
        RESTMetadata << "ADC Gain : " << fCalibrationGain << RESTendl;
        RESTMetadata << "ADC Offset : " << fCalibrationOffset << RESTendl;
 
        if (IsShapingEnabled()) {
            RESTMetadata << "Shaping time : " << fShapingTime << " us" << RESTendl;
        }
 
        EndPrintProcess();
    }
 
    TRestEventProcess* Maker() { return new TRestDetectorSignalToRawSignalProcess; }
 
    const char* GetProcessName() const override { return "signalToRawSignal"; }
 
    // Constructor
    TRestDetectorSignalToRawSignalProcess();
    TRestDetectorSignalToRawSignalProcess(const char* configFilename);
 
    // Destructor
    ~TRestDetectorSignalToRawSignalProcess();
 
    ClassDefOverride(TRestDetectorSignalToRawSignalProcess, 3);
};
#endif