DEtaDPhiFiller

#include <MNjets.h>

Inheritance diagram for DEtaDPhiFiller:

Collaboration diagram for DEtaDPhiFiller:

Public Member Functions
	DEtaDPhiFiller (const DEtaDPhi &obs, TTreeReader &reader)
std::list< int >	fillRec (DistVariation &) override
void	match () override
void	fillMC (DistVariation &) override
Public Member Functions inherited from Filler
virtual	~Filler ()=default

Public Attributes
DEtaDPhi	obs
std::optional< TTreeReaderArray< GenJet > >	genJets
TTreeReaderArray< RecJet >	recJets
std::optional< TTreeReaderValue< GenEvent > >	gEv
TTreeReaderValue< RecEvent >	rEv
Di< const GenJet, const GenJet >	genMNjets
Di< const RecJet, const RecJet >	recMNjets
std::optional< bool >	matched

Constructor & Destructor Documentation

DEtaDPhiFiller()

DEtaDPhiFiller	(	const DEtaDPhi &	obs,
		TTreeReader &	reader
	)

Constructor.

    : obs(obs)
    , genJets(initOptionalBranch<decltype(genJets)>(reader, "genJets"))
    , recJets(reader, "recJets")
    , gEv(initOptionalBranch<decltype(gEv)>(reader, "genEvent"))
    , rEv(reader, "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 irec = getBinNumber(recMNjets, v, obs.recBinning);
    double recW = rEv->Weight(v);
    if (irec > 0) recW *= recMNjets.Weight(v);
 
    auto igen = getBinNumber(genMNjets, v, obs.genBinning);
    double genW = (*gEv)->Weight(v);
    if (igen > 0) genW *= genMNjets.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.

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

match()

void match ( )

overridevirtual

Check if the MN jets can be matched at both levels.

Todo:

is optional necessary? (taken from MjjYmax...)

Reimplemented from Filler.

{
    matched.reset(); 
 
    using ROOT::Math::VectorUtil::DeltaR;
 
    genMNjets = selection(*genJets);
 
    // matching (swapping leading and subleading is allowed)
    matched = genMNjets && recMNjets
          && DeltaR(genMNjets.first ->p4, recMNjets.first ->p4) < obs.maxDR
          && DeltaR(genMNjets.second->p4, recMNjets.second->p4) < obs.maxDR;
}

Member Data Documentation

genJets

std::optional<TTreeReaderArray<GenJet> > genJets

genMNjets

Di<const GenJet, const GenJet> genMNjets

gEv

std::optional<TTreeReaderValue<GenEvent> > gEv

matched

std::optional<bool> matched

obs

DEtaDPhi obs

Backreference to the observable.

recJets

TTreeReaderArray<RecJet> recJets

recMNjets

Di<const RecJet, const RecJet> recMNjets

rEv

TTreeReaderValue<RecEvent> rEv

---

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

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