#include <SubJetAlgorithm.h>

Public Member Functions
double	rcut_factor () const

void	run (const std::vector< fastjet::PseudoJet > &cell_particles, std::vector< CompoundPseudoJet > &hardjetsOutput)
	Find the ProtoJets from the collection of input Candidates. More...

void	set_rcut_factor (double r)

void	set_zcut (double z)

	SubJetAlgorithm (double ptMin, unsigned int subjets, double zcut, double rcut_factor, std::shared_ptr< fastjet::JetDefinition > fjJetDefinition, bool doAreaFastjet, std::shared_ptr< fastjet::GhostedAreaSpec > fjActiveArea, double voronoiRfact)

double	zcut () const

Private Attributes
bool	doAreaFastjet_

std::shared_ptr < fastjet::GhostedAreaSpec >	fjActiveArea_

std::shared_ptr < fastjet::JetDefinition >	fjJetDefinition_

int	nSubjets_

double	ptMin_

double	rcut_factor_

double	voronoiRfact_

double	zcut_

Detailed Description

Definition at line 14 of file SubJetAlgorithm.h.

Constructor & Destructor Documentation

SubJetAlgorithm::SubJetAlgorithm	(	double	ptMin,
		unsigned int	subjets,
		double	zcut,
		double	rcut_factor,
		std::shared_ptr< fastjet::JetDefinition >	fjJetDefinition,
		bool	doAreaFastjet,
		std::shared_ptr< fastjet::GhostedAreaSpec >	fjActiveArea,
		double	voronoiRfact
	)

inline

Definition at line 16 of file SubJetAlgorithm.h.

       : ptMin_(ptMin),
         nSubjets_(subjets),
         zcut_(zcut),
         rcut_factor_(rcut_factor),
         fjJetDefinition_(fjJetDefinition),
         doAreaFastjet_(doAreaFastjet),
         fjActiveArea_(fjActiveArea),
         voronoiRfact_(voronoiRfact) {}

Member Function Documentation

double SubJetAlgorithm::rcut_factor ( ) const

inline

Definition at line 36 of file SubJetAlgorithm.h.

References rcut_factor_.

36 { return rcut_factor_; }

SubJetAlgorithm::rcut_factor_

double rcut_factor_

Definition: SubJetAlgorithm.h:45

void SubJetAlgorithm::run	(	const std::vector< fastjet::PseudoJet > &	cell_particles,
		std::vector< CompoundPseudoJet > &	hardjetsOutput
	)

Find the ProtoJets from the collection of input Candidates.

Definition at line 17 of file SubJetAlgorithm.cc.

References dqmdumpme::indices, and ptMin_().

Referenced by cms::SubJetProducer::runAlgorithm().

                                                                                                                    {
   //for actual jet clustering, either the pruned or the original version is used.
   //For the pruned version, a new jet definition using the PrunedRecombPlugin is required:
   fastjet::FastPrunePlugin PRplugin(*fjJetDefinition_, *fjJetDefinition_, zcut_, rcut_factor_);
   fastjet::JetDefinition pjetdef(&PRplugin);
 
   // cluster the jets with the jet definition jetDef:
   // run algorithm
   std::shared_ptr<fastjet::ClusterSequence> fjClusterSeq;
   if (!doAreaFastjet_) {
     fjClusterSeq = std::make_shared<fastjet::ClusterSequence>(cell_particles, pjetdef);
   } else if (voronoiRfact_ <= 0) {
     fjClusterSeq = std::shared_ptr<fastjet::ClusterSequence>(
         new fastjet::ClusterSequenceActiveArea(cell_particles, pjetdef, *fjActiveArea_));
   } else {
     fjClusterSeq = std::shared_ptr<fastjet::ClusterSequence>(
         new fastjet::ClusterSequenceVoronoiArea(cell_particles, pjetdef, fastjet::VoronoiAreaSpec(voronoiRfact_)));
   }
 
   vector<fastjet::PseudoJet> inclusiveJets = fjClusterSeq->inclusive_jets(ptMin_);
 
   // These will store the indices of each subjet that
   // are present in each jet
   vector<vector<int> > indices(inclusiveJets.size());
   // Loop over inclusive jets, attempt to find substructure
   vector<fastjet::PseudoJet>::iterator jetIt = inclusiveJets.begin();
   for (; jetIt != inclusiveJets.end(); ++jetIt) {
     //decompose into requested number of subjets:
     vector<fastjet::PseudoJet> subjets = fjClusterSeq->exclusive_subjets(*jetIt, nSubjets_);
     //create the subjets objects to put into the "output" objects
     vector<CompoundPseudoSubJet> subjetsOutput;
     std::vector<fastjet::PseudoJet>::const_iterator itSubJetBegin = subjets.begin(), itSubJet = itSubJetBegin,
                                                     itSubJetEnd = subjets.end();
     for (; itSubJet != itSubJetEnd; ++itSubJet) {
       // Get the transient subjet constituents from fastjet
       vector<fastjet::PseudoJet> subjetFastjetConstituents = fjClusterSeq->constituents(*itSubJet);
       // Get the indices of the constituents:
       vector<int> constituents;
       vector<fastjet::PseudoJet>::const_iterator fastSubIt = subjetFastjetConstituents.begin(),
                                                  transConstEnd = subjetFastjetConstituents.end();
       for (; fastSubIt != transConstEnd; ++fastSubIt) {
         if (fastSubIt->user_index() >= 0) {
           constituents.push_back(fastSubIt->user_index());
         }
       }
 
       double subJetArea =
           (doAreaFastjet_) ? dynamic_cast<fastjet::ClusterSequenceActiveArea&>(*fjClusterSeq).area(*itSubJet) : 0.0;
 
       // Make a CompoundPseudoSubJet object to hold this subjet and the indices of its constituents
       subjetsOutput.push_back(CompoundPseudoSubJet(*itSubJet, subJetArea, constituents));
     }
 
     double fatJetArea =
         (doAreaFastjet_) ? dynamic_cast<fastjet::ClusterSequenceActiveArea&>(*fjClusterSeq).area(*jetIt) : 0.0;
 
     // Make a CompoundPseudoJet object to hold this hard jet, and the subjets that make it up
     hardjetsOutput.push_back(CompoundPseudoJet(*jetIt, fatJetArea, subjetsOutput));
   }
 }