CONEXhybrid.hpp

/*
 * (c) Copyright 2020 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/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/core/ParticleProperties.hpp>
#include <corsika/framework/geometry/Point.hpp>
#include <corsika/framework/geometry/Vector.hpp>
#include <corsika/framework/process/SecondariesProcess.hpp>
#include <corsika/media/ShowerAxis.hpp>

#include <corsika/modules/conex/CONEX_f.hpp>

namespace corsika {

  namespace conex {
    LengthType constexpr earthRadius{6371315 * meter};
  } // namespace conex

  class CONEXhybrid : public SecondariesProcess<CONEXhybrid> {

  public:
    CONEXhybrid(Point center, ShowerAxis const& showerAxis, LengthType groundDist,
                LengthType injectionHeight, HEPEnergyType primaryEnergy, PDGCode pdg);

    template <typename TStackView>
    void doSecondaries(TStackView&);

    void solveCE();

    bool addParticle(Code pid, HEPEnergyType energy, HEPEnergyType mass,
                     Point const& position, Vector<dimensionless_d> const& direction,
                     TimeType t);

    CoordinateSystemPtr const& getObserverCS() const { return conexObservationCS_; }

    HEPEnergyType getEnergyEM() const;
    void reset();

  private:
    // data members
    static std::array<std::pair<Code, int>, 3> constexpr egs_em_codes_{
        {{Code::Photon, 0}, {Code::Electron, -1}, {Code::Positron, -1}}};

    Point const center_; 
    ShowerAxis const& showerAxis_;
    LengthType groundDist_;  
    Point const showerCore_; 
    CoordinateSystemPtr const conexObservationCS_; 
    DirectionVector const x_sf_,
        y_sf_; 
    HEPEnergyType energy_em_;
  };
} // namespace corsika

#include <corsika/detail/modules/conex/CONEXhybrid.inl>