rarRelBreitWigner.cc

Go to the documentation of this file.

/*****************************************************************************
* Project: BaBar detector at the SLAC PEP-II B-factory
* Package: RooRarFit
 *    File: $Id: rarRelBreitWigner.cc,v 1.4 2011/08/26 17:54:18 fwilson Exp $
 * Authors: Fergus Wilson
 * History:
 * 
 * Copyright (C) 2005-2012, RAL
 *****************************************************************************/

// -- CLASS DESCRIPTION [RooRarFit] --
// This class provides RelBreitWigner Pdf class for RooRarFit
//
// BEGIN_HTML
// This class provides RelBreitWigner 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/RooRelBreitWigner.hh"
#include "RooRarFit/rarRelBreitWigner.hh"

ClassImp(rarRelBreitWigner);

rarRelBreitWigner::rarRelBreitWigner()
  : rarBasePdf(),
    _x(0), _mean(0), _width(0), _radius(0), _mass_a(0), _mass_b(0), _spin(0)
{
  init();
}

rarRelBreitWigner::rarRelBreitWigner(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),  _radius(0), _mass_a(0), _mass_b(0), _spin(0)
{
  init();
}

rarRelBreitWigner::~rarRelBreitWigner()
{
}

void rarRelBreitWigner::init()
{
  // read in _parSpin from config section
  //_parSpin=atoi(readConfStr("spin", Form("%d", _parSpin), getVarSec()));

  cout <<"init of rarRelBreitWigner for "<<GetName()<<endl;
  
  // first get its obs
  _x=createAbsReal("x", "observable"); assert(_x);

  // Config pdf params
  _mean=createAbsReal("mean", "mean", 0, -10, 10);
  _width=createAbsReal("width", "width", 0, 0, 100);
  _radius=createAbsReal("radius", "radius", 3.1, "RooConstVar");
  _mass_a=createAbsReal("mass_a", "mass_a", 0.493677, "RooConstVar");
  _mass_b=createAbsReal("mass_b", "mass_b", 0.139570, "RooConstVar");
  _spin=createAbsReal("spin", "spin", 1, "RooConstVar");

  _params.Print("v");
  
  // create pdf
  _thePdf=new RooRelBreitWigner(Form("the_%s", GetName()),_pdfType+" "+GetTitle(),
                                *_x, *_mean, *_width, *_radius, *_mass_a, *_mass_b, *_spin);
}

---