Deprecated/UpgradeASCIITests/SdPMTSimulatorASCII/SdPMTSimulator.h

Go to the documentation of this file.

#ifndef _SdPMTSimulatorASCII_SdPMTSimulator_h_
#define _SdPMTSimulatorASCII_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 <vector>

namespace utl {
  class RandomEngine;
}

class TH1D;

namespace SdPMTSimulatorASCII {

  class SdPMTSimulator : public fwk::VModule {
  public:
    SdPMTSimulator();
    virtual ~SdPMTSimulator();

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

  private:
    CLHEP::RandGeneral* fChargeDist;

    sevt::SEvent::StationIterator fStationIterator;
    bool fFirstCycle;
    bool fLimitStationsPerCycle;
    unsigned int fNStationsPerCycle;
    
    // points to maximum of the pulse
    utl::TabulatedFunction::Iterator fMaxPulseHeight;

    utl::TabulatedFunction* fAveragePulseShape;
    double fPMTGain_WCD;
    double fPMTGain_ASCII;

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

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

    double fMinCharge;
    double fMaxCharge;
    bool fWantRootHistos;
    std::string fRootHistoFilename;

    // Histograms for testing
    TH1D* fChargeHisto;
    TH1D* fAnodePulseHisto;

    utl::TabulatedFunction* fSaturationFunction;
    double fCurrent2VoltageMultiplier;

    utl::RandomEngine* fRandomEngine;

    std::vector<sevt::StationConstants::SignalComponent> fFundamentalComponents,fFundamentalComponents_ASCII;

    REGISTER_MODULE("SdPMTSimulatorASCII", SdPMTSimulator);

  };

} // namespace SdPMTSimulatorOG


#endif


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