fastjet::atlas::JetSplitMergeTool Class Reference

#include <JetSplitMergeTool.hh>

Public Types
typedef Jet::jet_list_t	jetcollection_t
Public Member Functions
	JetSplitMergeTool ()
virtual	~JetSplitMergeTool ()
virtual int	execute (jetcollection_t *theJets)
double	etaTrue (Jet::constit_vect_t::iterator pj)
double	phiTrue (Jet::constit_vect_t::iterator pj)
void	split_merge ()
Public Attributes
double	m_f
jetcollection_t	m_preJet
jetcollection_t	m_jet
jetcollection_t *	m_jetVec
int	m_ctr
int	m_dctr

Detailed Description

Definition at line 52 of file JetSplitMergeTool.hh.

Member Typedef Documentation

typedef Jet::jet_list_t fastjet::atlas::JetSplitMergeTool::jetcollection_t

Definition at line 55 of file JetSplitMergeTool.hh.

Constructor & Destructor Documentation

fastjet::atlas::JetSplitMergeTool::JetSplitMergeTool ( )

Definition at line 27 of file JetSplitMergeTool.cc.

00028   :  m_f( 0.5 )
00029 {}

fastjet::atlas::JetSplitMergeTool::~JetSplitMergeTool ( ) [virtual]

Definition at line 31 of file JetSplitMergeTool.cc.

00032 {}

Member Function Documentation

double fastjet::atlas::JetSplitMergeTool::etaTrue ( Jet::constit_vect_t::iterator pj )

Definition at line 181 of file JetSplitMergeTool.cc.

00182 {
00183   double s = ((*pj)->e() > 0) ? +1.0 : -1.0;
00184   double px = (*pj)->px;
00185   double py = (*pj)->py;
00186   double pz = (*pj)->pz;
00187   double theta = acos(pz*s/sqrt(px*px+py*py+pz*pz));
00188   return -log(tan(theta/2.0));
00189 }

int fastjet::atlas::JetSplitMergeTool::execute ( jetcollection_t * theJets ) [virtual]

Definition at line 37 of file JetSplitMergeTool.cc.

References fastjet::atlas::Jet::addJet(), fastjet::atlas::clear_list(), fastjet::atlas::LorentzVector::e(), m_ctr, m_dctr, m_jet, m_preJet, and split_merge().

Referenced by fastjet::ATLASConePlugin::run_clustering().

00038 {
00039   m_ctr = 0;
00040   m_dctr = 0;
00041 
00043   // From the input, collection create a list of Jet//
00045   m_preJet.clear();
00046   m_jet.clear();
00047 
00048   jetcollection_t::iterator itrB = theJets->begin();
00049   jetcollection_t::iterator itrE = theJets->end(); 
00050 
00051   double etot =0.;
00052   for (;itrB!=itrE;itrB++) {
00053     Jet* j = new Jet(); j->addJet(*itrB);
00054     m_ctr +=1;
00055     m_preJet.push_back(j);    
00056 
00057     etot += j->e();    
00058   }
00059 
00061   // Split Merge Jets//
00063   this->split_merge();
00064  
00066   // Empty and re-fill input jetcollection_t //
00068   clear_list(*theJets);
00069   jetcollection_t::iterator it = m_jet.begin();
00070   jetcollection_t::iterator itE = m_jet.end();
00071   for(; it!=itE; ++it){    
00072     theJets->push_back(*it);
00073   }
00074 
00075   return 1;
00076 }

double fastjet::atlas::JetSplitMergeTool::phiTrue ( Jet::constit_vect_t::iterator pj )

Definition at line 191 of file JetSplitMergeTool.cc.

00192 {
00193   double s = ((*pj)->e() > 0) ? +1.0 : -1.0;
00194   double px = (*pj)->px;
00195   double py = (*pj)->py;
00196   return atan2(py*s,px*s);
00197 }

void fastjet::atlas::JetSplitMergeTool::split_merge ( )

Definition at line 82 of file JetSplitMergeTool.cc.

References fastjet::atlas::JetDistances::deltaPhi(), fastjet::atlas::Jet::firstConstituent(), fastjet::atlas::Jet::getConstituentNum(), fastjet::atlas::jet_from_overlap(), fastjet::atlas::Jet::lastConstituent(), m_ctr, m_dctr, m_f, m_jet, m_preJet, fastjet::atlas::LorentzVector::px, fastjet::atlas::LorentzVector::py, fastjet::atlas::Jet::removeConstituent(), and fastjet::atlas::sort_list_et().

Referenced by execute().

00084 00085 00086 00087 00088 00089 00090 00091 00092 00093       ++itr; 00094 00095 00096 00097 00098 00099 00100 00101 00102 00103 00104 00105 00106 00107 00108 00109 00110 00111 00112 00113 00114 00115 00116 00117 00118 00119 00120 00121 00122 00123 00124 00125 00126 00127 00128 00129 00130 00131 00132 00133 00134 00135 00136 00137 00138 00139 00140 00141 00142 00143 00144           } 00145 00146 00147 00148 00149         } 00150 00151 00152 00153 00154 00155       } 00156 00157 00158 00159 00160 00161       } 00162 00163 00164   } 00165 00166 00167   } 00168 00169 }

class="fragment">00083 { if ( m_preJet.size() >= 2 ) { do { sort_list_et(m_preJet); jetcollection_t::iterator itr; jetcollection_t::iterator first = m_preJet.begin(); jetcollection_t::iterator last = m_preJet.end(); itr=first; bool overlap = false; for (;itr != last;++itr) { double etaF = (*first)->eta(); double phiF = (*first)->phi(); double etaS = (*itr)->eta(); double phiS = (*itr)->phi(); Jet* oJet = jet_from_overlap( (*first),*itr); m_ctr +=1; Jet::constit_vect_t::iterator itro = oJet->firstConstituent(); Jet::constit_vect_t::iterator itroE = oJet->lastConstituent(); if ( oJet->getConstituentNum() != 0 ) { overlap = true; // fraction double f = sqrt(pow(oJet->px,2)+pow(oJet->py,2))/ sqrt(pow((*itr)->px,2)+pow((*itr)->py,2)); // merge if ( f > m_f) { // we need to remove constituents ! Jet *j = (*first); for ( ;itro != itroE; ++itro ) j->removeConstituent(*itro); (*first)->addJet(*itr); //m_preJet.remove(*itr); delete *itr; m_preJet.erase(itr); m_dctr +=1; } // split if ( f <= m_f) { for ( ;itro != itroE; ++itro ) { // Distance of first jet from ProtoJet double deta1 = etaF - (*itro)->eta(); double dphi1 = fabs(JetDistances::deltaPhi(phiF,(*itro)->phi())); double dist1 = pow( deta1 , 2 ) + pow( dphi1 , 2 ); // Distance of second jet from ProtoJet double deta2 = etaS - (*itro)->eta(); double dphi2 = fabs(JetDistances::deltaPhi(phiS,(*itro)->phi())); double dist2 = pow( deta2 , 2 ) + pow( dphi2 , 2 ); // Remove protojet from farther Jet if ( dist1 > dist2 ) (*first)->removeConstituent(*itro); if ( dist1 <= dist2 ) (*itr)->removeConstituent(*itro); } // Delete overlap jet delete oJet; m_dctr +=1; break; else { // Delete overlap jet delete oJet; m_dctr +=1; } if ( overlap == false ) { m_jet.push_back(*first); //m_preJet.remove(*first); m_preJet.erase(first); } while ( m_preJet.size() != 0 ); else if ( m_preJet.size() == 1) { m_jet.push_back( *(m_preJet.begin()) );