40static const char*
NLL_NAME =
"nll_MarkovChain_local_";
45MarkovChain::MarkovChain() :
120 for(
int i=0; i < otherChain.
Size(); i++ ) {
123 if( counter > burnIn ) {
135 double counter = 0.0;
136 for(
int i=0; i < otherChain.
Size(); i++ ) {
138 counter += otherChain.
Weight();
139 if( counter > discardEntries ) {
158 if (whichVars ==
nullptr) {
163 args.
add(*whichVars);
185 if (whichVars ==
nullptr) {
190 args.
add(*whichVars);
207 hist =
data->binnedClone();
216 if (whichVars ==
nullptr)
219 axes.
add(*whichVars);
222 std::vector<double> min(dim);
223 std::vector<double> max(dim);
224 std::vector<Int_t> bins(dim);
225 std::vector<const char *> names(dim);
227 for (
auto const *var : static_range_cast<RooRealVar *>(axes)) {
228 names[i] = var->GetName();
229 min[i] = var->getMin();
230 max[i] = var->getMax();
231 bins[i] = var->numBins();
236 dim, &bins[0], &min[0], &max[0]);
246 double*
x =
new double[dim];
247 for ( i = 0; i <
size; i++) {
250 for (
Int_t ii = 0; ii < dim; ii++) {
static const char * DATASET_NAME
static const char * NLL_NAME
static const char * DEFAULT_TITLE
static const char * WEIGHT_NAME
static const char * DEFAULT_NAME
size_t size(const MatrixT &matrix)
retrieve the size of a square matrix
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void data
THnSparseT< TArrayF > THnSparseF
RooAbsCollection is an abstract container object that can hold multiple RooAbsArg objects.
Int_t getSize() const
Return the number of elements in the collection.
virtual bool add(const RooAbsArg &var, bool silent=false)
Add the specified argument to list.
virtual RooAbsArg * addClone(const RooAbsArg &var, bool silent=false)
Add a clone of the specified argument to list.
double getRealValue(const char *name, double defVal=0, bool verbose=false) const
Get value of a RooAbsReal stored in set with given name.
RooAbsArg * find(const char *name) const
Find object with given name in list.
virtual Int_t numEntries() const
Return number of entries in dataset, i.e., count unweighted entries.
RooAbsData * reduce(const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg(), const RooCmdArg &arg3=RooCmdArg(), const RooCmdArg &arg4=RooCmdArg(), const RooCmdArg &arg5=RooCmdArg(), const RooCmdArg &arg6=RooCmdArg(), const RooCmdArg &arg7=RooCmdArg(), const RooCmdArg &arg8=RooCmdArg())
Create a reduced copy of this dataset.
RooArgList is a container object that can hold multiple RooAbsArg objects.
RooArgSet is a container object that can hold multiple RooAbsArg objects.
RooCmdArg is a named container for two doubles, two integers two object points and three string point...
The RooDataHist is a container class to hold N-dimensional binned data.
RooDataSet is a container class to hold unbinned data.
const RooArgSet * get(Int_t index) const override
Return RooArgSet with coordinates of event 'index'.
void add(const RooArgSet &row, double weight=1.0, double weightError=0) override
Add one ore more rows of data.
virtual void addFast(const RooArgSet &row, double weight=1.0, double weightError=0)
Add a data point, with its coordinates specified in the 'data' argset, to the data set.
double weight() const override
Return event weight of current event.
RooRealVar represents a variable that can be changed from the outside.
void setVal(double value) override
Set value of variable to 'value'.
Stores the steps in a Markov Chain of points.
virtual RooDataHist * GetAsDataHist(RooArgSet *whichVars=nullptr) const
get this MarkovChain as a RooDataHist whose entries contain the values of whichVars.
virtual void AddFast(RooArgSet &entry, double nllValue, double weight=1.0)
add an entry to the chain ONLY IF you have constructed with parameters or called SetParameters
virtual double NLL() const
get the NLL value of the current (last indexed) entry
virtual void AddWithBurnIn(MarkovChain &otherChain, Int_t burnIn=0)
add another markov chain
virtual void Add(RooArgSet &entry, double nllValue, double weight=1.0)
safely add an entry to the chain
virtual double NLL(Int_t i) const
get the NLL value of entry at position i
virtual THnSparse * GetAsSparseHist(RooAbsCollection *whichVars=nullptr) const
Get a clone of the markov chain on which this interval is based as a sparse histogram.
virtual void SetParameters(RooArgSet ¶meters)
set which of your parameters this chain should store
virtual const RooArgSet * Get(Int_t i) const
get the entry at position i
virtual double Weight() const
get the weight of the current (last indexed) entry
virtual RooDataSet * GetAsDataSet(RooArgSet *whichVars=nullptr) const
get this MarkovChain as a RooDataSet whose entries contain the values of whichVars.
virtual Int_t Size() const
get the number of steps in the chain
Long64_t Fill(const Double_t *x, Double_t w=1.)
Templated implementation of the abstract base THnSparse.
Efficient multidimensional histogram.
void Sumw2() override
Enable calculation of errors.
The TNamed class is the base class for all named ROOT classes.
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
Namespace for the RooStats classes.
void SetParameters(const RooArgSet *desiredVals, RooArgSet *paramsToChange)