DrumPhotonGenerator.h

Go to the documentation of this file.

#ifndef _DrumPhotonGeneratorOG_DrumPhotonGenerator_h_
#define _DrumPhotonGeneratorOG_DrumPhotonGenerator_h_

#include <fwk/VModule.h>

#include <utl/CoordinateSystemPtr.h>
#include <utl/ShadowPtr.h>

#include <map>
#include <vector>
#include <fstream>
#include <string>


namespace utl {
  class Point;
  class Vector;
  class RandomEngine;
}

namespace fdet {
  class Telescope;
}

namespace evt {
  class Event;
}

namespace fevt {
  class Telescope;
}


namespace DrumPhotonGeneratorOG {

  class DrumPhotonGenerator : public fwk::VModule {

    enum Status {
      eGeneratePhotons,
      eCalibrate,
      eCalibrated
    };

    struct CalibResult {
      double calib;
      double calibErr;
      int nPhotGen;
      int nPhotPixel;
      double electronicsFactor;
    };

  public:
    DrumPhotonGenerator();
    ~DrumPhotonGenerator();

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

  private:
    int Index(const int eyeId, const int telId) const { return eyeId*100 + telId; }

    bool DoDrum(evt::Event& event);
    bool CalculateCalibrationConstants(evt::Event& event);
    void GenerateDrumPhotons(fevt::Telescope& tel);

    std::map<int, CalibResult> CalibrateTelescope(fevt::Telescope& tel);

  private:
    // drum mode
    double fMinAccuracy;
    int fNDrumPhotons;
    //double fWeightDrumPhotons; // unused. LN.
    double fDrumNCos;
    double fThetaMinDrumPhotons;
    double fThetaMaxDrumPhotons;
    double fPhiMinDrumPhotons;
    double fPhiMaxDrumPhotons;
    unsigned int fMinBin;
    unsigned int fMaxBin;
    double fRDiaMin;
    std::string fDrumDataOutName;

    utl::RandomEngine* fRandomEngine;

    std::map<int, unsigned int> fTotalGeneratedPhotonsPerTelescope;
    bool fInit;
    Status fStatus;

    REGISTER_MODULE("DrumPhotonGeneratorOG", DrumPhotonGenerator);

  };

}


#endif

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