pmcmc code up and running