code: def _get_stddevs(self, C, stddev_types, num_sites, mag, c1_rrup,
                       log_phi_ss, mean_phi_ss):
          
          phi_ss = _compute_phi_ss(C, mag, c1_rrup, log_phi_ss, mean_phi_ss)
  
          stddevs = []
          for stddev_type in stddev_types:
              assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
  
              if stddev_type == const.StdDev.TOTAL:
                  stddevs.append(np.sqrt(
                                 C['tau'] * C['tau'] +
                                 phi_ss * phi_ss) +
                                 np.zeros(num_sites))
  
              elif stddev_type == const.StdDev.INTRA_EVENT:
                  stddevs.append(phi_ss + np.zeros(num_sites))
  
              elif stddev_type == const.StdDev.INTER_EVENT:
                  stddevs.append(C['tau'] + np.zeros(num_sites))
          return stddevs
docs: Return standard deviations