Rescale x-axis of phase space graph

The unit of the x-axis of the phase space integral for a wave was MeV, while
for all other graphs it is GeV. Corrected this.
Removed larger portion of commented code

Removed some commented code.
Removed not required argument to 'getPhaseSpace'

Removed one argument that was not used in 'getPhaseSpace', and modified the
calls accordingly.
Removed a member variable from class 'massDepFitLikeli'

Removed the '_finalStatePS' member variable from the class 'massDepFitLikeli'
as it was not used there. The final state phase space is nowadays introduced
via the 'pwacompset'.
Removed a few more variables

Removed a few more unused variables.
MassDepFit: replace hard-coded phase-space parametrization

Replaced the so-far hard-coded phase-space parametrization by the possibility
to define an overall final-state mass dependence of the intensities. This
final-state mass dependence can also be omitted. Some renaming had to be done
across the 'pwacompset' class and the main routine of the executable itself.
Correct handling of zero length arrays from config

Only check whether the elements of an array from the configuration file are
numbers, if there is at least one element.
Wave-specific mass-range for mass-dependent fit

It is now possible to set the fitted mass-range individually for each wave in
the mass-dependent fit. This required a number of changes to basically all
parts of the mass-dependent fitter.

In the configuration file one now has to specify the waves to be used in the
fit. At this place the mass-range can be specified with the options 'massLower'
and 'massUpper'. The example '5piMassDep.conf' has been adjusted accordingly
without defining any mass-range.

The default mass-range is now extracted from the mass-range covered by the bins
of the fit result, and is no longer hard-coded from 0 to 5000 MeV/c^2. Some
additional output has been added, so the result and correctness of this can
easily be checked.

The likelihood/chi^2 calculation has been adjusted to no take into account
mass-bins that are not in the range of each pair of waves, the highest lower
bound and the lowest higher bound define which mass-bins are used. A new
command-line option "-R" has been introduced controlling whether to plot the
result only in the specified mass-range or the total.

There is one issue if the plotting is not done of the complete mass-range. So
far the graphs are always filled at the index of the mass-bin, so that entries
with x- and y-values of 0 might exist in the beginning and end of the graph.
Combined with the drawing option this might create some strange loops in the
plots. Those are purely artifacts of the plotting, not of the fit itself.
Changed application of final-state mass-dependence

If the spin-density matrix is written as
                                 rho = T x T*
then the final-stat mass-dependence so far has been applied directly to the
product
                              rho = T x T* x fsmd
This has now been changed to
                         rho = (T x fsmd) x (T* x fsmd)
so the final-state mass-dependence now works on the "production amplitudes"
directly, and not on the spin-density matrix. This should make it simpler to in
the future experiment with directly fitting the production amplitudes, and not
the spin-density matrix.
Limited plot-range also for individual components

Fixed a bug that prevented the individual components from being drawn in the
specified mass-range only.
Cleaning command-line options

Moved the specification of the fit result from the mass-independent fit and the
specification of the files used for the plotting of the systematic errors from
the command-line to the configuration file. Removed the command-line switch
'-i', the name of the configuration file is now taken from the single
not-processed option.
Correct number of data points for limited mass-ranges

The calculation of the number of data points still assumed that the fit was
done over the complete mass range. Now the limits are correctly taken into
account.
A bit of code cleaning

Removed a couple of unused class member variables in 'massDepFitLikeli', and
also removed some commented code no longer required.
Re-structure massDepFit.cc

Added a new class 'massDepFit' which is suppossed to do everything that is
currently done in the 'main' function. So far most of what the class does is
not yet used, so for the moment the same things are done twice. This will of
course be removed later on.
Change source of input for 'massDepFitLikeli'

The spin-density matrices used as an input for 'massDepFitLikeli' are now
created in the 'massDepFit' class. Cleaned the source of the 'massDepFitLikeli'
class, all calculations during the initialisation are now done outside, the
'init' method only sets some member variables.

The 'overlap' method of 'pwacompset' can now take an additional argument, the
bin number of the current mass bin. This might be used in the future to use the
value of the phase-space from an array, and not to have to calculate it from a
TGraph.
Added reverse mapping of wave name to wave index

A map was introduced to allow mapping of a wave name to the wave index.
Removed 'getPhaseSpace' function

Replaced the 'getPhaseSpace' function by reading the phase-space integrals in
the 'massDepFit' class. The new methods no longer provide a TGraph, but a plain
vector of doubles instead, so the 'pwachannel' class had to be adopted.
Replace vectors of vectors (of vectors ...) by boost::multi_array

Part I of this work.
Replace vectors of vectors (of vectors ...) by boost::multi_array

Part II of this work.
Create graphs from content of 'massDepFit' class

The plots are now created from the content of the 'massDepFit' class, and no
longer from the content of the TTree itself. The next step is to also achieve
this for the systematic plots.