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