#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 boost::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 20 of file ReclusterJets.h.

Member Typedef Documentation

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

private

fastjet outputs

Definition at line 60 of file ReclusterJets.h.

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

Definition at line 23 of file ReclusterJets.h.

Constructor & Destructor Documentation

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

Definition at line 12 of file ReclusterJets.cc.

References Exception, fjClusterSeq_, fjInputs_, cmsHarvester::index, j, fwrapper::jet_def, ktpower_, connectstrParser::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_ = ClusterSequencePtr( new fastjet::ClusterSequence( fjInputs_, jet_def)); 
 }

Member Function Documentation

std::vector< math::XYZTLorentzVector > heppy::ReclusterJets::getGrouping ( double ptMin = 0.0 )

get grouping (inclusive jets)

Definition at line 52 of file ReclusterJets.cc.

References fjClusterSeq_, inclusiveJets_, and makeP4s().

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

std::vector< math::XYZTLorentzVector > heppy::ReclusterJets::getGroupingExclusive ( int njets )

get grouping (exclusive jets, until n are left)

Definition at line 66 of file ReclusterJets.cc.

References exclusiveJets_, fjClusterSeq_, and makeP4s().

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

std::vector< math::XYZTLorentzVector > heppy::ReclusterJets::getGroupingExclusive ( double dcut )

get grouping (exclusive jets, up to cut dcut)

Definition at line 59 of file ReclusterJets.cc.

References exclusiveJets_, fjClusterSeq_, and makeP4s().

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

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

get pruned 4-vector

Definition at line 73 of file ReclusterJets.cc.

References exclusiveJets_, fjClusterSeq_, and getPrunedSubjetExclusive().

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);
 }

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

private

Definition at line 93 of file ReclusterJets.cc.

                                                                                                            {
  // 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() );
 }

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 80 of file ReclusterJets.cc.

References Exception, exclusiveJets_, and getPruned().

Referenced by getPruned().

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

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 86 of file ReclusterJets.cc.

References Exception, getPruned(), and inclusiveJets_.

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

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

private

Definition at line 45 of file ReclusterJets.cc.

Referenced by getGrouping(), and getGroupingExclusive().

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