rarOsipDisc.cc

Go to the documentation of this file.

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

// -- CLASS DESCRIPTION [RooRarFit] --
// This class provides the Osipenkov Discriminator Pdf class for RooRarFit
//
// BEGIN_HTML
// This class provides Osipenkov Discriminator 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 "RooRarFit/RooOsipDisc.hh"
#include "RooRarFit/rarOsipDisc.hh"

ClassImp(rarOsipDisc)
  ;

rarOsipDisc::rarOsipDisc()
  : rarBasePdf(),
    _x(0), _r(0), _b(0), _a(0), _h(0), _s(0)
{
  init();
}

rarOsipDisc::rarOsipDisc(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), _r(0), _b(0), _a(0), _h(0), _s(0)
{
  init();
}

rarOsipDisc::~rarOsipDisc()
{
}

void rarOsipDisc::init()
{
  cout<<"init of rarOsipDisc 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
  
  _r=createAbsReal("r", "#r", .1, 0., 1., _x->getUnit());
  _b=createAbsReal("b", "#b", .1, 0., 1., _x->getUnit());
  _a=createAbsReal("a", "#a", .1, 0., 1., _x->getUnit());
  _h=createAbsReal("h", "#h", 1.0, -10, 10, _x->getUnit());
  _s=createAbsReal("s", "#s", .1, 0., 1., _x->getUnit());

  _params.Print("v");
  
  // create pdf
  if("OsipDisc"==_pdfType) {
    _thePdf=new RooOsipDisc(Form("the_%s", GetName()),
                               _pdfType+" "+GetTitle(), *_x, *_r, *_b, *_a, *_h, *_s);
  }
}

---