Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
RooRecursiveFraction.cxx
Go to the documentation of this file.
1/*****************************************************************************
2 * Project: RooFit *
3 * Package: RooFitCore *
4 * @(#)root/roofitcore:$Id$
5 * Authors: *
6 * WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu *
7 * DK, David Kirkby, UC Irvine, dkirkby@uci.edu *
8 * *
9 * Copyright (c) 2000-2005, Regents of the University of California *
10 * and Stanford University. All rights reserved. *
11 * *
12 * Redistribution and use in source and binary forms, *
13 * with or without modification, are permitted according to the terms *
14 * listed in LICENSE (http://roofit.sourceforge.net/license.txt) *
15 *****************************************************************************/
16
17/**
18\file RooRecursiveFraction.cxx
19\class RooRecursiveFraction
20\ingroup Roofitcore
21
22A RooAbsReal implementation that
23calculates the plain fraction of sum of RooAddPdf components
24from a set of recursive fractions: for a given set of input fractions
25\f$ {a_i} \f$, it returns \f$ a_n * \prod_{i=0}^{n-1} (1 - a_i) \f$.
26**/
27
28#include "Riostream.h"
29#include <cmath>
30
32#include "RooAbsReal.h"
33#include "RooAbsPdf.h"
34#include "RooErrorHandler.h"
35#include "RooArgSet.h"
36#include "RooMsgService.h"
37
39
40////////////////////////////////////////////////////////////////////////////////
41/// Constructor of plain RooAddPdf fraction from list of recursive fractions
42
43RooRecursiveFraction::RooRecursiveFraction(const char* name, const char* title, const RooArgList& fracList) :
44 RooAbsReal(name, title),
45 _list("list","First set of components",this)
46{
47 for (Int_t ifrac=fracList.size()-1 ; ifrac>=0 ; ifrac--) {
48 RooAbsArg* comp = fracList.at(ifrac) ;
49 if (!dynamic_cast<RooAbsReal*>(comp)) {
50 std::stringstream errorMsg;
51 errorMsg << "RooRecursiveFraction::ctor(" << GetName() << ") ERROR: component " << comp->GetName()
52 << " is not of type RooAbsReal" << std::endl;
53 coutE(InputArguments) << errorMsg.str();
54 throw std::invalid_argument(errorMsg.str());
55 }
56
57 _list.add(*comp) ;
58 }
59}
60
61
62
63////////////////////////////////////////////////////////////////////////////////
64/// Copy constructor
65
67 RooAbsReal(other, name),
68 _list("list",this,other._list)
69{
70
71}
72
73////////////////////////////////////////////////////////////////////////////////
74/// Calculate and return value of \f$ a_n * \prod_{i=0}^{n-1} (1 - a_i) \f$.
76{
77 const RooArgSet* nset = _list.nset() ;
78
79 // Note that input coefficients are saved in reverse in this list.
80 double prod = static_cast<RooAbsReal&>(_list[0]).getVal(nset);
81
82 for (unsigned int i=1; i < _list.size(); ++i) {
83 prod *= (1 - static_cast<RooAbsReal&>(_list[i]).getVal(nset));
84 }
85
86 return prod ;
87}
#define coutE(a)
#define ClassImp(name)
Definition Rtypes.h:377
char name[80]
Definition TGX11.cxx:110
Common abstract base class for objects that represent a value and a "shape" in RooFit.
Definition RooAbsArg.h:77
Storage_t::size_type size() const
const RooArgSet * nset() const
Definition RooAbsProxy.h:52
Abstract base class for objects that represent a real value and implements functionality common to al...
Definition RooAbsReal.h:59
double getVal(const RooArgSet *normalisationSet=nullptr) const
Evaluate object.
Definition RooAbsReal.h:103
RooArgList is a container object that can hold multiple RooAbsArg objects.
Definition RooArgList.h:22
RooAbsArg * at(Int_t idx) const
Return object at given index, or nullptr if index is out of range.
Definition RooArgList.h:110
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition RooArgSet.h:55
bool add(const RooAbsArg &var, bool valueServer, bool shapeServer, bool silent)
Overloaded RooCollection_t::add() method insert object into set and registers object as server to own...
A RooAbsReal implementation that calculates the plain fraction of sum of RooAddPdf components from a ...
double evaluate() const override
Calculate and return value of .
RooRecursiveFraction()=default
const char * GetName() const override
Returns name of object.
Definition TNamed.h:47