Added test for statistics + fix bug variance

This is the basic example of the paper.
Also fixed a strange problem in computation of the variance:

I have identified cases where:

$$ 1/N sum_i (x_i- m)^2  \neq 1/N sum_i x_i^2 - m^2 $$

with $ m = 1/N sum_i x_i $

In particular, the RHS (which corresponds to previous implementation) may give **negative** variances!! LHS seems to be more accurate (especially in cases where x_i is very close to m) --> changed implementation to LHS
[doc] update to preliminary introduction, prettification
[doc] Fixed example #99
[doc] stated faqs.rst with text file printing info
rm parameter class, fix #124

- now replaced by simpler and better code generation solution
- TODO: write the documentation, implemented for cthyb solver
[wrapper] Split the converters code, fix #121 and hello_world

- The wrapper_tools.hpp is split into several converters/xxx.hpp
- to include only on demand.
- This fixes the hello_world from the paper