Sudbury Neutrino Observatory

Sudbury Neutrino Observatory
Carleton University Site

Run It! (Ensemble Test)

QSigEx has the ability to execute many fits quickly.
It is a nice feature which allow the user the ability to run
ensemble tests and to store the results in a ROOT file.

The following instructions assume QSIGEX is installed on your system.

Compile

Here is an example on how to run ensembles tests with qSigEx.

1) Go into the directory QSIGEXDIR/runme/mcgen

> cd $QSIGEXDIR/runme/mcgen

2) Compile the executable (bias)

> make

3) Read about it! [highly recommanded]

The ensemble test study is posted on manhattan (click here for the results).

4) Run It!

Run the ensemble test with no correlation

> bias bias_xxx_pdfs.dat bias_xxx_pdfs.root

Here xxx=ucor, corr, or 2dim. The results are in the root
file bias_xxx_pdfs.root which is produce at run time.

Run the ensemble test with correlation (rho=0.5)

> bias bias_xxx_pdfs.dat bias_xxx_pdfs.root 0.5

Here xxx=ucor, corr, or 2dim. The results are in the root
file bias_xxx_pdfs.root which is produce at run time.

The card files for qsigex are:

bias_ucor_pdfs.dat (1D fit)
bias_corr_pdfs.dat (PCA fit)
bias_2dim_pdfs.dat (2D fit)

The user MUST provide the dat file and the program bias will produce
the ROOT file.

Description of the code

    bias.cxx: Main program
    datagen.cxx: Produce the template PDFs
    generator.cxx: Generate the Monte Carlo using the acceptance-rejection method
    mcgen.cxx: Produce the samples for the ensemble test

    bias_2dim_pdfs.dat: Card file for the multi-D fits
    bias_corr_pdfs.dat: Card file for the PCA fits
    bias_ucor_pdfs.dat: Card file for the 1D fits

    local: Contains local libraries needed for the C routine rnmngn. Look in local/datan/datsrc.

    bias.C: An example of a ROOT macro to plot the results.

Results

Here's the results: POSTSCRIPT (gzip) or HTML