IDApplier

Description

Guesses the selection from the histogram name in the maps, and applies it by setting the weights to 0 to muons that don't pass the selection.

Note: tuned for Run 2 UL

Inheritance diagram for IDApplier:

Collaboration diagram for IDApplier:

Public Member Functions
	IDApplier (const fs::path &filePath, TString histname, bool correction, bool uncertainties)
Public Member Functions inherited from GenericSFApplier< DAS::RecMuon >
	GenericSFApplier (const std::filesystem::path &filePath, bool correction, bool uncertainties)
virtual	~GenericSFApplier () noexcept=default
void	operator() (DAS::RecMuon &object, const Context &... ctx) const
void	operator() (std::vector< DAS::RecMuon > &objects, const Context &... ctx) const
std::vector< std::string >	weightNames () const

Public Attributes
TString	name
TString	NUM
TString	DEN

Protected Member Functions
bool	passes (const RecMuon &muon) const override
int	binIndex (const RecMuon &muon, const std::unique_ptr< TH1 > &hist) const override
Protected Member Functions inherited from GenericSFApplier< DAS::RecMuon >
void	loadNominal (const std::string &histPath)
void	loadBinWiseUnc (const std::string &name, const std::string &histPath, Interpretation intp=UseBinContent)
void	loadGlobalUnc (const std::string &name, const std::string &histPath, Interpretation intp=UseBinContent)
virtual bool	passes (const DAS::RecMuon &obj, const Context &... ctx) const=0
virtual int	binIndex (const DAS::RecMuon &obj, const Context &... ctx, const std::unique_ptr< TH1 > &hist) const=0

Static Private Member Functions
static Selector	getMask (const string &s)

Private Attributes
const bool	abseta
const bool	charge
const bool	IPCut
Selector	num
Selector	den

Additional Inherited Members
Protected Types inherited from GenericSFApplier< DAS::RecMuon >
enum	Interpretation

Constructor & Destructor Documentation

IDApplier()

IDApplier ( const fs::path &  filePath, TString  histname, bool  correction, bool  uncertainties  )

inline

        : GenericSFApplier(filePath, correction, uncertainties)
        , abseta(!histname.Contains("_eta_pt" )) // default is to take $|\eta|$
        , charge( histname.Contains("_charge" ))
        , IPCut ( histname.Contains("andIPCut"))
        , name  ( histname                     )
    {
        loadNominal(histname.Data());
        auto var = histname + "_stat";
        loadBinWiseUnc(var.Data(), var.Data(), UseBinError);
        var = histname + "_syst";
        loadGlobalUnc(var.Data(), var.Data(), UseBinError);
 
        histname.ReplaceAll("_eta_pt", "");
        histname.ReplaceAll("_abseta_pt", "");
        histname.ReplaceAll("_charge", "");
        histname.ReplaceAll("andIPCut", "");
 
        DEN = histname(TRegexp("_DEN_.*"));
        histname.ReplaceAll(DEN, "");
        DEN.ReplaceAll("_DEN_", "");
        den = getMask(DEN.Data());
 
        NUM = histname;
        NUM.ReplaceAll("NUM_", "");
        num = getMask(NUM.Data());
 
        cout << "Initialisation of " << __func__ << " done" << endl;
    }

Member Function Documentation

binIndex()

int binIndex ( const RecMuon &  muon, const std::unique_ptr< TH1 > &  hist  ) const

inlineoverrideprotected

    {
        const float eta = abseta ? fabsf(muon.p4.Eta()) : muon.p4.Eta();
        if (charge) {
            // 3D: Q eta pT
            const float maxpt = hist->GetZaxis()->GetXmax() - 1;
            const float pt = min(muon.p4.Pt(), maxpt);
            return hist->FindBin(muon.Q, eta, pt);
        } else {
            // 2D: eta pT
            const float maxpt = hist->GetYaxis()->GetXmax() - 1;
            const float pt = min(muon.p4.Pt(), maxpt);
            return hist->FindBin(eta, pt);
        }
    }

getMask()

static Selector getMask ( const string &  s )

inlinestaticprivate

Translate table name from Muon POG ROOT file into selector.

    {
        // RECO
        if (s == "TrackerMuons") return Selector::None;
        if (s == "genTracks"   ) return Selector::None;
 
        // Muon ID
        if (s ==        "LooseID") return Selector::CutBasedIdLoose       ;
        if (s ==       "MediumID") return Selector::CutBasedIdMedium      ;
        if (s == "MediumPromptID") return Selector::CutBasedIdMediumPrompt;
        if (s ==         "SoftID") return Selector::SoftCutBasedId        ;
        if (s ==        "TightID") return Selector::CutBasedIdTight       ; // andIPCut?
        if (s ==       "HighPtID") return Selector::CutBasedIdGlobalHighPt; // andIPCut?
        if (s ==    "TrkHighPtID") return Selector::CutBasedIdTrkHighPt   ; // andIPCut?
 
        // Muon ISO
        if (s == "VeryLooseRelIso") return Selector::PFIsoVeryLoose;
        if (s ==     "LooseRelIso") return Selector::PFIsoLoose    ;
        if (s ==    "MediumRelIso") return Selector::PFIsoMedium   ;
        if (s ==     "TightRelIso") return Selector::PFIsoTight    ;
        if (s == "VeryTightRelIso") return Selector::PFIsoVeryTight;
        if (s ==   "LooseRelTkIso") return Selector::TkIsoLoose    ;
        if (s ==   "TightRelTkIso") return Selector::TkIsoTight    ;
 
        BOOST_THROW_EXCEPTION( invalid_argument(s + " is not recognised.") );
    }

passes()

bool passes ( const RecMuon &  muon ) const

inlineoverrideprotected

Todo:

IPCut? (unless it's already included in IPCut?)

    {
        return num == 0 // Non standard, tracker SF
            || (num & muon.selectors) == num;
    }

Member Data Documentation

abseta

const bool abseta

private

charge

const bool charge

private

den

Selector den

private

DEN

TString DEN

IPCut

const bool IPCut

private

name

TString name

num

Selector num

private

NUM

TString NUM

---

The documentation for this class was generated from the following file:

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Muons/bin/applyMuonSelection.cc