ROOT   Reference Guide
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 *
11 * *
12 * Redistribution and use in source and binary forms, *
13 * with or without modification, are permitted according to the terms *
15 *****************************************************************************/
16
17/**
18\file RooRecursiveFraction.cxx
19\class RooRecursiveFraction
20\ingroup Roofitcore
21
22Class RooRecursiveFraction is a 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
29#include "RooFit.h"
30
31#include "Riostream.h"
32#include <math.h>
33
35#include "RooAbsReal.h"
36#include "RooAbsPdf.h"
37#include "RooErrorHandler.h"
38#include "RooArgSet.h"
39#include "RooNLLVar.h"
40#include "RooChi2Var.h"
41#include "RooMsgService.h"
42
43using namespace std;
44
46
47
48
49////////////////////////////////////////////////////////////////////////////////
50/// Default constructor
51
53{
54
55}
56
57
58
59////////////////////////////////////////////////////////////////////////////////
60/// Constructor of plain RooAddPdf fraction from list of recursive fractions
61
62RooRecursiveFraction::RooRecursiveFraction(const char* name, const char* title, const RooArgList& fracList) :
63 RooAbsReal(name, title),
64 _list("list","First set of components",this)
65{
66 for (Int_t ifrac=fracList.getSize()-1 ; ifrac>=0 ; ifrac--) {
67 RooAbsArg* comp = fracList.at(ifrac) ;
68 if (!dynamic_cast<RooAbsReal*>(comp)) {
69 std::stringstream errorMsg;
70 errorMsg << "RooRecursiveFraction::ctor(" << GetName() << ") ERROR: component " << comp->GetName()
71 << " is not of type RooAbsReal" << endl ;
72 coutE(InputArguments) << errorMsg.str();
73 throw std::invalid_argument(errorMsg.str());
74 }
75
77 }
78}
79
80
81
82////////////////////////////////////////////////////////////////////////////////
83/// Copy constructor
84
86 RooAbsReal(other, name),
87 _list("list",this,other._list)
88{
89
90}
91
92
93
94////////////////////////////////////////////////////////////////////////////////
95/// Destructor
96
98{
99
100}
101
102
103
104////////////////////////////////////////////////////////////////////////////////
105/// Calculate and return value of \f$a_n * \prod_{i=0}^{n-1} (1 - a_i) \f$.
107{
108 const RooArgSet* nset = _list.nset() ;
109
110 // Note that input coefficients are saved in reverse in this list.
111 Double_t prod = static_cast<RooAbsReal&>(_list[0]).getVal(nset);
112
113 for (unsigned int i=1; i < _list.size(); ++i) {
114 prod *= (1 - static_cast<RooAbsReal&>(_list[i]).getVal(nset));
115 }
116
117 return prod ;
118}
119
#define coutE(a)
Definition: RooMsgService.h:33
#define ClassImp(name)
Definition: Rtypes.h:361
char name[80]
Definition: TGX11.cxx:109
RooAbsArg is the common abstract base class for objects that represent a value (of arbitrary type) an...
Definition: RooAbsArg.h:73
Int_t getSize() const
Storage_t::size_type size() const
const RooArgSet * nset() const
Definition: RooAbsProxy.h:46
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
Definition: RooAbsReal.h:60
Double_t getVal(const RooArgSet *normalisationSet=nullptr) const
Evaluate object.
Definition: RooAbsReal.h:90
RooArgList is a container object that can hold multiple RooAbsArg objects.
Definition: RooArgList.h:21
RooAbsArg * at(Int_t idx) const
Return object at given index, or nullptr if index is out of range.
Definition: RooArgList.h:74
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition: RooArgSet.h:28
virtual Bool_t add(const RooAbsArg &var, Bool_t silent=kFALSE)