Gep::ConeJetMaker Class Reference

#include <ConeJetMaker.h>

Inheritance diagram for Gep::ConeJetMaker:

Collaboration diagram for Gep::ConeJetMaker:

Public Member Functions
	ConeJetMaker (float jetR, const xAOD::jFexSRJetRoIContainer &seeds, float seedEtThreshold=5.e3, const std::string &recombScheme="EScheme")
std::string	toString () const override
virtual std::vector< Gep::Jet >	makeJets (const std::vector< Gep::Cluster > &clusters) const override
float	getJetR () const
float	getSeedEtThreshold () const
std::string	recombSchemeAsString () const

Private Types
enum class	RecombScheme { EScheme , SeedScheme }

Private Member Functions
RecombScheme	string2RecombScheme (const std::string &) const

Private Attributes
float	m_jetR
const xAOD::jFexSRJetRoIContainer &	m_seeds
float	m_seedEtThreshold
std::vector< std::pair< std::string, RecombScheme > >	m_knownSchemes
RecombScheme	m_recombScheme {RecombScheme::EScheme}

Detailed Description

Definition at line 22 of file ConeJetMaker.h.

Member Enumeration Documentation

RecombScheme

enum class Gep::ConeJetMaker::RecombScheme

strongprivate

Enumerator
EScheme
SeedScheme

Definition at line 44 of file ConeJetMaker.h.

{EScheme, SeedScheme};

Constructor & Destructor Documentation

ConeJetMaker()

Gep::ConeJetMaker::ConeJetMaker	(	float	jetR,
		const xAOD::jFexSRJetRoIContainer &	seeds,
		float	seedEtThreshold = 5.e3,
		const std::string &	recombScheme = "EScheme" )

Definition at line 13 of file ConeJetMaker.cxx.

                                         :
  m_jetR{jetR},
  m_seeds{seeds},
  m_seedEtThreshold{seedEtThreshold},
  m_recombScheme{string2RecombScheme(scheme)}{
}

Member Function Documentation

getJetR()

float Gep::ConeJetMaker::getJetR ( ) const

inline

Definition at line 34 of file ConeJetMaker.h.

{ return m_jetR; }

getSeedEtThreshold()

float Gep::ConeJetMaker::getSeedEtThreshold ( ) const

inline

Definition at line 35 of file ConeJetMaker.h.

{ return m_seedEtThreshold; }

makeJets()

std::vector< Gep::Jet > Gep::ConeJetMaker::makeJets ( const std::vector< Gep::Cluster > & clusters ) const

overridevirtual

Implements Gep::IJetMaker.

Definition at line 62 of file ConeJetMaker.cxx.

{
  std::vector<Gep::Jet> jets;
  for (const auto seed: m_seeds) {
 
    float seedEt = seed->et();
 
    // skip seeds with Et below threshold
    if(seedEt < m_seedEtThreshold) continue; 
    float seedEta = seed->eta();
    float seedPhi = seed->phi();
 
    Gep::Jet jet;
 
    jet.radius = m_jetR;
    jet.seedEta = seedEta;
    jet.seedPhi = seedPhi;
    jet.seedEt = seedEt;
 
    TLorentzVector jetVec;
    float px{0}, py{0};
    int clusterIndex {0};
    
    //build jet with clusters within dR from seed
    for (const auto &cl: clusters) {
      float dR_seed_cl = deltaR(seedEta, cl.vec.Eta(), seedPhi, cl.vec.Phi());
 
      if (dR_seed_cl < m_jetR) {
        jetVec += cl.vec;
        px += cl.vec.Px();
        py += cl.vec.Py();
    
        jet.constituentsIndices.push_back(clusterIndex);
      }
      clusterIndex++;
    }
 
 
    // skip cone jets with 0 constituents
    if (jet.constituentsIndices.empty()) {continue;}
    
 
    // recombination scheme
    if (m_recombScheme == RecombScheme::EScheme) {
      // default option: add four-vectors of constituents
      jet.vec = jetVec;
    } else if (m_recombScheme == RecombScheme::SeedScheme) {
      // massless jet, correct pt, re-using seed (eta,phi)
      float m = 0;
      float pt = std::sqrt(px*px + py*py);
      jet.vec.SetPtEtaPhiM(pt, seedEta, seedPhi, m);
    } else {
      throw std::runtime_error("GEP ConeJetMaker recombinaton scheme error");
    }
    jets.emplace_back(std::move(jet));
  }
 
  return jets;
}

recombSchemeAsString()

std::string Gep::ConeJetMaker::recombSchemeAsString

(

)

const

Definition at line 42 of file ConeJetMaker.cxx.

                                            {
  auto schemeMatches =
    [&scheme=m_recombScheme](const std::pair<std::string, RecombScheme>& p){
      return (p.second) == scheme;};
  
  auto scheme_it = std::find_if(m_knownSchemes.cbegin(),
                m_knownSchemes.cend(),
                schemeMatches);
  
  if (scheme_it == m_knownSchemes.end()){
    throw std::runtime_error("GEP ConeJetMaker unknown recombinaton scheme " +
                 recombSchemeAsString());
  }
  
  return scheme_it->first;
}

string2RecombScheme()

Gep::ConeJetMaker::RecombScheme Gep::ConeJetMaker::string2RecombScheme ( const std::string & scheme ) const

private

Definition at line 25 of file ConeJetMaker.cxx.

                                                                  {
  auto schemeMatches =
    [&scheme](const std::pair<std::string, RecombScheme>& p){return (p.first) == scheme;};
  
  auto scheme_it = std::find_if(m_knownSchemes.cbegin(),
                m_knownSchemes.cend(),
                schemeMatches);
  
  if (scheme_it == m_knownSchemes.end()){
    throw std::runtime_error("GEP ConeJetMaker unknown recombinaton scheme " +
                 scheme);
  }
  
  return scheme_it->second;
}

toString()

std::string Gep::ConeJetMaker::toString ( ) const

overridevirtual

Implements Gep::IJetMaker.

Definition at line 123 of file ConeJetMaker.cxx.

                                          {
  std::stringstream ss;
  ss << "ConeJetMaker. seed thresholds: " << m_seedEtThreshold
     << " rad: " << m_jetR << " recombinaton scheme: " << recombSchemeAsString();
  return ss.str();
}

Member Data Documentation

m_jetR

float Gep::ConeJetMaker::m_jetR

private

Definition at line 40 of file ConeJetMaker.h.

m_knownSchemes

std::vector<std::pair<std::string, RecombScheme> > Gep::ConeJetMaker::m_knownSchemes

private

Initial value:

{
    {"EScheme", RecombScheme::EScheme},
      {"SeedScheme", RecombScheme::SeedScheme}
  }

Definition at line 46 of file ConeJetMaker.h.

                                                               {
    {"EScheme", RecombScheme::EScheme},
      {"SeedScheme", RecombScheme::SeedScheme}
  };

m_recombScheme

RecombScheme Gep::ConeJetMaker::m_recombScheme {RecombScheme::EScheme}

private

Definition at line 51 of file ConeJetMaker.h.

{RecombScheme::EScheme};

m_seedEtThreshold

float Gep::ConeJetMaker::m_seedEtThreshold

private

Definition at line 42 of file ConeJetMaker.h.

m_seeds

const xAOD::jFexSRJetRoIContainer& Gep::ConeJetMaker::m_seeds

private

Definition at line 41 of file ConeJetMaker.h.

---

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

• ConeJetMaker.h
• ConeJetMaker.cxx