SdPMTSimulatorOG/SdPMTSimulator.h

Go to the documentation of this file.

#ifndef _SdPMTSimulatorOG_SdPMTSimulator_h_
#define _SdPMTSimulatorOG_SdPMTSimulator_h_

#include <fwk/VModule.h>
#include <utl/TabulatedFunction.h>
#include <utl/TimeDistribution-fwd.h>
#include <utl/config.h>

#include <sevt/SEvent.h>
#include <sevt/StationConstants.h>
#include <CLHEP/Random/Randomize.h>

#include <set>


namespace utl {
  class RandomEngine;
}

class TH1D;

namespace SdPMTSimulatorOG {

  class SdPMTSimulator : public fwk::VModule {

  public:
    virtual ~SdPMTSimulator() = default;

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

  private:
    void FillPulseShape(const utl::TabulatedFunction& rawPulse);

    void ConvertPEToBaseSignal(const double scaling,
                               const utl::TimeDistributionI& pe,
                               utl::TimeDistributionD& baseSignal);

    void SimulateSaturation(utl::TimeDistributionD& baseSignal,
                            const utl::TimeDistributionD& totalBaseSignal);

    // points to maximum of the pulse
    double fPulseMaxTime = 0;

    utl::TabulatedFunction fMeanPulseShape;
    CLHEP::RandGeneral* fChargeDistribution = nullptr;

    double fPMTGain = 0;
    double fKFactor = 0;
    double fMinCharge = 0;
    double fMaxCharge = 0;

    utl::TabulatedFunction* fSaturationFunction = nullptr;
    double fCurrent2VoltageMultiplier = 0;

    utl::RandomEngine* fRandomEngine = nullptr;

    typedef std::set<sevt::StationConstants::SignalComponent> FundamentalComponents;
    FundamentalComponents fFundamentalComponents;

    REGISTER_MODULE("SdPMTSimulatorOG", SdPMTSimulator);

  };

}


#endif