docs: def compute_Wp(self, Epmin=None, Epmax=None):
          """ Total energy in protons between energies Epmin and Epmax
  
          Parameters
          ----------
          Epmin : :class:`~astropy.units.Quantity` float, optional
              Minimum proton energy for energy content calculation.
  
          Epmax : :class:`~astropy.units.Quantity` float, optional
              Maximum proton energy for energy content calculation.
          """
          if Epmin is None and Epmax is None:
              Wp = self.Wp
          else:
              if Epmax is None:
                  Epmax = self.Epmax
              if Epmin is None:
                  Epmin = self.Epmin
  
              log10Epmin = np.log10(Epmin.to("GeV").value)
              log10Epmax = np.log10(Epmax.to("GeV").value)
              Ep = (
                  np.logspace(
                      log10Epmin,
                      log10Epmax,
                      int(self.nEpd * (log10Epmax - log10Epmin)),
                  )
                  * u.GeV
              )
              pdist = self.particle_distribution(Ep)
              Wp = trapz_loglog(Ep * pdist, Ep).to("erg")
  
          return Wp