class TUnfoldSys: public TUnfold



TUnfoldSys, an extension of the class TUnfold to correct for
migration effects. It provides methods for background subtraction
and propagation of systematic uncertainties

Citation: S.Schmitt, JINST 7 (2012) T10003 [arXiv:1205.6201]

Function Members (Methods)

public:

virtual	~TUnfoldSys()
void	TObject::AbstractMethod(const char* method) const
void	AddSysError(const TH2* sysError, const char* name, TUnfold::EHistMap histmap, TUnfoldSys::ESysErrMode mode)
virtual void	TObject::AppendPad(Option_t* option = "")
virtual void	TObject::Browse(TBrowser* b)
static TClass*	Class()
virtual const char*	TObject::ClassName() const
virtual void	TObject::Clear(Option_t* = "")
virtual TObject*	TObject::Clone(const char* newname = "") const
virtual Int_t	TObject::Compare(const TObject* obj) const
virtual void	TObject::Copy(TObject& object) const
virtual void	TObject::Delete(Option_t* option = "")MENU
virtual Int_t	TObject::DistancetoPrimitive(Int_t px, Int_t py)
virtual Double_t	TUnfold::DoUnfold(Double_t tau)
Double_t	TUnfold::DoUnfold(Double_t tau, const TH1* hist_y, Double_t scaleBias = 0.)
virtual void	TObject::Draw(Option_t* option = "")
virtual void	TObject::DrawClass() constMENU
virtual TObject*	TObject::DrawClone(Option_t* option = "") constMENU
virtual void	TObject::Dump() constMENU
virtual void	TObject::Error(const char* method, const char* msgfmt) const
virtual void	TObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual void	TObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual void	TObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
virtual void	TObject::Fatal(const char* method, const char* msgfmt) const
virtual TObject*	TObject::FindObject(const char* name) const
virtual TObject*	TObject::FindObject(const TObject* obj) const
void	GetBackground(TH1* bgr, const char* bgrSource = 0, const Int_t* binMap = 0, Int_t includeError = 3, Bool_t clearHist = kTRUE) const
void	TUnfold::GetBias(TH1* bias, const Int_t* binMap = 0) const
Double_t	TUnfold::GetChi2A() const
Double_t	TUnfold::GetChi2L() const
Double_t	GetChi2Sys()
Bool_t	GetDeltaSysBackgroundScale(TH1* delta, const char* source, const Int_t* binMap = 0)
Bool_t	GetDeltaSysSource(TH1* hist_delta, const char* source, const Int_t* binMap = 0)
Bool_t	GetDeltaSysTau(TH1* delta, const Int_t* binMap = 0)
virtual Option_t*	TObject::GetDrawOption() const
static Long_t	TObject::GetDtorOnly()
void	TUnfold::GetEmatrix(TH2* ematrix, const Int_t* binMap = 0) const
void	GetEmatrixInput(TH2* ematrix, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)
void	GetEmatrixSysBackgroundScale(TH2* ematrix, const char* source, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)
void	GetEmatrixSysBackgroundUncorr(TH2* ematrix, const char* source, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)
void	GetEmatrixSysSource(TH2* ematrix, const char* source, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)
void	GetEmatrixSysTau(TH2* ematrix, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)
void	GetEmatrixSysUncorr(TH2* ematrix, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)
void	GetEmatrixTotal(TH2* ematrix, const Int_t* binMap = 0)
Double_t	TUnfold::GetEpsMatrix() const
void	TUnfold::GetFoldedOutput(TH1* folded, const Int_t* binMap = 0) const
virtual const char*	TObject::GetIconName() const
void	TUnfold::GetInput(TH1* inputData, const Int_t* binMap = 0) const
void	TUnfold::GetInputInverseEmatrix(TH2* ematrix)
void	TUnfold::GetL(TH2* l) const
virtual Double_t	TUnfold::GetLcurveX() const
virtual Double_t	TUnfold::GetLcurveY() const
void	TUnfold::GetLsquared(TH2* lsquared) const
virtual const char*	TObject::GetName() const
Int_t	TUnfold::GetNdf() const
void	TUnfold::GetNormalisationVector(TH1* s, const Int_t* binMap = 0) const
Int_t	TUnfold::GetNpar() const
Int_t	TUnfold::GetNr() const
virtual char*	TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_t	TObject::GetObjectStat()
virtual Option_t*	TObject::GetOption() const
void	TUnfold::GetOutput(TH1* output, const Int_t* binMap = 0) const
void	TUnfold::GetProbabilityMatrix(TH2* A, TUnfold::EHistMap histmap) const
Double_t	TUnfold::GetRhoAvg() const
Double_t	TUnfold::GetRhoI(TH1* rhoi, const Int_t* binMap = 0, TH2* invEmat = 0) const
void	TUnfold::GetRhoIJ(TH2* rhoij, const Int_t* binMap = 0) const
void	GetRhoItotal(TH1* rhoi, const Int_t* binMap = 0, TH2* invEmat = 0)
Double_t	TUnfold::GetRhoMax() const
Double_t	TUnfold::GetTau() const
virtual const char*	TObject::GetTitle() const
static const char*	TUnfold::GetTUnfoldVersion()
virtual UInt_t	TObject::GetUniqueID() const
virtual Bool_t	TObject::HandleTimer(TTimer* timer)
virtual ULong_t	TObject::Hash() const
virtual void	TObject::Info(const char* method, const char* msgfmt) const
virtual Bool_t	TObject::InheritsFrom(const char* classname) const
virtual Bool_t	TObject::InheritsFrom(const TClass* cl) const
virtual void	TObject::Inspect() constMENU
void	TObject::InvertBit(UInt_t f)
virtual TClass*	IsA() const
virtual Bool_t	TObject::IsEqual(const TObject* obj) const
virtual Bool_t	TObject::IsFolder() const
Bool_t	TObject::IsOnHeap() const
virtual Bool_t	TObject::IsSortable() const
Bool_t	TObject::IsZombie() const
virtual void	TObject::ls(Option_t* option = "") const
void	TObject::MayNotUse(const char* method) const
virtual Bool_t	TObject::Notify()
void	TObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const
void	TObject::operator delete(void* ptr)
void	TObject::operator delete(void* ptr, void* vp)
void	TObject::operator delete[](void* ptr)
void	TObject::operator delete[](void* ptr, void* vp)
void*	TObject::operator new(size_t sz)
void*	TObject::operator new(size_t sz, void* vp)
void*	TObject::operator new[](size_t sz)
void*	TObject::operator new[](size_t sz, void* vp)
TUnfoldSys&	operator=(const TUnfoldSys&)
virtual void	TObject::Paint(Option_t* option = "")
virtual void	TObject::Pop()
virtual void	TObject::Print(Option_t* option = "") const
virtual Int_t	TObject::Read(const char* name)
virtual void	TObject::RecursiveRemove(TObject* obj)
Int_t	TUnfold::RegularizeBins(int start, int step, int nbin, TUnfold::ERegMode regmode)
Int_t	TUnfold::RegularizeBins2D(int start_bin, int step1, int nbin1, int step2, int nbin2, TUnfold::ERegMode regmode)
Int_t	TUnfold::RegularizeCurvature(int left_bin, int center_bin, int right_bin, Double_t scale_left = 1., Double_t scale_right = 1.)
Int_t	TUnfold::RegularizeDerivative(int left_bin, int right_bin, Double_t scale = 1.)
Int_t	TUnfold::RegularizeSize(int bin, Double_t scale = 1.)
void	TObject::ResetBit(UInt_t f)
virtual void	TObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU
virtual void	TObject::SavePrimitive(ostream& out, Option_t* option = "")
virtual Int_t	TUnfold::ScanLcurve(Int_t nPoint, Double_t tauMin, Double_t tauMax, TGraph lCurve, TSpline logTauX = 0, TSpline** logTauY = 0)
void	TUnfold::SetBias(const TH1* bias)
void	TObject::SetBit(UInt_t f)
void	TObject::SetBit(UInt_t f, Bool_t set)
void	TUnfold::SetConstraint(TUnfold::EConstraint constraint)
virtual void	TObject::SetDrawOption(Option_t* option = "")MENU
static void	TObject::SetDtorOnly(void* obj)
void	TUnfold::SetEpsMatrix(Double_t eps)
virtual Int_t	SetInput(const TH1* hist_y, Double_t scaleBias = 0., Double_t oneOverZeroError = 0., const TH2* hist_vyy = 0, const TH2* hist_vyy_inv = 0)
static void	TObject::SetObjectStat(Bool_t stat)
void	SetTauError(Double_t delta_tau)
virtual void	TObject::SetUniqueID(UInt_t uid)
virtual void	ShowMembers(TMemberInspector& insp) const
virtual void	Streamer(TBuffer&)
void	StreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
void	SubtractBackground(const TH1* hist_bgr, const char* name, Double_t scale = 1., Double_t scale_error = 0.)
virtual void	TObject::SysError(const char* method, const char* msgfmt) const
Bool_t	TObject::TestBit(UInt_t f) const
Int_t	TObject::TestBits(UInt_t f) const
	TUnfoldSys(const TUnfoldSys&)
	TUnfoldSys(const TH2* hist_A, TUnfold::EHistMap histmap, TUnfold::ERegMode regmode = kRegModeSize, TUnfold::EConstraint constraint = kEConstraintArea)
virtual void	TObject::UseCurrentStyle()
virtual void	TObject::Warning(const char* method, const char* msgfmt) const
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const

protected:

void	TUnfold::AddMSparse(TMatrixDSparse* dest, Double_t f, const TMatrixDSparse* src) const
Bool_t	TUnfold::AddRegularisationCondition(Int_t nEle, const Int_t* indices, const Double_t* rowData)
Bool_t	TUnfold::AddRegularisationCondition(Int_t i0, Double_t f0, Int_t i1 = -1, Double_t f1 = 0., Int_t i2 = -1, Double_t f2 = 0.)
void	TUnfold::ClearHistogram(TH1* h, Double_t x = 0.) const
virtual void	ClearResults()
TMatrixDSparse*	TUnfold::CreateSparseMatrix(Int_t nrow, Int_t ncol, Int_t nele, Int_t* row, Int_t* col, Double_t* data) const
static void	TUnfold::DeleteMatrix(TMatrixD** m)
static void	TUnfold::DeleteMatrix(TMatrixDSparse** m)
void	DoBackgroundSubtraction()
virtual void	TObject::DoError(int level, const char* location, const char* fmt, va_list va) const
virtual Double_t	TUnfold::DoUnfold()
void	TUnfold::ErrorMatrixToHist(TH2* ematrix, const TMatrixDSparse* emat, const Int_t* binMap, Bool_t doClear) const
const TMatrixDSparse*	TUnfold::GetAx() const
Int_t	TUnfold::GetBinFromRow(int ix) const
const TMatrixDSparse*	TUnfold::GetDXDAM(int i) const
const TMatrixDSparse*	TUnfold::GetDXDAZ(int i) const
const TMatrixDSparse*	TUnfold::GetDXDtauSquared() const
const TMatrixDSparse*	TUnfold::GetDXDY() const
const TMatrixDSparse*	TUnfold::GetE() const
const TMatrixDSparse*	TUnfold::GetEinv() const
void	GetEmatrixFromVyy(const TMatrixDSparse* vyy, TH2* ematrix, const Int_t* binMap, Bool_t clearEmat)
Int_t	TUnfold::GetNx() const
Int_t	TUnfold::GetNy() const
virtual TString	TUnfold::GetOutputBinName(Int_t iBinX) const
Double_t	TUnfold::GetRhoIFromMatrix(TH1* rhoi, const TMatrixDSparse* eOrig, const Int_t* binMap, TH2* invEmat) const
Int_t	TUnfold::GetRowFromBin(int ix) const
TMatrixDSparse*	GetSummedErrorMatrixXX()
TMatrixDSparse*	GetSummedErrorMatrixYY()
const TMatrixDSparse*	TUnfold::GetVxx() const
const TMatrixDSparse*	TUnfold::GetVxxInv() const
const TMatrixDSparse*	TUnfold::GetVyyInv() const
const TMatrixD*	TUnfold::GetX() const
TMatrixDSparse*	TUnfold::InvertMSparseSymmPos(const TMatrixDSparse* A, Int_t* rank) const
void	TObject::MakeZombie()
TMatrixDSparse*	TUnfold::MultiplyMSparseM(const TMatrixDSparse* a, const TMatrixD* b) const
TMatrixDSparse*	TUnfold::MultiplyMSparseMSparse(const TMatrixDSparse* a, const TMatrixDSparse* b) const
TMatrixDSparse*	TUnfold::MultiplyMSparseMSparseTranspVector(const TMatrixDSparse* m1, const TMatrixDSparse* m2, const TMatrixTBase<Double_t>* v) const
TMatrixDSparse*	TUnfold::MultiplyMSparseTranspMSparse(const TMatrixDSparse* a, const TMatrixDSparse* b) const
virtual TMatrixDSparse*	PrepareCorrEmat(const TMatrixDSparse* m1, const TMatrixDSparse* m2, const TMatrixDSparse* dsys)
virtual void	PrepareSysError()
virtual TMatrixDSparse*	PrepareUncorrEmat(const TMatrixDSparse* m1, const TMatrixDSparse* m2)
void	ScaleColumnsByVector(TMatrixDSparse* m, const TMatrixTBase<Double_t>* v) const
	TUnfoldSys()
void	VectorMapToHist(TH1* hist_delta, const TMatrixDSparse* delta, const Int_t* binMap)

private:

void	InitTUnfoldSys()

Data Members

public:

static TObject::(anonymous)	TObject::kBitMask
static TObject::EStatusBits	TObject::kCanDelete
static TObject::EStatusBits	TObject::kCannotPick
static TUnfold::EConstraint	TUnfold::kEConstraintArea
static TUnfold::EConstraint	TUnfold::kEConstraintNone
static TObject::EStatusBits	TObject::kHasUUID
static TUnfold::EHistMap	TUnfold::kHistMapOutputHoriz
static TUnfold::EHistMap	TUnfold::kHistMapOutputVert
static TObject::EStatusBits	TObject::kInvalidObject
static TObject::(anonymous)	TObject::kIsOnHeap
static TObject::EStatusBits	TObject::kIsReferenced
static TObject::EStatusBits	TObject::kMustCleanup
static TObject::EStatusBits	TObject::kNoContextMenu
static TObject::(anonymous)	TObject::kNotDeleted
static TObject::EStatusBits	TObject::kObjInCanvas
static TObject::(anonymous)	TObject::kOverwrite
static TUnfold::ERegMode	TUnfold::kRegModeCurvature
static TUnfold::ERegMode	TUnfold::kRegModeDerivative
static TUnfold::ERegMode	TUnfold::kRegModeMixed
static TUnfold::ERegMode	TUnfold::kRegModeNone
static TUnfold::ERegMode	TUnfold::kRegModeSize
static TObject::(anonymous)	TObject::kSingleKey
static TUnfoldSys::ESysErrMode	kSysErrModeMatrix
static TUnfoldSys::ESysErrMode	kSysErrModeRelative
static TUnfoldSys::ESysErrMode	kSysErrModeShift
static TObject::(anonymous)	TObject::kWriteDelete
static TObject::(anonymous)	TObject::kZombie

protected:

TMatrixDSparse*	TUnfold::fA	Input: matrix
TMatrixD*	fAoutside	Input: underflow/overflow bins
TMap*	fBgrErrScaleIn	Input: background sources correlated error
TMap*	fBgrErrUncorrInSq	Input: uncorr error squared from bgr sources
TMap*	fBgrIn	Input: size of background sources
Double_t	TUnfold::fBiasScale	Input: scale factor for the bias
TUnfold::EConstraint	TUnfold::fConstraint	Input: type of constraint to use
TMatrixD*	fDAinColRelSq	Input: normalized column err.sq. (inp.matr.)
TMatrixDSparse*	fDAinRelSq	Input: normalized errors from input matrix
TMap*	fDeltaCorrAx	Result: syst.shift from fSysIn on fAx
TMap*	fDeltaCorrX	Result: syst.shift from fSysIn on fX
TMatrixDSparse*	fDeltaSysTau	Result: systematic shift from tau
Double_t	fDtau	Input: error on tau
TMatrixDSparse*	fEmatUncorrAx	Result: syst.error from fDA2 on fAx
TMatrixDSparse*	fEmatUncorrX	Result: syst.error from fDA2 on fX
TArrayI	TUnfold::fHistToX	Input: histogram bins -> matrix indices
TMatrixDSparse*	TUnfold::fL	Input: regularisation conditions
TUnfold::ERegMode	TUnfold::fRegMode	Input: type of regularisation
TArrayD	TUnfold::fSumOverY	Input: sum of all columns
TMap*	fSysIn	Input: correlated errors
Double_t	TUnfold::fTauSquared	Input: regularisation parameter
TMatrixDSparse*	TUnfold::fVyy	Input: covariance matrix for y
TMatrixDSparse*	fVyyData	Input: error on fY prior to bgr subtraction
TMatrixD*	TUnfold::fX0	Input: x0
TArrayI	TUnfold::fXToHist	Input: matrix indices -> histogram bins
TMatrixD*	TUnfold::fY	Input: y
TMatrixD*	fYData	Input: fY prior to bgr subtraction

Class Charts

Inheritance Chart:

TObject

←

TUnfold

←

TUnfoldSys

←

TUnfoldDensity

Function documentation

void InitTUnfoldSys(void)

TUnfoldSys(const TUnfoldSys& )

void ClearResults(void)

void PrepareSysError(void)

TMatrixDSparse * PrepareUncorrEmat(const TMatrixDSparse* m1, const TMatrixDSparse* m2)

TMatrixDSparse * PrepareCorrEmat(const TMatrixDSparse* m1, const TMatrixDSparse* m2, const TMatrixDSparse* dsys)

void ScaleColumnsByVector(TMatrixDSparse* m, const TMatrixTBase<Double_t>* v) const

void VectorMapToHist(TH1* hist_delta, const TMatrixDSparse* delta, const Int_t* binMap)

void GetEmatrixFromVyy(const TMatrixDSparse* vyy, TH2* ematrix, const Int_t* binMap, Bool_t clearEmat)

void DoBackgroundSubtraction(void)

TMatrixDSparse * GetSummedErrorMatrixYY(void)

TMatrixDSparse * GetSummedErrorMatrixXX(void)

TUnfoldSys(const TH2* hist_A, TUnfold::EHistMap histmap, TUnfold::ERegMode regmode = kRegModeSize, TUnfold::EConstraint constraint = kEConstraintArea)

void AddSysError(const TH2* sysError, const char* name, TUnfold::EHistMap histmap, TUnfoldSys::ESysErrMode mode)

void SubtractBackground(const TH1* hist_bgr, const char* name, Double_t scale = 1., Double_t scale_error = 0.)

Int_t SetInput(const TH1* hist_y, Double_t scaleBias = 0., Double_t oneOverZeroError = 0., const TH2* hist_vyy = 0, const TH2* hist_vyy_inv = 0)

Bool_t GetDeltaSysBackgroundScale(TH1* delta, const char* source, const Int_t* binMap = 0)

Bool_t GetDeltaSysTau(TH1* delta, const Int_t* binMap = 0)

void GetEmatrixSysUncorr(TH2* ematrix, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)

void GetEmatrixSysSource(TH2* ematrix, const char* source, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)

void GetEmatrixSysBackgroundScale(TH2* ematrix, const char* source, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)

void GetEmatrixInput(TH2* ematrix, const Int_t* binMap = 0, Bool_t clearEmat = kTRUE)

void GetEmatrixTotal(TH2* ematrix, const Int_t* binMap = 0)

void GetRhoItotal(TH1* rhoi, const Int_t* binMap = 0, TH2* invEmat = 0)

Double_t GetChi2Sys(void)