docs: def Wp(self):
          """Total energy in protons above 1.22 GeV threshold (erg).
          """
          from scipy.integrate import quad
  
          Eth = 1.22e-3
  
          with warnings.catch_warnings():
              warnings.simplefilter("ignore")
              Wp = quad(
                  lambda x: x * self._particle_distribution(x), Eth, np.Inf
              )[0]
  
          return (Wp * u.TeV).to("erg")