rf506_msgservice.py

## \ingroup tutorial_roofit_main

## \notebook -nodraw

## Organization and simultaneous fits: tuning and customizing the ROOT.RooFit message logging facility

##

## \macro_code

## \macro_output

##

## \date February 2018

## \authors Clemens Lange, Wouter Verkerke (C++ version)

import ROOT

# Create pdf

# --------------------

# Construct gauss(x,m,s)

x = ROOT.RooRealVar("x", "x", -10, 10)

m = ROOT.RooRealVar("m", "m", 0, -10, 10)

s = ROOT.RooRealVar("s", "s", 1, -10, 10)

gauss = ROOT.RooGaussian("g", "g", x, m, s)

# Construct poly(x,p0)

p0 = ROOT.RooRealVar("p0", "p0", 0.01, 0.0, 1.0)

poly = ROOT.RooPolynomial("p", "p", x, [p0])

# model = f*gauss(x) + (1-f)*poly(x)

f = ROOT.RooRealVar("f", "f", 0.5, 0.0, 1.0)

model = ROOT.RooAddPdf("model", "model", [gauss, poly], [f])

data = model.generate({x}, 10)

# Print configuration of message service

# ------------------------------------------

# Print streams configuration

ROOT.RooMsgService.instance().Print()

# Adding integration topic to existing INFO stream

# ---------------------------------------------------

# Print streams configuration

ROOT.RooMsgService.instance().Print()

# Add Integration topic to existing INFO stream

ROOT.RooMsgService.instance().getStream(1).addTopic(ROOT.RooFit.Integration)

# Construct integral over gauss to demonstrate message stream

igauss = gauss.createIntegral({x})

igauss.Print()

# Print streams configuration in verbose, also shows inactive streams

ROOT.RooMsgService.instance().Print()

# Remove stream

ROOT.RooMsgService.instance().getStream(1).removeTopic(ROOT.RooFit.Integration)

# Examples of pdf value tracing

# -----------------------------------------------------------------------

# Show DEBUG level message on function tracing, ROOT.RooGaussian only

ROOT.RooMsgService.instance().addStream(ROOT.RooFit.DEBUG, Topic=ROOT.RooFit.Tracing, ClassName="RooGaussian")

# Perform a fit to generate some tracing messages

model.fitTo(data, Verbose=True)

# Reset message service to default stream configuration

ROOT.RooMsgService.instance().reset()

# Show DEBUG level message on function tracing on all objects, output to

# file

ROOT.RooMsgService.instance().addStream(ROOT.RooFit.DEBUG, Topic=ROOT.RooFit.Tracing, OutputFile="rf506_debug.log")

# Perform a fit to generate some tracing messages

model.fitTo(data, Verbose=True)

# Reset message service to default stream configuration

ROOT.RooMsgService.instance().reset()

# Example of another debugging stream

# ---------------------------------------------------------------------

# Show DEBUG level messages on client/server link state management

ROOT.RooMsgService.instance().addStream(ROOT.RooFit.DEBUG, Topic=ROOT.RooFit.LinkStateMgmt)

ROOT.RooMsgService.instance().Print("v")

# Clone composite pdf g to trigger some link state management activity

gprime = gauss.cloneTree()

gprime.Print()

# Reset message service to default stream configuration

ROOT.RooMsgService.instance().reset()