rarSimPdf.cc

Go to the documentation of this file.

/*****************************************************************************
* Project: BaBar detector at the SLAC PEP-II B-factory
* Package: RooRarFit
 *    File: $Id: rarSimPdf.cc,v 1.17 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 Simultaneous Pdf class for RooRarFit
//
// BEGIN_HTML
// This class provides Simultaneous Pdf class for RooRarFit
// END_HTML
//

#include "RooRarFit/rarVersion.hh"

#include "Riostream.h"
//#include <sstream>

#include "RooFitCore/RooAbsPdf.hh"
#include "RooFitCore/RooArgList.hh"
#include "RooFitCore/RooDataSet.hh"
#include "RooFitCore/RooFitResult.hh"
#include "RooFitCore/RooProdPdf.hh"
#include "RooFitCore/RooRealVar.hh"
#include "RooFitCore/RooSimPdfBuilder.hh"
#include "RooFitCore/RooStringVar.hh"

#include "RooRarFit/rarMLFitter.hh"
#include "RooRarFit/rarSimPdf.hh"

ClassImp(rarSimPdf)
  ;

rarSimPdf::rarSimPdf()
  : rarCompBase(),
    _simBuilder(0), _simConfig(0)
{
  init();
}

rarSimPdf::rarSimPdf(const char *configFile, const char *configSec,
                     const char *configStr,
                     rarDatasets *theDatasets, RooDataSet *theData,
                     const char *name, const char *title)
  : rarCompBase(configFile, configSec, configStr,
                theDatasets, theData, name, title, kFALSE),
    _simBuilder(0), _simConfig(0)
{
  init();
}

rarSimPdf::~rarSimPdf()
{
}

void rarSimPdf::init()
{
  // it should always have fitData
  if (!_theData) {
    cout<<" No dataset for Sim Pdf"<<endl;
    exit(-1);
  }
  // create SimPdfBuilder
  _simBuilder=new RooSimPdfBuilder(_subPdfs);
  _simConfig=_simBuilder->createProtoBuildConfig();
  _simConfig->readFromFile(_configFile, 0, getVarSec());
  cout<<"simPdfBuilder configs in the config file:"<<endl;
  _simConfig->Print("v");
  // do we need to use _simBuilder->addSpecializations?
  // Yes, let's do it
  if (getFitter()->getSimBuilder()) {
    RooSimPdfBuilder *theFitterBuilder=getFitter()->getSimBuilder();
    _simBuilder->addSpecializations(theFitterBuilder->splitLeafList());
  }
  // build SimPdf
  _thePdf=(RooAbsPdf*)_simBuilder->buildPdf(*_simConfig, _theData, 0, kTRUE);
  _thePdf->SetName(Form("the_%s", GetName()));
  // fit the total pdf for pdfFit as well
  setControlBits("UseAlsoBasePdfFit");
  
  //_thePdf->Print();
  //_thePdf->Print("v");
  //_thePdf->getComponents()->Print("v");
  cout<<"done init of rarSimPdf for "<<GetName()<<endl<<endl;
}

RooAbsPdf *rarSimPdf::getPdfWOvar(RooArgList ignoredObs)
{
  RooAbsPdf *thePdf=_thePdf;
  cout<<"This pdf can not be used in model"<<endl;
  exit(-1);
  ignoredObs.getHashTableSize() ; // avoid "unused" warning
  return thePdf;
}

RooAbsPdf *rarSimPdf::getProtGen()
{
  cout<<"This pdf can not be used in model"<<endl;
  exit(-1);
  return _theProtGen;
}

---