CONEXhybrid.hpp

Go to the documentation of this file.

/*
 * (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/process/SecondariesProcess.hpp>
#include <corsika/framework/process/CascadeEquationsProcess.hpp>

#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/media/ShowerAxis.hpp>

#include <corsika/modules/writers/SubWriter.hpp>
#include <corsika/modules/writers/EnergyLossWriter.hpp>
#include <corsika/modules/writers/LongitudinalWriter.hpp>

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

namespace corsika {

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

  template <typename TOutputE, typename TOutputN>
  class CONEXhybrid : public CascadeEquationsProcess<CONEXhybrid<TOutputE, TOutputN>>,
                      public SecondariesProcess<CONEXhybrid<TOutputE, TOutputN>>,
                      public SubWriter<TOutputE>,
                      public SubWriter<TOutputN> {

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

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

    void initCascadeEquations();

    template <typename TStack>
    void doCascadeEquations(TStack& stack);

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

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

    YAML::Node getConfig() const final override;

  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_;       
    LengthType injectionHeight_;  
    HEPEnergyType primaryEnergy_; 
    PDGCode primaryPDG_;          
    Point const showerCore_;      
    CoordinateSystemPtr const conexObservationCS_; 
    DirectionVector const x_sf_,
        y_sf_; 
  };
} // namespace corsika

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