code: def get_stddevs(self, C, nsites, mag, stddev_types):
          
          tau = self._get_tau(C, mag) + np.zeros(nsites)
          phi = self._get_phi(C, mag) + np.zeros(nsites)
          stddevs = []
          for stddev in stddev_types:
              assert stddev in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
              if stddev == const.StdDev.TOTAL:
                  stddevs.append(np.sqrt(tau ** 2. + phi ** 2.))
              elif stddev == const.StdDev.INTER_EVENT:
                  stddevs.append(tau)
              elif stddev == const.StdDev.INTRA_EVENT:
                  stddevs.append(phi)
          return stddevs
docs: Returns the standard deviations