def calc_theoretical_noise_dr(imax, theoretical_noise, nsigma):
    """
    Calculates the theoretical noise dynamic range
    theoretical_noise (float): the theoretical noise of the observation
    returns (float): the theoretical noise dynamic range
    """
    theoretical_noise_dr = imax / (nsigma * theoretical_noise)
    return theoretical_noise_dr
