RdStationInterpolatorStarShape.h

Go to the documentation of this file.

#ifndef RdStationInterpolatorStarShape_H
#define RdStationInterpolatorStarShape_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>

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

namespace evt {
  class Event;
  class ShowerSimData;
  }
namespace revt {
  class Channel;
  }
namespace utl {
  class Point;
  }
namespace rdet {
  class RDetector;
}


namespace RdStationInterpolatorStarShape {
class SimPulseEntry {
  public:
    double fDistanceFromShowerAxis;
    double fObserverAzimuthAngle;
    const evt::SimRadioPulse* fRadioPulse;
    
    SimPulseEntry(double distance, double azimuth, const evt::SimRadioPulse* radiopulse)
    : fDistanceFromShowerAxis(distance),
      fObserverAzimuthAngle(azimuth),
      fRadioPulse(radiopulse)
    {}

    bool operator< (const SimPulseEntry& rhs) const {   // sort by radius in the shower plane, then azimuth
      if (fabs(fDistanceFromShowerAxis - rhs.fDistanceFromShowerAxis)/fDistanceFromShowerAxis > 1.e-3)
        return (fDistanceFromShowerAxis < rhs.fDistanceFromShowerAxis);
      else
        return (fObserverAzimuthAngle < rhs.fObserverAzimuthAngle);
    }
};


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

    private:

      revt::StationFFTDataContainer Interpolate( const revt::StationFFTDataContainer& simData1, const revt::StationFFTDataContainer& simData2, double x1, double x2, double x, double designlowerfreq, double designupperfreq);

      revt::StationTimeSeries PadTimeSeries(const evt::SimRadioPulse& simtimeseries);

      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("RdStationInterpolatorStarShape", RdStationInterpolatorStarShape);
  };

} // namespace
#endif