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

StJetMaker: add idTruth to StJetTrack #152

Merged
merged 3 commits into from
Sep 24, 2021

Conversation

veprbl
Copy link
Member

@veprbl veprbl commented Sep 21, 2021

This should be useful for studies of simulation (e.g. select real tracks vs pile-up), to match tracks to the Pythia particles. For matching to secondary particles an external geant.root or minimc.root would be required.

Copy link
Member

@plexoos plexoos left a comment

Choose a reason for hiding this comment

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

Looks like a straightforward addition to me

Raise the version number by 1 for StjTrack
@@ -57,7 +57,7 @@ class StjTrack : public TObject {
double nSigmaTofProton;
double nSigmaTofElectron;

ClassDef(StjTrack,5);
ClassDef(StjTrack,6);
Copy link
Member Author

Choose a reason for hiding this comment

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

This should never appear on disk, so really should be a version=0

Copy link
Member Author

Choose a reason for hiding this comment

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

but changing that is not the point of this PR

@zlchang
Copy link
Member

zlchang commented Sep 22, 2021 via email

@veprbl
Copy link
Member Author

veprbl commented Sep 22, 2021

This is a very old code maybe from Tai Sakuma. This has never been used since 2009 analysis, but I remembered before the 2009 analysis tracks are saved in this StjTrack format in the jet tree.

It is old, but looks like the previous tree format didn't use streamers
https://github.com/star-bnl/star-sw/blob/e0f02802d3ac24161117286395bf64f007cae923/StRoot/StJetMaker/tree/StjTrackListWriter.cxx

@zlchang
Copy link
Member

zlchang commented Sep 22, 2021 via email

@zlchang
Copy link
Member

zlchang commented Sep 22, 2021

@veprbl do you want to send an email to spin pwg to let people know that you have updated the code or do you want me to do it? My concern is that now people start using code from dev, today they submit job with one code and the next day, they would submit the same job without realizing using a new code. Once you sent the email, what about we wait for one or two days on the list and see if people have objections? If none, then we will merge the branch.

@veprbl
Copy link
Member Author

veprbl commented Sep 22, 2021

Regarding merging this, I'm fine with whatever you decide, @zlchang

@zlchang
Copy link
Member

zlchang commented Sep 22, 2021

Regarding merging this, I'm fine with whatever you decide, @zlchang

What don't you write an email to the pwg, just explain what the two variables are (idtrue and qatrue) and their potential impact to your analysis? If no objections, we will merge the code in 24 hours.

@BassamAboona
Copy link
Contributor

Hello All,

I was wondering if these new track properties need to be added to StUEMaker2009.cxx to make sure that the tracks from the UE region trees have access to these properties and are consistent with the off-axis cone and jet trees?

@zlchang
Copy link
Member

zlchang commented Sep 23, 2021 via email

@veprbl
Copy link
Member Author

veprbl commented Sep 24, 2021

The region code was never used and contains an unfixed bug. The PicoDst trees are better than it in many respects. I will pass on implementing this suggestion.

@zlchang
Copy link
Member

zlchang commented Sep 24, 2021 via email

@BassamAboona
Copy link
Contributor

BassamAboona commented Sep 24, 2021 via email

@veprbl
Copy link
Member Author

veprbl commented Sep 24, 2021

For my analysis, I do produce and use the use region trees.

I advise that you don't. It's a waste of disk space to have those.

I am curious, what is the bug you are referring to in the UE region codes, Dmitry?

Actually, there are two bugs, both unresolved:
https://drupal.star.bnl.gov/STAR/system/files/jet_meeting_2017_02_01_v2.pdf
There was also bug in the only analysis that used those trees:
https://drupal.star.bnl.gov/STAR/blog/veprbl/run12-pp200-region-ue-plots-qa
hence, in the end, there is no physics yield from this software

@zlchang
Copy link
Member

zlchang commented Sep 24, 2021 via email

@veprbl
Copy link
Member Author

veprbl commented Sep 24, 2021

Pass

@zlchang
Copy link
Member

zlchang commented Sep 24, 2021 via email

Comment on lines +111 to +112
int mIdTruth;
int mQaTruth;
Copy link
Member Author

Choose a reason for hiding this comment

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

Suggested change
int mIdTruth;
int mQaTruth;
unsigned short mIdTruth;
unsigned short mQaTruth;

https://www.star.bnl.gov/webdata/dox/html/classStMuTrack.html
StMuTrack::mIdTruth is UShort_t while StMuTrack::idTruth() is Int_t

Comment on lines +238 to +240

track->mIdTruth = t.idTruth;
track->mQaTruth = t.qaTruth;
Copy link
Member Author

Choose a reason for hiding this comment

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

If you knew that this was what you wanted, just post a patch next time. Please.

@zlchang
Copy link
Member

zlchang commented Sep 24, 2021 via email

@zlchang
Copy link
Member

zlchang commented Sep 24, 2021 via email

@zlchang
Copy link
Member

zlchang commented Sep 24, 2021

Are we good with recent commits? Keeping the idtruth and qatruth as int is fine to me because int has as twice or equal much storage bits as short.

@veprbl
Copy link
Member Author

veprbl commented Sep 24, 2021

You also need to modify stjtracklist

I don't think we need to have narrow type for stjtrack, it's not written to disk.

@veprbl
Copy link
Member Author

veprbl commented Sep 24, 2021

Keeping int for stjettrack should be also fine.

@zlchang zlchang merged commit 0151294 into star-bnl:main Sep 24, 2021
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.

None yet

4 participants