-
Notifications
You must be signed in to change notification settings - Fork 62
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
Add calculation dX in TPC using local curvature #383
Conversation
Yuri, my understanding is that this request is to add a user function for dX calculation. Seem to be straightforward to me. Jason/Dmitri, any further comments? |
To first order this looks fine. My only concern is whether there are other codes which write the dX value in StTpcHit? If so, what is the reason for replacing the calculation? (Or is the calculation just being moved into StiStEventFiller?) |
As it was explained in message this method used local curvature for each hit instead of using curvature from the first and last hit. The only place where we can access local curvature (StKalanTrackNode parameters) is StiStEventFiller. This is essential for low pT tracks. |
If no other objections, I would suggest to merge this PR. Thanks |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Diff LGTM
Is the time required for this extra pair of pathlength calculations for every TPC hit on a track negligible? |
node->getHelicity()); | ||
StThreeVectorD upper(tpcHit->positionU().x(),tpcHit->positionU().y(),tpcHit->positionU().z()); | ||
StThreeVectorD lower(tpcHit->positionL().x(),tpcHit->positionL().y(),tpcHit->positionL().z()); | ||
StThreeVectorD middle(tpcHit->position().x(),tpcHit->position().y(),tpcHit->position().z()); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This variable is unused
StThreeVectorD middle(tpcHit->position().x(),tpcHit->position().y(),tpcHit->position().z()); |
A quick grep shows that StTpcHit::setdX is called from StdEdxY2Maker. I presume that dEdxY2 runs after StiStEventFiller. So the dX values will be overwritten when it runs. Don't we need to be able to shut this new calculation off in order to maintain the ability to reproduce old results? |
Yuri, do you have any response to Jason's question? And regarding the time consumption raised by Gene, have you tested it also? Thanks |
Was there any response to these questions? There was also an unaddressed comment about a defined but unused variable in the submitted code. Why would we want to have it in the code? Doesn't it create a warning at compile time? |
No description provided.