rarGaussian.cc

Go to the documentation of this file.

/*****************************************************************************
* Project: BaBar detector at the SLAC PEP-II B-factory
* Package: RooRarFit
 *    File: $Id: rarGaussian.cc,v 1.7 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 Gaussian/BreitWigner Pdf class for RooRarFit
//
// BEGIN_HTML
// This class provides Gaussian/BreitWigner 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/RooBreitWigner.hh"
#include "RooFitModels/RooGaussian.hh"
#include "RooFitModels/RooLandau.hh"

#include "RooRarFit/rarGaussian.hh"

ClassImp(rarGaussian)
  ;

rarGaussian::rarGaussian()
  : rarBasePdf(),
    _x(0), _mean(0), _sigma(0),
    _scale(0), _shift(0)
{
  init();
}

rarGaussian::rarGaussian(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),
    _scale(0), _shift(0)
{
  init();
}

rarGaussian::~rarGaussian()
{
}

void rarGaussian::init()
{
  cout<<"init of rarGaussian 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());
  if ("notSet"!=readConfStr("shift", "notSet", getVarSec())) { // use shift
    _shift=createAbsReal("shift", "d#mu", 0.);
    _mean=(RooAbsReal*)
      createAbsVar(Form("%s %s", "meanS RooFormulaVar","@0+@1 mean shift"));
  }
  _sigma=createAbsReal("sigma", "#sigma", .1, 0., 1., _x->getUnit());
  if ("notSet"!=readConfStr("scale", "notSet", getVarSec())) { // use scale
    _scale=createAbsReal("scale", "S#sigma_{C}", 1.);
    _sigma=(RooAbsReal*)
      createAbsVar(Form("%s %s","sigmaS RooFormulaVar","@0*@1 sigma scale"));
  }
  _params.Print("v");
  
  // create pdf
  if("BreitWigner"==_pdfType) {
    _thePdf=new RooBreitWigner(Form("the_%s", GetName()),
                               _pdfType+" "+GetTitle(), *_x, *_mean, *_sigma);
  } else if ("Landau"==_pdfType) {
    _thePdf=new RooLandau(Form("the_%s", GetName()),
                          _pdfType+" "+GetTitle(), *_x, *_mean, *_sigma);
  } else { // default
    _thePdf=new RooGaussian(Form("the_%s", GetName()), _pdfType+" "+GetTitle(),
                            *_x, *_mean, *_sigma);
  }
}

---