RdScintSignalReconstructor.h

Go to the documentation of this file.

#ifndef _RdScintSignalReconstructor_h_
#define _RdScintSignalReconstructor_h_

#include <fwk/VModule.h>
#include <revt/Station.h>
#include <revt/Channel.h>
#include <utl/Point.h>
#include <string>
#include <vector>
#include <map>


namespace evt {
  class Event;
}

namespace revt {
  class REvent;
}

namespace RdScintSignalReconstructor {

  class RdScintSignalReconstructor : public fwk::VModule {

  public:
    RdScintSignalReconstructor();
    virtual ~RdScintSignalReconstructor();

    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event& event);
    fwk::VModule::ResultFlag Finish();

  private:
    enum InfoLevel {
      eNone = 0,
      eFinal = 1,
      eIntermediate = 2,
      eObscure = 3,
      eDebug = 4
    };

    // xml settings

    double fScintSignalSearchWindowStart;
    double fScintSignalSearchWindowStop;
    double fNoiseWindowStart;
    unsigned int fSamplesToCalculateOffset;
    unsigned int fMinNumberOfScintTop;
    unsigned int fMinNumberOfScintBottom;
    unsigned int fMinNumberOfScint;
    double fMinSignal;
    bool fWeightedBaryValues;
    double fSimSingleMuonEnergyDepositTop;
    double fSimSingleMuonEnergyDepositTopError;
    double fSimSingleMuonEnergyDepositBottom;
    double fSimSingleMuonEnergyDepositBottomError;
    std::string fFilenameMuonCalibrationTop;
    std::string fFilenameMuonCalibrationBottom;

    struct CalibrationData {
      double signalMuonPeak;
      double signalMuonPeakError;

      CalibrationData() : signalMuonPeak(0), signalMuonPeakError(0) { }
    };
    std::map<int, CalibrationData> fCalibrationDataTop;
    std::map<int, CalibrationData> fCalibrationDataBottom;

    //std::vector<int> fExcludedStationsId;
    std::vector<std::string> fExcludedStationsName;

    int fInfoLevel;

    // calculates Mean and RMS in selected noise window
    void Noisefinder(const revt::ChannelTimeSeries& channeltrace, double NoiseWindowStart,
                     double NoiseWindowSize, double& SignalError, unsigned int samples_offset,
                     double& RMSNoise) const;

    // calculates the amplitude and the position of the maximum in the signal search window
    void Pulsefinder(const revt::ChannelTimeSeries& channeltrace,
                     double& PeakAmplitude, double& PeakTime, double& PeakTimeError,
                     double SignalWindowStart, double SignalWindowStop, unsigned int& sample,
                     unsigned int samples_offset) const;

    // calculates the integrated signal in a fixed window
    void PulseFixedWindowIntegrator(const revt::ChannelTimeSeries& channeltrace, double IntegrationTime,
                                    double& IntegratedSignal, double SignalWindowStart,
                                    double SignalWindowStop, unsigned int samples_offset) const;

    void Signalwindowfinder(const revt::ChannelTimeSeries channeltrace, double ScintSearchWindowStart,
                            double ScintSearchWindowStop, double& SignalWindowStart, double& SignalWindowStop,
                            unsigned int sample, unsigned int samples_offset) const;

    void ComputeBaryCenterTop(const revt::REvent& rEvent, utl::Point& BaryCenterTop) const;

    void ComputeBaryCenterBottom(const revt::REvent& rEvent, utl::Point& BaryCenterBottom) const;

    void ComputeBaryCenter(const revt::REvent& rEvent, utl::Point& BaryCenter) const;

    void ComputeBaryTimeTop(const revt::REvent& rEvent, double& BaryTimeTop) const;

    void ComputeBaryTimeBottom(const revt::REvent& rEvent, double& BaryTimeBottom) const;

    void ComputeBaryTime(const revt::REvent& rEvent, double& BaryTime) const;

    REGISTER_MODULE("RdScintSignalReconstructor", RdScintSignalReconstructor);

  };

}


#endif

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -k"
// End: