rarVoigtian.cc

Go to the documentation of this file.

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

// -- CLASS DESCRIPTION [RooRarFit] --
// This class provides Voigtian Pdf class for RooRarFit
//
// BEGIN_HTML
// This class provides Voigtian 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/RooVoigtian.hh"

#include "RooRarFit/rarVoigtian.hh"

ClassImp(rarVoigtian)
  ;

rarVoigtian::rarVoigtian()
  : rarBasePdf(),
    _x(0), _mean(0), _width(0), _sigma(0)
{
  init();
}

rarVoigtian::rarVoigtian(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), _width(0), _sigma(0)
{
  init();
}

rarVoigtian::~rarVoigtian()
{
}

void rarVoigtian::init()
{
  cout<<"init of rarVoigtian 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());

  _width=createAbsReal("width", "#mu", (x->getMin()+x->getMax())/8,
                      x->getMin(), x->getMax(), _x->getUnit());

  _sigma=createAbsReal("sigma", "#sigma", 0.0, 0., 1., _x->getUnit());

  _params.Print("v");
  
  // create pdf
  Bool_t doFast = kTRUE;
  _thePdf=new RooVoigtian(Form("the_%s", GetName()), _pdfType+" "+GetTitle(),
                          *_x, *_mean, *_width, *_sigma, doFast);
}

---