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