RadioSimulation.h

Go to the documentation of this file.

#ifndef _evt_RadioSimulation_h_
#define _evt_RadioSimulation_h_

#include <evt/SimRadioPulse.h>
#include <utl/LameShadowPtr.h>
#include <utl/Point.h>
#include <utl/TimeStamp.h>
#include <utl/Validated.h>
#include <fwk/LocalCoordinateSystem.h>
#include <vector>


namespace evt {
  // forward declaration for friendship
  class ShowerSimData;
}

namespace io {
  class RadioSimulation_ROOT;
}

namespace evt {

  class RadioSimulation {

  public:
    typedef std::vector<SimRadioPulse>::const_iterator PulseIterator;

    // PulseIterator FirstPulse() { return fPulseIter.Get().begin(); }
    // PulseIterator LastPulse() { return fPulseIter.Get().end(); }
    const SimRadioPulse& GetNextSimRadioPulse(bool& ok);
    // const SimRadioPulse& GetSimRadioPulse();

    const SimRadioPulse& GetSimPulseByIndex(const int index);

    bool HasSimPulseByIndex(const int index) { return -1 < index && index < int(fSimRadioPulseVector.size()); }

    long GetNumPulses() const { return fSimRadioPulseVector.size(); }

    bool HasCorePosition() const { return fCorePosition.IsValid(); }

    utl::CoordinateSystemPtr GetLocalCoordinateSystem() const;

//#warning this interface is implicitly assuming x,y,z is in detector cs
    void SetCoreCoordinates(const double x, const double y, const double z);

    void SetCorePosition(const utl::Point& core);

    utl::Point GetCorePosition() const;

    int GetRunNumber() const { return fRunNumber; }

    void SetRunNumber(const int runnum) { fRunNumber = runnum; }

    int GetEventNumber() const { return fEventNumber; }

    void SetEventNumber(const int eventnum) { fEventNumber = eventnum; }

    const utl::TimeStamp& GetEventTime() const { return fEventTime; }

    void SetEventTime(const utl::TimeStamp& t) { fEventTime = t; }

    void AddSimRadioPulse(const SimRadioPulse& rp) { fSimRadioPulseVector.push_back(rp); }

    void ClearSimRadioPulses() { fSimRadioPulseVector.clear(); }

    bool GoToFirstSimRadioPulse();

    // Set refractive index at sea level specified in *.reas file
    void SetRefractiveIndexAtSeaLevel(const double n) { fRefractiveIndexAtSeaLevel = n; }

    // Get refractive index at sea level specified in *.reas file
    double GetRefractiveIndexAtSeaLevel() const { return fRefractiveIndexAtSeaLevel; }

    void SetRadiationEnergy(const double radiationEnergy) { fRadiationEnergy = radiationEnergy; }

    double GetRadiationEnergy() const { return fRadiationEnergy; }

  private:
    RadioSimulation() { }
    ~RadioSimulation() { }

    utl::Validated<utl::Point> fCorePosition;
    int fRunNumber = -1;
    int fEventNumber = -1;
    utl::TimeStamp fEventTime;

    double fRefractiveIndexAtSeaLevel = 0;
    double fRadiationEnergy = 0;

    utl::Validated<utl::CoordinateSystemPtr> fLocalCoordinateSystem;

    std::vector<SimRadioPulse> fSimRadioPulseVector;
    utl::Validated<PulseIterator> fPulseIter;

    friend class evt::ShowerSimData;
    friend class utl::LameShadowPtr<RadioSimulation>;
    friend class io::RadioSimulation_ROOT;
  };

}


#endif