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:

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


Public Types
enum	SplitMergeScale { SM_pt, SM_Et, SM_mt }
	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