JetSubStructureUtils::ZCut Class Reference

#include <ZCut.h>

Inheritance diagram for JetSubStructureUtils::ZCut:

Collaboration diagram for JetSubStructureUtils::ZCut:

Public Member Functions
	ZCut (unsigned int nSubJets)
virtual double	result (const fastjet::PseudoJet &jet) const
virtual TOut	result (const xAOD::Jet &jet) const

Private Attributes
unsigned int	m_nSubJets

Detailed Description

Definition at line 11 of file ZCut.h.

Constructor & Destructor Documentation

ZCut()

JetSubStructureUtils::ZCut::ZCut ( unsigned int nSubJets )

inline

Definition at line 13 of file ZCut.h.

: m_nSubJets(nSubJets) {};

Member Function Documentation

result() [1/2]

double ZCut::result ( const fastjet::PseudoJet & jet ) const

virtual

Definition at line 13 of file ZCut.cxx.

{
  vector<fastjet::PseudoJet> constit_pseudojets = jet.constituents();
  if(constit_pseudojets.empty()) return -1;
 
  fastjet::JetDefinition jet_def = fastjet::JetDefinition(fastjet::kt_algorithm, 1.5,
                                                          fastjet::E_scheme, fastjet::Best);
  fastjet::ClusterSequence kt_clust_seq(constit_pseudojets, jet_def);
 
  if(constit_pseudojets.size() < m_nSubJets) {
    //We were asked to calculate zCut, but there are not enough constituents
    return 0.0;
  }
  vector<fastjet::PseudoJet> subjets = kt_clust_seq.exclusive_jets((int)m_nSubJets);
 
  // Find last split jet (cluster_hist_index should be highest for the last created jet)
  fastjet::PseudoJet *lastSplitSubjet = nullptr;
  int max_cluster_hist_index = -1;
  for(size_t iSubjet=0; iSubjet < subjets.size(); iSubjet++) {
    fastjet::PseudoJet parent1, parent2;
    if(kt_clust_seq.has_parents(subjets[iSubjet], parent1, parent2) &&
       subjets[iSubjet].cluster_hist_index() > max_cluster_hist_index) {
      max_cluster_hist_index = subjets[iSubjet].cluster_hist_index();
      lastSplitSubjet = &subjets[iSubjet];
    }
  }
 
  if(lastSplitSubjet == nullptr) {
    //None of the subjets were split
    return 0.0;
  }
 
  KtSplittingScale split(m_nSubJets);
  double dmin = pow(split.result(jet), 2.0);
 
  double zcut = -1;
  if(dmin == 0)
    zcut = 0;
  else
    zcut = dmin / (dmin + lastSplitSubjet->m2());
 
  return zcut;
}

result() [2/2]

virtual TOut JetSubStructureUtils::SubstructureCalculator< TOut >::result ( const xAOD::Jet & jet ) const

inlinevirtual

Reimplemented from JetSubStructureUtils::SubstructureCalculator< double >.

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_nSubJets

unsigned int JetSubStructureUtils::ZCut::m_nSubJets

private

Definition at line 19 of file ZCut.h.

---

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

• ZCut.h
• ZCut.cxx