Write an equivalent of genPDBs.py for map to map fitting