| 170 |
|
double deltarmin = double_infinity(); |
| 171 |
|
Jet* nextjet=0; |
| 172 |
|
for(unsigned int i=0; i<jets->size(); ++i) { |
| 173 |
+ |
Jet ji = jets->at(i); |
| 174 |
+ |
if (fabs(p->pt() - ji.pt())<1e-8) continue; // skip identical particle |
| 175 |
|
if(jets->at(i).deltaR(*p) < deltarmin) { |
| 176 |
|
deltarmin = jets->at(i).deltaR(*p); |
| 177 |
|
nextjet = &jets->at(i); |
| 219 |
|
{ |
| 220 |
|
|
| 221 |
|
double ptrel=0; |
| 220 |
– |
|
| 222 |
|
Jet* nextjet = nextJet(p,jets); |
| 223 |
+ |
if (!nextjet) return ptrel; |
| 224 |
|
|
| 225 |
|
TVector3 p3(p->v4().Px(),p->v4().Py(),p->v4().Pz()); |
| 226 |
|
TVector3 jet3(nextjet->v4().Px(),nextjet->v4().Py(),nextjet->v4().Pz()); |
| 237 |
|
|
| 238 |
|
double deltaRmin(const Particle *p, std::vector<Jet> *jets) |
| 239 |
|
{ |
| 240 |
< |
return nextJet(p,jets)->deltaR(*p); |
| 240 |
> |
Jet* j = nextJet(p,jets); |
| 241 |
> |
double dr = 999.; |
| 242 |
> |
if (j) dr = j->deltaR(*p); |
| 243 |
> |
return dr; |
| 244 |
|
} |
| 245 |
|
|
| 246 |
|
TVector3 toVector(LorentzVector v4) |
