MjjYmax

#include <MjjYmax.h>

Inheritance diagram for MjjYmax:

Collaboration diagram for MjjYmax:

Public Member Functions
	MjjYmax ()
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

Public Attributes
int	nGenMjjBinsFwd
int	nGenMjjBinsCnt
Public Attributes inherited from Observable
TUnfoldBinning *	recBinning
TUnfoldBinning *	genBinning

Additional Inherited Members
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

MjjYmax()

MjjYmax

(

)

Constructor.

                  :
    Observable(__FUNCTION__, "Dijet mass double differential cross section")
{
    vector<double> recMjj_edges;
    recMjj_edges.reserve(Mjj_edges.size()*2);
 
    // rec binning with twice more bins
    for (size_t j = 0; j < Mjj_edges.size()-1; ++j) {
        double m = Mjj_edges.at(j),
               M = Mjj_edges.at(j+1);
        recMjj_edges.push_back(m);
        recMjj_edges.push_back((m+M)/2);
    }
    recMjj_edges.push_back(Mjj_edges.back());
 
    int nRecMjjBins = recMjj_edges.size()-1;
    int nGenMjjBins =    Mjj_edges.size()-1;
 
    recBinning->AddAxis("Mjj", nRecMjjBins, recMjj_edges.data(),false,false);
    recBinning->AddAxis("y"  ,   nYmaxBins,   ymax_edges.data(),false,false);
    genBinning->AddAxis("Mjj", nGenMjjBins,    Mjj_edges.data(),false,false);
    genBinning->AddAxis("y"  ,   nYmaxBins,   ymax_edges.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<MjjYmaxFiller>(*this, reader);
}

setLmatrix()

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

overridevirtual

See Observable::getLmatrix

Todo:

fix

Reimplemented from Observable.

{
    for (size_t iy = 1; iy <= nYmaxBins; ++iy) {
        const TVectorD * Mjj_edges = genBinning->GetDistributionBinning(0);
        int nMjjBins = Mjj_edges->GetNrows()-1;
        for (int iM = 1; iM <= nMjjBins; ++iM) {
 
            double Ymax = (  y_edges.at(iy-1) +   y_edges.at(iy)) / 2,
                   Mjj  = ((*Mjj_edges)[iM-1] + (*Mjj_edges)[iM]) / 2;
            int i = genBinning->GetGlobalBinNumber(Mjj, Ymax);
            if (i == 0)
                BOOST_THROW_EXCEPTION( runtime_error("Mjj = "s + Mjj + " and Ymax = "s + Ymax
                                        + " do not correspond to any bin index"s) );
 
            // b(in)
            int                 bUp     = i-nMjjBins,
                 bLeft   = i-1, bCenter = i         , bRight  = i+1,
                                bDown   = i+nMjjBins;
 
            // values (curvature regularisation)
            auto get = [&bias](int i) {
                double content = bias->GetBinContent(i);
                if (content < 0)
                    BOOST_THROW_EXCEPTION( DE::BadInput("Expecting only positive entries", bias) );
                return content > 0 ? 1./content : 0;
            };
            double cUp     = get(bUp   ),
                   cLeft   = get(bLeft ),
                   cRight  = get(bRight),
                   cDown   = get(bDown );
 
            cout << setw(3) << iy << setw(3) << iM
                 << 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;
}

Member Data Documentation

nGenMjjBinsCnt

int nGenMjjBinsCnt

nGenMjjBinsFwd

int nGenMjjBinsFwd

---

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

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