 |
 |
Hi All, me and Sho found it after a long debug session last night. It was a bug in lcsim, in the org.lcsim.fit.helicaltrack.HelixParamCalculator. Another one of those famous tracking vs. detector coordinates issues where the helix z0 got the origin from the MC particle in detector x instead of detector z. Since there was no offset from 0 (with beamspot equal to zero) in x it meant that the helix always started out from 0 in the non-bend direction. I've committed the fix to lcsim. Looking at it more carefully we think it's used in only one more case where it can make a difference (Matt can find it). I will change the rest of the suspicious cases and commit. Attached are the same residuals for a case where there is 0 beamspot width, 0 offset. It looks a lot better. Thanks for the suggestions. /Pelle [cid:1E80CD58-AE00-405B-8428-ADB8C73370C8] [cid:CF4265AF-165F-4473-B116-497E04888176] On Aug 26, 2013, at 6:22 PM, Hansson Adrian, Per Ola wrote: aida file attached. /Pelle On Aug 26, 2013, at 6:05 PM, Hansson Adrian, Per Ola wrote: > Dear Experts, > > I'm doing, what I thought would be, a simple cross-check in my tracking developments that the first SimTrackerHit in the active silicon is at the correct place as predicted by the truth particle momentum . > > This is what I'm doing: > - Find SimTrackerHit of first sensor layer associated with a reconstructed track > - Find the MC particle associated with the SimTrackerHit > - Convert the MC particle momentum, and position, vectors into a helix > - Use standard tools to find the path length along the truth particle helix to the z-position of the SimTracker hit > - Use standard tools to find the 3D point on the truth particle helix that corresponds to that path length. > > The plots attached shows the residual in bend plane and non-bendplane for each layer (starting from top left). The aida file is also attached. Didn't put proper axis labels but the main point is that the residual for layer 1 (top left is close to zero in the bend plane but has a shift of about 30um and also RMS~40um in the non-bend plane. > > Since there is no material in between, and I verified that the truth particle is the correct A' daughter, we would expect this to be zero in the first layer. > > Anyone has any idea on why I get this result? > > Normal explained to me that the position of the SimTrackerHit is located in the "middle" of the active silicon. Since I'm propagating to the SimTrackerHit 3D position this should not matter but even if I got that wrong (if the SimTrackerHit z-position was somehow wrong) the mean slope is about 25mrad which means that I would get corrections of about 0.025/0.32mm=8um; not enough, I think, to be responsible for this? > > It's also strange that it seems to be ok in the bend plane. > > /Pelle > > > <Screen Shot 2013-08-26 at 5.44.28 PM.png><Screen Shot 2013-08-26 at 5.44.01 PM.png> <gbldriverplots_aprime40MeV_2.2GeV.aida> ######################################################################## Use REPLY-ALL to reply to list To unsubscribe from the HPS-SOFTWARE list, click the following link: https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1