def METHOD_NAME(self):
    theta = 5  # initial guess
    delta_inverse = 1 / self.delta
    while ((theta - 3) / (theta - 2)) * math.factorial(theta) < delta_inverse:
        theta += 1
    while theta < np.e ** (self.epsilon / self.MAX_L) - 1:
        theta += 1
    print(f'Theta used by TrieHH: {theta}')
    return theta