- genevb's home page
- Posts
- 2024
- 2023
- 2022
- September (1)
- 2021
- 2020
- 2019
- December (1)
- October (4)
- September (2)
- August (6)
- July (1)
- June (2)
- May (4)
- April (2)
- March (3)
- February (3)
- 2018
- 2017
- December (1)
- October (3)
- September (1)
- August (1)
- July (2)
- June (2)
- April (2)
- March (2)
- February (1)
- 2016
- November (2)
- September (1)
- August (2)
- July (1)
- June (2)
- May (2)
- April (1)
- March (5)
- February (2)
- January (1)
- 2015
- December (1)
- October (1)
- September (2)
- June (1)
- May (2)
- April (2)
- March (3)
- February (1)
- January (3)
- 2014
- December (2)
- October (2)
- September (2)
- August (3)
- July (2)
- June (2)
- May (2)
- April (9)
- March (2)
- February (2)
- January (1)
- 2013
- December (5)
- October (3)
- September (3)
- August (1)
- July (1)
- May (4)
- April (4)
- March (7)
- February (1)
- January (2)
- 2012
- December (2)
- November (6)
- October (2)
- September (3)
- August (7)
- July (2)
- June (1)
- May (3)
- April (1)
- March (2)
- February (1)
- 2011
- November (1)
- October (1)
- September (4)
- August (2)
- July (4)
- June (3)
- May (4)
- April (9)
- March (5)
- February (6)
- January (3)
- 2010
- December (3)
- November (6)
- October (3)
- September (1)
- August (5)
- July (1)
- June (4)
- May (1)
- April (2)
- March (2)
- February (4)
- January (2)
- 2009
- November (1)
- October (2)
- September (6)
- August (4)
- July (4)
- June (3)
- May (5)
- April (5)
- March (3)
- February (1)
- 2008
- 2005
- October (1)
- My blog
- Post new blog entry
- All blogs
NHitsPoss systematics
Updated on Thu, 2017-03-23 13:18. Originally created by genevb on 2017-03-23 13:13.
Victor noted that NHitsPoss (number of hits possible on a track) is incorrectly calculated. I've observed that he is correct in that possible hits at radii larger than the last hit included on a track are ignored for primary tracks and some global tracks. However, possible hits at radii smaller than the first hit included on a track appear to be counted as there is no obvious correlation. This is demonstrated in plots made from tracks whose first hit is at |z|<50 cm, |η|<0.5, and 0.2 < pT < 5.0 GeV/c below.
(note: first hit on primary tracks is always at the primary vertex, so it isn't helpful to make a plot with it)
While the plot for primary tracks is quite clear in its correlation, the plot using last hit of global tracks appears to have multiple components:
The above plots can be generated from any MuDst by the following bit of code:
-Gene
(note: first hit on primary tracks is always at the primary vertex, so it isn't helpful to make a plot with it)
While the plot for primary tracks is quite clear in its correlation, the plot using last hit of global tracks appears to have multiple components:
- A correlated band (similar to the primary tracks but shifted down by 1 hit accounting for the primary vertex)
- A band at large radius but smaller number of possible hits
- A group with up to 45 possible hits despite the last hit not being all the way out
NHitsPoss vs. radius of the first hit included on global tracks |
NHitsPoss vs. radius of the last hit included on global tracks |
NHitsPoss vs. radius of the last hit included on primary tracks |
---|---|---|
The above plots can be generated from any MuDst by the following bit of code:
MuDst.SetAlias("firstHitRadius","sqrt(sq(GlobalTracks.mFirstPoint.mX1)+sq(GlobalTracks.mFirstPoint.mX2))"); MuDst.SetAlias("lastHitRadius","sqrt(sq(GlobalTracks.mLastPoint.mX1)+sq(GlobalTracks.mLastPoint.mX2))"); MuDst.SetAlias("plastHitRadius","sqrt(sq(PrimaryTracks.mLastPoint.mX1)+sq(PrimaryTracks.mLastPoint.mX2))"); TCut cg1 = "abs(GlobalTracks.mFirstPoint.mX3)<50&&abs(GlobalTracks.mEta)<0.5&&GlobalTracks.mPt>0.2&&GlobalTracks.mPt<5"; TCut cp1 = "abs(PrimaryTracks.mFirstPoint.mX3)<50&&abs(PrimaryTracks.mEta)<0.5&&PrimaryTracks.mPt>0.2&&PrimaryTracks.mPt<5"; MuDst.Draw("GlobalTracks.mNHitsPoss:firstHitRadius>>fhg(50,50,200,50,0.5,50.5)",cg1,"zcol",100); gPad->SetLogz(); MuDst.Draw("GlobalTracks.mNHitsPoss:lastHitRadius>>lhg(50,50,200,50,0.5,50.5)",cg1,"zcol",100); MuDst.Draw("PrimaryTracks.mNHitsPoss:plastHitRadius>>lhp(50,50,200,50,0.5,50.5)",cp1,"zcol",100);
-Gene
»
- genevb's blog
- Login or register to post comments