/usr/include/root/RooMultiVarGaussian.h is in libroot-roofit-dev 5.34.00-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 | /*****************************************************************************
* Project: RooFit *
* Package: RooFitModels *
* File: $Id: RooMultiVarGaussian.h 44368 2012-05-30 15:38:44Z axel $
* Authors: *
* WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu *
* DK, David Kirkby, UC Irvine, dkirkby@uci.edu *
* *
* Copyright (c) 2000-2005, Regents of the University of California *
* and Stanford University. All rights reserved. *
* *
* Redistribution and use in source and binary forms, *
* with or without modification, are permitted according to the terms *
* listed in LICENSE (http://roofit.sourceforge.net/license.txt) *
*****************************************************************************/
#ifndef ROO_MULTI_VAR_GAUSSIAN
#define ROO_MULTI_VAR_GAUSSIAN
#include "RooAbsPdf.h"
#include "RooListProxy.h"
#include "TMatrixDSym.h"
#include "TMatrixD.h"
#include "TVectorD.h"
class RooRealVar;
class RooFitResult ;
#include <map>
#include <vector>
class RooMultiVarGaussian : public RooAbsPdf {
public:
RooMultiVarGaussian() {} ;
RooMultiVarGaussian(const char *name, const char *title, const RooArgList& xvec, const RooFitResult& fr, Bool_t reduceToConditional=kTRUE) ;
RooMultiVarGaussian(const char *name, const char *title, const RooArgList& xvec, const RooArgList& mu, const TMatrixDSym& covMatrix) ;
RooMultiVarGaussian(const char *name, const char *title, const RooArgList& xvec, const TVectorD& mu, const TMatrixDSym& covMatrix) ;
RooMultiVarGaussian(const char *name, const char *title, const RooArgList& xvec,const TMatrixDSym& covMatrix) ;
void setAnaIntZ(Double_t z) { _z = z ; }
RooMultiVarGaussian(const RooMultiVarGaussian& other, const char* name=0) ;
virtual TObject* clone(const char* newname) const { return new RooMultiVarGaussian(*this,newname); }
inline virtual ~RooMultiVarGaussian() { }
Int_t getAnalyticalIntegral(RooArgSet& allVars, RooArgSet& analVars, const char* rangeName=0) const ;
Double_t analyticalIntegral(Int_t code, const char* rangeName=0) const ;
Int_t getGenerator(const RooArgSet& directVars, RooArgSet &generateVars, Bool_t staticInitOK=kTRUE) const;
void initGenerator(Int_t code) ;
void generateEvent(Int_t code);
const TMatrixDSym& covarianceMatrix() const { return _cov ; }
class AnaIntData {
public:
TMatrixD S22bar ;
Double_t S22det ;
std::vector<int> pmap ;
Int_t nint ;
} ;
class GenData {
public:
TMatrixD UT ;
std::vector<int> omap ;
std::vector<int> pmap ;
TVectorD mu1 ;
TVectorD mu2 ;
TMatrixD S12S22I ;
} ;
class BitBlock {
public:
BitBlock() : b0(0), b1(0), b2(0), b3(0) {} ;
void setBit(Int_t ibit) ;
Bool_t getBit(Int_t ibit) ;
Bool_t operator==(const BitBlock& other) ;
Int_t b0 ;
Int_t b1 ;
Int_t b2 ;
Int_t b3 ;
} ;
static void blockDecompose(const TMatrixD& input, const std::vector<int>& map1, const std::vector<int>& map2, TMatrixDSym& S11, TMatrixD& S12, TMatrixD& S21, TMatrixDSym& S22) ;
protected:
void decodeCode(Int_t code, std::vector<int>& map1, std::vector<int>& map2) const;
AnaIntData& anaIntData(Int_t code) const ;
GenData& genData(Int_t code) const ;
mutable std::map<int,AnaIntData> _anaIntCache ; //!
mutable std::map<int,GenData> _genCache ; //!
mutable std::vector<BitBlock> _aicMap ; //!
RooListProxy _x ;
RooListProxy _mu ;
TMatrixDSym _cov ;
TMatrixDSym _covI ;
Double_t _det ;
Double_t _z ;
void syncMuVec() const ;
mutable TVectorD _muVec ; //! Do not persist
Double_t evaluate() const ;
private:
ClassDef(RooMultiVarGaussian,1) // Multivariate Gaussian PDF with correlations
};
#endif
|