removed thermal mass from particle.py and inert doublet model

Models/InertDoubletModel/inertDoubletModel.py

@@ -125,15 +125,11 @@ def topMsqVacuum(fields: Fields) -> Fields:
         def topMsqDerivative(fields: Fields) -> Fields:
             return self.modelParameters["yt"] ** 2 * fields.getField(0)
 
-        def topMsqThermal(T: float) -> float:
-            return self.modelParameters["g3"] ** 2 * T**2 / 6.0
-
         topQuarkL = Particle(
             name="TopL",
             index=0,
             msqVacuum=topMsqVacuum,
             msqDerivative=topMsqDerivative,
-            msqThermal=topMsqThermal,
             statistics="Fermion",
             totalDOFs=6,
         )
@@ -144,7 +140,6 @@ def topMsqThermal(T: float) -> float:
             index=1,
             msqVacuum=topMsqVacuum,
             msqDerivative=topMsqDerivative,
-            msqThermal=topMsqThermal,
             statistics="Fermion",
             totalDOFs=6,
         )
@@ -157,15 +152,11 @@ def WMsqVacuum(fields: Fields) -> Fields:  # pylint: disable=invalid-name
         def WMsqDerivative(fields: Fields) -> Fields:  # pylint: disable=invalid-name
             return self.modelParameters["g2"] ** 2 * fields.getField(0) / 2
 
-        def WMsqThermal(T: float) -> float:  # pylint: disable=invalid-name
-            return self.modelParameters["g2"] ** 2 * T**2 * 11.0 / 6.0
-
         wBoson = Particle(
             name="W",
             index=4,
             msqVacuum=WMsqVacuum,
             msqDerivative=WMsqDerivative,
-            msqThermal=WMsqThermal,
             statistics="Boson",
             totalDOFs=9,
         )
@@ -177,16 +168,12 @@ def heavyScalarMsqVacuum(fields: Fields) -> Fields:
 
         def heavyScalarMsqDerivative(fields: Fields) -> Fields:
             return self.modelParameters["lambda3"]*fields.getField(0)
-
-        def heavyScalarThermal(T: float) -> float:
-            return self.modelParameters["lambda3"] * T**2 / 24.0
 
         heavyScalar = Particle(
             name="A",
             index=6,
             msqVacuum=heavyScalarMsqVacuum,
             msqDerivative=heavyScalarMsqDerivative,
-            msqThermal=heavyScalarThermal,
             statistics="Boson",
             totalDOFs=3,
         )

src/WallGo/particle.py

@@ -22,7 +22,6 @@ def __init__(
         index: int,
         msqVacuum: typing.Callable[[Fields | FieldPoint], np.ndarray],
         msqDerivative: typing.Callable[[Fields | FieldPoint], np.ndarray],
-        msqThermal: typing.Callable[[float], float],
         statistics: str,
         totalDOFs: int,
     ) -> None:
@@ -42,10 +41,6 @@ def __init__(
         msqDerivative : function
             Function :math:`d(m_0^2)/d(\phi)`, should take a Fields or FieldPoints
             object and return an array of shape Fields.shape.
-        msqThermal : function
-            Function :math:`m^2_T(T)`, should take a float and return one. The
-            temperature dependent but background field independent part of the
-            effective mass squared.
         statistics : {\"Fermion\", \"Boson\"}
             Particle statistics.
         totalDOFs : int
@@ -63,15 +58,13 @@ def __init__(
             index,
             msqVacuum,
             msqDerivative,
-            msqThermal,
             statistics,
             totalDOFs,
         )
         self.name = name
         self.index = index
         self.msqVacuum = msqVacuum
         self.msqDerivative = msqDerivative
-        self.msqThermal = msqThermal
         self.statistics = statistics
         self.totalDOFs = totalDOFs
 
@@ -81,7 +74,6 @@ def _validateInput(  # pylint: disable=unused-argument
         index: int,
         msqVacuum: typing.Callable[[Fields], np.ndarray],
         msqDerivative: typing.Callable[[Fields], np.ndarray],
-        msqThermal: typing.Callable[[float], float],
         statistics: str,
         totalDOFs: int,
     ) -> None:
@@ -92,13 +84,6 @@ def _validateInput(  # pylint: disable=unused-argument
         # assert isinstance(msqVacuum(fields), float), \
         #    f"msqVacuum({fields}) must return float"
 
-        # LN: comment mass check out to prevent errors at model creation time if no valid params have yet been passed
-        """
-        temperature = 100
-        assert isinstance(
-            msqThermal(temperature), float
-        ), f"msqThermal({temperature}) must return float"
-        """
         if statistics not in Particle.STATISTICS_OPTIONS:
             raise ValueError(f"{statistics=} not in {Particle.STATISTICS_OPTIONS}")
         assert isinstance(totalDOFs, int), "totalDOFs must be an integer"