From e2e070289864f7138b1ab223c81dfdbf267b4abb Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jo=C3=A3o=20Viana?=
 <57032457+vollous@users.noreply.github.com>
Date: Thu, 25 Apr 2024 17:53:19 +0100
Subject: [PATCH] Remove unused function

---
 .../BSMPT/bounce_solution/bounce_solution.h   |  7 -----
 src/bounce_solution/bounce_solution.cpp       | 29 -------------------
 2 files changed, 36 deletions(-)

diff --git a/include/BSMPT/bounce_solution/bounce_solution.h b/include/BSMPT/bounce_solution/bounce_solution.h
index 9098a73a..eea0c393 100644
--- a/include/BSMPT/bounce_solution/bounce_solution.h
+++ b/include/BSMPT/bounce_solution/bounce_solution.h
@@ -442,13 +442,6 @@ class BounceSolution
    */
   double GetPTStrength() const;
 
-  /**
-   * @brief CalcGstar Calculate the number of effective degrees of freedom as a
-   * function of temperature, higgs masses are evaluated at global minimum
-   * @param Temp temperature
-   */
-  void CalcGstar(const double &Temp);
-
   /**
    * @brief CalcGstarPureRad Calculate the number of effective degrees of
    * freedom assuming a purely radiative universe
diff --git a/src/bounce_solution/bounce_solution.cpp b/src/bounce_solution/bounce_solution.cpp
index 6e64c97e..5871667d 100644
--- a/src/bounce_solution/bounce_solution.cpp
+++ b/src/bounce_solution/bounce_solution.cpp
@@ -577,35 +577,6 @@ double BounceSolution::HubbleRate(const double &Temp)
          modelPointer->SMConstants.MPl; // radiation dominated universe
 }
 
-void BounceSolution::CalcGstar(const double &Temp)
-{
-  std::size_t NHiggs = this->modelPointer->get_NHiggs();
-
-  const std::vector<double> HiggsMassesSquared =
-      this->modelPointer->HiggsMassesSquared(
-          this->modelPointer->MinimizeOrderVEV(
-              this->modelPointer->get_vevTreeMin()),
-          Temp); // todo: calculate HiggsMassesSquared not in global, but in
-                 // minimum in which we are still in (e.g. local one until
-                 // T < Tperc)
-
-  std::size_t HiggsAboveThreshold = 0;
-
-  for (double mass : HiggsMassesSquared)
-  {
-    if (mass > Temp)
-    {
-      HiggsAboveThreshold++;
-    }
-  }
-
-  double gb   = 8 * 2 + 3 * 3 + 2 + NHiggs - HiggsAboveThreshold;
-  double gf   = 6 * 3 * 2 * 2 + 3 * 2 * 2 + 3 * 2;
-  double geff = gb + 7. / 8 * gf;
-
-  this->SetGstar(geff);
-}
-
 void BounceSolution::CalcGstarPureRad()
 {
   std::size_t NHiggs = this->modelPointer->get_NHiggs();