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

Mott scattering #66

Merged
merged 24 commits into from
Oct 15, 2024
Merged

Conversation

juniorpena
Copy link
Contributor

This change is implementing Mott scattering of high energy electrons off of Helium nuclei. I uploaded the files from my Kassiopeia setup on my local machine into a repo I forked, and I made changes to the appropriate CMakeLists.txt files.

Copy link
Member

@zykure zykure left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for the PR and for providing documentation with the code! I have some comments that need to be addressed before merging.

Kassiopeia/Interactions/Source/KSIntCalculatorMott.cxx Outdated Show resolved Hide resolved
Kassiopeia/Interactions/Source/KSIntCalculatorMott.cxx Outdated Show resolved Hide resolved
@zykure
Copy link
Member

zykure commented Jan 9, 2023

Is there any update on this PR @juniorpena ?

@juniorpena
Copy link
Contributor Author

juniorpena commented Feb 1, 2023 via email

@zykure zykure requested a review from 2xB February 2, 2023 08:19
@juniorpena
Copy link
Contributor Author

Hello, I've eliminated the use of ROOT in this scattering calculator. Now the GetTheta function uses inverse transform sampling and rejection sampling. I've also updated it to allow for both He and Ne Mott Scattering

@2xB
Copy link
Member

2xB commented May 7, 2024

Hey, thank you very much! Could you add documentation of the added objects to https://github.com/KATRIN-Experiment/Kassiopeia/blob/main/Kassiopeia/XML/Complete/Interactions.xml ? That is helpful both for users and for the review.

(Related, but independent issue: #106 )

@juniorpena
Copy link
Contributor Author

I've added the documentation

Copy link
Member

@2xB 2xB left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

After looking through the code, I had three remaining comments after which this looks good to merge!

Btw. if there will be a thesis describing a use case of this, one can also advertise it in the top comment of KSIntCalculatorMott.h, since it is always useful to see what the original use case for components of Kassiopeia was.

juniorpena and others added 2 commits June 10, 2024 11:23
…ilder.h

Co-authored-by: 2xB <31772910+2xB@users.noreply.github.com>
I've added code and comments for the choice of size of theta_arr responding to previous comment. I've also added comments wherever there are hard coded values, for where those values come from
@2xB
Copy link
Member

2xB commented Jul 4, 2024

Our Ubuntu 20.04 CI check failed with

/__w/Kassiopeia/Kassiopeia/Kassiopeia/Interactions/Source/KSIntCalculatorMott.cxx:10:7: error: use of undeclared identifier 'KGeoBag'
using KGeoBag::KThreeVector;
      ^
/__w/Kassiopeia/Kassiopeia/Kassiopeia/Interactions/Source/KSIntCalculatorMott.cxx:43:5: error: unknown type name 'KThreeVector'; did you mean 'katrin::KThreeVector'?
    KThreeVector tPosition = anInitialParticle.GetPosition();
    ^~~~~~~~~~~~
    katrin::KThreeVector

, as can be seen by clicking on the check details...

@2xB
Copy link
Member

2xB commented Jul 4, 2024

Btw. additionally to fixing the above error, I believe that it would make sense to update your branch by merging our main branch into it, that would fix the fedora_37 test case.

Merged main branch, and fixed bug causing Ubuntu 20.04 CI check failure
Found a bug after running simulations, the momentum of the final state particle was being changed twice when executing the interaction
@2xB
Copy link
Member

2xB commented Aug 26, 2024

@juniorpena What is the status here? I have added one question to resolve one of the two open threads here, the other I think is not that important but I leave it open just in case you want to comment on it as well. I think once at least the thread regarding the provided isotopes is solved, we should finally be able to merge this from my side. (CC @richeldichel )

@juniorpena
Copy link
Contributor Author

I've replied to the isotopes thread. Thank you for the review!

Co-authored-by: 2xB <31772910+2xB@users.noreply.github.com>
@2xB
Copy link
Member

2xB commented Sep 12, 2024

Independent on the last nit-picky comments from my side, this looks good to merge! @richeldichel

@richeldichel
Copy link
Contributor

Since all of the remaining conversations are resolved, I will merge this now. Thanks @juniorpena for your work and @2xB for the careful review!

@richeldichel richeldichel merged commit 1870d94 into KATRIN-Experiment:main Oct 15, 2024
4 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants