ClassPotentialN2HDM.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_CLASSPOTENTIALN2HDM_H_
#define SRC_CLASSPOTENTIALN2HDM_H_
 
#include <BSMPT/models/ClassPotentialOrigin.h>
 
namespace BSMPT
{
namespace Models
{
 
class Class_Potential_N2HDM : public Class_Potential_Origin
{
public:
  Class_Potential_N2HDM(const ISMConstants &smConstants);
  virtual ~Class_Potential_N2HDM();
 
  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;
  double TanBeta = 0, C_CosBeta = 0, C_SinBeta = 0, C_CosBetaSquared = 0,
         C_SinBetaSquared = 0;
  double beta             = 0;
  int Type                = 0;
  double CTempC1 = 0, CTempC2 = 0, CTempCS = 0;
  double alpha1 = 0, alpha2 = 0, alpha3 = 0;
  double MSM = 0, MhUp = 0, MhDown = 0;
 
  double Nus = 0, NL6 = 0, NL7 = 0, NL8 = 0, Nvs = 0;
  double NDus = 0, NDL6 = 0, NDL7 = 0, NDL8 = 0, NDvs = 0, NDTS = 0;
  double DTCharged = 0;
 
  std::size_t pos_G0, pos_Gp, pos_Gm, pos_Hp, pos_Hm;
  std::size_t pos_h1, pos_h2, pos_h3, pos_A;
  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;
 
  /*
   * N2HDM interaction basis:
   * 0     1     2     3     4      5     6      7     8
   * rho1, rho2, eta1, eta2, psi1, psi2, zeta1, zeta2, rhoS
   */
  const std::size_t pos_rho1 = 0, pos_rho2 = 1, pos_eta1 = 2, pos_eta2 = 3,
                    pos_psi1 = 4, pos_psi2 = 5, pos_zeta1 = 6, pos_zeta2 = 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
#endif /* SRC_CLASSPOTENTIALN2HDM_H_ */