rarHistPdf.cc

Go to the documentation of this file.

/*****************************************************************************
* Project: BaBar detector at the SLAC PEP-II B-factory
* Package: RooRarFit
 *    File: $Id: rarHistPdf.cc,v 1.4 2011/06/16 13:18:50 fwilson Exp $
 * Authors: Lei Zhang
 * History:
 * 
 * Copyright (C) 2005-2012,, University of California, Riverside
 *****************************************************************************/

// -- CLASS DESCRIPTION [RooRarFit] --
// This class provides HistPdf class for RooRarFit
//
// BEGIN_HTML
// This class provides HistPdf class for RooRarFit
// END_HTML
//

#include "RooRarFit/rarVersion.hh"

#include "Riostream.h"

#include "RooFitCore/RooArgList.hh"
#include "RooFitCore/RooDataSet.hh"
#include "RooFitCore/RooDataHist.hh"
#include "RooFitCore/RooProdPdf.hh"
#include "RooFitCore/RooRealVar.hh"
#include "RooFitCore/RooStringVar.hh"

#include "RooFitCore/RooHistPdf.hh"
#include "RooRarFit/rarHistPdf.hh"

ClassImp(rarHistPdf)
  ;

rarHistPdf::rarHistPdf()
  : rarBasePdf(),
    _theHist(0)
{
  init();
}

rarHistPdf::rarHistPdf(const char *configFile, const char *configSec,
                       const char *configStr,
                       rarDatasets *theDatasets, RooDataSet *theData,
                       const char *name, const char *title)
  : rarBasePdf(configFile, configSec, configStr,
               theDatasets, theData, name, title),
    _theHist(0)
{
  init();
}

rarHistPdf::~rarHistPdf()
{
}

void rarHistPdf::init()
{
  cout<<"init of rarHistPdf for "<<GetName()<<":"<<endl;
  
  // no need for pdf fit
  setControlBits("noPdfFit");
  
  // read the obs string
  rarStrParser obsStrParser=readConfStr("obs", "", getVarSec());
  // get obs list
  getFormulaArgs(obsStrParser);
  
  cout<<" Obs in pdf:"<<endl;
  _obsSet.Print("v");
  
  // it should always have fitData
  if (!_theData) {
    cout<<" No dataset for HistPdf"<<endl;
    exit(-1);
  }
  // create the RooDataHist
  _theHist=new RooDataHist(Form("the_%s_Hist", GetName()),
                           "pdf histogram", _obsSet, *_theData);
  // create the generic pdf
  _thePdf=new RooHistPdf(Form("the_%s", GetName()), _pdfType+" "+GetTitle(),
                         _obsSet, *_theHist);
}

---