-
Notifications
You must be signed in to change notification settings - Fork 52
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
feat: Implement Hans Drevermanns z-finder and V-plot track finder #27
Comments
Out of curiosity, do you have any written sources I can use to learn more about this? 😄 |
Yes, I do - I also have the original fortran code, let me add the links here for the notes
I am happy to discuss these ones, I am also sure Hans would be delighted to explain the algorithm |
I see, so if I understand correctly the idea of the z-finder is to estimate a number of possible primary vertices along the z-axis, and then you can perform a much cheaper seed finding, because seed finding is much easier if you know the z-position of the vertex. Did I grok that correctly? For the average event what is the ratio (very roughly) of primary vertices compared to the number of spacepoints? |
Yes, although it has it's shortcomings, if you don't find the primary vertex (e.g. in Higgs-> gamma gamma) you miss the entire thing. A possible extension would be e.g. a scan of the vertex region, if the filter algorithm is fast enough and you can parallelise, then you can omit the z-finder and run e.g. on 1000 z-vertex bins in parallel. For Run-3: ~60 primary vertices, ~40-60k space points |
Alternatively to (re-)implementing the seed finder in
traccc
I suggest to add an implementation of Hans Drevermann's z-finder and v-plot. I have both the original Fortran code and a transcribed version in C++, but it would need serious adaption to fit into thetraccc
repository.However, it would offer a different pattern recognition algorithm with different container access and would thus build a great alternative reconstruction path and test of
vecmem
.The text was updated successfully, but these errors were encountered: