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, boost::shared_ptr< fastjet::JetDefinition > fjJetDefinition, bool doAreaFastjet, boost::shared_ptr< fastjet::ActiveAreaSpec > fjActiveArea, double voronoiRfact)
double	zcut () const

Private Attributes
bool	doAreaFastjet_
boost::shared_ptr < fastjet::ActiveAreaSpec >	fjActiveArea_
boost::shared_ptr < fastjet::JetDefinition >	fjJetDefinition_
int	nSubjets_
double	ptMin_
double	rcut_factor_
double	voronoiRfact_
double	zcut_

Detailed Description

Definition at line 15 of file SubJetAlgorithm.h.

Constructor & Destructor Documentation

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

inline

Definition at line 17 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 41 of file SubJetAlgorithm.h.

References rcut_factor_.

{ return rcut_factor_; }

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 20 of file SubJetAlgorithm.cc.

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
  boost::shared_ptr<fastjet::ClusterSequence> fjClusterSeq;
  if ( !doAreaFastjet_ ) {
    fjClusterSeq = 
      boost::shared_ptr<fastjet::ClusterSequence>( new fastjet::ClusterSequence( cell_particles, 
                                         pjetdef ) );
  } else if (voronoiRfact_ <= 0) {
    fjClusterSeq = 
      boost::shared_ptr<fastjet::ClusterSequence>( new fastjet::ClusterSequenceActiveArea( cell_particles, 
                                               pjetdef , 
                                               *fjActiveArea_ ) );
  } else {
    fjClusterSeq = 
      boost::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(),
    transConstBegin = 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));
  }
}

void SubJetAlgorithm::set_rcut_factor

(

double

r

)

Definition at line 13 of file SubJetAlgorithm.cc.

References csvReporter::r.

                                             {
    rcut_factor_ = r;
}

void SubJetAlgorithm::set_zcut

(

double

z

)

Definition at line 9 of file SubJetAlgorithm.cc.

References detailsBasic3DVector::z.

                                      {
    zcut_ = z;
}

double SubJetAlgorithm::zcut ( ) const

inline

Definition at line 40 of file SubJetAlgorithm.h.

References zcut_.

{ return zcut_;}

Member Data Documentation

bool SubJetAlgorithm::doAreaFastjet_

private

Definition at line 55 of file SubJetAlgorithm.h.

boost::shared_ptr<fastjet::ActiveAreaSpec> SubJetAlgorithm::fjActiveArea_

private

Definition at line 56 of file SubJetAlgorithm.h.

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

private

Definition at line 54 of file SubJetAlgorithm.h.

int SubJetAlgorithm::nSubjets_

private

Definition at line 51 of file SubJetAlgorithm.h.

double SubJetAlgorithm::ptMin_

private

Definition at line 50 of file SubJetAlgorithm.h.

double SubJetAlgorithm::rcut_factor_

private

Definition at line 53 of file SubJetAlgorithm.h.

Referenced by rcut_factor().

double SubJetAlgorithm::voronoiRfact_

private

Definition at line 57 of file SubJetAlgorithm.h.

double SubJetAlgorithm::zcut_

private

Definition at line 52 of file SubJetAlgorithm.h.

Referenced by zcut().