Essentially you need to run a counterfactual using simult_. See e.g. Counterfactual simulations based on smoothed shocks
simult_