ObsMiniJets

Description

MN_Helper::ObsMiniJets is a structure that implements observables based on the jets that are emmited in-between the Mueller-Navelet jets. These jets are called "mini jets" Definition of these observavbles ($N$ is the number of jets in the event, and the jets are labeled $i \in (0,N-1)$, including the two Mueller-Navelet jets, the jets are assumed ordred in rapidity from highest rapidity to lowest rapidity (signed rapidity).

• PtAveMini --> $ \left\langle P_{T,Mini} \right\rangle = \frac{1}{N-2} \sum_{i=1}^{N-2} p_{T,i} $
• DEtaAveMini --> & \left\langle \Delta \eta_{Mini} \right\rangle = \frac{1}{N-3} \sum_{i=1}^{N-3} \left\vert \eta_{i} - \eta_{i+1} \right\vert$
• REtaAveMini --> $ R_{k} = \frac{1}{N-3} \sum_{i=1}^{N-3} \frac{\eta_{i}}{\eta_{i+1}} $

  Todo:

  Check this and compare with the standard definition.

#include <MuellerNavelet.h>

Collaboration diagram for ObsMiniJets:

Public Member Functions
	ObsMiniJets (const std::vector< DAS::RecJet > &recJets)
	ObsMiniJets (const std::vector< DAS::GenJet > &genJets)
float	PtAveMini () const
float	DEtaAveMini () const
float	REtaAveMini () const
float	RPtExpDEta () const
size_t	size () const

Public Attributes
const float	weight

Private Attributes
std::vector< DAS::FourVector >	MiniJets

Constructor & Destructor Documentation

ObsMiniJets() [1/2]

ObsMiniJets

(

const std::vector< DAS::RecJet > &

recJets

)

                                                          :
     weight( (recJets.size()<1) ? 0 :
             accumulate( recJets.begin(), recJets.end(), 1., 
                [] ( double w, const RecJet& j ) { return w * j.weights.front(); } ) )
{
    for (auto it = recJets.begin() ; it != recJets.end() ; ++it) {
 
        FourVector p4 = it->p4;
        p4 *= it->scales.front();
        MiniJets.push_back( p4 );
    }
}

ObsMiniJets() [2/2]

ObsMiniJets

(

const std::vector< DAS::GenJet > &

genJets

)

                                                          :
    weight( (genJets.size() < 1) ? 0 :
            accumulate( genJets.begin(), genJets.end(), 1., 
                [] ( double w, const GenJet& j ) { return w * j.weights.front(); } ) )
{
    for (auto it = genJets.begin() ; it != genJets.end() ; ++it)
        MiniJets.push_back( it->p4 );
}

Member Function Documentation

DEtaAveMini()

float DEtaAveMini

(

)

const

{
    assert(MiniJets.size() >= 2);
    float detaavemini = 0;
    for (auto it=MiniJets.begin() ; it != prev(MiniJets.end()) ; ++it)
        detaavemini += it->Eta() - next(it)->Eta();
    return detaavemini/MiniJets.size();
}

PtAveMini()

float PtAveMini

(

)

const

{
    assert(MiniJets.size()>0);
    auto AddPt = [] ( float sumpt, const FourVector& v) { return sumpt + v.Pt(); };
    float ptavemini = accumulate(MiniJets.begin(), MiniJets.end(), 0., AddPt);
    return ptavemini/MiniJets.size();
}

REtaAveMini()

float REtaAveMini

(

)

const

{
    assert(MiniJets.size() >= 2); 
    float retaavemini = 0;
    for (auto it=MiniJets.begin() ; it != prev(MiniJets.end()) ; ++it)
        retaavemini += it->Eta()/next(it)->Eta();
    return retaavemini/MiniJets.size();
}

RPtExpDEta()

float RPtExpDEta

(

)

const

{
    assert(MiniJets.size() >= 2);
    float rptexpdeta = 0;
    for (auto it=MiniJets.begin() ; it != prev(MiniJets.end()) ; ++it) {
        auto it2 = next(it);
        auto Eta = it->Eta(), Eta2 = it2->Eta();
        auto expo = exp( Eta - Eta2 );
        auto Pt = it->Pt(), Pt2 = it2->Pt();
        assert( Pt2 > 0 );
        rptexpdeta += Pt / Pt2 * expo;
    }
    return rptexpdeta/MiniJets.size();
}

size()

size_t size

(

)

const

{ return MiniJets.size(); }

Member Data Documentation

MiniJets

std::vector<DAS::FourVector> MiniJets

private

The size should be >1.

The jets that are between the Mueller-navelet Jets.

weight

const float weight

---

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

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/JetObservables/interface/MuellerNavelet.h
• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/JetObservables/src/MuellerNavelet.cc