Line data    Source code

       1           0 : // Distributed under the MIT License.
       2             : // See LICENSE.txt for details.
       3             : 
       4             : #pragma once
       5             : 
       6             : #include <random>
       7             : #include <vector>
       8             : 
       9             : #include "DataStructures/DataBox/Tag.hpp"
      10             : #include "DataStructures/Tensor/TypeAliases.hpp"
      11             : #include "Evolution/Particles/MonteCarlo/MonteCarloOptions.hpp"
      12             : #include "Evolution/Particles/MonteCarlo/NeutrinoInteractionTable.hpp"
      13             : #include "Evolution/Particles/MonteCarlo/Packet.hpp"
      14             : 
      15             : /// Items related to the evolution of particles
      16             : /// Items related to Monte-Carlo radiation transport
      17             : /// Tags for MC
      18             : namespace Particles::MonteCarlo::Tags {
      19             : 
      20             : /// Simple tag containing the vector of Monte-Carlo
      21             : /// packets belonging to an element.
      22           1 : struct PacketsOnElement : db::SimpleTag {
      23           0 :   using type = std::vector<Particles::MonteCarlo::Packet>;
      24             : };
      25             : 
      26             : /// Simple tag containing an approximation of the light
      27             : /// crossing time for each cell (the shortest time among
      28             : /// all coordinate axis directions).
      29             : template <typename DataType>
      30           1 : struct CellLightCrossingTime : db::SimpleTag {
      31           0 :   using type = Scalar<DataType>;
      32             : };
      33             : 
      34             : /// Simple tag storing the coupling term between
      35             : /// MC and tilde_Tau (i.e. the energy variable)
      36             : template <typename DataType>
      37           1 : struct CouplingTildeTau : db::SimpleTag {
      38           0 :   using type = Scalar<DataType>;
      39             : };
      40             : 
      41             : /// Simple tag storing the coupling term between
      42             : /// MC and tilde_RhoYe (i.e. the composition variable)
      43             : template <typename DataType>
      44           1 : struct CouplingTildeRhoYe : db::SimpleTag {
      45           0 :   using type = Scalar<DataType>;
      46             : };
      47             : 
      48             : /// Simple tag storing the coupling term between
      49             : /// MC and tilde_S (i.e. the momentum variable)
      50             : template <typename DataType, size_t Dim>
      51           1 : struct CouplingTildeS : db::SimpleTag {
      52           0 :   using type = tnsr::i<DataType, Dim, Frame::Inertial>;
      53             : };
      54             : 
      55             : /// Simple tag storing the random number generator
      56             : /// used by Monte-Carlo
      57           1 : struct RandomNumberGenerator : db::SimpleTag {
      58           0 :   using type = std::mt19937;
      59             : };
      60             : 
      61             : /// Simple tag containing the desired energy of
      62             : /// packets in low-density regions. The energy
      63             : /// can be different for each neutrino species.
      64             : template <size_t NeutrinoSpecies>
      65           1 : struct DesiredPacketEnergyAtEmission : db::SimpleTag {
      66           0 :   using type = std::array<DataVector, NeutrinoSpecies>;
      67             : };
      68             : 
      69             : /// Simple tag for the table of neutrino-matter interaction
      70             : /// rates (emission, absorption and scattering for each
      71             : /// energy bin and neutrino species).
      72             : template <size_t EnergyBins, size_t NeutrinoSpecies>
      73           1 : struct InteractionRatesTable : db::SimpleTag {
      74           0 :   using type =
      75             :       std::unique_ptr<NeutrinoInteractionTable<EnergyBins, NeutrinoSpecies>>;
      76           0 :   static constexpr bool pass_metavariables = false;
      77           0 :   using option_tags =
      78             :       typename NeutrinoInteractionTable<EnergyBins, NeutrinoSpecies>::options;
      79           0 :   static type create_from_options(const std::string filename) {
      80             :     std::unique_ptr<Particles::MonteCarlo::NeutrinoInteractionTable<
      81             :         EnergyBins, NeutrinoSpecies>>
      82             :         interaction_table_ptr =
      83             :             std::make_unique<Particles::MonteCarlo::NeutrinoInteractionTable<
      84             :                 EnergyBins, NeutrinoSpecies>>(filename);
      85             :     return interaction_table_ptr;
      86             :     ;
      87             :   }
      88             : };
      89             : 
      90             : template <size_t NeutrinoSpecies>
      91           0 : struct MonteCarloOptions : db::SimpleTag {
      92           0 :   using type = std::unique_ptr<
      93             :       Particles::MonteCarlo::MonteCarloOptions<NeutrinoSpecies>>;
      94           0 :   static constexpr bool pass_metavariables = false;
      95           0 :   using option_tags = typename Particles::MonteCarlo::MonteCarloOptions<
      96             :       NeutrinoSpecies>::options;
      97           0 :   static type create_from_options(
      98             :       const std::array<double, NeutrinoSpecies> initial_packet_energy) {
      99             :     std::unique_ptr<Particles::MonteCarlo::MonteCarloOptions<NeutrinoSpecies>>
     100             :         mc_options_ptr = std::make_unique<
     101             :             Particles::MonteCarlo::MonteCarloOptions<NeutrinoSpecies>>(
     102             :             initial_packet_energy);
     103             :     return mc_options_ptr;
     104             :     ;
     105             :   }
     106             : };
     107             : 
     108             : }  // namespace Particles::MonteCarlo::Tags