rf512_wsfactory_oper.C File Reference

Detailed Description

'ORGANIZATION AND SIMULTANEOUS FITS' RooFit tutorial macro #512

Illustration of operator expressions and expression-based basic p.d.f.s in the workspace factory syntax

0.576235055923
12.7450568676

Processing /mnt/vdb/lsf/workspace/root-makedoc-v608/rootspi/rdoc/src/v6-08-00-patches/tutorials/roofit/rf512_wsfactory_oper.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::autoImportClass(w) importing code of class RooCFAuto000Pdf from /mnt/vdb/lsf/workspace/root-makedoc-v608/rootspi/rdoc/src/v6-08-00-patches/documentation/doxygen/RooCFAuto000Pdf.cxx and RooCFAuto000Pdf.h

RooWorkspace(w) w contents

variables
---------
(Nbkg,Nsig,a,c,f,m,mg,s,u,v,x,y,z)

p.d.f.s
-------
RooGenericPdf::G[ actualVars=(x) formula="x*x+1" ] = 1
RooCFAuto000Pdf::GC[ x=x a=a ] = 1
RooChebychev::ch[ x=x coefList=(0.1,0.2,-0.3) ] = 0.8
RooAddPdf::extsummodel[ Nsig * gx + Nbkg * ch ] = 0.9
RooGaussian::g[ x=x mean=mg sigma=1 ] = 1
RooGaussian::gs[ x=x mean=m sigma=s ] = 1
RooGaussian::gx[ x=x mean=m sigma=1 ] = 1
RooProdPdf::gxycond[ gx * gy|x ] = 1
RooProdPdf::gxz[ gx * gz ] = 1
RooGaussian::gy[ x=y mean=x sigma=1 ] = 1
RooGaussian::gz[ x=z mean=0 sigma=1 ] = 1
RooLandau::lc[ x=x mean=0 sigma=1 ] = 0.178854
RooLandau::ls[ x=x mean=0 sigma=1 ] = 0.178854
RooFFTConvPdf::lxg[ g(x) (*) lc(x) ] = 0.172055
RooSimultaneous::smodel[ indexCat=c A=gs B=ls ] = 1
RooAddPdf::summodel[ f * gx + [%] * ch ] = 0.9

functions
--------
RooProduct::uv[ u * v ] = 100
RooAddition::uv2[ u + v ] = 20

embedded class code
-------------------
RooCFAuto000Pdf

embedded precalculated expensive components
-------------------------------------------
uid = 0 key=lxg_g_CONV_lc_CACHEHIST_Obs[x]_BufFrac0.1_BufStrat0 value=RooDataHist::lxg_g_CONV_lc_CACHEHIST_Obs[x]_BufFrac0.1_BufStrat0 parameters=( mg=0 )

#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooGaussian.h"
#include "RooConstVar.h"
#include "RooChebychev.h"
#include "RooAddPdf.h"
#include "RooWorkspace.h"
#include "RooPlot.h"
#include "TCanvas.h"
#include "TAxis.h"
using namespace RooFit ;


void rf512_wsfactory_oper()
{
   RooWorkspace* w = new RooWorkspace("w") ;

   // You can define typedefs for even shorter construction semantics
   w->factory("$Typedef(Gaussian,Gaus)") ;
   w->factory("$Typedef(Chebychev,Cheby)") ;


   // O p e r a t o r   p . d . f .   e x a m p l e s 
   // ------------------------------------------------

   // PDF addition is done with SUM (coef1*pdf1,pdf2)
   w->factory("SUM::summodel( f[0,1]*Gaussian::gx(x[-10,10],m[0],1.0), Chebychev::ch(x,{0.1,0.2,-0.3}) )") ;
   

   // Extended PDF addition is done with SUM (yield1*pdf1,yield2*pdf2)
   w->factory("SUM::extsummodel( Nsig[0,1000]*gx, Nbkg[0,1000]*ch )") ;
   

   // PDF multiplication is done with PROD ( pdf1, pdf2 ) 
   w->factory("PROD::gxz( gx, Gaussian::gz(z[-10,10],0,1) )") ;


   // Conditional p.d.f multiplication is done with PROD ( pdf1|obs, pdf2 )
   w->factory("Gaussian::gy( y[-10,10], x, 1.0 )") ;
   w->factory("PROD::gxycond( gy|x, gx )") ;


   // Convolution (numeric/ fft) is done with NCONV/FCONV (obs,pdf1,pdf2)
   w->factory("FCONV::lxg( x, Gaussian::g(x,mg[0],1), Landau::lc(x,0,1) )") ;


   // Simultaneous p.d.f.s are constructed with SIMUL( index, state1=pdf1, state2=pdf2,...)
   w->factory("SIMUL::smodel( c[A=0,B=1], A=Gaussian::gs(x,m,s[1]), B=Landau::ls(x,0,1) )") ; 




   // O p e r a t o r   f u n c t i o n   e x a m p l e s 
   // ---------------------------------------------------
   
   // Function multiplication is done with prod (func1, func2,...)
   w->factory("prod::uv(u[10],v[10])") ;


   // Function addition is done with sum(func1,func2)
   w->factory("sum::uv2(u,v)") ;



   // I n t e r p r e t e d   a n d   c o m p i l e d   e x p r e s s i o n   b a s e d   p . d . f . s .
   // ---------------------------------------------------------------------------------------------------

   // Create a RooGenericPdf interpreted p.d.f. You can use single quotes to pass the expression string argument
   w->factory("EXPR::G('x*x+1',x)") ;

   // Create a custom compiled p.d.f similar to the above interpreted p.d.f.
   // The code required to make this p.d.f. is automatically embedded in the workspace
   w->factory("CEXPR::GC('x*x+a',{x,a[1]})") ;


   // Compiled and interpreted functions (rather than p.d.f.s) can be made with the lower case 
   // 'expr' and 'cexpr' types

   // Print workspace contents
   w->Print() ;

   // Make workspace visible on command line
   gDirectory->Add(w) ;
}

Author

04/2009 - Wouter Verkerke

Definition in file rf512_wsfactory_oper.C.