added control_simulation, and refactored ctrl_exp
population breakdown is now recorded