Remove unused function

phbasler · Apr 25, 2024 · e2e0702 · e2e0702
1 parent 22564ec
commit e2e0702
Show file tree

Hide file tree

Showing 2 changed files with 0 additions and 36 deletions.
diff --git 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
@@ -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();