code: def prior_H0(self, H0, H0_min=0, H0_max=200):
          
          if H0 < H0_min or H0 > H0_max:
              penalty = -10**15
              return penalty, False
          else:
              return 0, True
docs: checks whether the parameter vector has left its bound, if so, adds a big number