CoREAS.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/RadioProcess.hpp>
#include <corsika/framework/geometry/QuantityVector.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>

namespace corsika {

  template <typename TRadioDetector, typename TPropagator>
  class CoREAS final
      : public RadioProcess<TRadioDetector, CoREAS<TRadioDetector, TPropagator>,
                            TPropagator> {

  public:
    // an identifier for which algorithm was used
    static constexpr auto algorithm = "CoREAS";

    template <typename... TArgs>
    CoREAS(TRadioDetector& detector, TArgs&&... args)
        : RadioProcess<TRadioDetector, CoREAS, TPropagator>(detector, args...) {}

    template <typename Particle>
    ProcessReturn simulate(Step<Particle> const& step);
    static constexpr auto emConstant_{1.0 / (4.0 * M_PI) / (constants::epsilonZero) /
                                      constants::c};

    using Base =
        RadioProcess<TRadioDetector, CoREAS<TRadioDetector, TPropagator>, TPropagator>;
    using Base::antennas_;

  }; // end of class CoREAS

} // namespace corsika

#include <corsika/detail/modules/radio/CoREAS.inl>