kreczko / roounfold

Clone from http://hepunx.rl.ac.uk/~adye/software/unfold/RooUnfold.html

RooUnfold

Fork from http://hepunx.rl.ac.uk/~adye/software/unfold/RooUnfold.html

RooUnfold is a framework for unfolding (AKA "deconvolution" or "unsmearing"). It currently implements several methods: iterative (Bayes), singular value decomposition (SVD), unregularised matrix inversion, bin-by-bin (simple correction factors), and an interface to ROOT's TUnfold. It can be used from the ROOT prompt, or linked against the ROOT libraries. RooUnfold requires ROOT 4 or later.

For the latest version and documentation, see here

http://hepunx.rl.ac.uk/~adye/software/unfold/RooUnfold.html

RooUnfold was written by Tim Adye, Kerstin Tackmann, and Fergus Wilson.

Building

Since version 2 RooUnfold is using cmake to find dependencies and create the MakeFile:

cmake CMakeLists.txt
make jobs=2

to build the RooUnfold shared library.

Tests

In order to run the boost tests for RooUnfold simply

make test

If you want to run specific tests you can do this with

# a single test suite
./RooUnfold_test --log_level=message --run_test=TestRooUnfoldSvd
# a single test case
./RooUnfold_test --log_level=message --run_test=TestRooUnfoldSvd/test_get_tau_from_constructor

Running

RooUnfoldExample.cxx makes a simple test of RooUnfold.

% root root [0] .L examples/RooUnfoldExample.cxx root [1] gSystem->Load("libRooUnfold.so"); root [2] RooUnfoldExample()

See the web page (http://hepunx.rl.ac.uk/~adye/software/unfold/RooUnfold.html) for more examples and documentation.

roounfold's Issues

Root 6 problems

• rlibmap is no longer provided
• TUnfold::GetOutput function changed:

RooUnfold/src/RooUnfoldTUnfold.cxx:202:60: error: no matching function for call to ‘TUnfold::GetOutput(const char [5], const char [19], int, int)’
   TH1* reco=_unf->GetOutput("_rec","reconstructed dist",0,0);
                                                            ^

Synchronise Eclipse project configuration with GNUmakefile

Some links are currently missing in the Eclipse project configuration and prevent the library from being loaded in ROOT.

Hreco(RooUnfold::kNoError) returns 0 errors

In its current implementation RooUnfoldSvd::Hreco(RooUnfold::kNoError) returns errors of value 0 for all bins. This is in conflict with the comments in RooUnfold.h:

enum ErrorTreatment {  // Error treatment:
    kNoError,            //   no error treatment: returns sqrt(N)
    kErrors,             //   bin-by-bin errors (diagonal covariance matrix)
    kCovariance,         //   covariance matrix from unfolding
    kCovToy              //   covariance matrix from toy MC
};