rarPoly.cc

Go to the documentation of this file.

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

#include "RooRarFit/rarVersion.hh"

#include "Riostream.h"

#include "RooFitCore/RooArgList.hh"
#include "RooFitCore/RooDataSet.hh"
#include "RooFitCore/RooFormulaVar.hh"
#include "RooFitCore/RooRealVar.hh"
#include "RooFitCore/RooStringVar.hh"

#include "RooFitModels/RooChebychev.hh"
#include "RooFitModels/RooPolynomial.hh"

#include "RooRarFit/rarPoly.hh"

ClassImp(rarPoly)
  ;

rarPoly::rarPoly()
  : rarBasePdf(),
    _x(0), _nOrder(1)
{
  init();
}

rarPoly::rarPoly(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),
    _x(0), _nOrder(1)
{
  init();
}

rarPoly::~rarPoly()
{
}

void rarPoly::init()
{
  cout<<"init of rarPoly for "<<GetName()<<":"<<endl;
  
  // first get its obs
  _x=createAbsReal("x", "observable"); assert(_x);
  // read in pdf params
  _nOrder=atoi(readConfStr("nOrder", "1", getVarSec()));
  if (_nOrder<=0) _nOrder=1;
  _nOrder++;
  for (Int_t i=1; i<_nOrder; i++) {
    RooAbsReal *P=createAbsReal(Form("P%02d",i),Form("p_{%d}", i),
                                0, -10000, +10000);
    _coeffs.add(*P);
  }
  _params.Print("v");
  
  // create pdf
  if ("Chebychev"==_pdfType) {
    _thePdf=new RooChebychev(Form("the_%s", GetName()),_pdfType+" "+GetTitle(),
                             *_x, _coeffs);
  } else { // default
    _thePdf=new RooPolynomial(Form("the_%s",GetName()),_pdfType+" "+GetTitle(),
                              *_x, _coeffs);
  }
}

---