ClassPotentialC2HDM.h

Go to the documentation of this file.

// 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_CLASSPOTENTIALC2HDM_H_
#define SRC_CLASSPOTENTIALC2HDM_H_
 
#include <BSMPT/models/ClassPotentialOrigin.h>
namespace BSMPT
{
namespace Models
{
 
class Class_Potential_C2HDM : public Class_Potential_Origin
{
public:
  Class_Potential_C2HDM(const ISMConstants &smConstants);
  virtual ~Class_Potential_C2HDM() override;
 
  bool UseHsmNotationInTripleHiggs = false;
 
  double L1 = 0, L2 = 0, L3 = 0, L4 = 0, RL5 = 0, RealMMix = 0, u1 = 0, u2 = 0;
  double IL5 = 0, Iu3 = 0;
  double DL1CT = 0, DL2CT = 0, DL3CT = 0, DL4CT = 0, DRL5CT = 0, Du2CT = 0,
         Du1CT = 0, DRu3CT = 0;
  double DIL5CT = 0, DIu3CT = 0;
  double DT1 = 0, DT2 = 0, DT3 = 0, DTCharged = 0;
  double DIL6CT  = 0;
  double TanBeta = 0, C_CosBeta = 0, C_SinBeta = 0, C_CosBetaSquared = 0,
         C_SinBetaSquared = 0;
  double beta             = 0;
  double M1 = 0, M2 = 0, M3 = 0, alpha1 = 0, alpha2 = 0, alpha3 = 0, MHC = 0;
  int Type       = 0;
  double CTempC1 = 0, CTempC2 = 0, CTempCS = 0;
  double R_Hh_1 = 0, R_Hh_2 = 0, R_Hh_3 = 0, R_Hl_1 = 0, R_Hl_2 = 0, R_Hl_3 = 0,
         R_Hsm_1 = 0, R_Hsm_2 = 0, R_Hsm_3 = 0;
 
  std::size_t pos_G0, pos_Gp, pos_Gm, pos_Hp, pos_Hm, pos_h1, pos_h2, pos_h3;
  std::size_t pos_h_SM, pos_h_l, 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;
 
  /*
   * C2HDM 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;
 
  /*
   * C2HDM semi-interaction basis (neutral Goldstone rotated out)
   * G^0, rho1, eta1, rho2, eta2, zeta1, zeta2, zeta3
   */
  const std::size_t pos_si_G0 = 0, pos_si_rho1 = 1, pos_si_eta1 = 2,
                    pos_si_rho2 = 3, pos_si_eta2 = 4, pos_si_zeta1 = 5,
                    pos_si_zeta2 = 6, pos_si_zeta3 = 7;
 
  void FindMassBasisIndices(const std::vector<double> &HiggsMasses,
                            const Eigen::MatrixXd &HiggsRot);
  std::tuple<std::size_t,
             std::size_t,
             std::size_t,
             std::size_t,
             std::size_t,
             std::size_t,
             std::size_t>
  FindSemiMassBasisIndices(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;
 
  bool IncludeChargeBreakingVEV = true;
  bool CTAlternative            = false;
  double DiffDelta              = 0;
 
  int PosSM = -1;
};
 
} // namespace Models
} // namespace BSMPT
#endif /* SRC_CLASSPOTENTIALC2HDM_H_ */