ClassPotentialCPintheDark.h

Go to the documentation of this file.

/*
 * ClassPotentialCPintheDark.h
 */
// Copyright (C) 2018  Philipp Basler and Margarete Mühlleitner
// SPDX-FileCopyrightText: 2021 Philipp Basler, Margarete Mühlleitner, Jonas
// Müller and Lisa Biermann
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#pragma once
 
#include <BSMPT/models/ClassPotentialOrigin.h>
namespace BSMPT
{
namespace Models
{
 
class Class_Potential_CPintheDark : public Class_Potential_Origin
{
public:
  Class_Potential_CPintheDark(const ISMConstants &smConstants);
  virtual ~Class_Potential_CPintheDark();
 
  // parameters of scalar potential
  double m11s, m22s, mSs, ReA, ImA, L1, L2, L3, L4, L5, L6, L7, L8;
  // counterterms
  double dm11s, dm22s, dmSs, dReA, dImA, dL1, dL2, dL3, dL4, dL5, dL6, dL7, dL8,
      dTCB, dT1, dT2, dTCP, dTS;
  double dImL5;
  // vev
  double v1;
 
  // dark sector mixing angles
  double alpha1 = 0, alpha2 = 0, alpha3 = 0;
 
  std::size_t pos_Gp, pos_Gm, pos_Hp, pos_Hm, pos_HSM;
  std::size_t pos_G0, pos_h1, pos_h2, pos_h3;
 
  void ReadAndSet(const std::string &linestr,
                  std::vector<double> &par) override;
  std::vector<std::string> addLegendCT() const override;
  std::vector<std::string> addLegendTemp() const override;
  std::vector<std::string> addLegendTripleCouplings() const override;
  std::vector<std::string> addLegendVEV() const override;
 
  void set_gen(const std::vector<double> &par) override;
  void set_CT_Pot_Par(const std::vector<double> &par) override;
  void write() const override;
 
  /*
   * CP in the Dark interaction basis:
   * 0     1     2     3     4      5     6      7     8
   * rho1, eta1, rho2, eta2, zeta1, psi1, zeta2, psi2, rhoS
   */
  const std::size_t pos_rho1 = 0, pos_eta1 = 1, pos_rho2 = 2, pos_eta2 = 3,
                    pos_zeta1 = 4, pos_psi1 = 5, pos_zeta2 = 6, pos_psi2 = 7,
                    pos_rhoS = 8;
 
  void FindMassBasisIndices(const std::vector<double> &HiggsMasses,
                            const Eigen::MatrixXd &HiggsRot);
 
  void AdjustRotationMatrix() override;
  void TripleHiggsCouplings() override;
  std::vector<double> calc_CT() const override;
 
  void SetCurvatureArrays() override;
  bool CalculateDebyeSimplified() override;
  bool CalculateDebyeGaugeSimplified() override;
  double VTreeSimplified(const std::vector<double> &v) const override;
  double VCounterSimplified(const std::vector<double> &v) const override;
  void Debugging(const std::vector<double> &input,
                 std::vector<double> &output) const override;
};
 
} // namespace Models
} // namespace BSMPT