rarArgusBG.cc

Go to the documentation of this file.

/*****************************************************************************
* Project: BaBar detector at the SLAC PEP-II B-factory
* Package: RooRarFit
 *    File: $Id: rarArgusBG.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 ArgusBG Pdf class for RooRarFit
//
// BEGIN_HTML
// This class provides ArgusBG 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/RooArgusBG.hh"

#include "RooRarFit/rarArgusBG.hh"

ClassImp(rarArgusBG)
  ;

rarArgusBG::rarArgusBG()
  : rarBasePdf(),
    _x(0), _max(0), _c(0), _pow(0)
{
  init();
}

rarArgusBG::rarArgusBG(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), _max(0), _c(0), _pow(0)
{
  init();
}

rarArgusBG::~rarArgusBG()
{
}

void rarArgusBG::init()
{
  cout<<"init of rarArgusBG 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
  _max=createAbsReal("max", "E_{end}", x->getMax(), _x->getUnit());
  _c=createAbsReal("c", "#xi", -23, -80, -1);
  _pow=createAbsReal("pow", "n", 0.5);
  _params.Print("v");
  
  // create pdf
  _thePdf=new RooArgusBG(Form("the_%s",GetName()), _pdfType+" "+GetTitle(),
                         *_x, *_max, *_c, *_pow);
}

---