ZPtY

Description

Following SMP-17-010 (https://doi.org/10.17182/hepdata.91215.v4)

#include <ZPtY.h>

Inheritance diagram for ZPtY:

Collaboration diagram for ZPtY:

Public Member Functions
	ZPtY ()
std::unique_ptr< Filler >	getFiller (TTreeReader &reader) const override
void	setLmatrix (const std::unique_ptr< TH1 > &, std::unique_ptr< TH2 > &) override
Public Member Functions inherited from Observable
virtual std::unique_ptr< Transformer >	getTransformer (TUnfoldBinning *) const
virtual	~Observable ()=default

Additional Inherited Members
Public Attributes inherited from Observable
TUnfoldBinning *	recBinning
TUnfoldBinning *	genBinning
Static Public Attributes inherited from Observable
static bool	isMC = false
static double	maxDR = 0.2
Protected Member Functions inherited from Observable
	Observable (const char *, const char *)

Constructor & Destructor Documentation

ZPtY()

ZPtY

(

)

Constructor.

            :
    Observable(__FUNCTION__, "Z boson dsigma/dpT dy")
{
    recBinning->AddAxis("pt", nRecPtBins, recPtBins.data(), false, false);
    recBinning->AddAxis("y" ,     nYbins,     yBins.data(), false, false);
    genBinning->AddAxis("pt", nGenPtBins, genPtBins.data(), false, false);
    genBinning->AddAxis("y" ,     nYbins,     yBins.data(), false, false);
}

Member Function Documentation

getFiller()

unique_ptr< DAS::Unfolding::Filler > getFiller ( TTreeReader &  reader ) const

overridevirtual

Constructs a filler for the observable.

Reimplemented from Observable.

{
    return make_unique<ZPtYFiller>(*this, reader);
}

setLmatrix()

void setLmatrix ( const std::unique_ptr< TH1 > &  bias, std::unique_ptr< TH2 > &  L  )

overridevirtual

See Observable::setLmatrix

Todo:

some docs... this is quite obscure

Todo:

Most existing L-matrices follow a similar structure, this should be implemented in a function instead of repeating the same code

Reimplemented from Observable.

{
    for (int iy = 1; iy <= nRecPtBins; ++iy)
    for (int ipt = 1; ipt <= nGenPtBins; ++ipt) {
 
        auto  y = (yBins.at(iy - 1) + yBins.at(iy)) / 2,
             pt = (genPtBins.at(ipt - 1) + genPtBins.at(ipt)) / 2;
        int i = genBinning->GetGlobalBinNumber(pt, y);
        if (i == 0)
            BOOST_THROW_EXCEPTION( logic_error(
                Form("pt = %f and y = %f do not correspond to any bin index", pt, y)) );
 
        // b(in)
        int               bUp     = i-nGenPtBins,
             bLeft = i-1, bCenter = i           , bRight  = i+1,
                          bDown   = i+nGenPtBins;
 
        // values (curvature regularisation)
        auto get = [&bias](int i) {
            auto content = bias->GetBinContent(i);
            if (content < 0)
                BOOST_THROW_EXCEPTION( DE::BadInput("Expecting only positive entries", bias) );
            return content > 0 ? 1./content : 0;
        };
 
        auto cUp    = get(bUp   ),
             cLeft  = get(bLeft ),
             cRight = get(bRight),
             cDown  = get(bDown );
 
        cout << setw(3) << iy << setw(3) << ipt
             << setw(5) << bCenter << setw(15) << (cUp+cLeft+cRight+cDown)
             << setw(5) << bUp     << setw(15) << -cUp
             << setw(5) << bLeft   << setw(15) <<     -cLeft
             << setw(5) << bRight  << setw(15) <<           -cRight
             << setw(5) << bDown   << setw(15) <<                  -cDown << '\n';
 
        // filling L-matrix
                        L->SetBinContent(i, bCenter,  cUp+cLeft+cRight+cDown );
        if (cUp    > 0) L->SetBinContent(i, bUp    , -cUp                    );
        if (cLeft  > 0) L->SetBinContent(i, bLeft  ,     -cLeft              );
        if (cRight > 0) L->SetBinContent(i, bRight ,           -cRight       );
        if (cDown  > 0) L->SetBinContent(i, bDown  ,                  -cDown );
    }
    cout << flush;
}

---

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

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