code: def _get_stddevs(self, C, stddev_types, stddev_shape):
          
          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(C["sigtot"] + np.zeros(stddev_shape))
              elif stddev_type == const.StdDev.INTRA_EVENT:
                  stddevs.append(C['sig2'] + np.zeros(stddev_shape))
              elif stddev_type == const.StdDev.INTER_EVENT:
                  stddevs.append(C['sig1'] + np.zeros(stddev_shape))
          return stddevs
docs: Returns the standard deviations given in Table 2