run_counterfactual function

Full counterfactual simulation run