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
41
42////////////////////////////////////////////////////////////////////////////////
43/// Constructor of plain RooAddPdf fraction from list of recursive fractions
44
45RooRecursiveFraction::RooRecursiveFraction(const char* name, const char* title, const RooArgList& fracList) :
46 RooAbsReal(name, title),
47 _list("list","First set of components",this)
48{
49 for (Int_t ifrac=fracList.size()-1 ; ifrac>=0 ; ifrac--) {
50 RooAbsArg* comp = fracList.at(ifrac) ;
51 if (!dynamic_cast<RooAbsReal*>(comp)) {
52 std::stringstream errorMsg;
53 errorMsg << "RooRecursiveFraction::ctor(" << GetName() << ") ERROR: component " << comp->GetName()
54 << " is not of type RooAbsReal" << std::endl;
55 coutE(InputArguments) << errorMsg.str();
56 throw std::invalid_argument(errorMsg.str());
57 }
58
59 _list.add(*comp) ;
60 }
61}
62
63
64
65////////////////////////////////////////////////////////////////////////////////
66/// Copy constructor
67
69 RooAbsReal(other, name),
70 _list("list",this,other._list)
71{
72
73}
74
75////////////////////////////////////////////////////////////////////////////////
76/// Calculate and return value of \f$ a_n * \prod_{i=0}^{n-1} (1 - a_i) \f$.
78{
79 const RooArgSet* nset = _list.nset() ;
80
81 // Note that input coefficients are saved in reverse in this list.
82 double prod = static_cast<RooAbsReal&>(_list[0]).getVal(nset);
83
84 for (unsigned int i=1; i < _list.size(); ++i) {
85 prod *= (1 - static_cast<RooAbsReal&>(_list[i]).getVal(nset));
86 }
87
88 return prod ;
89}
90
92{
93 ctx.addResult(this, ctx.buildCall("RooFit::Detail::MathFuncs::recursiveFraction", _list, _list.size()));
94}
#define coutE(a)
#define ClassImp(name)
Definition Rtypes.h:382
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:79
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:24
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 class to maintain the context for squashing of RooFit models into code.
std::string buildCall(std::string const &funcname, Args_t const &...args)
Build the code to call the function with name funcname, passing some arguments.
void addResult(RooAbsArg const *key, std::string const &value)
A function to save an expression that includes/depends on the result of the input node.
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
void translate(RooFit::Detail::CodeSquashContext &ctx) const override
This function defines a translation for each RooAbsReal based object that can be used to express the ...
const char * GetName() const override
Returns name of object.
Definition TNamed.h:47