double Qjets::d_ij(const fastjet::PseudoJet& v1,const fastjet::PseudoJet& v2) { double p1 = v1.perp(); double p2 = v2.perp(); // small fix here -- nhan // double ret = pow(min(p1,p2),_exp_min) * pow(max(p1,p2),_exp_max) * v1.squared_distance(v2); double ret = pow(10.,-5.); if(v1.squared_distance(v2) != 0.) { ret = pow(min(p1,p2),_exp_min) * pow(max(p1,p2),_exp_max) * v1.squared_distance(v2); } assert(!std::isnan(ret)); return ret; }
double Qjets::d_ij(const fastjet::PseudoJet& v1,const fastjet::PseudoJet& v2){ double p1 = v1.perp(); double p2 = v2.perp(); double ret = pow(min(p1,p2),_exp_min) * pow(max(p1,p2),_exp_max) * v1.squared_distance(v2); assert(!isnan((float)ret)); return ret; }