MNjets.cc File Reference

#include "Core/Unfolding/interface/MNjets.h"
#include "Core/JetObservables/interface/MuellerNavelet.h"
#include "Math/VectorUtil.h"
#include <exceptions.h>

Include dependency graph for MNjets.cc:

Functions
template<typename Jet >
Di< const Jet, const Jet >	selection (const TTreeReaderArray< Jet > &jets, const Uncertainties::Variation &v=Uncertainties::nominal)
template<typename Jet >
double	getBinNumber (const Di< Jet, Jet > &MNjets, const Uncertainties::Variation &v, TUnfoldBinning *bng)

Function Documentation

getBinNumber()

double getBinNumber	(	const Di< Jet, Jet > &	MNjets,
		const Uncertainties::Variation &	v,
		TUnfoldBinning *	bng
	)

Find the bin number.

Returns

bin index (0 means out of PS, like TUnfoldBinning::GetGlobalBinNumber)

{
    if (!MNjets) return 0;
 
    float DEta = std::abs(MNjets.DeltaEta(v)),
          DPhi = std::abs(MNjets.DeltaPhi(v));
 
    return bng->GetGlobalBinNumber(DEta, DPhi);
}

selection()

Di<const Jet, const Jet> selection	(	const TTreeReaderArray< Jet > &	jets,
		const Uncertainties::Variation &	v = Uncertainties::nominal
	)

Mueller-Navelet jet selection.

{
    Di<const Jet, const Jet> MNjets;
    for (size_t i = 0; i < jets.GetSize(); ++i) {
        const Jet& jet = jets.At(i);
        if (jet.CorrPt(v) < minpt) continue;
        if (jet.AbsRap() > maxy) continue;
 
        if (jet.p4.Eta() > (MNjets.first ? MNjets.first->p4.Eta() : 0.))
            MNjets.first = &jet;
        else if (jet.p4.Eta() < (MNjets.second ? MNjets.second->p4.Eta() : 0.))
            MNjets.second = &jet;
    }
 
    return MNjets;
}