ClassPotentialR2HDM.h

// Copyright (C) 2018  Philipp Basler and Margarete Mühlleitner
// SPDX-FileCopyrightText: 2021 Philipp Basler, Margarete Mühlleitner and Jonas
// Müller
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#ifndef SRC_CLASSPOTENTIALR2HDM_H_
#define SRC_CLASSPOTENTIALR2HDM_H_
 
#include <BSMPT/models/ClassPotentialOrigin.h>
using namespace Eigen;
 
namespace BSMPT
{
namespace Models
{
 
class Class_Potential_R2HDM : public Class_Potential_Origin
{
public:
  Class_Potential_R2HDM(const ISMConstants &smConstants);
  virtual ~Class_Potential_R2HDM();
 
  double L1 = 0, L2 = 0, L3 = 0, L4 = 0, RL5 = 0, RealMMix = 0, u1 = 0, u2 = 0;
  double DL1CT = 0, DL2CT = 0, DL3CT = 0, DL4CT = 0, DRL5CT = 0, Du2CT = 0,
         Du1CT = 0, DRu3CT = 0;
  double DT1 = 0, DT2 = 0, DT3 = 0, DTCharged = 0;
  double TanBeta = 0, C_CosBeta = 0, C_SinBeta = 0, C_CosBetaSquared = 0,
         C_SinBetaSquared = 0;
  double beta             = 0;
  long double Mh = 0, MH = 0, MA = 0, MHP = 0, alpha = 0;
  int Type       = 0;
  double CTempC1 = 0, CTempC2 = 0, CTempCS = 0;
 
  std::size_t pos_Gp, pos_Gm, pos_Hp, pos_Hm, pos_G0, pos_A, pos_H, pos_h;
  std::size_t pos_h_SM, pos_h_H;
 
  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;
 
  /*
   * R2HDM interaction basis:
   * 0     1     2     3     4      5     6      7
   * rho1, eta1, rho2, eta2, zeta1, psi1, zeta2, psi2
   */
  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;
 
  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
 
#endif /* SRC_CLASSPOTENTIALR2HDM_H_ */