Generating data from marginalized model with fixed value for the parameter of interest

Perhaps do is the tool for the job?

The functionality is already in PyMC

https://www.pymc.io/projects/docs/en/stable/api/model/generated/pymc.model.transform.conditioning.do.html