Added the final piece. Now everything should be present in order to run on gaussian fit bases. Fishy energies tho. needs debugging.