MjjYmaxFiller

#include <MjjYmax.h>

Inheritance diagram for MjjYmaxFiller:

Collaboration diagram for MjjYmaxFiller:

Public Member Functions
	MjjYmaxFiller (const MjjYmax &, Darwin::Tools::Flow &)
std::list< int >	fillRec (DistVariation &) override
void	match () override
void	fillMC (DistVariation &) override
Public Member Functions inherited from Filler
virtual	~Filler ()=default

Public Attributes
MjjYmax	obs
std::vector< GenJet > *	genJets
std::vector< RecJet > *	recJets
GenEvent *	gEv
RecEvent *	rEv
std::optional< bool >	matched

Constructor & Destructor Documentation

MjjYmaxFiller()

MjjYmaxFiller	(	const MjjYmax &	obs,
		Darwin::Tools::Flow &	flow
	)

Constructor.

    : obs(obs)
    , genJets(obs.isMC ? flow.GetBranchReadOnly<vector<GenJet>>("genJets") : nullptr)
    , recJets(flow.GetBranchReadOnly<vector<RecJet>>("recJets"))
    , gEv(obs.isMC ? flow.GetBranchReadOnly<GenEvent>("genEvent") : nullptr)
    , rEv(flow.GetBranchReadOnly<RecEvent>("recEvent"))
{
}

Member Function Documentation

fillMC()

void fillMC ( DistVariation &  v )

overridevirtual

See Observable::fillMC

Reimplemented from Filler.

{
    if (!obs.isMC)
        BOOST_THROW_EXCEPTION( runtime_error(__func__ + " should only be called for MC"s) );
 
    auto recdijet = selection(*recJets, v);
    auto irec = getBinNumber(recdijet, v, obs.recBinning);
    double recW = rEv->Weight(v);
    if (irec > 0) recW *= recdijet.Weight(v);
 
    auto gendijet = selection(*genJets, v);
    auto igen = getBinNumber(gendijet, v, obs.genBinning);
    double genW = gEv->Weight(v);
    if (igen > 0) genW *= gendijet.Weight(v);
 
    if (igen > 0) v.gen->Fill(igen, genW);
 
    if (*matched) {
        if      (irec > 0  && igen >  0) {    v.RM->Fill(igen, irec, genW *    recW );
                                     v.missNoMatch->Fill(igen,       genW * (1-recW)); }
        else if (irec == 0 && igen >  0) v.missOut->Fill(igen,       genW           );
        else if (irec > 0  && igen == 0) v.fakeOut->Fill(      irec, genW *    recW );
    }
    else {
        if (igen > 0) v.missNoMatch->Fill(igen, genW       );
        if (irec > 0) v.fakeNoMatch->Fill(irec, genW * recW);
    }
}

fillRec()

list< int > fillRec ( DistVariation &  v )

overridevirtual

See Observable::fillRec

Reimplemented from Filler.

{
    auto dijet = selection(*recJets, v);
    if (!dijet) return {};
 
    int i = getBinNumber(dijet, v, obs.recBinning);
    if (i == 0) return {};
 
    double w = rEv->Weight(v);
    if (obs.isMC) w *= gEv->Weight(v);
    w *= dijet.Weight(v);
 
    v.tmp->Fill(i, w);
    v.rec->Fill(i, w);
 
    return list<int>{i};
}

match()

void match ( )

overridevirtual

Match the two pairs of leading jets (if any) and set matched member.

Reimplemented from Filler.

{
    matched.reset();
 
    auto match = [this](size_t i, size_t j) {
        const FourVector& g = genJets->at(i).p4,
                          r = recJets->at(j).p4;
        using ROOT::Math::VectorUtil::DeltaR;
        auto DR = DeltaR(g, r);
        //cout << g << '\t' << r << '\t' << DR << '\t' << result << '\n';
        return DR < obs.maxDR;
    };
 
    // matching (swapping leading and subleading is allowed)
    matched = genJets->size() > 1 && recJets->size() > 1
              && (   (match(0,0) && match(1,1))
                  || (match(0,1) && match(1,0)) );
}

Member Data Documentation

genJets

std::vector<GenJet>* genJets

gEv

GenEvent* gEv

matched

std::optional<bool> matched

obs

MjjYmax obs

Backreference to the observable.

recJets

std::vector<RecJet>* recJets

rEv

RecEvent* rEv

---

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