05/02/23 I dont' remember what this code did anymore. Is this a replacement of the mass_binary_analysis code? Treat all stars as binaries, resolve them to estimate their secondary masses and make the selection a posteriori