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.

(cherry picked from commit f69a72b6fc33b20cab6986a96132a2dc513b4bc1)
Removed larger portion of commented code

Removed some commented code.

(cherry picked from commit 7106ffec0ba30fb71003a90749f89fcc15a2f03f)
Removed not required argument to 'getPhaseSpace'

Removed one argument that was not used in 'getPhaseSpace', and modified the
calls accordingly.

(cherry picked from commit 3466fa1d889cfa86210632d207a776201f755f60)
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'.

(cherry picked from commit 44e1477a2b5391b5cd1489b5da7298efd8d7cca9)
Removed a few more variables

Removed a few more unused variables.

(cherry picked from commit 310fda712ce3cb5590b9fb6432c0c850cbb9f030)
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.

(cherry picked from commit b20206fc071359dcb2eb554095b06ced4a678654)
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.

(cherry picked from commit 81420534f97081f95f9b756a42483bd6ba40af54)
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.

(cherry picked from commit fe2dcf3255c37610d903fdb56bdb8a9a82e3bc71)
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.

(cherry picked from commit b08f825f4f5f4dd3cd6c296150a4173f6b053164)
Limited plot-range also for individual components

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

(cherry picked from commit c8257bdf0b26619c754a09727d43e4f305b49351)
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.

(cherry picked from commit 7ee62d9727388835d936a18a5f2cf2e9b6244346)
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.

(cherry picked from commit 803b974cbbb4f288936c2e904cda496b8b007aa5)
A bit of code cleaning

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

(cherry picked from commit cb9c98b35a249073ebab45360f37c7f879410cb2)
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.

(cherry picked from commit aba4f52483398558c79a723d31da49cc2cfb2a53)
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.

(cherry picked from commit bd4544ab655090ad3171c323fb1fab412241e273)
Added reverse mapping of wave name to wave index

A map was introduced to allow mapping of a wave name to the wave index.

(cherry picked from commit 57ae516677fe1090b83a2a2b39decb5f31221bf2)
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.

(cherry picked from commit 39e58560b1a93ad657fdb025b2ce351928647259)
Replace vectors of vectors (of vectors ...) by boost::multi_array

Part I of this work.

(cherry picked from commit 205cfd3d75df3113a4c19a7c3a87a9cd80871f92)
Replace vectors of vectors (of vectors ...) by boost::multi_array

Part II of this work.

(cherry picked from commit 591fe09f5060e8a35e04469a29385a64f552cd44)
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.

(cherry picked from commit 92ef8a9182b14d1a980b8400c62368f941d8267d)