num_gms = getOptimalModuleNumber_HC(
                        readcount_logscaled, hc, 1e10, numcores,
                        display=display,
                        pdf_plot_path=basename,
                        old_heuristic=FALSE
                        )