 |
 |
Another option is turn off the magnetic field. Takashi -----Original Message----- From: [log in to unmask] [mailto:[log in to unmask]] On Behalf Of Maurik Holtrop Sent: Tuesday, August 27, 2013 8:06 AM To: Hansson Adrian, Per Ola Cc: hps-software Subject: Re: Matching truth particle and sim tracker hit Hi Pelle, I agree with your last statement, which is confirmed by the width of the peaks for layer 2, the width in x is smaller than that for y in layer 1, so there is some other process going on. One thing that may help would be to turn off all multiple scattering in the simulation. I know this possible, just don't know how this is done for SLIC. Perhaps Jeremy or Norman know the trick. That would take one of the variations out of the picture. The other thing you could do is to have a point source beam rather than the spread out beam. If there is a mismatch between the beam source and the expected beam source in the analysis, this could skew the distributions. Best, Maurik On Aug 26, 2013, at 9:05 PM, "Hansson Adrian, Per Ola" <[log in to unmask]> 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 > > > [cid:[log in to unmask]][cid:[log in to unmask]] > > ######################################################################## > 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 > <winmail.dat> ######################################################################## 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 ######################################################################## 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