RdStationInterpolator.h

Go to the documentation of this file.

#ifndef _RdStationInterpolator_RdStationInterpolator_h_
#define _RdStationInterpolator_RdStationInterpolator_h_

#include <fwk/VModule.h>
#include <utl/CoordinateSystem.h>
#include <utl/TimeStamp.h>
#include <vector>

#include <utl/FFTDataContainer.h>
#include <evt/SimRadioPulse.h>
#include <evt/RadioSimulation.h>
#include <revt/Station.h>
#include <complex>
#include <map>
#include <utl/Point.h>

#include <TRandom2.h>
#include <string>

namespace evt {
  class Event;
  class ShowerSimData;
}

namespace revt {
  class Channel;
}

namespace utl {
  class Point;
}

namespace rdet {
  class RDetector;
}


namespace RdStationInterpolator {

  class RdStationInterpolator : public fwk::VModule {

  public:
    RdStationInterpolator();
    ~RdStationInterpolator();
    VModule::ResultFlag Init();
    VModule::ResultFlag Run(evt::Event& theevent);
    fwk::VModule::ResultFlag Finish();

  private:
    typedef std::multimap<double,const evt::SimRadioPulse&> PulseMMap;

    PulseMMap FindClosestSimPulsesToStation(const utl::Point& pt, const utl::CoordinateSystemPtr& coord,
                                            const PulseMMap& mm_srpulse_x, const PulseMMap& mm_srpulse_y) const;
    PulseMMap FindSurroundingNN(const PulseMMap& mm_simradiopulse, const utl::Point& pt, const utl::CoordinateSystemPtr& coord, const bool inverse) const;
    std::vector<double> GetAreaOfTriangles(const utl::Point& stationpos, const std::vector<utl::Point>& v_simloc, const utl::CoordinateSystemPtr& coord ) const;
    revt::StationTimeSeries PadTimeSeries(const evt::SimRadioPulse& simtimeseries) const;
    revt::StationFFTDataContainer InterpolateInTwoD(const std::vector<revt::StationFFTDataContainer>& simData, const std::vector<double>& area,
                                                    const double designLowerFreq, const double designUpperFreq) const;
    void BuildSimShower(evt::ShowerSimData& theshower, const utl::CoordinateSystemPtr& coord);
    double fMaximumDistance;
    double fTracePaddingAtFront;
    double fMinimumTraceLength;
    double fEventRate;
    utl::TimeInterval fEventTimeOffset;
    utl::TimeStamp fFirstEventTime;
    bool fUseRateTiming;
    std::vector<int> fIncludedStationIds;
    std::vector<int> fExcludedStationIds;
    int fInfoLevel;
    unsigned int fEventcounter;

    enum InfoLevel {
      eNone = 0,
      eFinal = 1,
      eIntermediate = 2,
      eObscure = 3,
      eMinuit = 4
    };

    REGISTER_MODULE("RdStationInterpolator", RdStationInterpolator);

  };

}


#endif