OPT: clean a bit, add relevant plots, optimize fitting of mean PD in CSF