BalanceHists< Jet >

Description

template<typename Jet>
struct DAS::JetEnergy::BalanceHists< Jet >

Dijet balance histograms for all variations.

Collaboration diagram for BalanceHists< Jet >:

Public Member Functions
	BalanceHists (const char *name, const vector< TString > vars)
void	operator() (const Jet &tag, const Jet &probe, const double ev_w)
void	Write (TDirectory *f)

Public Attributes
const TString	n
vector< TH3 * >	hs

Constructor & Destructor Documentation

BalanceHists()

BalanceHists ( const char *  name, const vector< TString >  vars  )

inline

Parameters

name	general name (will be used for the directory)
vars	variation names

                   : n(name)
    {
        int nBins = 60;
        auto balBins = getBinning(nBins, -0.6, 0.6);
        assertValidBinning(balBins);
 
        for (const auto& var: vars) {
            TString name = n + '_' + var;
            hs.push_back(new TH3D(name, var, nBins, balBins.data(),
                                             nPtBins, pt_edges.data(),
                                             nYbins, y_edges.data()));
        }
    }

Member Function Documentation

operator()()

void operator() ( const Jet &  tag, const Jet &  probe, const double  ev_w  )

inline

    {
        if (std::abs(tag.p4.Eta()) > 1.3) return; // barrel acceptance
 
        auto absy_probe = probe.AbsRap();
        for (size_t i = 0; i < hs.size(); ++i ) {
 
            double pt_tag, pt_probe;
            if constexpr (is_same<Jet,RecJet>()) {
                pt_tag   = tag  .CorrPt(i);
                pt_probe = probe.CorrPt(i);
            }
            else {
                pt_tag   = tag  .p4.Pt();
                pt_probe = probe.p4.Pt();
            }
 
            if (pt_probe < 28) continue;
 
            auto ptMean = (pt_tag+pt_probe)/2.0;
            auto bal    = (pt_probe-pt_tag)/ptMean;
 
            auto jw = tag  .weights.front()
                    * probe.weights.front();
            hs[i]->Fill(bal, ptMean, absy_probe, ev_w*jw);
        } // end of loop on variations
    }

Write()

void Write ( TDirectory *  f )

inline

    {
        f->cd();
        TDirectory * d = f->mkdir(n); // member of the class
        d->cd();
        for (TH3* h: hs) {
            h->SetDirectory(d);
            TString hname = h->GetName();
            hname.ReplaceAll(n + '_', ""); // n is a member of the class
            h->Write(hname);
        }
    }

Member Data Documentation

hs

vector<TH3*> hs

variations

n

const TString n

name to indicate the selection

---

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

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/JEC/bin/getPtBalance.cc