1) Adds check if the cached data will be used to avoid loading the
search in the grid search
2) Adds better handling of errors in the par read
3) Adds ability to add 2d plots to an existing figure

Adds a switch to calculate (or optimize over) BSGL rather than 2F. This
should help to quiten noise events which exist only in one detector.

- Adds ability to truncate data to a specified period
- Removes the write-after option from the grid search
- Adds option to just return the data array from the grid search

It seems that , due to the choice of transient windows, occasionally
the calculated Fstat is +/- inf or a nan. This behaviour is not
understood. To address this, we apply a simple check and if it is an
undesirable value, return 0. This will inherently bias the chains away
from regions of parameter space in which the output breaks.

1) Improves the output
2) Fixes a bug in which the while loop got stcuk attempting to rename
the glitch parameters (when ng > 1)

If the initial p0 values are -np.inf due to the prior, it will iterate
through the other p0 values to replace the duff values. The replacement
takes another randomly selected p0 value and shifts it round by a
fractional gaussian with std=1e-10

1) Add rcparams to prevent all use of offsets
2) Fixes error when nglitch > 1 in which the summary did not show all
values
3) Fixes error in which the max_twoF used incorrect indexes

Also various improvements to logging commands

1) Use calculate_thetas to precalculate all thetas and then iterate
through to calculate twoF values. This will be important when we want to
allow theta to run both from and two theta0
2) Add diagnostic information on the mean acceptance fraction

1) Adds option to force the use of cached data
2) Orders the output .par files
3) Adds maxtwoF to output par files

Gives a more general fully-coherent grid search along with the
specialised grid glitch search

- Also adds examples to readme
- Fixes useOffset when ndim>1

Note: this also adds a check for transient, if not called, then the
whole data is searched. This does not however mean only transient, but
also glitch like searches which use the transient tools.

Additionally fixes the plot_walkers command to allow it to plot walkers
when ndim =1 and to switch off the use of offsets

This makes the MCMCGlitchSearch a subclass of the more general
MCMCSearch

1) Fixes grid search and adds example
2) Automates min/max Cover Freq to simplify input

Create a separate class for pure ComputeFstat wrapper, then gives this
as a parent to the other searches. Ultimately, this will allow the use
of the ComputeFstat class without needing the glitch part, for example
for a transient search.

Removes the fully-coherent search which was just a wrapper to CFS_v2 and
all of the `slow' methods which used command line arguments rather than
the python interface.