Main Page | Namespace List | Class Hierarchy | Class List | Directories | File List | Namespace Members | Class Members | File Members

fastjet::CDFJetCluPlugin Class Reference

a plugin for fastjet-v2.1 that provides an interface to the CDF jetclu algorithm More...

#include <CDFJetCluPlugin.hh>

Inheritance diagram for fastjet::CDFJetCluPlugin:

Inheritance graph
[legend]
Collaboration diagram for fastjet::CDFJetCluPlugin:

Collaboration graph
[legend]
List of all members.

Public Member Functions

 CDFJetCluPlugin (double cone_radius, double overlap_threshold=0.5, double seed_threshold=1.0)
 a compact constructor
 CDFJetCluPlugin (double seed_threshold, double cone_radius, int adjacency_cut, int max_iterations, int iratch, double overlap_threshold)
 a constructor that looks like the one provided by CDF
double seed_threshold () const
double cone_radius () const
int adjacency_cut () const
int max_iterations () const
int iratch () const
double overlap_threshold () const
virtual std::string description () const
 return a textual description of the jet-definition implemented in this plugin
virtual void run_clustering (ClusterSequence &) const
 given a ClusterSequence that has been filled up with initial particles, the following function should fill up the rest of the ClusterSequence, using the following member functions of ClusterSequence:
  • plugin_do_ij_recombination(.


Private Member Functions

void _insert_unique (PseudoJet &jet, std::map< double, int > &jetmap) const
 given a jet try inserting its energy into the map -- if that energy entry already exists, modify the jet infinitesimally so as ensure that the jet energy is unique

Private Attributes

double _seed_threshold
double _cone_radius
int _adjacency_cut
int _max_iterations
int _iratch
double _overlap_threshold

Detailed Description

a plugin for fastjet-v2.1 that provides an interface to the CDF jetclu algorithm

Definition at line 44 of file CDFJetCluPlugin.hh.


Constructor & Destructor Documentation

fastjet::CDFJetCluPlugin::CDFJetCluPlugin double  cone_radius,
double  overlap_threshold = 0.5,
double  seed_threshold = 1.0
[inline]
 

a compact constructor

Definition at line 47 of file CDFJetCluPlugin.hh.

00049                                                   : 
00050     _seed_threshold    ( seed_threshold    ),    
00051     _cone_radius       ( cone_radius       ),
00052     _adjacency_cut     (   2               ),
00053     _max_iterations    ( 100               ),
00054     _iratch            (   1               ),
00055     _overlap_threshold ( overlap_threshold )  {}

fastjet::CDFJetCluPlugin::CDFJetCluPlugin double  seed_threshold,
double  cone_radius,
int  adjacency_cut,
int  max_iterations,
int  iratch,
double  overlap_threshold
[inline]
 

a constructor that looks like the one provided by CDF

Definition at line 58 of file CDFJetCluPlugin.hh.


Member Function Documentation

void fastjet::CDFJetCluPlugin::_insert_unique PseudoJet jet,
std::map< double, int > &  jetmap
const [private]
 

given a jet try inserting its energy into the map -- if that energy entry already exists, modify the jet infinitesimally so as ensure that the jet energy is unique

Definition at line 151 of file CDFJetCluPlugin.cc.

Referenced by run_clustering().

00152                                                                      {
00153   while (jetmap.find(jet.E()) != jetmap.end()) {
00154     // deal with cases where something else has the same energy, and
00155     // also with situation where that energy is zero.
00156     if (jet.E() != 0.0) {
00157       jet *= 1.0+1e-12;
00158     } else {
00159       jet += PseudoJet(0.0,0.0,0.0,1e-300);
00160     }
00161   }
00162   jetmap[jet.E()] = jet.cluster_hist_index();
00163 }

int fastjet::CDFJetCluPlugin::adjacency_cut  )  const [inline]
 

Definition at line 75 of file CDFJetCluPlugin.hh.

Referenced by description(), and run_clustering().

00075 {return _adjacency_cut     ;}

double fastjet::CDFJetCluPlugin::cone_radius  )  const [inline]
 

Definition at line 74 of file CDFJetCluPlugin.hh.

Referenced by description(), and run_clustering().

00074 {return _cone_radius       ;}

string fastjet::CDFJetCluPlugin::description  )  const [virtual]
 

return a textual description of the jet-definition implemented in this plugin

Implements fastjet::JetDefinition::Plugin.

Definition at line 44 of file CDFJetCluPlugin.cc.

References adjacency_cut(), cone_radius(), iratch(), max_iterations(), overlap_threshold(), and seed_threshold().

00044                                            {
00045   ostringstream desc;
00046   
00047   desc << "CDF JetClu jet algorithm with " 
00048        << "seed_threshold = "     << seed_threshold    () << ", "
00049        << "cone_radius = "        << cone_radius       () << ", "
00050        << "adjacency_cut = "      << adjacency_cut     () << ", " 
00051        << "max_iterations = "     << max_iterations    () << ", "
00052        << "iratch = "             << iratch            () << ", "
00053        << "overlap_threshold = "  << overlap_threshold () ;
00054 
00055   return desc.str();
00056 }

int fastjet::CDFJetCluPlugin::iratch  )  const [inline]
 

Definition at line 77 of file CDFJetCluPlugin.hh.

Referenced by description(), and run_clustering().

00077 {return _iratch            ;}

int fastjet::CDFJetCluPlugin::max_iterations  )  const [inline]
 

Definition at line 76 of file CDFJetCluPlugin.hh.

Referenced by description(), and run_clustering().

00076 {return _max_iterations    ;}

double fastjet::CDFJetCluPlugin::overlap_threshold  )  const [inline]
 

Definition at line 78 of file CDFJetCluPlugin.hh.

Referenced by description(), and run_clustering().

00078 {return _overlap_threshold ;}

void fastjet::CDFJetCluPlugin::run_clustering ClusterSequence  )  const [virtual]
 

given a ClusterSequence that has been filled up with initial particles, the following function should fill up the rest of the ClusterSequence, using the following member functions of ClusterSequence:

  • plugin_do_ij_recombination(.

..)

  • plugin_do_iB_recombination(...)

Implements fastjet::JetDefinition::Plugin.

Definition at line 59 of file CDFJetCluPlugin.cc.

References _insert_unique(), adjacency_cut(), cone_radius(), fastjet::PseudoJet::E(), iratch(), fastjet::ClusterSequence::jets(), max_iterations(), overlap_threshold(), fastjet::ClusterSequence::plugin_record_iB_recombination(), fastjet::ClusterSequence::plugin_record_ij_recombination(), fastjet::PseudoJet::px(), fastjet::PseudoJet::py(), fastjet::PseudoJet::pz(), and seed_threshold().

00059                                                                       {
00060  
00061   // create the physics towers needed by the CDF code
00062   vector<PhysicsTower> towers;
00063   towers.reserve(clust_seq.jets().size());
00064 
00065   // create a map to identify jets (actually just the input particles)...
00066   map<double,int> jetmap;
00067 
00068   for (unsigned i = 0; i < clust_seq.jets().size(); i++) {
00069     PseudoJet particle(clust_seq.jets()[i]);
00070     _insert_unique(particle, jetmap);
00071     LorentzVector fourvect(particle.px(), particle.py(),
00072                            particle.pz(), particle.E());
00073     PhysicsTower tower(fourvect);
00074     // cannot use MidPoint trick of misusing one of the indices for
00075     // tracking, since the JetClu implementation _does_ seem to make
00076     // use of these indices 
00077     //tower.calTower.iEta = i;
00078     towers.push_back(tower);
00079   }
00080 
00081   // prepare the CDF algorithm
00082   JetCluAlgorithm j(seed_threshold(), cone_radius(), adjacency_cut(),
00083                     max_iterations(), iratch(), overlap_threshold());
00084     
00085   // run the CDF algorithm
00086   std::vector<Cluster> jets;
00087   j.run(towers,jets);
00088 
00089 
00090   // now transfer the jets back into our own structure -- we will
00091   // mimic the cone code with a sequential recombination sequence in
00092   // which the jets are built up by adding one particle at a time
00093 
00094   // NB: with g++-4.0, the reverse iterator code gave problems, so switch
00095   //     to indices instead
00096   //for(vector<Cluster>::const_reverse_iterator jetIter = jets.rbegin(); 
00097   //                                    jetIter != jets.rend(); jetIter++) {
00098   //  const vector<PhysicsTower> & tower_list = jetIter->towerList;
00099   //  int jet_k = jetmap[tower_list[0].fourVector.E];
00100   //
00101   //  int ntow = int(jetIter->towerList.size());
00102 
00103   for(int iCDFjets = jets.size()-1; iCDFjets >= 0; iCDFjets--) {
00104     const vector<PhysicsTower> & tower_list = jets[iCDFjets].towerList;
00105     int jet_k = jetmap[tower_list[0].fourVector.E];
00106   
00107     int ntow = int(tower_list.size());
00108     for (int itow = 1; itow < ntow; itow++) {
00109       int jet_i = jet_k;
00110       // retrieve our misappropriated index for the jet
00111       //int jet_j = tower_list[itow].calTower.iEta;
00112       int jet_j = jetmap[tower_list[itow].fourVector.E];
00113       // do a fake recombination step with dij=0
00114       double dij = 0.0;
00115       // JetClu does E-scheme recombination so we can stick with the
00116       // simple option
00117       clust_seq.plugin_record_ij_recombination(jet_i, jet_j, dij, jet_k);
00118 
00119       //if (itow != ntow) {
00120       //  clust_seq.plugin_record_ij_recombination(jet_i, jet_j, dij, jet_k);
00121       //} else {
00122       //  clust_seq.plugin_record_ij_recombination(jet_i, jet_j, dij, 
00123       //       PseudoJet(jetIter->fourVector.px,jetIter->fourVector.py,
00124       //                 jetIter->fourVector.pz,jetIter->fourVector.E),
00125       //                                           jet_k);
00126       //}
00127     }
00128   
00129     // NB: put a sensible looking d_iB just to be nice...
00130     double d_iB = clust_seq.jets()[jet_k].perp2();
00131     clust_seq.plugin_record_iB_recombination(jet_k, d_iB);
00132   }
00133 
00134 
00135   // following code is for testing only
00136   //cout << endl;
00137   //for(vector<Cluster>::const_iterator jetIter = jets.begin(); 
00138   //                                    jetIter != jets.end(); jetIter++) {
00139   //  cout << jetIter->fourVector.pt() << " " << jetIter->fourVector.y() << endl;
00140   //}
00141   //cout << "-----------------------------------------------------\n";
00142   //vector<PseudoJet> ourjets(clust_seq.inclusive_jets());
00143   //for (vector<PseudoJet>::const_iterator ourjet = ourjets.begin();
00144   //     ourjet != ourjets.end(); ourjet++) {
00145   //  cout << ourjet->perp() << " " << ourjet->rap() << endl;
00146   //}
00147   //cout << endl;
00148 }

double fastjet::CDFJetCluPlugin::seed_threshold  )  const [inline]
 

Definition at line 73 of file CDFJetCluPlugin.hh.

Referenced by description(), and run_clustering().

00073 {return _seed_threshold    ;}


Member Data Documentation

int fastjet::CDFJetCluPlugin::_adjacency_cut [private]
 

Definition at line 91 of file CDFJetCluPlugin.hh.

double fastjet::CDFJetCluPlugin::_cone_radius [private]
 

Definition at line 90 of file CDFJetCluPlugin.hh.

int fastjet::CDFJetCluPlugin::_iratch [private]
 

Definition at line 93 of file CDFJetCluPlugin.hh.

int fastjet::CDFJetCluPlugin::_max_iterations [private]
 

Definition at line 92 of file CDFJetCluPlugin.hh.

double fastjet::CDFJetCluPlugin::_overlap_threshold [private]
 

Definition at line 94 of file CDFJetCluPlugin.hh.

double fastjet::CDFJetCluPlugin::_seed_threshold [private]
 

Definition at line 89 of file CDFJetCluPlugin.hh.


The documentation for this class was generated from the following files:
Generated on Mon Apr 2 20:58:24 2007 for fastjet by  doxygen 1.4.2