For ab-initio, you should be relaxing the resulting structures after abinitio in a full-atom scorefunction (score12). Then you want to cluster the structures + rescore them using the full-atom scorefunction, if you don’t have the scores from the relax step. For clustering, I would recommend calibur. Its quick and easy http://sourceforge.net/projects/calibur/
In the C++ application, this is known as Ab-Relax. To do this in PyRosetta, create + run a fastrelax instance after you bring the model to all-atom + repack:
rel = FastRelax(scorefxn)
rel.apply(pose)
Hope this helps