TimeDomainAntenna.hpp

/*
 * (c) Copyright 2022 CORSIKA Project, corsika-project@lists.kit.edu
 *
 * This software is distributed under the terms of the GNU General Public
 * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
 * the license.
 */
#pragma once

#include <corsika/modules/radio/antennas/Antenna.hpp>
#include <yaml-cpp/yaml.h>
#include <vector>

namespace corsika {

  class TimeDomainAntenna : public Antenna<TimeDomainAntenna> {

    TimeType const start_time_;         
    TimeType const duration_;           
    InverseTimeType const sample_rate_; 
    TimeType const ground_hit_time_; 
    uint64_t const num_bins_;        
    std::vector<double> waveformEX_; 
    std::vector<double> waveformEY_; 
    std::vector<double> waveformEZ_; 
    std::vector<long double> const
        time_axis_; 

  public:
    // import the methods from the antenna
    using Antenna<TimeDomainAntenna>::getName;
    using Antenna<TimeDomainAntenna>::getLocation;

    TimeDomainAntenna(std::string const& name, Point const& location,
                      CoordinateSystemPtr coordinateSystem, TimeType const& start_time,
                      TimeType const& duration, InverseTimeType const& sample_rate,
                      TimeType const ground_hit_time);

    // TODO: rethink this method a bit. If the endpoint is at the end of the antenna
    // resolution then you get the startpoint signal but you lose the endpoint signal!
    void receive(TimeType const time, Vector<dimensionless_d> const& receive_vector,
                 ElectricFieldVector const& efield);

    void receive(TimeType const time, Vector<dimensionless_d> const& receive_vector,
                 VectorPotential const& vectorP);

    auto const& getWaveformX() const;

    auto const& getWaveformY() const;

    auto const& getWaveformZ() const;

    std::string const getDomainLabel();

    std::vector<long double> createTimeAxis() const;

    auto const getAxis() const;

    InverseTimeType const& getSampleRate() const;

    TimeType const& getStartTime() const;

    void reset();

    YAML::Node getConfig() const;

  }; // END: class TimeDomainAntenna

} // namespace corsika

#include <corsika/detail/modules/radio/antennas/TimeDomainAntenna.inl>