QPMonteCarlo.hpp

#ifndef __QPMONTECARLO__
#define __QPMONTECARLO__ 1
 
#include <string>
#include <vector>
 
class QPDigitizerBase;
class QPRandomBase;
 
class QPMonteCarlo{
    public:
    struct Record{
        double fAmp;   
        double fX;     
        double fY;     
        double fSigma; 
        int    fNpix;  
    };
 
    QPMonteCarlo(QPDigitizerBase* digitizer, QPRandomBase& random);
 
    void SetPositionRange(const double& xmin, const double& xmax,
                          const double& ymin, const double& ymax);
    void SetAmplitude(const double& amp);                        // fixed amplitude
    void SetAmplitudeRange(const double& min, const double& max);// uniform per event
    void SetSigma(const double& sigma);
    void SetThreshold(const double& threshold);
 
    Record Shoot();                      // one event (see draw order above)
    std::vector<Record> Run(int ntrial); // Shoot() ntrial times, accumulated
 
    static bool WriteCSV(const std::vector<Record>& records,
                         const std::string& filename);
 
    QPDigitizerBase* GetDigitizer() const {return fDigitizer;}
 
    protected:
    private:
    QPDigitizerBase* fDigitizer = nullptr;
    QPRandomBase&    fRandom;
 
    double fXmin = 0.;
    double fXmax = 0.;
    double fYmin = 0.;
    double fYmax = 0.;
    double fAmp  = 0.;
    double fAmpMin = 0.;
    double fAmpMax = 0.;
    bool   fRandomAmp = false;
    double fSigma = 1.;
    double fThreshold = 0.;
};
 
#endif