fastjet: fastjet::CDFMidPointPlugin Class Reference

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:

00086                                                                         {
00087  
00088   // create the physics towers needed by the CDF code
00089   vector<PhysicsTower> towers;
00090   towers.reserve(clust_seq.jets().size());
00091   for (unsigned i = 0; i < clust_seq.jets().size(); i++) {
00092     LorentzVector fourvect(clust_seq.jets()[i].px(),
00093                            clust_seq.jets()[i].py(),
00094                            clust_seq.jets()[i].pz(),
00095                            clust_seq.jets()[i].E());
00096     PhysicsTower tower(fourvect);
00097     // misuse one of the indices for tracking, since the MidPoint
00098     // implementation doesn't seem to make use of these indices
00099     tower.calTower.iEta = i;
00100     towers.push_back(tower);
00101   }
00102 
00103   // prepare the CDF algorithm
00104   MidPointAlgorithm m(_seed_threshold,_cone_radius,_cone_area_fraction,
00105                       _max_pair_size,_max_iterations,_overlap_threshold,
00106                       MidPointAlgorithm::SplitMergeScale(_sm_scale));
00107     
00108   // run the CDF algorithm
00109   std::vector<Cluster> jets;
00110   m.run(towers,jets);
00111 
00112 
00113   // now transfer the jets back into our own structure -- we will
00114   // mimic the cone code with a sequential recombination sequence in
00115   // which the jets are built up by adding one particle at a time
00116   for(vector<Cluster>::const_iterator jetIter = jets.begin(); 
00117                                       jetIter != jets.end(); jetIter++) {
00118     const vector<PhysicsTower> & tower_list = jetIter->towerList;
00119     int jet_k = tower_list[0].calTower.iEta;
00120   
00121     int ntow = int(jetIter->towerList.size());
00122     for (int itow = 1; itow < ntow; itow++) {
00123       int jet_i = jet_k;
00124       // retrieve our misappropriated index for the jet
00125       int jet_j = tower_list[itow].calTower.iEta;
00126       // do a fake recombination step with dij=0
00127       double dij = 0.0;
00128       clust_seq.plugin_record_ij_recombination(jet_i, jet_j, dij, jet_k);
00129     }
00130   
00131     // NB: put a sensible looking d_iB just to be nice...
00132     double d_iB = clust_seq.jets()[jet_k].perp2();
00133     clust_seq.plugin_record_iB_recombination(jet_k, d_iB);
00134   }
00135 
00136 
00137   // following code is for testing only
00138   //cout << endl;
00139   //for(vector<Cluster>::const_iterator jetIter = jets.begin(); 
00140   //                                    jetIter != jets.end(); jetIter++) {
00141   //  cout << jetIter->fourVector.pt() << " " << jetIter->fourVector.y() << endl;
00142   //}
00143   //cout << "-----------------------------------------------------\n";
00144   //vector<PseudoJet> ourjets(clust_seq.inclusive_jets());
00145   //for (vector<PseudoJet>::const_reverse_iterator ourjet = ourjets.rbegin();
00146   //     ourjet != ourjets.rend(); ourjet++) {
00147   //  cout << ourjet->perp() << " " << ourjet->rap() << endl;
00148   //}
00149   //cout << endl;
00150 }


Public Types
enum	SplitMergeScale { SM_pt, SM_Et, SM_mt, SM_pttilde }
	the choice of scale to be used in the split-merge step More...
Public Member Functions
	CDFMidPointPlugin (double seed_threshold, double cone_radius, double cone_area_fraction, int max_pair_size, int max_iterations, double overlap_threshold, SplitMergeScale sm_scale=SM_pt)
	A CDFMidPointPlugin constructor that looks like the one provided by CDF.
	CDFMidPointPlugin (double cone_radius, double overlap_threshold=0.5, double seed_threshold=1.0, double cone_area_fraction=1.0)
	a compact constructor
double	seed_threshold () const
double	cone_radius () const
double	cone_area_fraction () const
int	max_pair_size () const
int	max_iterations () 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 Attributes
double	_seed_threshold
double	_cone_radius
double	_cone_area_fraction
int	_max_pair_size
int	_max_iterations
double	_overlap_threshold
SplitMergeScale	_sm_scale

fastjet::CDFMidPointPlugin Class Reference

Public Types

Public Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation