int HiForest::getMatchedHBHEAllowReuse(int j) { int best=-1; float bestDR2 = 1e10; float eta = track.trkEta[j]; float phi = track.trkPhi[j]; for (int i=0; i<hbhe.n; i++) { float dR2 = deltaR2(eta,phi,hbhe.eta[i],hbhe.phi[i]); if (dR2<bestDR2) { bestDR2 = dR2; best = i; } } return best; }
bool operator() (const pair <lorentzVector *, lorentzVector *> & x, const pair <lorentzVector *, lorentzVector *> & y) { return (deltaR2 (x.first->Eta (), x.first->Phi (), x.second->Eta (), x.second->Phi()) < deltaR2 (y.first->Eta (), y.first->Phi (), y.second->Eta (), y.second->Phi())) ; }