reconstruction seems to be working; first simulation conditional on one family