ROOT logo
// @(#)root/roostats:$Id: LikelihoodIntervalPlot.cxx 31876 2009-12-14 11:11:27Z brun $

/*************************************************************************
 * Project: RooStats                                                     *
 * Package: RooFit/RooStats                                              *
 * Authors:                                                              *
 *   Kyle Cranmer, Lorenzo Moneta, Gregory Schott, Wouter Verkerke       *
 *************************************************************************
 * Copyright (C) 1995-2008, Rene Brun and Fons Rademakers.               *
 * All rights reserved.                                                  *
 *                                                                       *
 * For the licensing terms see $ROOTSYS/LICENSE.                         *
 * For the list of contributors see $ROOTSYS/README/CREDITS.             *
 *************************************************************************/

//____________________________________________________________________
/*
LikelihoodIntervalPlot : 

This class provides simple and straightforward utilities to plot a LikelihoodInterval
object.
*/

#include "RooStats/LikelihoodIntervalPlot.h"

#include <algorithm>
#include <iostream>

#include "TROOT.h"
#include "TMath.h"
#include "TLine.h"
#include "TObjArray.h"
#include "TList.h"
#include "TGraph.h"
#include "TPad.h"


#include "RooRealVar.h"
#include "RooPlot.h"
//#include "RooProfileLL.h"
#include "TF1.h"

/// ClassImp for building the THtml documentation of the class 
ClassImp(RooStats::LikelihoodIntervalPlot);

using namespace RooStats;

//_______________________________________________________
LikelihoodIntervalPlot::LikelihoodIntervalPlot()
{
  // LikelihoodIntervalPlot default constructor
  // with default parameters
  fInterval = 0;
  fNdimPlot = 0;
  fParamsPlot = 0;
  fColor = 0;
  fFillStyle = 4050; // half transparent
  fLineColor = 0;
  fMaximum = 2.;
  fNPoints = 0;  // default depends if 1D or 2D 
  // default is variable range
  fXmin = 0;
  fXmax = -1;
  fYmin = 0;
  fYmax = -1;
  fPrecision = -1; // use default 
}

//_______________________________________________________
LikelihoodIntervalPlot::LikelihoodIntervalPlot(LikelihoodInterval* theInterval)
{
  // LikelihoodIntervalPlot copy constructor
  fInterval = theInterval;
  fParamsPlot = fInterval->GetParameters();
  fNdimPlot = fParamsPlot->getSize();
  fColor = kBlue;
  fLineColor = kGreen;
  fFillStyle = 4050; // half transparent
  fMaximum = 2.;
  fNPoints = 0;  // default depends if 1D or 2D 
  // default is variable range
  fXmin = 0;
  fXmax = -1;
  fYmin = 0;
  fYmax = -1;
  fPrecision = -1; // use default 
}

//_______________________________________________________
LikelihoodIntervalPlot::~LikelihoodIntervalPlot()
{
  // LikelihoodIntervalPlot destructor
}

//_____________________________________________________________________________
void LikelihoodIntervalPlot::SetLikelihoodInterval(LikelihoodInterval* theInterval)
{
  fInterval = theInterval;
  fParamsPlot = fInterval->GetParameters();
  fNdimPlot = fParamsPlot->getSize();

  return;
}

//_____________________________________________________________________________
void LikelihoodIntervalPlot::SetPlotParameters(const RooArgSet *params) 
{
  fNdimPlot = params->getSize();
  fParamsPlot = (RooArgSet*) params->clone((std::string(params->GetName())+"_clone").c_str());

  return;
}


//_____________________________________________________________________________
void LikelihoodIntervalPlot::Draw(const Option_t *options) 
{
   // draw the Likelihood interval or contour plot
   // For 1D problem draw the log profile likelihood function ratio and its interval 
   // The curve is draws in a RooPLot by default (i.e as a RooCurve) 
   // The plotting range (default is the full parameter range) and the precision of the RooCurve
   // can be specified by using SetRange(x1,x2) and SetPrecision(eps). 
   // SetNPoints(npoints) can also be used  (default is npoints=100) 
   // Optionally the function can be drawn as a TF1 (option="tf1") obtained by sampling the npoints
   // For 2D case, a contour is drawn. The number of contour points is controlled by 
   // SetNPoints(npoints) (default is npoints=40)
   // In case of problems finding the contour with Minuit, the option "nominuit" can be used. 
   // In this case the profile likelihood function is sampled in the npoints x npoints values and then 
   // an approximate contour is obtained. 

   if(fNdimPlot > 2){
      std::cout << "LikelihoodIntervalPlot::Draw(" << GetName() 
                << ") ERROR: contours for more than 2 dimensions not implemented!" << std::endl;
      return;
   }
   
   TIter it = fParamsPlot->createIterator();
   RooRealVar *myparam = (RooRealVar*)it.Next();
      
   RooAbsReal* newProfile = fInterval->GetLikelihoodRatio(); 

   // analyze options 
   TString opt = options; 
   opt.ToLower(); 
   // use RooPLot for drawing the 1D PL
   // if option is TF1 use TF1 for drawing
   bool useRooPlot = opt.Contains("rooplot") ||  ! (opt.Contains("tf1"));
   opt.ReplaceAll("rooplot","");
   opt.ReplaceAll("tf1","");
   // use Minuit for drawing the contours of the PL 
   bool useMinuit = !opt.Contains("nominuit");
   opt.ReplaceAll("nominuit","");

   RooPlot * frame = 0; 

   TString title = GetTitle(); 
   int nPoints = fNPoints; 
   
   if(fNdimPlot == 1){

      if (title.Length() == 0) 
         title = "- log profile likelihood ratio";

      if (nPoints <=0) nPoints = 100; // default in 1D

      const Double_t xcont_min = fInterval->LowerLimit(*myparam);
      const Double_t xcont_max = fInterval->UpperLimit(*myparam);

      RooRealVar* myarg = (RooRealVar *) newProfile->getVariables()->find(myparam->GetName());
      double x1 = myarg->getMin(); 
      double x2 = myarg->getMax(); 


      // use TF1 for drawing the function
      if (!useRooPlot) { 

         // set a first estimate of range including 2 times upper and lower limit
         double xmin = std::max( x1, 2*xcont_min - xcont_max); 
         double xmax = std::min( x2, 2*xcont_max - xcont_min); 
         if (fXmin < fXmax) { xmin = fXmin; xmax = fXmax; }
         
         TF1 * tmp = newProfile->asTF(*myarg); 
         tmp->SetRange(xmin, xmax);      
         tmp->SetNpx(nPoints);

         // clone the function to avoid later to sample it
         TF1 * f1 = (TF1*) tmp->Clone(); 
         delete tmp;
         
         f1->SetTitle(title);
         TString name = TString(GetName()) + TString("_PLL_") + TString(myarg->GetName());
         f1->SetName(name);
         
         // set range for displaying x values where function <=  fMaximum
         // if no range is set amd 
         // if no reasanable value found mantain first estimate
         x1 = xmin; x2 = xmax;  
         if (fMaximum > 0 && fXmin >= fXmax ) { 
            double x0 = f1->GetX(0, xmin, xmax);
            // check that minimum is between xmin and xmax
            if ( x0 > x1 && x0 < x2) { 
               x1 = f1->GetX(fMaximum, xmin, x0); 
               x2 = f1->GetX(fMaximum, x0, xmax); 
               f1->SetMaximum(fMaximum);
            //std::cout << "setting range to " << x1 << " , " << x2 << " x0 = " << x0 << std::endl;
            }
         }
         
         f1->SetRange(x1,x2);
         
         
         f1->SetLineColor(kBlue);
         f1->GetXaxis()->SetTitle(myarg->GetName());
         f1->GetYaxis()->SetTitle("- log #lambda");
         f1->Draw(opt);

      } 
      else { 
         // use a RooPlot for drawing the PL function        
         double xmin = myparam->getMin(); double xmax =  myparam->getMax();
         if (fXmin < fXmax) { xmin = fXmin; xmax = fXmax; }  

         // want to set range on frame not function
         frame = myarg->frame(xmin,xmax,nPoints);
         frame->SetTitle(title);
         frame->GetYaxis()->SetTitle("- log #lambda");
         //    frame->GetYaxis()->SetTitle("- log profile likelihood ratio");
         
         // plot 
         RooCmdArg cmd; 
         if (fPrecision > 0) cmd = RooFit::Precision(fPrecision); 
         newProfile->plotOn(frame,cmd); 
         
         frame->SetMaximum(fMaximum);
         frame->SetMinimum(0.);
      }

      
      myarg->setVal(xcont_max);
      const Double_t Yat_Xmax = newProfile->getVal();
         
      TLine *Yline_cutoff = new TLine(x1,Yat_Xmax,x2,Yat_Xmax);
      TLine *Yline_min = new TLine(xcont_min,0.,xcont_min,Yat_Xmax);
      TLine *Yline_max = new TLine(xcont_max,0.,xcont_max,Yat_Xmax);
      
      Yline_cutoff->SetLineColor(fLineColor);
      Yline_min->SetLineColor(fLineColor);
      Yline_max->SetLineColor(fLineColor);
      
      if (!useRooPlot) { 
         // need to draw the line 
         Yline_cutoff->Draw();
         Yline_min->Draw();
         Yline_max->Draw();
      } 
      else { 
         // add line in the RooPlot
         frame->addObject(Yline_min);
         frame->addObject(Yline_max);
         frame->addObject(Yline_cutoff);
         frame->Draw(opt);
      }
      

      return;
   }
   else if(fNdimPlot == 2){

      RooRealVar *myparamY = (RooRealVar*)it.Next();

      Double_t cont_level = TMath::ChisquareQuantile(fInterval->ConfidenceLevel(),fNdimPlot); // level for -2log LR
      cont_level = cont_level/2; // since we are plotting -log LR

      RooArgList params(*newProfile->getVariables());
      // set values and error for the POI to the best fit values 
      for (int i = 0; i < params.getSize(); ++i) { 
         RooRealVar & par =  (RooRealVar &) params[i];
         RooRealVar * fitPar =  (RooRealVar *) (fInterval->GetBestFitParameters()->find(par.GetName() ) );
         if (fitPar) {
            par.setVal( fitPar->getVal() );
         }
      }
      // do a profile evaluation to start from the best fit values of parameters 
      newProfile->getVal(); 

      if (title.Length() == 0)
         title = TString("Contour of ") + TString(myparamY->GetName() ) + TString(" vs ") + TString(myparam->GetName() ); 

      if (nPoints <=0) nPoints = 40; // default in 2D

      if (!useMinuit) { 
      
         // draw directly the TH2 from the profile LL
         TH2F* hist2D = (TH2F*)newProfile->createHistogram("_hist2D",*myparam,RooFit::YVar(*myparamY),RooFit::Binning(nPoints),RooFit::Scaling(kFALSE));


         hist2D->SetTitle(title);
         hist2D->SetStats(kFALSE);

         hist2D->SetContour(1,&cont_level);

         hist2D->SetFillColor(fColor); 
         hist2D->SetFillStyle(fFillStyle); 
         hist2D->SetLineColor(fLineColor);

         TString tmpOpt(options);

         if(!tmpOpt.Contains("CONT")) tmpOpt.Append("CONT");
         if(!tmpOpt.Contains("LIST")) tmpOpt.Append("LIST"); // if you want the contour TGraphs

         hist2D->Draw(tmpOpt.Data());
         //    hist2D->Draw("cont2,list,same");

         gPad->Update();  // needed for get list of specials 

         // get TGraphs and add them
         //    gROOT->GetListOfSpecials()->Print();
         TObjArray *contours = (TObjArray*) gROOT->GetListOfSpecials()->FindObject("contours"); 
         if(contours){
            TList *list = (TList*)contours->At(0); 
            TGraph *gr1 = (TGraph*)list->First();
            gr1->SetLineColor(kBlack);
            gr1->SetLineStyle(kDashed);
            gr1->Draw("same");
         } else{
            std::cout << "no countours found in ListOfSpecials" << std::endl;
         }
      }
      else { 

         // find contours  using Minuit       
         TGraph * gr = new TGraph(nPoints+1); 
         
         int ncp = fInterval->GetContourPoints(*myparam, *myparamY, gr->GetX(), gr->GetY(),nPoints); 

         if (int(ncp) < nPoints) {
            std::cout << "Warning - Less points calculated in contours np = " << ncp << " / " << nPoints << std::endl;
            for (int i = ncp; i < nPoints; ++i) gr->RemovePoint(i);
         }
         // add last point to same as first one to close the contour
         gr->SetPoint(ncp, gr->GetX()[0], gr->GetY()[0] );
         opt.Append("LF");
         // draw first a dummy 2d histogram gfor the axis 
         if (!opt.Contains("same")) { 

            double xmin = myparam->getMin(); double xmax =  myparam->getMax();
            double ymin = myparamY->getMin(); double ymax =  myparamY->getMax();
            if (fXmin < fXmax) { xmin = fXmin; xmax = fXmax; }  
            if (fYmin < fYmax) { ymin = fYmin; ymax = fYmax; }  

            TH2F* hist2D = new TH2F("_hist2D",title, nPoints, xmin, xmax, nPoints, ymin, ymax );
            hist2D->GetXaxis()->SetTitle(myparamY->GetName());
            hist2D->GetYaxis()->SetTitle(myparam->GetName());
            hist2D->SetBit(TH1::kNoStats); // do not draw statistics
            hist2D->SetFillStyle(fFillStyle); 
            hist2D->SetMaximum(1);  // to avoid problem with subsequents draws
            hist2D->Draw("AXIS");
         }
         gr->SetFillColor(fColor); 
         //gr->SetFillStyle(fFillStyle); 
         gr->SetLineColor(kBlack); 
         if (opt.Contains("same"))  gr->SetFillStyle(fFillStyle); // put transparent
         gr->Draw(opt);
         TString name = TString("Graph_of_") + TString(fInterval->GetName());
         gr->SetName(name);
         
      }

   }

   return;
}
 LikelihoodIntervalPlot.cxx:1
 LikelihoodIntervalPlot.cxx:2
 LikelihoodIntervalPlot.cxx:3
 LikelihoodIntervalPlot.cxx:4
 LikelihoodIntervalPlot.cxx:5
 LikelihoodIntervalPlot.cxx:6
 LikelihoodIntervalPlot.cxx:7
 LikelihoodIntervalPlot.cxx:8
 LikelihoodIntervalPlot.cxx:9
 LikelihoodIntervalPlot.cxx:10
 LikelihoodIntervalPlot.cxx:11
 LikelihoodIntervalPlot.cxx:12
 LikelihoodIntervalPlot.cxx:13
 LikelihoodIntervalPlot.cxx:14
 LikelihoodIntervalPlot.cxx:15
 LikelihoodIntervalPlot.cxx:16
 LikelihoodIntervalPlot.cxx:17
 LikelihoodIntervalPlot.cxx:18
 LikelihoodIntervalPlot.cxx:19
 LikelihoodIntervalPlot.cxx:20
 LikelihoodIntervalPlot.cxx:21
 LikelihoodIntervalPlot.cxx:22
 LikelihoodIntervalPlot.cxx:23
 LikelihoodIntervalPlot.cxx:24
 LikelihoodIntervalPlot.cxx:25
 LikelihoodIntervalPlot.cxx:26
 LikelihoodIntervalPlot.cxx:27
 LikelihoodIntervalPlot.cxx:28
 LikelihoodIntervalPlot.cxx:29
 LikelihoodIntervalPlot.cxx:30
 LikelihoodIntervalPlot.cxx:31
 LikelihoodIntervalPlot.cxx:32
 LikelihoodIntervalPlot.cxx:33
 LikelihoodIntervalPlot.cxx:34
 LikelihoodIntervalPlot.cxx:35
 LikelihoodIntervalPlot.cxx:36
 LikelihoodIntervalPlot.cxx:37
 LikelihoodIntervalPlot.cxx:38
 LikelihoodIntervalPlot.cxx:39
 LikelihoodIntervalPlot.cxx:40
 LikelihoodIntervalPlot.cxx:41
 LikelihoodIntervalPlot.cxx:42
 LikelihoodIntervalPlot.cxx:43
 LikelihoodIntervalPlot.cxx:44
 LikelihoodIntervalPlot.cxx:45
 LikelihoodIntervalPlot.cxx:46
 LikelihoodIntervalPlot.cxx:47
 LikelihoodIntervalPlot.cxx:48
 LikelihoodIntervalPlot.cxx:49
 LikelihoodIntervalPlot.cxx:50
 LikelihoodIntervalPlot.cxx:51
 LikelihoodIntervalPlot.cxx:52
 LikelihoodIntervalPlot.cxx:53
 LikelihoodIntervalPlot.cxx:54
 LikelihoodIntervalPlot.cxx:55
 LikelihoodIntervalPlot.cxx:56
 LikelihoodIntervalPlot.cxx:57
 LikelihoodIntervalPlot.cxx:58
 LikelihoodIntervalPlot.cxx:59
 LikelihoodIntervalPlot.cxx:60
 LikelihoodIntervalPlot.cxx:61
 LikelihoodIntervalPlot.cxx:62
 LikelihoodIntervalPlot.cxx:63
 LikelihoodIntervalPlot.cxx:64
 LikelihoodIntervalPlot.cxx:65
 LikelihoodIntervalPlot.cxx:66
 LikelihoodIntervalPlot.cxx:67
 LikelihoodIntervalPlot.cxx:68
 LikelihoodIntervalPlot.cxx:69
 LikelihoodIntervalPlot.cxx:70
 LikelihoodIntervalPlot.cxx:71
 LikelihoodIntervalPlot.cxx:72
 LikelihoodIntervalPlot.cxx:73
 LikelihoodIntervalPlot.cxx:74
 LikelihoodIntervalPlot.cxx:75
 LikelihoodIntervalPlot.cxx:76
 LikelihoodIntervalPlot.cxx:77
 LikelihoodIntervalPlot.cxx:78
 LikelihoodIntervalPlot.cxx:79
 LikelihoodIntervalPlot.cxx:80
 LikelihoodIntervalPlot.cxx:81
 LikelihoodIntervalPlot.cxx:82
 LikelihoodIntervalPlot.cxx:83
 LikelihoodIntervalPlot.cxx:84
 LikelihoodIntervalPlot.cxx:85
 LikelihoodIntervalPlot.cxx:86
 LikelihoodIntervalPlot.cxx:87
 LikelihoodIntervalPlot.cxx:88
 LikelihoodIntervalPlot.cxx:89
 LikelihoodIntervalPlot.cxx:90
 LikelihoodIntervalPlot.cxx:91
 LikelihoodIntervalPlot.cxx:92
 LikelihoodIntervalPlot.cxx:93
 LikelihoodIntervalPlot.cxx:94
 LikelihoodIntervalPlot.cxx:95
 LikelihoodIntervalPlot.cxx:96
 LikelihoodIntervalPlot.cxx:97
 LikelihoodIntervalPlot.cxx:98
 LikelihoodIntervalPlot.cxx:99
 LikelihoodIntervalPlot.cxx:100
 LikelihoodIntervalPlot.cxx:101
 LikelihoodIntervalPlot.cxx:102
 LikelihoodIntervalPlot.cxx:103
 LikelihoodIntervalPlot.cxx:104
 LikelihoodIntervalPlot.cxx:105
 LikelihoodIntervalPlot.cxx:106
 LikelihoodIntervalPlot.cxx:107
 LikelihoodIntervalPlot.cxx:108
 LikelihoodIntervalPlot.cxx:109
 LikelihoodIntervalPlot.cxx:110
 LikelihoodIntervalPlot.cxx:111
 LikelihoodIntervalPlot.cxx:112
 LikelihoodIntervalPlot.cxx:113
 LikelihoodIntervalPlot.cxx:114
 LikelihoodIntervalPlot.cxx:115
 LikelihoodIntervalPlot.cxx:116
 LikelihoodIntervalPlot.cxx:117
 LikelihoodIntervalPlot.cxx:118
 LikelihoodIntervalPlot.cxx:119
 LikelihoodIntervalPlot.cxx:120
 LikelihoodIntervalPlot.cxx:121
 LikelihoodIntervalPlot.cxx:122
 LikelihoodIntervalPlot.cxx:123
 LikelihoodIntervalPlot.cxx:124
 LikelihoodIntervalPlot.cxx:125
 LikelihoodIntervalPlot.cxx:126
 LikelihoodIntervalPlot.cxx:127
 LikelihoodIntervalPlot.cxx:128
 LikelihoodIntervalPlot.cxx:129
 LikelihoodIntervalPlot.cxx:130
 LikelihoodIntervalPlot.cxx:131
 LikelihoodIntervalPlot.cxx:132
 LikelihoodIntervalPlot.cxx:133
 LikelihoodIntervalPlot.cxx:134
 LikelihoodIntervalPlot.cxx:135
 LikelihoodIntervalPlot.cxx:136
 LikelihoodIntervalPlot.cxx:137
 LikelihoodIntervalPlot.cxx:138
 LikelihoodIntervalPlot.cxx:139
 LikelihoodIntervalPlot.cxx:140
 LikelihoodIntervalPlot.cxx:141
 LikelihoodIntervalPlot.cxx:142
 LikelihoodIntervalPlot.cxx:143
 LikelihoodIntervalPlot.cxx:144
 LikelihoodIntervalPlot.cxx:145
 LikelihoodIntervalPlot.cxx:146
 LikelihoodIntervalPlot.cxx:147
 LikelihoodIntervalPlot.cxx:148
 LikelihoodIntervalPlot.cxx:149
 LikelihoodIntervalPlot.cxx:150
 LikelihoodIntervalPlot.cxx:151
 LikelihoodIntervalPlot.cxx:152
 LikelihoodIntervalPlot.cxx:153
 LikelihoodIntervalPlot.cxx:154
 LikelihoodIntervalPlot.cxx:155
 LikelihoodIntervalPlot.cxx:156
 LikelihoodIntervalPlot.cxx:157
 LikelihoodIntervalPlot.cxx:158
 LikelihoodIntervalPlot.cxx:159
 LikelihoodIntervalPlot.cxx:160
 LikelihoodIntervalPlot.cxx:161
 LikelihoodIntervalPlot.cxx:162
 LikelihoodIntervalPlot.cxx:163
 LikelihoodIntervalPlot.cxx:164
 LikelihoodIntervalPlot.cxx:165
 LikelihoodIntervalPlot.cxx:166
 LikelihoodIntervalPlot.cxx:167
 LikelihoodIntervalPlot.cxx:168
 LikelihoodIntervalPlot.cxx:169
 LikelihoodIntervalPlot.cxx:170
 LikelihoodIntervalPlot.cxx:171
 LikelihoodIntervalPlot.cxx:172
 LikelihoodIntervalPlot.cxx:173
 LikelihoodIntervalPlot.cxx:174
 LikelihoodIntervalPlot.cxx:175
 LikelihoodIntervalPlot.cxx:176
 LikelihoodIntervalPlot.cxx:177
 LikelihoodIntervalPlot.cxx:178
 LikelihoodIntervalPlot.cxx:179
 LikelihoodIntervalPlot.cxx:180
 LikelihoodIntervalPlot.cxx:181
 LikelihoodIntervalPlot.cxx:182
 LikelihoodIntervalPlot.cxx:183
 LikelihoodIntervalPlot.cxx:184
 LikelihoodIntervalPlot.cxx:185
 LikelihoodIntervalPlot.cxx:186
 LikelihoodIntervalPlot.cxx:187
 LikelihoodIntervalPlot.cxx:188
 LikelihoodIntervalPlot.cxx:189
 LikelihoodIntervalPlot.cxx:190
 LikelihoodIntervalPlot.cxx:191
 LikelihoodIntervalPlot.cxx:192
 LikelihoodIntervalPlot.cxx:193
 LikelihoodIntervalPlot.cxx:194
 LikelihoodIntervalPlot.cxx:195
 LikelihoodIntervalPlot.cxx:196
 LikelihoodIntervalPlot.cxx:197
 LikelihoodIntervalPlot.cxx:198
 LikelihoodIntervalPlot.cxx:199
 LikelihoodIntervalPlot.cxx:200
 LikelihoodIntervalPlot.cxx:201
 LikelihoodIntervalPlot.cxx:202
 LikelihoodIntervalPlot.cxx:203
 LikelihoodIntervalPlot.cxx:204
 LikelihoodIntervalPlot.cxx:205
 LikelihoodIntervalPlot.cxx:206
 LikelihoodIntervalPlot.cxx:207
 LikelihoodIntervalPlot.cxx:208
 LikelihoodIntervalPlot.cxx:209
 LikelihoodIntervalPlot.cxx:210
 LikelihoodIntervalPlot.cxx:211
 LikelihoodIntervalPlot.cxx:212
 LikelihoodIntervalPlot.cxx:213
 LikelihoodIntervalPlot.cxx:214
 LikelihoodIntervalPlot.cxx:215
 LikelihoodIntervalPlot.cxx:216
 LikelihoodIntervalPlot.cxx:217
 LikelihoodIntervalPlot.cxx:218
 LikelihoodIntervalPlot.cxx:219
 LikelihoodIntervalPlot.cxx:220
 LikelihoodIntervalPlot.cxx:221
 LikelihoodIntervalPlot.cxx:222
 LikelihoodIntervalPlot.cxx:223
 LikelihoodIntervalPlot.cxx:224
 LikelihoodIntervalPlot.cxx:225
 LikelihoodIntervalPlot.cxx:226
 LikelihoodIntervalPlot.cxx:227
 LikelihoodIntervalPlot.cxx:228
 LikelihoodIntervalPlot.cxx:229
 LikelihoodIntervalPlot.cxx:230
 LikelihoodIntervalPlot.cxx:231
 LikelihoodIntervalPlot.cxx:232
 LikelihoodIntervalPlot.cxx:233
 LikelihoodIntervalPlot.cxx:234
 LikelihoodIntervalPlot.cxx:235
 LikelihoodIntervalPlot.cxx:236
 LikelihoodIntervalPlot.cxx:237
 LikelihoodIntervalPlot.cxx:238
 LikelihoodIntervalPlot.cxx:239
 LikelihoodIntervalPlot.cxx:240
 LikelihoodIntervalPlot.cxx:241
 LikelihoodIntervalPlot.cxx:242
 LikelihoodIntervalPlot.cxx:243
 LikelihoodIntervalPlot.cxx:244
 LikelihoodIntervalPlot.cxx:245
 LikelihoodIntervalPlot.cxx:246
 LikelihoodIntervalPlot.cxx:247
 LikelihoodIntervalPlot.cxx:248
 LikelihoodIntervalPlot.cxx:249
 LikelihoodIntervalPlot.cxx:250
 LikelihoodIntervalPlot.cxx:251
 LikelihoodIntervalPlot.cxx:252
 LikelihoodIntervalPlot.cxx:253
 LikelihoodIntervalPlot.cxx:254
 LikelihoodIntervalPlot.cxx:255
 LikelihoodIntervalPlot.cxx:256
 LikelihoodIntervalPlot.cxx:257
 LikelihoodIntervalPlot.cxx:258
 LikelihoodIntervalPlot.cxx:259
 LikelihoodIntervalPlot.cxx:260
 LikelihoodIntervalPlot.cxx:261
 LikelihoodIntervalPlot.cxx:262
 LikelihoodIntervalPlot.cxx:263
 LikelihoodIntervalPlot.cxx:264
 LikelihoodIntervalPlot.cxx:265
 LikelihoodIntervalPlot.cxx:266
 LikelihoodIntervalPlot.cxx:267
 LikelihoodIntervalPlot.cxx:268
 LikelihoodIntervalPlot.cxx:269
 LikelihoodIntervalPlot.cxx:270
 LikelihoodIntervalPlot.cxx:271
 LikelihoodIntervalPlot.cxx:272
 LikelihoodIntervalPlot.cxx:273
 LikelihoodIntervalPlot.cxx:274
 LikelihoodIntervalPlot.cxx:275
 LikelihoodIntervalPlot.cxx:276
 LikelihoodIntervalPlot.cxx:277
 LikelihoodIntervalPlot.cxx:278
 LikelihoodIntervalPlot.cxx:279
 LikelihoodIntervalPlot.cxx:280
 LikelihoodIntervalPlot.cxx:281
 LikelihoodIntervalPlot.cxx:282
 LikelihoodIntervalPlot.cxx:283
 LikelihoodIntervalPlot.cxx:284
 LikelihoodIntervalPlot.cxx:285
 LikelihoodIntervalPlot.cxx:286
 LikelihoodIntervalPlot.cxx:287
 LikelihoodIntervalPlot.cxx:288
 LikelihoodIntervalPlot.cxx:289
 LikelihoodIntervalPlot.cxx:290
 LikelihoodIntervalPlot.cxx:291
 LikelihoodIntervalPlot.cxx:292
 LikelihoodIntervalPlot.cxx:293
 LikelihoodIntervalPlot.cxx:294
 LikelihoodIntervalPlot.cxx:295
 LikelihoodIntervalPlot.cxx:296
 LikelihoodIntervalPlot.cxx:297
 LikelihoodIntervalPlot.cxx:298
 LikelihoodIntervalPlot.cxx:299
 LikelihoodIntervalPlot.cxx:300
 LikelihoodIntervalPlot.cxx:301
 LikelihoodIntervalPlot.cxx:302
 LikelihoodIntervalPlot.cxx:303
 LikelihoodIntervalPlot.cxx:304
 LikelihoodIntervalPlot.cxx:305
 LikelihoodIntervalPlot.cxx:306
 LikelihoodIntervalPlot.cxx:307
 LikelihoodIntervalPlot.cxx:308
 LikelihoodIntervalPlot.cxx:309
 LikelihoodIntervalPlot.cxx:310
 LikelihoodIntervalPlot.cxx:311
 LikelihoodIntervalPlot.cxx:312
 LikelihoodIntervalPlot.cxx:313
 LikelihoodIntervalPlot.cxx:314
 LikelihoodIntervalPlot.cxx:315
 LikelihoodIntervalPlot.cxx:316
 LikelihoodIntervalPlot.cxx:317
 LikelihoodIntervalPlot.cxx:318
 LikelihoodIntervalPlot.cxx:319
 LikelihoodIntervalPlot.cxx:320
 LikelihoodIntervalPlot.cxx:321
 LikelihoodIntervalPlot.cxx:322
 LikelihoodIntervalPlot.cxx:323
 LikelihoodIntervalPlot.cxx:324
 LikelihoodIntervalPlot.cxx:325
 LikelihoodIntervalPlot.cxx:326
 LikelihoodIntervalPlot.cxx:327
 LikelihoodIntervalPlot.cxx:328
 LikelihoodIntervalPlot.cxx:329
 LikelihoodIntervalPlot.cxx:330
 LikelihoodIntervalPlot.cxx:331
 LikelihoodIntervalPlot.cxx:332
 LikelihoodIntervalPlot.cxx:333
 LikelihoodIntervalPlot.cxx:334
 LikelihoodIntervalPlot.cxx:335
 LikelihoodIntervalPlot.cxx:336
 LikelihoodIntervalPlot.cxx:337
 LikelihoodIntervalPlot.cxx:338
 LikelihoodIntervalPlot.cxx:339
 LikelihoodIntervalPlot.cxx:340
 LikelihoodIntervalPlot.cxx:341
 LikelihoodIntervalPlot.cxx:342
 LikelihoodIntervalPlot.cxx:343
 LikelihoodIntervalPlot.cxx:344
 LikelihoodIntervalPlot.cxx:345
 LikelihoodIntervalPlot.cxx:346
 LikelihoodIntervalPlot.cxx:347
 LikelihoodIntervalPlot.cxx:348
 LikelihoodIntervalPlot.cxx:349
 LikelihoodIntervalPlot.cxx:350
 LikelihoodIntervalPlot.cxx:351
 LikelihoodIntervalPlot.cxx:352
 LikelihoodIntervalPlot.cxx:353
 LikelihoodIntervalPlot.cxx:354
 LikelihoodIntervalPlot.cxx:355
 LikelihoodIntervalPlot.cxx:356
 LikelihoodIntervalPlot.cxx:357
 LikelihoodIntervalPlot.cxx:358
 LikelihoodIntervalPlot.cxx:359
 LikelihoodIntervalPlot.cxx:360
 LikelihoodIntervalPlot.cxx:361
 LikelihoodIntervalPlot.cxx:362
 LikelihoodIntervalPlot.cxx:363
 LikelihoodIntervalPlot.cxx:364
 LikelihoodIntervalPlot.cxx:365
 LikelihoodIntervalPlot.cxx:366
 LikelihoodIntervalPlot.cxx:367
 LikelihoodIntervalPlot.cxx:368
 LikelihoodIntervalPlot.cxx:369
 LikelihoodIntervalPlot.cxx:370
 LikelihoodIntervalPlot.cxx:371
 LikelihoodIntervalPlot.cxx:372