RdChannelGalacticConstantsGenerator.h

Go to the documentation of this file.

#ifndef _RdChannelGalacticConstantsGenerator_RdChannelGalacticConstantsGenerator_h_
#define _RdChannelGalacticConstantsGenerator_RdChannelGalacticConstantsGenerator_h_

#include <fwk/VModule.h>
#include "TProfile3D.h"
#include <string>
#include <fstream>


namespace evt {
  class Event;
}

namespace revt {
  class Station;
}


namespace RdChannelGalacticConstantsGenerator {

  class RdChannelGalacticConstantsGenerator : public fwk::VModule {

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

  private:
    std::vector<int> fIgnoreStations;
    double fMaxIntensity = 0;
    std::string fOutputFileName;
    bool fSaveROOTFile = true;
    std::string fFileNameROOTFile;

    std::ofstream fOutfile;
    bool fFirstEvent = true;

    // TProfiles of NS and EW spectra: [station id][LST][frequency][spectral amplitude]
    TProfile3D* fLSTSpecNS = nullptr;
    TProfile3D* fLSTSpecEW = nullptr;

    // ROOT quantities and statistics
    unsigned long fLowestGPSSec = std::numeric_limits<unsigned long>::max();
    unsigned long fHighestGPSSec = 0;

    static double GPSSecToLSTHourAuger(const int gpsSec);

    REGISTER_MODULE("RdChannelGalacticConstantsGenerator", RdChannelGalacticConstantsGenerator);
  };

}


#endif