HTnFiller

#include <HTn.h>

Inheritance diagram for HTnFiller:

Collaboration diagram for HTnFiller:

Public Member Functions
	HTnFiller (const HTn &, 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
HTn	obs
std::vector< GenJet > *	genJets
std::vector< RecJet > *	recJets
GenEvent *	gEv
RecEvent *	rEv
std::vector< std::pair< GenJet, RecJet > >	matches
std::vector< GenJet >	misses
std::vector< RecJet >	fakes
std::optional< bool >	matched

Constructor & Destructor Documentation

HTnFiller()

HTnFiller	(	const HTn &	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 Filler::fillMC

Reimplemented from Filler.

{
    if (!obs.isMC)
        BOOST_THROW_EXCEPTION( runtime_error(__func__ + " should only be called for MC"s) );
 
    auto recdijet = selection<RecJet>(*recJets, v);
    auto [recn,recW] = getMultiplicity<RecJet>(*recJets, v);
    auto irec = getBinNumber<RecJet>(recdijet, recn, v, obs.recBinning);
    recW *= rEv->Weight(v);
    if (irec > 0) recW *= recdijet.Weight(v);
 
    auto gendijet = selection<GenJet>(*genJets, v);
    auto [genn,genW] = getMultiplicity<GenJet>(*genJets, v);
    auto igen = getBinNumber<GenJet>(gendijet, genn, v, obs.genBinning);
    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 Filler::fillRec

Reimplemented from Filler.

{
    auto dijet = selection<RecJet>(*recJets, v);
    if (!dijet) return {};
 
    auto [n,w] = getMultiplicity(*recJets, v);
    int i = getBinNumber(dijet, n, v, obs.recBinning);
    if (i == 0) return {};
 
    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

fakes

std::vector<RecJet> fakes

genJets

std::vector<GenJet>* genJets

gEv

GenEvent* gEv

matched

std::optional<bool> matched

matches

std::vector<std::pair<GenJet,RecJet> > matches

misses

std::vector<GenJet> misses

obs

HTn 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/HTn.h
• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Unfolding/src/HTn.cc