start to split the likelihood into bitesize chunks
The current likelihood function does a lot. I think we could split it into neater chunks that will make it easier to use and adapt if we want to start analytically marginalising some parameters.
The current version shouldn't break anything but I feel like this should probably get some discussion before being merged.