transition_tracer.h

// Copyright (C) 2024 Lisa Biermann, Margarete Mühlleitner, Rui Santos, João
// Viana
//
// SPDX-FileCopyrightText: 2024 Lisa Biermann, Margarete Mühlleitner, Rui
// Santos, João Viana
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#pragma once
 
#include "BSMPT/bounce_solution/bounce_solution.h" // BounceSolution
#include "BSMPT/gravitational_waves/gw.h"          // GravitationalWaves
#include "BSMPT/minimum_tracer/minimum_tracer.h"   // MinimumTracer
 
namespace BSMPT
{
 
struct user_input
{
  std::shared_ptr<Class_Potential_Origin> modelPointer;
  double T_low            = 0;
  double T_high           = 300;
  double vwall            = 0.95;
  double perc_prbl        = 0.71;
  double compl_prbl       = 0.01;
  double epsturb          = 0.1;
  int maxpathintegrations = 7;
  int multistepmode       = -1;
  int num_points          = 10;
  int ewsr_check          = 0;
  int nlo_check           = 1;
 
  int which_minimizer     = Minimizer::WhichMinimizerDefault;
  bool use_multithreading = false;
 
  bool gw_calculation = false;
  TransitionTemperature which_transition_temp =
      TransitionTemperature::Percolation;
  int PNLO_scaling                           = 1;
  size_t number_of_initial_scan_temperatures = 25;
};
 
struct status_codes
{
  StatusNLOStability status_nlo_stability = StatusNLOStability::NotSet;
  StatusEWSR status_ewsr                  = StatusEWSR::NotSet;
  StatusTracing status_tracing            = StatusTracing::NotSet;
  StatusCoexPair status_coex_pairs        = StatusCoexPair::NotSet;
  // index of vectors is coex_phase_id
  std::vector<StatusCrit> status_crit;
  std::vector<StatusGW> status_bounce_sol;
  std::vector<StatusTemperature> status_nucl_approx;
  std::vector<StatusTemperature> status_nucl;
  std::vector<StatusTemperature> status_perc;
  std::vector<StatusTemperature> status_compl;
};
 
struct transition_data
{
  std::optional<double> crit_temp;
  std::optional<double> nucl_approx_temp;
  std::optional<double> nucl_temp;
  std::optional<double> perc_temp;
  std::optional<double> compl_temp;
 
  std::vector<double> crit_true_vev;
  std::vector<double> crit_false_vev;
  std::vector<double> nucl_approx_true_vev;
  std::vector<double> nucl_approx_false_vev;
  std::vector<double> nucl_true_vev;
  std::vector<double> nucl_false_vev;
  std::vector<double> perc_true_vev;
  std::vector<double> perc_false_vev;
  std::vector<double> compl_true_vev;
  std::vector<double> compl_false_vev;
};
 
struct gw_data
{
  std::optional<double> vwall;
 
  std::optional<double> alpha;
  std::optional<double> beta_over_H;
 
  std::optional<double> kappa_col;
  std::optional<double> kappa_sw;
  std::optional<double> Epsilon_Turb;
  std::optional<double> cs_f;
  std::optional<double> cs_t;
 
  std::optional<double> fb_col;
  std::optional<double> omegab_col;
 
  std::optional<double> f1_sw;
  std::optional<double> f2_sw;
  std::optional<double> omega_2_sw;
 
  std::optional<double> f1_turb;
  std::optional<double> f2_turb;
  std::optional<double> omega_2_turb;
 
  std::optional<double> SNR_col;
  std::optional<double> SNR_sw;
  std::optional<double> SNR_turb;
  std::optional<double> SNR;
 
  StatusGW status_gw = StatusGW::NotSet;
  std::optional<double> trans_temp;
  std::optional<double> reh_temp;
};
 
struct output
{
  std::vector<std::string> legend;
  status_codes status;
  std::vector<transition_data> vec_trans_data;
  std::vector<gw_data> vec_gw_data;
  std::size_t num_coex_phase_pairs = 0;
 
  // stores string-identifier of phase id's throughout the transition history of
  // the universe
  std::string transition_history = "not_set";
};
 
class TransitionTracer
{
protected:
  std::shared_ptr<Class_Potential_Origin> modelPointer;
 
private:
  std::size_t num_vev;
 
public:
  std::vector<CoexPhases> vec_coex;
 
  TransitionTracer(user_input &input);
  ~TransitionTracer();
 
  std::vector<BounceSolution> ListBounceSolution;
 
  output output_store;
 
  double CheckMassRatio(const user_input &input,
                        const std::vector<double> &vec,
                        const double &temp) const;
};
 
} // namespace BSMPT