rf506_msgservice.C

Go to the documentation of this file.

/// \file
/// \ingroup tutorial_roofit
/// \notebook -nodraw
///  'ORGANIZATION AND SIMULTANEOUS FITS' RooFit tutorial macro #506
///
///  Tuning and customizing the RooFit message logging facility
///
/// \macro_output
/// \macro_code
/// \author 07/2008 - Wouter Verkerke 
 
 
#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooGaussian.h"
#include "RooConstVar.h"
#include "RooPolynomial.h"
#include "RooAddPdf.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "RooPlot.h"
#include "RooMsgService.h"
 
using namespace RooFit ;
 
 
void rf506_msgservice()
{
   // C r e a t e   p d f 
   // --------------------
 
   // Construct gauss(x,m,s)
   RooRealVar x("x","x",-10,10) ;
   RooRealVar m("m","m",0,-10,10) ;
   RooRealVar s("s","s",1,-10,10) ;
   RooGaussian gauss("g","g",x,m,s) ;
 
   // Construct poly(x,p0)
   RooRealVar p0("p0","p0",0.01,0.,1.) ;
   RooPolynomial poly("p","p",x,p0) ;      
 
   // Construct model = f*gauss(x) + (1-f)*poly(x)
   RooRealVar f("f","f",0.5,0.,1.) ;
   RooAddPdf model("model","model",RooArgSet(gauss,poly),f) ;
 
   RooDataSet* data = model.generate(x,10) ;
 
 
 
   // P r i n t   c o n f i g u r a t i o n   o f   m e s s a g e   s e r v i c e
   // ---------------------------------------------------------------------------
 
   // Print streams configuration
   RooMsgService::instance().Print() ;
   cout << endl ;
 
 
 
   // A d d i n g   I n t e g r a t i o n   t o p i c   t o   e x i s t i n g   I N F O   s t r e a m
   // -----------------------------------------------------------------------------------------------
 
   // Print streams configuration
   RooMsgService::instance().Print() ;
   cout << endl ;
 
   // Add Integration topic to existing INFO stream
   RooMsgService::instance().getStream(1).addTopic(Integration) ;
 
   // Construct integral over gauss to demonstrate new message stream
   RooAbsReal* igauss = gauss.createIntegral(x) ;
   igauss->Print() ;
 
   // Print streams configuration in verbose, which also shows inactive streams
   cout << endl ;
   RooMsgService::instance().Print() ;
   cout << endl ;
 
   // Remove stream
   RooMsgService::instance().getStream(1).removeTopic(Integration) ;
 
 
 
   // E x a m p l e s   o f   p d f   v a l u e   t r a c i n g   s t r e a m
   // -----------------------------------------------------------------------
   
   // Show DEBUG level message on function tracing, trace RooGaussian only
   RooMsgService::instance().addStream(DEBUG,Topic(Tracing),ClassName("RooGaussian")) ;
 
   // Perform a fit to generate some tracing messages
   model.fitTo(*data,Verbose(kTRUE)) ;
 
   // Reset message service to default stream configuration
   RooMsgService::instance().reset() ;
 
 
 
   // Show DEBUG level message on function tracing on all objects, redirect output to file
   RooMsgService::instance().addStream(DEBUG,Topic(Tracing),OutputFile("rf506_debug.log")) ;
 
   // Perform a fit to generate some tracing messages
   model.fitTo(*data,Verbose(kTRUE)) ;
 
   // Reset message service to default stream configuration
   RooMsgService::instance().reset() ;
 
 
 
   // E x a m p l e   o f   a n o t h e r   d e b u g g i n g   s t r e a m
   // ---------------------------------------------------------------------
 
   // Show DEBUG level messages on client/server link state management
   RooMsgService::instance().addStream(DEBUG,Topic(LinkStateMgmt)) ;
   RooMsgService::instance().Print("v") ;
 
   // Clone composite pdf g to trigger some link state management activity
   RooAbsArg* gprime = gauss.cloneTree() ;
   gprime->Print() ;
 
   // Reset message service to default stream configuration
   RooMsgService::instance().reset() ;
 
 
 
}