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())) ;
}
