SubJetAlgorithm Class Reference

#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::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

rcut_factor()

double SubJetAlgorithm::rcut_factor ( ) const

inline

Definition at line 36 of file SubJetAlgorithm.h.

References rcut_factor_.

{ return rcut_factor_; }

run()

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.

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

set_rcut_factor()

void SubJetAlgorithm::set_rcut_factor

(

double

r

)

Definition at line 13 of file SubJetAlgorithm.cc.

{ rcut_factor_ = r; }

set_zcut()

void SubJetAlgorithm::set_zcut

(

double

z

)

Definition at line 11 of file SubJetAlgorithm.cc.

{ zcut_ = z; }

zcut()

double SubJetAlgorithm::zcut ( ) const

inline

Definition at line 35 of file SubJetAlgorithm.h.

References zcut_.

{ return zcut_; }

Member Data Documentation

doAreaFastjet_

bool SubJetAlgorithm::doAreaFastjet_

private

Definition at line 47 of file SubJetAlgorithm.h.

fjActiveArea_

std::shared_ptr<fastjet::GhostedAreaSpec> SubJetAlgorithm::fjActiveArea_

private

Definition at line 48 of file SubJetAlgorithm.h.

fjJetDefinition_

std::shared_ptr<fastjet::JetDefinition> SubJetAlgorithm::fjJetDefinition_

private

Definition at line 46 of file SubJetAlgorithm.h.

nSubjets_

int SubJetAlgorithm::nSubjets_

private

Definition at line 43 of file SubJetAlgorithm.h.

ptMin_

double SubJetAlgorithm::ptMin_

private

Definition at line 42 of file SubJetAlgorithm.h.

rcut_factor_

double SubJetAlgorithm::rcut_factor_

private

Definition at line 45 of file SubJetAlgorithm.h.

Referenced by rcut_factor().

voronoiRfact_

double SubJetAlgorithm::voronoiRfact_

private

Definition at line 49 of file SubJetAlgorithm.h.

zcut_

double SubJetAlgorithm::zcut_

private

Definition at line 44 of file SubJetAlgorithm.h.

Referenced by zcut().