Logo ROOT   6.12/07
Reference Guide
rs601_HLFactoryexample.C File Reference
Tutorials » RooStats Tutorials

Detailed Description

View in nbviewer Open in SWAN 'High Level Factory Example' RooStats tutorial macro #601

This tutorial shows an example of creating a simple model using the High Level model Factory.

pict1_rs601_HLFactoryexample.C.png
Processing /mnt/build/workspace/root-makedoc-v612/rootspi/rdoc/src/v6-12-00-patches/tutorials/roostats/rs601_HLFactoryexample.C...
RooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby
Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University
All rights reserved, please read http://roofit.sourceforge.net/license.txt
[#1] INFO:ObjectHandling -- RooWorkspace::exportToCint(HLFavtoryexample_ws) INFO: references to all objects in this workspace will be created in CINT in 'namespace HLFavtoryexample_ws'
[#1] INFO:Minization -- p.d.f. provides expected number of events, including extended term in likelihood.
[#1] INFO:Minization -- createNLL: caching constraint set under name CONSTR_OF_PDF_sum_FOR_OBS_mes with 0 entries
[#1] INFO:Minization -- RooMinimizer::optimizeConst: activating const optimization
[#1] INFO:Minization -- The following expressions will be evaluated in cache-and-track mode: (gauss,argus)
**********
** 1 **SET PRINT 1
**********
**********
** 2 **SET NOGRAD
**********
PARAMETER DEFINITIONS:
NO. NAME VALUE STEP SIZE LIMITS
1 argpar -2.00000e+01 9.50000e+00 -1.00000e+02 -1.00000e+00
2 mean 5.28000e+00 1.00000e-02 5.20000e+00 5.30000e+00
3 nbkg 8.00000e+02 4.00000e+02 0.00000e+00 1.00000e+04
4 nsig 2.00000e+02 1.00000e+02 0.00000e+00 1.00000e+04
5 width 2.70000e-03 8.50000e-04 1.00000e-03 1.00000e+00
**********
** 3 **SET ERR 0.5
**********
**********
** 4 **SET PRINT 1
**********
**********
** 5 **SET STR 1
**********
NOW USING STRATEGY 1: TRY TO BALANCE SPEED AGAINST RELIABILITY
**********
** 6 **MIGRAD 2500 1
**********
FIRST CALL TO USER FUNCTION AT NEW START POINT, WITH IFLAG=4.
START MIGRAD MINIMIZATION. STRATEGY 1. CONVERGENCE WHEN EDM .LT. 1.00e-03
FCN=-17788.8 FROM MIGRAD STATUS=INITIATE 20 CALLS 21 TOTAL
EDM= unknown STRATEGY= 1 NO ERROR MATRIX
EXT PARAMETER CURRENT GUESS STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 argpar -2.00000e+01 9.50000e+00 2.51381e-01 1.48295e+01
2 mean 5.28000e+00 1.00000e-02 2.57889e-01 1.26241e+02
3 nbkg 8.00000e+02 4.00000e+02 1.52384e-01 -2.75719e+03
4 nsig 2.00000e+02 1.00000e+02 7.39156e-02 -1.30866e+03
5 width 2.70000e-03 8.50000e-04 2.13708e-02 -5.96944e+02
ERR DEF= 0.5
MIGRAD MINIMIZATION HAS CONVERGED.
MIGRAD WILL VERIFY CONVERGENCE AND ERROR MATRIX.
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=-18178.2 FROM MIGRAD STATUS=CONVERGED 127 CALLS 128 TOTAL
EDM=8.23524e-07 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER STEP FIRST
NO. NAME VALUE ERROR SIZE DERIVATIVE
1 argpar -2.28748e+01 3.42311e+00 6.81552e-03 -1.05902e-04
2 mean 5.27987e+00 2.15784e-04 4.99568e-04 1.14694e-01
3 nbkg 1.61271e+03 4.46744e+01 1.05575e-03 -2.82924e-02
4 nsig 3.87281e+02 2.77659e+01 1.13759e-03 3.18941e-02
5 width 3.01017e-03 1.99068e-04 3.69588e-04 2.63544e-02
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5
1.174e+01 -9.561e-06 -4.054e+01 4.054e+01 2.300e-04
-9.561e-06 4.656e-08 3.241e-04 -3.248e-04 -2.699e-09
-4.054e+01 3.241e-04 1.996e+03 -3.831e+02 -2.181e-03
4.054e+01 -3.248e-04 -3.831e+02 7.710e+02 2.184e-03
2.300e-04 -2.699e-09 -2.181e-03 2.184e-03 3.963e-08
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4 5
1 0.47848 1.000 -0.013 -0.265 0.426 0.337
2 0.07490 -0.013 1.000 0.034 -0.054 -0.063
3 0.35856 -0.265 0.034 1.000 -0.309 -0.245
4 0.52842 0.426 -0.054 -0.309 1.000 0.395
5 0.45058 0.337 -0.063 -0.245 0.395 1.000
**********
** 7 **SET ERR 0.5
**********
**********
** 8 **SET PRINT 1
**********
**********
** 9 **HESSE 2500
**********
COVARIANCE MATRIX CALCULATED SUCCESSFULLY
FCN=-18178.2 FROM HESSE STATUS=OK 31 CALLS 159 TOTAL
EDM=8.27129e-07 STRATEGY= 1 ERROR MATRIX ACCURATE
EXT PARAMETER INTERNAL INTERNAL
NO. NAME VALUE ERROR STEP SIZE VALUE
1 argpar -2.28748e+01 3.42862e+00 2.72621e-04 5.92076e-01
2 mean 5.27987e+00 2.15791e-04 9.99136e-05 6.40145e-01
3 nbkg 1.61271e+03 4.47114e+01 4.22301e-05 -7.44300e-01
4 nsig 3.87281e+02 2.78173e+01 2.27518e-04 -1.17462e+00
5 width 3.01017e-03 1.99415e-04 1.47835e-05 -1.48105e+00
ERR DEF= 0.5
EXTERNAL ERROR MATRIX. NDIM= 25 NPAR= 5 ERR DEF=0.5
1.178e+01 -9.115e-06 -4.092e+01 4.092e+01 2.327e-04
-9.115e-06 4.657e-08 3.206e-04 -3.207e-04 -2.752e-09
-4.092e+01 3.206e-04 1.999e+03 -3.865e+02 -2.207e-03
4.092e+01 -3.207e-04 -3.865e+02 7.739e+02 2.207e-03
2.327e-04 -2.752e-09 -2.207e-03 2.207e-03 3.977e-08
PARAMETER CORRELATION COEFFICIENTS
NO. GLOBAL 1 2 3 4 5
1 0.48108 1.000 -0.012 -0.267 0.429 0.340
2 0.07539 -0.012 1.000 0.033 -0.053 -0.064
3 0.36057 -0.267 0.033 1.000 -0.311 -0.248
4 0.53093 0.429 -0.053 -0.311 1.000 0.398
5 0.45364 0.340 -0.064 -0.248 0.398 1.000
[#1] INFO:Minization -- RooMinimizer::optimizeConst: deactivating const optimization
[#1] INFO:Plotting -- RooAbsPdf::plotOn(sum) directly selected PDF components: (argus)
[#1] INFO:Plotting -- RooAbsPdf::plotOn(sum) indirectly selected PDF components: ()
#include <fstream>
#include "TString.h"
#include "TROOT.h"
#include "RooGlobalFunc.h"
#include "RooWorkspace.h"
#include "RooRealVar.h"
#include "RooAbsPdf.h"
#include "RooDataSet.h"
#include "RooPlot.h"
#include "RooStats/HLFactory.h"
// use this order for safety on library loading
using namespace RooFit;
using namespace RooStats;
using namespace std;
void rs601_HLFactoryexample() {
// --- Build the datacard and dump to file---
TString card_name("HLFavtoryexample.rs");
ofstream ofile(card_name);
ofile << "// The simplest card\n\n"
<< "gauss = Gaussian(mes[5.20,5.30],mean[5.28,5.2,5.3],width[0.0027,0.001,1]);\n"
<< "argus = ArgusBG(mes,5.291,argpar[-20,-100,-1]);\n"
<< "sum = SUM(nsig[200,0,10000]*gauss,nbkg[800,0,10000]*argus);\n\n";
ofile.close();
HLFactory hlf("HLFavtoryexample",
card_name,
false);
// --- Take elements out of the internal workspace ---
RooWorkspace* w = hlf.GetWs();
RooRealVar* mes = dynamic_cast<RooRealVar*>(w->arg("mes"));
RooAbsPdf* sum = dynamic_cast<RooAbsPdf*>(w->pdf("sum"));
RooAbsPdf* argus = dynamic_cast<RooAbsPdf*>(w->pdf("argus"));
// RooRealVar* mean = dynamic_cast<RooRealVar*>(w->arg("mean"));
// RooRealVar* argpar = dynamic_cast<RooRealVar*>(w->arg("argpar"));
// --- Generate a toyMC sample from composite PDF ---
RooDataSet *data = sum->generate(*mes,2000) ;
// --- Perform extended ML fit of composite PDF to toy data ---
sum->fitTo(*data) ;
// --- Plot toy data and composite PDF overlaid ---
RooPlot* mesframe = mes->frame() ;
data->plotOn(mesframe) ;
sum->plotOn(mesframe) ;
sum->plotOn(mesframe,Components(*argus),LineStyle(kDashed)) ;
gROOT->SetStyle("Plain");
mesframe->Draw() ;
}
Author
Danilo Piparo

Definition in file rs601_HLFactoryexample.C.