code for tau analysis at UCL