rarMinuit Class Reference

#include <rarMinuit.hh>

List of all members.

Public Member Functions
	rarMinuit (RooAbsReal &function)
RooPlot *	contour (RooRealVar &var1, RooRealVar &var2, Double_t n1=1, Double_t n2=2, Double_t n3=0, Double_t n4=0, Double_t n5=0, Double_t n6=0)
void	fixGraph (TGraph *graph, Int_t lineStyle=1)
	make it close graph (fix the graph, looks a bug in Root)
Private Attributes
RooArgList *	_floatParamList
RooAbsReal *	_func

---

Detailed Description

Definition at line 27 of file rarMinuit.hh.

---

Constructor & Destructor Documentation

rarMinuit::rarMinuit

(

RooAbsReal &

function

)

Definition at line 40 of file rarMinuit.cc.

References _floatParamList, and _func.

                                         : RooMinuit(function)
{
  // The following are private in the base class so we have to set them up here
  _func = &function ;
  // Examine parameter list
  RooArgSet* paramSet = function.getParameters(RooArgSet()) ;
  RooArgList paramList(*paramSet) ;
  delete paramSet ;

  _floatParamList = (RooArgList*) paramList.selectByAttrib("Constant",kFALSE) ;
  if (_floatParamList->getSize()>1) {
    _floatParamList->sort() ;
  }
  _floatParamList->setName("floatParamList") ;
  
}

---

Member Function Documentation

RooPlot * rarMinuit::contour	(	RooRealVar &	var1,
		RooRealVar &	var2,
		Double_t	n1 = 1,
		Double_t	n2 = 2,
		Double_t	n3 = 0,
		Double_t	n4 = 0,
		Double_t	n5 = 0,
		Double_t	n6 = 0
	)

Definition at line 57 of file rarMinuit.cc.

References _floatParamList, _func, and fixGraph().

Referenced by rarMLFitter::doContourPlot().

{
  // Verify that both variables are floating parameters of PDF
  Int_t index1= _floatParamList->index(&var1);
  if(index1 < 0) {
    cout << "rarMinuit::contour(" << GetName() 
         << ") ERROR: " << var1.GetName()
         << " is not a floating parameter of " << _func->GetName() << endl ;
    return 0;
  }
  
  Int_t index2= _floatParamList->index(&var2);
  if(index2 < 0) {
    cout << "rarMinuit::contour(" << GetName() 
         << ") ERROR: " << var2.GetName()
         << " is not a floating parameter of PDF "
         << _func->GetName() << endl ;
    return 0;
  }
  
  // create and draw a frame
  RooPlot *frame = new RooPlot(var1, var2) ;
  frame->SetStats(kFALSE);
  
  // draw a point at the current parameter values
  TMarker *point= new TMarker(var1.getVal(), var2.getVal(), 8);
  
  // remember our original value of ERRDEF
  Double_t errdef= gMinuit->fUp;
  
  TGraph* graph1 = 0;
  if(n1 > 0) {
    // set the value corresponding to an n1-sigma contour
    gMinuit->SetErrorDef(n1*n1*errdef);
    // calculate and draw the contour
    graph1= (TGraph*)gMinuit->Contour(25, index1, index2);
    fixGraph(graph1);
  }
  
  TGraph* graph2 = 0;
  if(n2 > 0) {
    // set the value corresponding to an n2-sigma contour
    gMinuit->SetErrorDef(n2*n2*errdef);
    // calculate and draw the contour
    graph2= (TGraph*)gMinuit->Contour(25, index1, index2);
    fixGraph(graph2,2);
  }
  
  TGraph* graph3 = 0;
  if(n3 > 0) {
    // set the value corresponding to an n3-sigma contour
    gMinuit->SetErrorDef(n3*n3*errdef);
    // calculate and draw the contour
    graph3= (TGraph*)gMinuit->Contour(25, index1, index2);
    fixGraph(graph3,3);
  }
  
  TGraph* graph4 = 0;
  if(n4 > 0) {
    // set the value corresponding to an n4-sigma contour
    gMinuit->SetErrorDef(n4*n4*errdef);
    // calculate and draw the contour
    graph4= (TGraph*)gMinuit->Contour(25, index1, index2);
    fixGraph(graph4,4);
  }
  
  TGraph* graph5 = 0;
  if(n5 > 0) {
    // set the value corresponding to an n5-sigma contour
    gMinuit->SetErrorDef(n5*n5*errdef);
    // calculate and draw the contour
    graph5= (TGraph*)gMinuit->Contour(25, index1, index2);
    fixGraph(graph5,5);
  }
  
  TGraph* graph6 = 0;
  if(n6 > 0) {
    // set the value corresponding to an n6-sigma contour
    gMinuit->SetErrorDef(n6*n6*errdef);
    // calculate and draw the contour
    graph6= (TGraph*)gMinuit->Contour(25, index1, index2);
    fixGraph(graph6,6);
  }
  // restore the original ERRDEF
  gMinuit->SetErrorDef(errdef);
  
  
  // Add all objects to the frame
  frame->addObject(point);
  if (graph1) frame->addObject(graph1, "C");
  if (graph2) frame->addObject(graph2, "C");
  if (graph3) frame->addObject(graph3, "C");
  if (graph4) frame->addObject(graph4, "C");
  if (graph5) frame->addObject(graph5, "C");
  if (graph6) frame->addObject(graph6, "C");
  
  return frame;
}

void rarMinuit::fixGraph	(	TGraph *	graph,
		Int_t	lineStyle = 1
	)

make it close graph (fix the graph, looks a bug in Root)

Definition at line 159 of file rarMinuit.cc.

Referenced by contour().

                                                       {
  if (!graph) return;
  // set line style
  graph->SetLineStyle(lineStyle);
  // add a point to form closure
  Int_t nPoint=graph->GetN();
  graph->Set(nPoint+1);
  Double_t x, y;
  graph->GetPoint(0, x, y);
  graph->SetPoint(nPoint, x, y);
  return;
}

---

Member Data Documentation

RooArgList* rarMinuit::_floatParamList [private]

Definition at line 40 of file rarMinuit.hh.

Referenced by contour(), and rarMinuit().

RooAbsReal* rarMinuit::_func [private]

Definition at line 41 of file rarMinuit.hh.

Referenced by contour(), and rarMinuit().

---

The documentation for this class was generated from the following files:

• rarMinuit.hh
• rarMinuit.cc

---