heppy::ReclusterJets Class Reference

#include <ReclusterJets.h>

Public Types
typedef math::XYZTLorentzVector	LorentzVector

Public Member Functions
std::vector< LorentzVector >	getGrouping (double ptMin=0.0)
	get grouping (inclusive jets) More...
std::vector< LorentzVector >	getGroupingExclusive (int njets)
	get grouping (exclusive jets, until n are left) More...
std::vector< LorentzVector >	getGroupingExclusive (double dcut)
	get grouping (exclusive jets, up to cut dcut) More...
LorentzVector	getPruned (double zcut, double rcutFactor)
	get pruned 4-vector More...
LorentzVector	getPrunedSubjetExclusive (unsigned int isubjet, double zcut, double rcutFactor)
	get pruned 4-vector for a given subject (must be called after getGroupingExclusive) More...
LorentzVector	getPrunedSubjetInclusive (unsigned int isubjet, double zcut, double rcutFactor)
	get pruned 4-vector for a given subject (must be called after getGroupingInclusive) More...
	ReclusterJets (const std::vector< LorentzVector > &objects, double ktpower, double rparam)

Private Types
typedef std::shared_ptr< fastjet::ClusterSequence >	ClusterSequencePtr
	fastjet outputs More...

Private Member Functions
LorentzVector	getPruned (const fastjet::PseudoJet &jet, double zcut, double rcutFactor)
std::vector< LorentzVector >	makeP4s (const std::vector< fastjet::PseudoJet > &jets)

Private Attributes
std::vector< fastjet::PseudoJet >	exclusiveJets_
ClusterSequencePtr	fjClusterSeq_
std::vector< fastjet::PseudoJet >	fjInputs_
std::vector< fastjet::PseudoJet >	inclusiveJets_
double	ktpower_
double	rparam_

Detailed Description

Definition at line 19 of file ReclusterJets.h.

Member Typedef Documentation

ClusterSequencePtr

typedef std::shared_ptr<fastjet::ClusterSequence> heppy::ReclusterJets::ClusterSequencePtr

private

fastjet outputs

Definition at line 57 of file ReclusterJets.h.

LorentzVector

typedef math::XYZTLorentzVector heppy::ReclusterJets::LorentzVector

Definition at line 21 of file ReclusterJets.h.

Constructor & Destructor Documentation

ReclusterJets()

heppy::ReclusterJets::ReclusterJets	(	const std::vector< LorentzVector > &	objects,
		double	ktpower,
		double	rparam
	)

Definition at line 14 of file ReclusterJets.cc.

References Exception, fjClusterSeq_, fjInputs_, dqmiolumiharvest::j, ktpower_, EcalTangentSkim_cfg::o, and rparam_.

      : ktpower_(ktpower), rparam_(rparam) {
    // define jet inputs
    fjInputs_.clear();
    int index = 0;
    for (const LorentzVector &o : objects) {
      fastjet::PseudoJet j(o.Px(), o.Py(), o.Pz(), o.E());
      j.set_user_index(index);
      index++;  // in case we want to know which piece ended where
      fjInputs_.push_back(j);
    }

    // choose a jet definition
    fastjet::JetDefinition jet_def;

    // prepare jet def
    if (ktpower_ == 1.0) {
      jet_def = JetDefinition(kt_algorithm, rparam_);
    } else if (ktpower_ == 0.0) {
      jet_def = JetDefinition(cambridge_algorithm, rparam_);
    } else if (ktpower_ == -1.0) {
      jet_def = JetDefinition(antikt_algorithm, rparam_);
    } else {
      throw cms::Exception("InvalidArgument", "Unsupported ktpower value");
    }

    // print out some infos
    //  cout << "Clustering with " << jet_def.description() << endl;
    // define jet clustering sequence
    fjClusterSeq_ = std::make_shared<fastjet::ClusterSequence>(fjInputs_, jet_def);
  }

Member Function Documentation

getGrouping()

std::vector< math::XYZTLorentzVector > heppy::ReclusterJets::getGrouping

(

double

ptMin = 0.0

)

get grouping (inclusive jets)

Definition at line 55 of file ReclusterJets.cc.

References fjClusterSeq_, inclusiveJets_, makeP4s(), and ptMin.

                                                                          {
    // recluster jet
    inclusiveJets_ = fastjet::sorted_by_pt(fjClusterSeq_->inclusive_jets(ptMin));
    // return
    return makeP4s(inclusiveJets_);
  }

getGroupingExclusive() [1/2]

std::vector< math::XYZTLorentzVector > heppy::ReclusterJets::getGroupingExclusive

(

int

njets

)

get grouping (exclusive jets, until n are left)

Definition at line 69 of file ReclusterJets.cc.

References exclusiveJets_, fjClusterSeq_, makeP4s(), and BTaggingMonitoring_cff::njets.

                                                                                {
    // recluster jet
    exclusiveJets_ = fastjet::sorted_by_pt(fjClusterSeq_->exclusive_jets(njets));
    // return
    return makeP4s(exclusiveJets_);
  }

getGroupingExclusive() [2/2]

std::vector< math::XYZTLorentzVector > heppy::ReclusterJets::getGroupingExclusive

(

double

dcut

)

get grouping (exclusive jets, up to cut dcut)

Definition at line 62 of file ReclusterJets.cc.

References exclusiveJets_, fjClusterSeq_, and makeP4s().

                                                                                  {
    // recluster jet
    exclusiveJets_ = fastjet::sorted_by_pt(fjClusterSeq_->exclusive_jets(dcut));
    // return
    return makeP4s(exclusiveJets_);
  }

getPruned() [1/2]

math::XYZTLorentzVector heppy::ReclusterJets::getPruned	(	double	zcut,
		double	rcutFactor
	)

get pruned 4-vector

Definition at line 76 of file ReclusterJets.cc.

References exclusiveJets_, fjClusterSeq_, getPrunedSubjetExclusive(), and HLT_2023v12_cff::zcut.

Referenced by getPrunedSubjetExclusive(), and getPrunedSubjetInclusive().

                                                                               {
    // cluster everything first
    exclusiveJets_ = fastjet::sorted_by_pt(fjClusterSeq_->exclusive_jets(1));
    // get pruned exclusive
    return getPrunedSubjetExclusive(0, zcut, rcutFactor);
  }

getPruned() [2/2]

math::XYZTLorentzVector heppy::ReclusterJets::getPruned ( const fastjet::PseudoJet &  jet, double  zcut, double  rcutFactor  )

private

Definition at line 100 of file ReclusterJets.cc.

References metsig::jet, and HLT_2023v12_cff::zcut.

                                                                                                            {
    // create pruner
    fastjet::Pruner pruner(fastjet::cambridge_algorithm, zcut, rcutFactor);
    // Prune the jet
    fastjet::PseudoJet pruned_jet = pruner(jet);
    return LorentzVector(pruned_jet.px(), pruned_jet.py(), pruned_jet.pz(), pruned_jet.e());
  }

getPrunedSubjetExclusive()

math::XYZTLorentzVector heppy::ReclusterJets::getPrunedSubjetExclusive	(	unsigned int	isubjet,
		double	zcut,
		double	rcutFactor
	)

get pruned 4-vector for a given subject (must be called after getGroupingExclusive)

Definition at line 83 of file ReclusterJets.cc.

References Exception, exclusiveJets_, getPruned(), and HLT_2023v12_cff::zcut.

Referenced by getPruned().

                                                                                     {
    if (isubjet >= exclusiveJets_.size()) {
      throw cms::Exception("InvalidArgument", "getPrunedSubjetExclusive called for non-existing exclusive subjet");
    }
    return getPruned(exclusiveJets_[isubjet], zcut, rcutFactor);
  }

getPrunedSubjetInclusive()

math::XYZTLorentzVector heppy::ReclusterJets::getPrunedSubjetInclusive	(	unsigned int	isubjet,
		double	zcut,
		double	rcutFactor
	)

get pruned 4-vector for a given subject (must be called after getGroupingInclusive)

Definition at line 91 of file ReclusterJets.cc.

References Exception, getPruned(), inclusiveJets_, and HLT_2023v12_cff::zcut.

                                                                                     {
    if (isubjet >= inclusiveJets_.size()) {
      throw cms::Exception("InvalidArgument", "getPrunedSubjetInclusive called for non-existing inclusive subjet");
    }
    return getPruned(inclusiveJets_[isubjet], zcut, rcutFactor);
  }

makeP4s()

std::vector< math::XYZTLorentzVector > heppy::ReclusterJets::makeP4s ( const std::vector< fastjet::PseudoJet > &  jets )

private

Definition at line 47 of file ReclusterJets.cc.

References PDWG_EXODelayedJetMET_cff::jets.

Referenced by getGrouping(), and getGroupingExclusive().

                                                                                                 {
    std::vector<math::XYZTLorentzVector> JetObjectsAll;
    JetObjectsAll.reserve(jets.size());
    for (const fastjet::PseudoJet &pj : jets) {
      JetObjectsAll.push_back(LorentzVector(pj.px(), pj.py(), pj.pz(), pj.e()));
    }
    return JetObjectsAll;
  }

Member Data Documentation

exclusiveJets_

std::vector<fastjet::PseudoJet> heppy::ReclusterJets::exclusiveJets_

private

Definition at line 60 of file ReclusterJets.h.

Referenced by getGroupingExclusive(), getPruned(), and getPrunedSubjetExclusive().

fjClusterSeq_

ClusterSequencePtr heppy::ReclusterJets::fjClusterSeq_

private

Definition at line 58 of file ReclusterJets.h.

Referenced by getGrouping(), getGroupingExclusive(), getPruned(), and ReclusterJets().

fjInputs_

std::vector<fastjet::PseudoJet> heppy::ReclusterJets::fjInputs_

private

Definition at line 51 of file ReclusterJets.h.

Referenced by ReclusterJets().

inclusiveJets_

std::vector<fastjet::PseudoJet> heppy::ReclusterJets::inclusiveJets_

private

Definition at line 59 of file ReclusterJets.h.

Referenced by getGrouping(), and getPrunedSubjetInclusive().

ktpower_

double heppy::ReclusterJets::ktpower_

private

Definition at line 53 of file ReclusterJets.h.

Referenced by ReclusterJets().

rparam_

double heppy::ReclusterJets::rparam_

private

Definition at line 54 of file ReclusterJets.h.

Referenced by ReclusterJets().