Hist

Description

Contains the histograms that are filled in getMNobservables. These histograms are filled through the different Fill() methods. Beware that the three Fill() methods are different.

The attributes of this structure are the following:

• name --> name of the instance to avoid memory leak
• njets --> histogram of the number of jets in the event
• dsigmadphi --> $(\pi - \Delta \phi)$ between the MN Jets
• ptleading --> $p_T$ of the leading jet (in term of pt)
• ptratmn --> ratio of pt of the mnjets (see the Obs2Jets)
• deta --> Delta Eta distribution between the MN Jets
• pt --> pt of all the jets in the event
• ptmnfwd --> pt of the most fwd MN jet
• ptmnbwd --> pt of the most bwd MN jet
• njetsptave --> 2D histogram jets multiplicity and average pt of the MN jets
• njetsdeta --> 2D histogram jets multiplicity and delat eta btw the MN jets
• ptavedeta --> 2D histogram binned in $p_{T,MN,ave}$ and $ \Delta \eta$
• dphideta --> 2D histogram binned in $\Delta\phi$ and $\Delta\eta$
• cos1,2,3 --> TProfile containing $ \left\langle \cos \left( n \left( \pi - \Delta \phi \right) \right) \right\rangle $ for $n=1,2,3$
• ptavemini --> See Obs2Jets
• detaavemini --> See Obs2Jets
• retaavemini --> See Obs2Jets
• rptexpdeta --> See Obs2Jets

Public Member Functions
	Hist (TString thename)
void	Fill (vector< RecJet > &recJets, double evweight)
void	Fill (vector< GenJet > &genJets, double evweight)
void	Fill (Obs2Jets &mn2Jets, ObsMiniJets &miniJets, double evweight)
void	Write (TDirectory *dir)

Public Attributes
TString	name
TH1F	njets
TH1F	dsigmadphi
TH1F	ptleading
TH1F	ptratmn
TH1F	pt
TH1F	ptmnfwd
TH1F	ptmnbwd
TH1F	deta
TH2F	njetsptave
TH2F	njetsdeta
TH2F	ptavedeta
TH2F	dphideta
TH2F	pteta
TProfile	cos1
TProfile	cos2
TProfile	cos3
TH1F	ptavemini
TH1F	detaavemini
TH1F	retaavemini
TH1F	rptexpdeta

Private Member Functions
template<typename HIST >
void	Loop (vector< HIST * > hists, TDirectory *dir)

Constructor & Destructor Documentation

Hist()

Hist ( TString  thename )

inline

                           :
        name(thename),
        njets(name+"njets",";N_{jets};N_{evts}",20,0,20),
        dsigmadphi(name+"dsigmadphi",";\\Delta \\phi;N_{jets}",binsPhi.size()-1, binsPhi.data()),
        ptleading(name+"ptlead","pt of the leading jet;p_{T} (GeV);N_{jets}", pt_edges.size()-1, pt_edges.data()),
        ptratmn(name+"ptratmn","Pt Balance;\\frac{p_{T,1}}{p_{T,2}};N_{jets}",30,0,1),
        pt(name+"pt","INclusive jets;p_{T} (GeV);N_{jets}",pt_edges.size()-1,pt_edges.data()),
        ptmnfwd(name+"ptmnfwd","Pt Backward MN jet;p_{T}^{fwd} (GeV);N_{jets}",pt_edges.size()-1,pt_edges.data()),
        ptmnbwd(name+"ptmnbwd","Pt Foward MN jet;p_{T}^{fwd} (GeV);N_{jets}",pt_edges.size()-1,pt_edges.data()),
        deta(name+"delta_eta",";\\Delta \\eta;\\frac{d \\sigma}{d p_{T} d \\eta};", binsY.size()-1, binsY.data()),
        njetsptave(name+"njetsptave","Jet multiplicity;p_{T}^{ave};N_{jets}", pt_edges.size()-1, pt_edges.data(),20,0,20),
        njetsdeta(name+"njetsdeta","Jet multiplicity;\\Delta \\eta;N_{jets}", binsY.size()-1, binsY.data(),10 ,0 ,10),
        ptavedeta(name+"ptavedeta","ptavedeta",pt_edges.size()-1, pt_edges.data(),binsY.size()-1, binsY.data()),
        dphideta(name+"dphideta","dphideta",30, 0, M_PI, binsY.size()-1, binsY.data()),
        pteta(name+"pteta",";p_{T};\\eta", pt_edges.size()-1, pt_edges.data(), y_edges.size()-1, y_edges.data()),
        cos1(name+"cos1","cos1",30,0,10,0,1),
        cos2(name+"cos2","cos2",30,0,10,0,1),
        cos3(name+"cos3","cos3",30,0,10,0,1),
        ptavemini(name+"ptavemini","ptavemini",50,0.,80.),
        detaavemini(name+"detaavemini","detaavemini",50,0.,5.),
        retaavemini(name+"retaavemini","retaavemini",50,0.,1.),
        rptexpdeta(name+"rptexpdeta","rptexpdeta",50,0.,2.)
    { }

Member Function Documentation

Fill() [1/3]

void Fill ( Obs2Jets &  mn2Jets, ObsMiniJets &  miniJets, double  evweight  )

inline

                                                                         {
        double w1 = evweight * mn2Jets.weight;
        double w2 = evweight * miniJets.weight;
 
        dsigmadphi.Fill( mn2Jets.DPhi() , w1 );
        deta.Fill( mn2Jets.DEta(), w1 );
        njetsptave.Fill( mn2Jets.PtAve(), miniJets.size()+2, w1 );
        njetsdeta.Fill( mn2Jets.DEta(), miniJets.size()+2, w1 );
        ptratmn.Fill( mn2Jets.PtRatMN(), w1 );
        ptavedeta.Fill( mn2Jets.PtAve(), mn2Jets.DEta(), w1 );
        dphideta.Fill( mn2Jets.DPhi(), mn2Jets.DEta(), w1 );
 
        ptmnfwd.Fill( mn2Jets.PtFWD(), evweight * mn2Jets.weightFWD );
        ptmnbwd.Fill( mn2Jets.PtBWD(), evweight * mn2Jets.weightBWD );
 
        cos1.Fill( mn2Jets.DEta(), mn2Jets.CosDPhi(1), w1 );
        cos2.Fill( mn2Jets.DEta(), mn2Jets.CosDPhi(2), w1 );
        cos3.Fill( mn2Jets.DEta(), mn2Jets.CosDPhi(3), w1 );
 
        if (miniJets.size()>0) ptavemini.Fill(miniJets.PtAveMini(), w2 );
        if (miniJets.size()>1) {
            detaavemini.Fill(miniJets.DEtaAveMini(), w2 );
            retaavemini.Fill(miniJets.REtaAveMini(), w2 );
            rptexpdeta.Fill(miniJets.RPtExpDEta(), w2 );
        }
    }

Fill() [2/3]

void Fill ( vector< GenJet > &  genJets, double  evweight  )

inline

                                                        {
        double jetswgt = (genJets.front()).weights.front();
        for (auto it = next(genJets.begin()) ; it != genJets.end() ; ++it)
            jetswgt *= (it->weights).front();
        for (auto it = genJets.begin() ; it != genJets.end() ; ++it) {
            pt.Fill( (it->p4).Pt(), evweight * (it->weights).front() );
            pteta.Fill( (it->p4).Pt(), (it->p4).Eta(), evweight * (it->weights).front());
        }
        njets.Fill(genJets.size(), evweight * jetswgt);
        ptleading.Fill( (genJets.front()).p4.Pt(), evweight * (genJets.front()).weights.front() );
    }

Fill() [3/3]

void Fill ( vector< RecJet > &  recJets, double  evweight  )

inline

                                                        {
        double jetswgt = (recJets.front()).weights.front();
        for (auto it = next(recJets.begin()) ; it != recJets.end() ; ++it)
            jetswgt *= (it->weights).front();
        for (auto it = recJets.begin() ; it != recJets.end() ; ++it) {
            pt.Fill(it->CorrPt(), evweight * (it->weights).front());
            pteta.Fill(it->CorrPt(), it->p4.Eta(), evweight * (it->weights).front());
        }
        njets.Fill( recJets.size(), evweight*jetswgt );
        ptleading.Fill( (recJets.front()).CorrPt(), evweight * (recJets.front()).weights.front() );
    }

Loop()

void Loop ( vector< HIST * >  hists, TDirectory *  dir  )

inlineprivate

    {
        for (HIST * h: hists){
            h->SetDirectory(dir);
            TString histoname(h->GetName());
            histoname.ReplaceAll(name,"");
            h->Write(histoname);
        }
    }

Write()

void Write ( TDirectory *  dir )

inline

    {
        dir->cd();
        auto subdir = dir->mkdir(name);
        subdir->cd();
        Loop<TH1F>({&njets, &dsigmadphi, &ptleading, &ptratmn, &ptmnfwd, &ptmnbwd, &pt,
                    &ptavemini, &detaavemini, &retaavemini, &rptexpdeta, &deta}, subdir);
        Loop<TH2F>({&njetsptave, &njetsdeta, &ptavedeta, &dphideta, &pteta}, subdir);
        Loop<TProfile>({&cos1, &cos2, &cos3}, subdir);
    }

Member Data Documentation

cos1

TProfile cos1

cos2

TProfile cos2

cos3

TProfile cos3

deta

TH1F deta

detaavemini

TH1F detaavemini

dphideta

TH2F dphideta

dsigmadphi

TH1F dsigmadphi

name

TString name

njets

TH1F njets

njetsdeta

TH2F njetsdeta

njetsptave

TH2F njetsptave

pt

TH1F pt

ptavedeta

TH2F ptavedeta

ptavemini

TH1F ptavemini

pteta

TH2F pteta

ptleading

TH1F ptleading

ptmnbwd

TH1F ptmnbwd

ptmnfwd

TH1F ptmnfwd

ptratmn

TH1F ptratmn

retaavemini

TH1F retaavemini

rptexpdeta

TH1F rptexpdeta

---

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

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/JetObservables/bin/getMNobservables.cc