Balance

Description

Balance histogram for a given variation.

#include <Balance.h>

Collaboration diagram for Balance:

Public Member Functions
	Balance (TDirectory *d, const DAS::Uncertainties::Variation &v, const vector< double > &pt_edges, const vector< double > &y_edges)
void	operator() (const DAS::AbstractPhysicsObject &tag, const DAS::AbstractPhysicsObject &probe, const DAS::MET *met, const double ev_w)
void	Write ()

Public Attributes
DAS::Uncertainties::Variation	v

Private Attributes
TDirectory *	d
unique_ptr< TH3D >	asymmetry
unique_ptr< TH2D >	response

Constructor & Destructor Documentation

Balance()

Balance ( TDirectory *  d, const DAS::Uncertainties::Variation &  v, const vector< double > &  pt_edges, const vector< double > &  y_edges  )

inline

                                            :
        d(d), v(v)
    {
        int nBins = 120;
        vector balBins = MakeBinning<double>({-1.2,1.2}, nBins);
 
        int nPtBins = pt_edges.size()-1,
            nYbins  =  y_edges.size()-1;
 
        TString dirname = d->GetName();
        TString hname = dirname + '_' + v.group + '_' + v.name;
        asymmetry = make_unique<TH3D>(hname, v.name, nBins, balBins.data(),
                                                     nPtBins, pt_edges.data(),
                                                     nYbins, y_edges.data());
 
        balBins = MakeBinning<double>({0., 1.2}, nBins);
        response = make_unique<TH2D>("response_" + hname, v.name, nBins, balBins.data(),
                                                                  nYbins, y_edges.data());
    }

Member Function Documentation

operator()()

void operator() ( const DAS::AbstractPhysicsObject &  tag, const DAS::AbstractPhysicsObject &  probe, const DAS::MET *  met, const double  ev_w  )

inline

Calculates balance for two objects following the definition of the Run 1 JetMET reference publication for dijet systems.

    {
        float absy_probe = probe.AbsRap(v);
 
        float pt_tag   = tag  .CorrPt(v),
              pt_probe = probe.CorrPt(v);
 
        double jw = tag  .Weight(v)
                  * probe.Weight(v);
 
        response->Fill(pt_probe/pt_tag, absy_probe, ev_w*jw);
 
        // \todo does this also apply to Z+jet balance?
        //       (e.g. not in 1804.05252)
        float pt_dif = (pt_tag - pt_probe),
              pt_ave = (pt_tag + pt_probe) * 0.5;
 
        float asym = pt_dif / pt_ave; // pT balance method
 
        if (met != nullptr) // MPF method ref:1607.03663 $\mathcal{A} =\frac{\vec{p}_T^{miss}\cdot\left(\vec{p}_{T,tag}/p_{T,tag}\right)}{2 p_{T,ave}}$
            asym = -met->Pt * cos(met->Phi - tag.CorrP4(v).Phi()) / pt_ave;
 
        asymmetry->Fill(asym, pt_ave, absy_probe, ev_w*jw);
    }

Write()

void Write ( )

inline

    {
        d->cd();
        TString dirname = d->GetName();
 
        asymmetry->SetDirectory(d);
        TString hname = asymmetry->GetName();
        hname.ReplaceAll(dirname + '_', "");
        asymmetry->Write(hname);
 
        response->SetDirectory(d);
        hname = response->GetName();
        hname.ReplaceAll(dirname + '_', "");
        response->Write(hname);
    }

Member Data Documentation

asymmetry

unique_ptr<TH3D> asymmetry

private

asymmetry histogram

d

TDirectory* d

private

directory in output ROOT file

response

unique_ptr<TH2D> response

private

response

v

DAS::Uncertainties::Variation v

systematic variation

---

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

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/JEC/bin/Balance.h