Matrix

Description

Functor to extract a 2D distribution with physical axes from a 2ND distribution with bin IDs.

Intended to be used together with axisLoop and nodeLoop.

Collaboration diagram for Matrix:

Public Member Functions
	Matrix (TH2 *h, TUnfoldBinning *bng2)
void	operator() (TDirectory *dOut, TString name2, TUnfoldBinning *node2, vector< double > x2, TString title2)

Public Attributes
TH2 *	hIn
TUnfoldBinning *	binning2
TString	name1
TUnfoldBinning *	node1
vector< int >	ibins1
vector< double >	x1
TString	title1

Constructor & Destructor Documentation

Matrix()

Matrix ( TH2 *  h, TUnfoldBinning *  bng2  )

inline

                                            :
        hIn(h), binning2(bng2), node1(nullptr)
    { }

Member Function Documentation

operator()()

void operator() ( TDirectory *  dOut, TString  name2, TUnfoldBinning *  node2, vector< double >  x2, TString  title2  )

inline

Parameters

dOut	output directory
name2	partial name (w.r.t. 2nd axis)
node2	local node (2nd axis)
x2	physical values (2nd axis)
title2	simplistic title, based on bin edges

    {
        const TVectorD * edges1 = node1->GetDistributionBinning(0),
                       * edges2 = node2->GetDistributionBinning(0);
        int nbins1 = edges1->GetNrows()-1,
            nbins2 = edges2->GetNrows()-1;
        TString name = name1 + "__" + name2,
                title = title1 + "-- " + title2;
 
//    const char * name = Form("%s_%s"           , X.h->GetName()             , Y.h->GetName()             ),
//               * title = Form("%s, %s;%s;%s;%s", X.h->GetTitle()            , Y.h->GetTitle()            ,
//                                                 X.h->GetXaxis()->GetTitle(), Y.h->GetXaxis()->GetTitle(),
//                                                 h->GetYaxis()->GetTitle()                               );
 
        TH2F hOut(name, title, nbins1, edges1->GetMatrixArray(),
                               nbins2, edges2->GetMatrixArray());
        for (int ibin1 = 1; ibin1 <= nbins1; ++ibin1)
        for (int ibin2 = 1; ibin2 <= nbins2; ++ibin2) {
            x1.front() = hOut.GetXaxis()->GetBinLowEdge(ibin1);
            x2.front() = hOut.GetYaxis()->GetBinLowEdge(ibin2);
            int Ibin1 = node1->GetGlobalBinNumber(x1.data()),
                Ibin2 = node2->GetGlobalBinNumber(x2.data());
            auto content = hIn->GetBinContent(Ibin1, Ibin2);
            hOut.SetBinContent(ibin1, ibin2, content);
        }
        hOut.SetDirectory(dOut);
        hOut.SetMinimum(hIn->GetMinimum());
        hOut.SetMaximum(hIn->GetMaximum());
        hOut.Write();
    }

Member Data Documentation

binning2

TUnfoldBinning* binning2

global binning (2nd axis)

hIn

TH2* hIn

input histogram with bin IDs

ibins1

vector<int> ibins1

global bin IDs of present iteration

name1

TString name1

partial name (w.r.t. 1st axis)

node1

TUnfoldBinning* node1

local node (1st axis)

title1

TString title1

simplistic title (1st axis)

x1

vector<double> x1

physical values (1st axis)

---

The documentation for this struct was generated from the following file:

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Unfolding/bin/getUnfPhysDist.cc