rarCBShape.cc

Go to the documentation of this file.

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

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

#include "RooRarFit/rarVersion.hh"

#include "Riostream.h"

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

#include "RooFitModels/RooCBShape.hh"

#include "RooRarFit/rarCBShape.hh"

ClassImp(rarCBShape)
  ;

rarCBShape::rarCBShape()
  : rarBasePdf(),
    _x(0), _mean(0), _sigma(0), _alpha(0), _n(0)
{
  init();
}

rarCBShape::rarCBShape(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), _mean(0), _sigma(0), _alpha(0), _n(0)
{
  init();
}

rarCBShape::~rarCBShape()
{
}

void rarCBShape::init()
{
  cout<<"init of rarCBShape for "<<GetName()<<":"<<endl;
  
  // first get its dependent/observable
  _x=createAbsReal("x", "observable"); assert(_x);
  RooRealVar *x=(RooRealVar *)RooArgList(_obsSet).at(0); assert(x);
  // Config pdf params
  _mean=createAbsReal("mean", "#mu", (x->getMin()+x->getMax())/2,
                      x->getMin(), x->getMax(), _x->getUnit());
  _sigma=createAbsReal("sigma", "#sigma", .1, 0., 1., _x->getUnit());
  _alpha=createAbsReal("alpha", "#alpha", .9);
  _n=createAbsReal("n", "n", 10, 0.1, 200.);  
  _params.Print("v");
  
  // create pdf
  _thePdf=new RooCBShape(Form("the_%s", GetName()), _pdfType+" "+GetTitle(),
                         *_x, *_mean, *_sigma, *_alpha, *_n);
}

---