Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Phi wraparound bug in PFDisplacedVertexCandidateFinder? #38068

Closed
rappoccio opened this issue May 24, 2022 · 13 comments · Fixed by #38080
Closed

Phi wraparound bug in PFDisplacedVertexCandidateFinder? #38068

rappoccio opened this issue May 24, 2022 · 13 comments · Fixed by #38080

Comments

@rappoccio
Copy link
Contributor

Hi,

I think there is a phi wraparound bug in this line of PFDisplacedVertexCandidateFinder:

if (pt1 > 2 && pt2 > 2 && std::abs(phi1 - phi2) > 1) {

std::abs(phi1 - phi2) > 1

should be

math::deltaPhi(phi1,phi2) > 1

This may be related to the issue reported here. Can someone confirm? We can fix it (and hopefully the other issue as well) if this is wrong as I expect.

@jpata @clacaputo @laurenhay

@cmsbuild
Copy link
Contributor

A new Issue was created by @rappoccio .

@Dr15Jones, @perrotta, @dpiparo, @makortel, @smuzaffar, @qliphy can you please review it and eventually sign/assign? Thanks.

cms-bot commands are listed here

@makortel
Copy link
Contributor

assign reconstruction

@cmsbuild
Copy link
Contributor

New categories assigned: reconstruction

@jpata,@slava77,@clacaputo you have been requested to review this Pull request/Issue and eventually sign? Thanks

@jpata
Copy link
Contributor

jpata commented May 25, 2022

cc @cms-sw/pf-l2

@jpata
Copy link
Contributor

jpata commented May 25, 2022

type pf

@cmsbuild cmsbuild added the pf label May 25, 2022
@jpata
Copy link
Contributor

jpata commented May 25, 2022

Digging through the code, it looks like Maxime was the original author. I don't know if he's on github, but I shot him an email.

@juska
Copy link
Contributor

juska commented May 25, 2022

Isn't v1.DeltaPhi(v2) equivalent to abs(v1.Phi() - v2.Phi()) ?

@perrotta
Copy link
Contributor

Isn't v1.DeltaPhi(v2) equivalent to abs(v1.Phi() - v2.Phi()) ?

Not if phi1 is 359 degrees and phi2 is 1 degree (e.g.)

@mmusich
Copy link
Contributor

mmusich commented May 25, 2022

This may be related to the issue reported #31158. Can someone confirm? We can fix it (and hopefully the other issue as well) if this is wrong as I expect.

I'd propose to make the PR and check the performance....

@mmusich
Copy link
Contributor

mmusich commented May 25, 2022

I'd propose to make the PR and check the performance....

here it is: #38080

@rappoccio
Copy link
Contributor Author

Thanks @mmusich, we were really just asking if we were missing something, @laurenhay has the fix. It sounds like this is indeed really a bug so we will go ahead with the tests.

@jpata
Copy link
Contributor

jpata commented May 30, 2022

+reconstruction

@cmsbuild
Copy link
Contributor

This issue is fully signed and ready to be closed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging a pull request may close this issue.

7 participants