JetSubStructureUtils::SubjetFinder Class Reference

#include <SubjetFinder.h>

Inheritance diagram for JetSubStructureUtils::SubjetFinder:

Collaboration diagram for JetSubStructureUtils::SubjetFinder:

Public Member Functions
	SubjetFinder (fastjet::JetAlgorithm fj_jetalg=fastjet::kt_algorithm, float jet_radius=0.3, float pt_min=5000, int exclusive_njets=-1)
virtual std::vector< fastjet::PseudoJet >	result (const fastjet::PseudoJet &jet) const
virtual TOut	result (const xAOD::Jet &jet) const
virtual std::vector< fastjet::PseudoJet >	result (const xAOD::Jet &jet) const

Private Attributes
fastjet::JetAlgorithm	m_fj_jetalg
float	m_jetrad
float	m_ptmin
int	m_exclusivenjets

Detailed Description

Definition at line 13 of file SubjetFinder.h.

Constructor & Destructor Documentation

SubjetFinder()

SubjetFinder::SubjetFinder	(	fastjet::JetAlgorithm	fj_jetalg = fastjet::kt_algorithm,
		float	jet_radius = 0.3,
		float	pt_min = 5000,
		int	exclusive_njets = -1
	)

Definition at line 15 of file SubjetFinder.cxx.

                                                                                                             :
  m_fj_jetalg(fj_jetalg), m_jetrad(jet_radius), m_ptmin(pt_min), m_exclusivenjets(exclusive_njets)
{
}

Member Function Documentation

result() [1/3]

vector< fastjet::PseudoJet > SubjetFinder::result ( const fastjet::PseudoJet &  jet ) const

virtual

Definition at line 20 of file SubjetFinder.cxx.

{
  vector<fastjet::PseudoJet> constit_pseudojets = jet.constituents();
  vector<fastjet::PseudoJet> subjets;
  if(constit_pseudojets.empty()) { 
    cout << "Warning in SubjetFinder: jet has no constituents" << endl;
    return subjets;
  }
 
  fastjet::ClusterSequence *clust_seq = nullptr;
 
  if (m_fj_jetalg==fastjet::ee_kt_algorithm) {
    fastjet::JetDefinition jet_def = fastjet::JetDefinition(fastjet::ee_kt_algorithm);
    clust_seq =  new fastjet::ClusterSequence(constit_pseudojets, jet_def);
  }else if (m_fj_jetalg==fastjet::plugin_algorithm) {
    fastjet::JetDefinition jet_def = fastjet::JetDefinition(new fastjet::EECambridgePlugin(m_jetrad));
    clust_seq =  new fastjet::ClusterSequence(constit_pseudojets, jet_def);
  }else{
    fastjet::JetDefinition jet_def = fastjet::JetDefinition(m_fj_jetalg, m_jetrad, fastjet::E_scheme, fastjet::Best);
    clust_seq =  new fastjet::ClusterSequence(constit_pseudojets, jet_def);
  }
 
  if(m_exclusivenjets < 0) { // Inclusive
    subjets = fastjet::sorted_by_pt(clust_seq->inclusive_jets(m_ptmin));
  }
  else {
    subjets = fastjet::sorted_by_pt(clust_seq->exclusive_jets_up_to(m_exclusivenjets));
  }
 
  if(subjets.empty()) {
    delete clust_seq;
  }
  else {
    clust_seq->delete_self_when_unused();
  }
 
  return subjets;
}

result() [2/3]

virtual TOut JetSubStructureUtils::SubstructureCalculator< TOut >::result

inline

Definition at line 25 of file SubstructureCalculator.h.

                                                    {
        // PS 4/18 master developent
        // std::vector<fastjet::PseudoJet> constit_pseudojets =
        //   jet::JetConstituentFiller::constituentPseudoJets(jet);
 
        std::vector<fastjet::PseudoJet> constit_pseudojets;
        std::transform(jet.getConstituents().begin(),
                       jet.getConstituents().end(),
                       std::back_inserter(constit_pseudojets),
                       [](const auto& c){
                         const xAOD::IParticle* ip = c->rawConstituent();
                         return
                           // fastjet::PseudoJet((c->rawConstituent())->p4());
                           fastjet::PseudoJet(ip->p4());
                       });
 
        fastjet::PseudoJet pjet = fastjet::join(constit_pseudojets);
 
        return result(pjet);
      }

result() [3/3]

virtual std::vector< fastjet::PseudoJet > JetSubStructureUtils::SubstructureCalculator< std::vector< fastjet::PseudoJet > >::result ( const xAOD::Jet &  jet ) const

inlinevirtualinherited

Definition at line 25 of file SubstructureCalculator.h.

                                                    {
        // PS 4/18 master developent
        // std::vector<fastjet::PseudoJet> constit_pseudojets =
        //   jet::JetConstituentFiller::constituentPseudoJets(jet);
 
        std::vector<fastjet::PseudoJet> constit_pseudojets;
        std::transform(jet.getConstituents().begin(),
                       jet.getConstituents().end(),
                       std::back_inserter(constit_pseudojets),
                       [](const auto& c){
                         const xAOD::IParticle* ip = c->rawConstituent();
                         return
                           // fastjet::PseudoJet((c->rawConstituent())->p4());
                           fastjet::PseudoJet(ip->p4());
                       });
 
        fastjet::PseudoJet pjet = fastjet::join(constit_pseudojets);
 
        return result(pjet);
      }

Member Data Documentation

m_exclusivenjets

int JetSubStructureUtils::SubjetFinder::m_exclusivenjets

private

Definition at line 24 of file SubjetFinder.h.

m_fj_jetalg

fastjet::JetAlgorithm JetSubStructureUtils::SubjetFinder::m_fj_jetalg

private

Definition at line 21 of file SubjetFinder.h.

m_jetrad

float JetSubStructureUtils::SubjetFinder::m_jetrad

private

Definition at line 22 of file SubjetFinder.h.

m_ptmin

float JetSubStructureUtils::SubjetFinder::m_ptmin

private

Definition at line 23 of file SubjetFinder.h.

---

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

• SubjetFinder.h
• SubjetFinder.cxx