rarProd Class Reference

ProdPdf builder. More...

#include <rarProd.hh>

Inheritance diagram for rarProd:

List of all members.

Public Member Functions
	rarProd ()
	Trivial ctor.
	rarProd (const char *configFile, const char *configSec, const char *configStr, rarDatasets *theDatasets, RooDataSet *theData, const char *name, const char *title)
	Default ctor.
virtual	~rarProd ()
virtual RooAbsPdf *	getPdfWOvar (RooArgList ignoredObs)
	Get a PDF not depending on specified observables for sPlot.
virtual RooAbsPdf *	getDPdfWvar (RooRealVar *theVar)
	Return the pdf if it depends on the var.
virtual RooAbsPdf *	getProtGen ()
	Return prototype var generator for toy study.
Protected Member Functions
void	init ()
	Initial function called by ctor.
Protected Attributes
TList	_condPdfList
	RooRarFit Pdf list.
RooArgList	_condPdfs
	condPdf ArgList
Private Member Functions
	rarProd (const rarProd &)
	ClassDef (rarProd, 0)

---

Detailed Description

ProdPdf builder.

Build composite pdfs through RooProdPdf.

Config Directives:

See doc for ProdPdf configs.

Definition at line 26 of file rarProd.hh.

---

Constructor & Destructor Documentation

rarProd::rarProd

(

)

Trivial ctor.

Usually the objects should be created using other ctors.

Definition at line 40 of file rarProd.cc.

References init().

  : rarCompBase()
{
  init();
}

rarProd::rarProd	(	const char *	configFile,
		const char *	configSec,
		const char *	configStr,
		rarDatasets *	theDatasets,
		RooDataSet *	theData,
		const char *	name,
		const char *	title
	)

Default ctor.

Parameters:

	configFile	The config file
	configSec	The config section
	configStr	The config string
	theDatasets	Available datasets
	theData	Default dataset for this PDF
	name	The name
	title	The title

The default ctor first initializes data members, and then calls init.

Definition at line 58 of file rarProd.cc.

References init().

  : rarCompBase(configFile, configSec, configStr,
                theDatasets, theData, name, title, kFALSE)
{
  init();
}

rarProd::~rarProd

(

)

[virtual]

Definition at line 68 of file rarProd.cc.

{
}

rarProd::rarProd

(

const rarProd &

)

[private]

---

Member Function Documentation

rarProd::ClassDef	(	rarProd	,
		0
	)			[private]

RooAbsPdf * rarProd::getDPdfWvar

(

RooRealVar *

theVar

)

[virtual]

Return the pdf if it depends on the var.

Parameters:

theVar

The var to check

Returns:

This (sim)Pdf if it depends on the var

It checks if the pdf depends on the obs; if yes, it returns this (sim)Pdf

Reimplemented from rarBasePdf.

Definition at line 175 of file rarProd.cc.

References rarCompBase::_nComp, and rarCompBase::_pdfList.

{
  RooAbsPdf *retVal(0);
  for (Int_t i=0; i<_nComp; i++) {
    rarBasePdf *thePdf=(rarBasePdf*)_pdfList.At(i);
    retVal=thePdf->getDPdfWvar(theVar);
    if (retVal) return retVal;
  }
  return retVal;
}

RooAbsPdf * rarProd::getPdfWOvar

(

RooArgList

ignoredObs

)

[virtual]

Get a PDF not depending on specified observables for sPlot.

Parameters:

ignoredObs

The observables to check

Returns:

The Pdf required

It checks of the pdf depends on the observables. if no, it returns the default pdf; if yes, it builds the non-dependent pdf.

Reimplemented from rarBasePdf.

Definition at line 147 of file rarProd.cc.

References rarConfig::_fullObs, rarBasePdf::_thePdf, and rarCompBase::getPdfsWOvar().

{
  RooAbsPdf *thePdf=_thePdf;
  // check if thePdf depends on theVar
  if (!(thePdf->dependsOn(ignoredObs))) return thePdf;
  // create individual Pdf
  TString theName=Form("the_%s_%s", ignoredObs[0].GetName(), GetName());
  TString theTitle=Form("%s w/o %s", GetTitle(), ignoredObs[0].GetName());
  RooArgList myPdfList=getPdfsWOvar(ignoredObs);
  thePdf=new RooProdPdf(theName, theTitle, myPdfList);
  //thePdf->Print("v");
  //thePdf->Print();
  if (!(thePdf->dependsOn(ignoredObs))) return thePdf;
  delete thePdf;
  // create a dummy for it
  thePdf=new RooGenericPdf(theName, theTitle, "1",
                           *_thePdf->getParameters(_fullObs));
  //thePdf->Print("v");
  //thePdf->Print();
  return thePdf;
}

RooAbsPdf * rarProd::getProtGen

(

)

[virtual]

Return prototype var generator for toy study.

Returns:

The generator

It constrcuts prototype var generator for toy study

Reimplemented from rarCompBase.

Definition at line 190 of file rarProd.cc.

References rarBasePdf::_myDummyPdf, rarBasePdf::_protGenPdfs, rarCompBase::_subProtGenPdfs, rarBasePdf::_theProtGen, rarBasePdf::getControlBit(), and rarCompBase::getProtGen().

{
  if (_theProtGen) return _theProtGen;
  // call super class' getProtGen
  rarCompBase::getProtGen();
  _theProtGen=_myDummyPdf;
  RooArgList pdfList(_protGenPdfs);
  if (!getControlBit("dummySubGenPdfs")) pdfList.add(_subProtGenPdfs);
  else pdfList.add(*_theProtGen);
  if ((pdfList.getSize()>1)||(!getControlBit("dummySubGenPdfs")))
    _theProtGen = new RooProdPdf (Form("protGen_%s",GetName()),
                                  Form("Prot Gen for %s", GetName()), pdfList);
  
  return _theProtGen;
}

void rarProd::init

(

)

[protected, virtual]

Initial function called by ctor.

init is called by the ctor. The initialization work to read in config info is done with rarCompBase::init(), and init here just builds the RooProdPdf with configured components.

Reimplemented from rarCompBase.

Definition at line 78 of file rarProd.cc.

References _condPdfList, _condPdfs, rarCompBase::_pdfList, rarBasePdf::_pdfType, rarBasePdf::_subPdfs, rarBasePdf::_thePdf, rarCompBase::getArgSet(), rarStrParser::nArgs(), rarConfig::readConfStr(), and rarBasePdf::setControlBit().

Referenced by rarProd().

{
  // Config pdfs done by parent Pdf
  // list configed pdfs done by parent pdf
  
  // set control bits
  setControlBit("noUseBasePdfFit", "ndFit");
  
  // list of cmd
  RooLinkedList cmdList;

  // do we have conditional obseervables?
  RooArgSet condObsSet(getArgSet(readConfStr("CondObss",""),kFALSE));
  // do we have conditional Pdf ?
  rarStrParser condPdfsStr=readConfStr("CondPdfs","");
  while (condPdfsStr.nArgs()>0) {
    rarBasePdf *condPdf=(rarBasePdf*)_pdfList.FindObject(condPdfsStr[0]);
    condPdfsStr.Remove();
    if (condPdf) {
      _condPdfList.Add(condPdf);
      _condPdfs.add(*condPdf->getPdf());
      // now set CondObss to this pdf
      //RooArgSet *thisPdfSet=new RooArgSet;
      //thisPdfSet->add(*condPdf->getPdf());
      RooArgSet *thisPdfSet=new RooArgSet(*condPdf->getPdf());
      RooArgSet *thisCondObsSet=new RooArgSet(condObsSet);
      thisCondObsSet->
        add(getArgSet(readConfStr
                      (Form("CondObss_%s",condPdf->GetName()),""),kFALSE));
      condPdf->setCondObss(*thisCondObsSet);
      cmdList.Add(Conditional(*thisPdfSet,*thisCondObsSet).Clone());
      cout<<" "<<condPdf->GetName()<<" is conditional PDF"<<endl;
      cout<<" The normalization of this PDF is on "<<endl;
      thisCondObsSet->Print();
    }
  }
  // Pdfs without conditional
  RooArgList subPdfs(_subPdfs);
  subPdfs.remove(_condPdfs);
  if (_condPdfs.getSize()>0) {
    // do not use ndFit by default
    //setControlBit("UseBasePdfFit", "ndFit");
    // cout...
    cout<<"The following PDFs are conditional PDFs:"<<endl;
    _condPdfs.Print();
    cout<<endl;
    //cout<<"The normalization of those PDFs are on observables:"<<endl;
    //condObsSet.Print();
    //cout<<endl;
    cout<<"The remaining PDFs in the production are:"<<endl;
    subPdfs.Print();
    cout<<endl;
  }
  
  // build ProdPdf
  //cmdList.Print("v");
  _thePdf=new RooProdPdf(Form("the_%s", GetName()), _pdfType+" "+GetTitle(),
                         subPdfs, cmdList);
  
  cout<<"done init of rarProd for "<<GetName()<<endl<<endl;
}

---

Member Data Documentation

TList rarProd::_condPdfList [protected]

RooRarFit Pdf list.

Definition at line 42 of file rarProd.hh.

Referenced by init().

RooArgList rarProd::_condPdfs [protected]

condPdf ArgList

Definition at line 43 of file rarProd.hh.

Referenced by init().

---

The documentation for this class was generated from the following files:

• rarProd.hh
• rarProd.cc

---