Difference: HiggsAnalysisAtATLASUsingRooStats (14 vs. 15)

Example code: examining PDFs and creating graphical representations of them

// Create variables and a PDF using those variables.
	RooRealVar mu("mu", "mu", 150);
	RooRealVar sigma("sigma", "sigma", 5, 0, 20);
	RooGaussian myGaussianPDF("myGaussianPDF", "Gaussian PDF", x, mu, sigma);
// Create a Graphviz DOT file with a representation of the object tree.
	myGaussianPDF.graphVizTree("myGaussianPDFTree.dot");
	// This produced DOT file can be converted to some graphical representation:
		// Convert the DOT file to a 'top-to-bottom graph' using UNIX commands:
			// dot -Tgif -o myGaussianPDF_top-to-bottom_graph.gif myGaussianPDFTree.dot
		// Convert the DOT file to a 'spring-model graph' using UNIX commands:
			// fdp -Tgif -o myGaussianPDF_spring-model_graph.gif myGaussianPDFTree.dot
// Print the PDF contents.
	myGaussianPDF.Print();
	// Example output:
		// RooGaussian::G[ x=x mean=mu sigma=sigma ] = 1
		// sigma.Print()
		// RooRealVar::sigma = 5  L(0 - 20) 
// Print the PDF contents to stdout.
	myGaussianPDF.Print("t");
	// Example output:
		// 0x166eab0 RooGaussian::G = 1 [Auto] 
		// 0x15f7fe0/V- RooRealVar::x = 150
		// 0x1487090/V- RooRealVar::mu = 150
		// 0x1487bc0/V- RooRealVar::sigma = 5
// Print the PDF contents to a file.
	myGaussianPDF.printCompactTree("", "myGaussianPDFTree.txt")
	// Example output file contents:
		// 0x166eab0 RooGaussian::G = 1 [Auto] 
		// 0x15f7fe0/V- RooRealVar::x = 150
		// 0x1487090/V- RooRealVar::mu = 150
		// 0x1487bc0/V- RooRealVar::sigma = 5

Example code: Using both data and PDF from file

// Note that the following code is independent of actual PDF in the file. So, for example, a full Higgs combination could work with identical code.

// Open a file and import the workspace.
	TFile myFile("myResults.root") ;
	RooWorkspace* myWorkspace = f.Get("myWorkspace") ;
// Plot the data and PDF
	RooPlot* xframe = w->var("x")->frame() ;
	w->data("d")->plotOn(xframe) ;
	w->pdf("g")->plotOn(xframe) ;
// Construct a likelihood and profile likelihood
	RooNLLVar nll("nll","nll",*myWorkspace->pdf("g"),*w->data("d")) ;
	RooProfileLL pll("pll","pll", nll,*myWorkspace->var("m")) ;
	RooPlot* myFrame = w->var("m")->frame(-1,1) ;
	pll.plotOn(myFrame) ;
	myFrame->Draw()