HTnTransformer

#include <HTn.h>

Inheritance diagram for HTnTransformer:

Collaboration diagram for HTnTransformer:

Public Member Functions
	HTnTransformer (TUnfoldBinning *)
	~HTnTransformer ()=default
void	Transform (const Eigen::VectorXd &) const override
Public Member Functions inherited from Transformer
virtual	~Transformer ()=default
virtual void	RemoveBadInputBins (TH1 *, TH2 *)
	Transformer (TUnfoldBinning *, bool=true)

Public Attributes
const size_t	nHtBins
Public Attributes inherited from Transformer
TUnfoldBinning *	preBinning
TUnfoldBinning *	postBinning

Additional Inherited Members
Static Public Attributes inherited from Transformer
static Eigen::VectorXd	y
Static Protected Member Functions inherited from Transformer
static void	AddAxis (TUnfoldBinning *, TUnfoldBinning *, int)

Constructor & Destructor Documentation

HTnTransformer()

HTnTransformer

(

TUnfoldBinning *

pre

)

Constructor.

Parameters

pre	a priori, gen or rec binning

                                                    :
    Transformer(pre, false), nHtBins(pre->GetDistributionBinning(0)->GetNoElements()-1)
{
    AddAxis(preBinning, postBinning, 0);
 
    postBinning->AddAxis("j", nJetsBins-1, &n_jets_edges[1] /* 3+ */, false, false);
    postBinning->AddAxis("i", nJetsBins-1, &n_jets_edges[0] /* 2+ */, false, false);
}

~HTnTransformer()

~HTnTransformer ( )

default

Destructor.

Member Function Documentation

Transform()

void Transform ( const Eigen::VectorXd &  ) const

overridevirtual

Calculate R32, R43, R42, etc. from the 2-jet, 3-jet, 4-jet, etc. inclusive jet cross sections in bins of HT.

Todo:

harmonise (either nJetsmin, or 2)

Reimplemented from Transformer.

{
    VectorXd z = x.segment(preBinning->GetStartBin()-1,
                           preBinning->GetDistributionNumberOfBins());
 
    // 1) from excl mult to incl mult
    for (int n = nJetsmax; n > nJetsmin; --n)
        for (long iHt = 0; iHt < static_cast<long>(nHtBins); ++iHt)
            z(nHtBins*(n-3) + iHt) += z(nHtBins*(n-2) + iHt);
 
    // 2) from incl mult to ratios
    for (int j = 3; j <= nJetsmax; ++j)
    for (int i = 2; i < j; ++i) {
        for (size_t iHt = 0; iHt < nHtBins; ++iHt) {
            double num = z(nHtBins*(j-2) + iHt),
                   den = z(nHtBins*(i-2) + iHt);
            if (std::abs(den) < deps) continue;
            double Ht = postBinning->GetDistributionBinCenter(0, iHt);
            int k = postBinning->GetGlobalBinNumber(Ht, j, i);
            y(k-1) = num / den;
        }
    }
}

Member Data Documentation

nHtBins

const size_t nHtBins

---

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

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Unfolding/interface/HTn.h
• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Unfolding/src/HTn.cc