<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Andrei,<div class=""><br class=""></div><div class="">Ok, that rules that out then. There are definitely edges in your plots at +- 2ns. This would point to some RF bunch selection issue. I’d follow Paul’s suggestion. If you want to quickly check the quality of the timing calibrations, run the plugin “HLDetectorTiming” along with yours. Of particular interest would be the histograms in the SC_Target_RF_Compare and the TRACKING directory. These seem OK in the last version of the offline monitoring, but if these have shifted, it could throw off the selection.</div><div class=""><br class=""></div><div class="">-Mike</div><div class=""><br class=""></div><div class=""><img height="430" width="640" apple-width="yes" apple-height="yes" apple-inline="yes" id="DB174E69-C412-42E8-9B08-441E91BEE44C" src="cid:F2F98F61-F1AC-4485-80BF-88800BF3C26E@phys.cmu.edu" class=""></div><div class=""><div class="">
<div class=""><div class="">--</div><div class="">Michael Staib</div><div class="">Graduate Student, Dept. of Physics</div><div class="">Carnegie Mellon University</div><div class=""><a href="mailto:mstaib@cmu.edu" class="">mstaib@cmu.edu</a></div><div class="">phone: 412-268-2983</div></div><div class=""><br class=""></div><br class="Apple-interchange-newline">
</div>
<br class=""><div><blockquote type="cite" class=""><div class="">On Apr 27, 2016, at 5:27 PM, <a href="mailto:semenov@jlab.org" class="">semenov@jlab.org</a> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Michael and Paul:<br class=""><br class="">In my plots, I have the only tracks that are matched with the showers in<br class="">BCAL; t1_detector values for the events I selected are always equal 4<br class="">(that corresponds to BCAL, I hope).<br class=""><br class="">Paul: I'll produce "beta-vs-momentum" 2-dim plots. (Probably,<br class="">"mass^2-vs-momentum" plots will be interesting also.)<br class=""><br class="">Thank you,<br class="">Andrei<br class=""><br class=""><br class=""><br class=""><br class=""><blockquote type="cite" class="">Right. You should probably make separate plots for whether m_t1_detector<br class="">is either SYS_BCAL, SYS_TOF, or SYS_FCAL.<br class=""><br class="">The t0 one should hopefully always be SYS_RF.<br class=""><br class=""> - Paul<br class=""><br class="">On Apr 27, 2016, at 3:58 PM, Michael Staib <<a href="mailto:mstaib@andrew.cmu.edu" class="">mstaib@andrew.cmu.edu</a>> wrote:<br class=""><br class=""><blockquote type="cite" class="">Hi Andrei,<br class=""><br class="">Have you checked that every track has the same m_t0_detector and<br class="">m_t1_detector? I wouldnât be surprised if there is some lingering<br class="">differences in the timing between the possible t0 and t1 detectors here.<br class=""><br class="">There is no plan to improve the calibrations in that old data, but feel<br class="">free if you find that it helps.<br class=""><br class="">--<br class="">Michael Staib<br class="">Graduate Student, Dept. of Physics<br class="">Carnegie Mellon University<br class=""><a href="mailto:mstaib@cmu.edu" class="">mstaib@cmu.edu</a><br class="">phone: 412-268-2983<br class=""><br class=""><br class=""><br class=""><blockquote type="cite" class="">On Apr 27, 2016, at 3:42 PM, semenov@jlab.org wrote:<br class=""><br class="">Paul:<br class=""><br class="">The requested 2-dim plots for pions+/-, protons and K+ are in the<br class="">attachment.<br class=""><br class="">Though the statistics is not very high, the 2-horn structure is still<br class="">visible for pions as well as a lot of "too-high-speed" events.<br class=""><br class="">Does it help?<br class=""><br class="">Thank you,<br class="">Andrei<br class=""><br class=""><br class=""><br class=""><blockquote type="cite" class="">Paul:<br class=""><br class="">You suggestion to use the google group is not good for me because I do<br class="">not<br class="">have google account, and I do not want create it because of the google<br class="">demand to share my cell phone number. I do believe that we are working<br class="">for<br class="">JLab (not for the Google corporation), so we should use JLab computer<br class="">resources. I would suggest to use "halld-offline" e-mail list. Sorry<br class="">and<br class="">thank you for understanding. (Alternatively, some shared "gluex"<br class="">account<br class="">for google might be a solution, but in this case it will be not<br class="">obvious<br class="">why's actually the author of the specific communication.)<br class=""><br class="">Talking of "t0 is set to the chosen RF time, propagated to the<br class="">vertex-z of<br class="">the track": Am I correct that in this case, "t0" should be exactly the<br class="">same as "T" value from the DVertex class? (If yes, it's strange<br class="">because<br class="">"t0" is close to the "T" from DVertex, but not exactly the same.)<br class=""><br class="">By the way, how the "t1" value is produced? Is it some estimation from<br class="">the<br class="">chambers timing, or it came from the "t1_detector" (which is in my<br class="">case<br class="">always =4; I hope that it corresponds to BCAL). If the last is<br class="">correct,<br class="">how good this "t1" value is? (Meaning time-walk corrections, channels<br class="">alignment etc.)<br class=""><br class="">My plot contains all PIDs.<br class=""><br class="">Talking of multiple hypotheses, I do believe that "PID" in this class<br class="">is<br class="">not a vector but one number. Does it mean that<br class="">"DChargedTrackHypothesis"<br class="">vector size is always bigger than the size of the correspondent<br class="">"DTrackTimeBased" vector (that contains one record per one "real"<br class="">track, I<br class="">hope)?<br class=""><br class="">I do not agree that the right-column plots do not help. In the case of<br class="">multiple hypotheses (viz., unknown particle mass), the right-column<br class="">plot<br class="">(which presents the TOF minus the time needed to travel the pathLength<br class="">distance with the speed-of-light) shows that we have a lot of tracks<br class="">with<br class="">TOF that suggests that the particle travels much faster than<br class="">speed-of-light (negative values in the histogram). If we do not<br class="">consider<br class="">tachyons seriously, it means that the "TOF" value has some problems.<br class=""><br class="">I'll produce the "(TOF-L/(beta*C)) vs momentum" 2-dim plots you<br class="">suggested,<br class="">but it will require a few hours on our computers; I hope that it will<br class="">be<br class="">ready tomorrow. (I'm not quite sure why you want to see "momentum<br class="">dimension": most probably, it will be just horizontal bands. But OK,<br class="">no<br class="">problem. Note that I can not just plot "TOF vs momentum" because<br class="">different<br class="">pathLengths.)<br class=""><br class="">Thank you,<br class="">Andrei<br class=""><br class=""><br class=""><br class=""><br class=""><blockquote type="cite" class=""><br class="">On Apr 25, 2016, at 8:26 PM, pmatt@jlab.org wrote:<br class=""><br class="">You may know this, but for the record, when you call<br class="">DKinematicData::TOF(), it is subtracting t1 - t0. In the<br class="">DChargedTrackHypothesis factory, t0 is set to the chosen RF time,<br class="">propagated to the vertex-z of the track.<br class=""><br class="">Are you histogramming this for all PID hypotheses? Or are you always<br class="">picking the ones with the same PID? (e.g. proton) Remember, each<br class="">physical track has multiple hypotheses, for different PIDs. If you<br class="">chose<br class="">all hypotheses, then only one can be correct, and you will see funny<br class="">structures. Although yes, they would not give rise to the<br class="">right-column<br class="">plots. However, the tracks certainly arenÂt going at beta = 1, so<br class="">the<br class="">right plots donÂt help much.<br class=""><br class="">You should have separate plots for each PID, and then you should make<br class="">this<br class="">a 2D plot vs. track momentum, and you will see a mass dependence that<br class="">will give rise to the structures. You should see distinct pion,<br class="">proton,<br class="">etc. bands.<br class=""><br class="">- Paul<br class=""><br class="">************************************************************************<br class=""><br class="">We now have a new method of asking software questions, a google<br class="">group:<br class="">gluex-software@googlegroups.com<br class=""><br class="">I am about to forward this message to there, and then will try to<br class="">answer<br class="">it there when I have time. If you are not part of the group yet,<br class="">please<br class="">sign up:<br class=""><br class="">https://groups.google.com/forum/#!forum/gluex-software<br class=""><br class="">- Paul<br class=""><br class="">On Apr 25, 2016, at 6:46 PM, semenov@jlab.org wrote:<br class=""><br class=""><blockquote type="cite" class="">Paul, Mark and Simon:<br class=""><br class="">I'm trying to use the information from DChargedTrackHypothesis, and<br class="">I<br class="">see<br class="">the things I do not understand. If I plot "TOF" (from the class)<br class="">minus<br class="">the<br class="">time-of-flight "pathLength*sqrt(mass*mass+pmag*pmag)/pmag/c" (that I<br class="">calculated from the variables in the class) for the charged tracks<br class="">that<br class="">are matched with the showers in BCAL, I see some interesting<br class="">structure<br class="">(see the top-left panel in the attached file) instead of expected<br class="">centered-at-zero peak with the RMS of the order of the uncertainty<br class="">on<br class="">TOF<br class="">(0.4-0.5ns). The bottom-left panel shows the same-variable histogram<br class="">with<br class="">the requirement to have 3 or more charged tracks in the event (to be<br class="">sure<br class="">that we have a good start time).<br class=""><br class="">For the "mass", I used the "PID" hypothesis in the class. To be sure<br class="">that<br class="">this structure is not affected much by PID hypothesis, I plotted<br class="">also<br class="">the<br class="">histograms for "TOF-pathLength/c" (see the right column in the<br class="">file);<br class="">everything that is visibly below zero are "tachyons" :) that means<br class="">that<br class="">the problem is probably with the "TOF" variable which is the<br class="">difference<br class="">between "t1" (stop time) and "t0" (start time) variables. "t0" is<br class="">pretty<br class="">close to the "vertex" time "T" from DVertex class.<br class=""><br class="">As you can see from the plot, the peaks in the structure are not<br class="">separated<br class="">with 4-ns intervals, so I would not blame "wrong-RF-bucket" reason.<br class=""><br class="">For this plot, I used run 3180, and my version of sim-recon is 1.10<br class="">.<br class=""><br class="">Most probably, I just using something incorrectly, and you already<br class="">have<br class="">a<br class="">good explanation for the plots I see. Could you please help me?<br class=""><br class="">Thank you,<br class="">Andrei<tdiff-160423.pdf><br class=""></blockquote><br class=""><br class=""></blockquote>_______________________________________________<br class="">Halld-offline mailing list<br class="">Halld-offline@jlab.org<br class="">https://mailman.jlab.org/mailman/listinfo/halld-offline<br class=""></blockquote><run3180-2d.pdf>_______________________________________________<br class="">Halld-offline mailing list<br class="">Halld-offline@jlab.org<br class="">https://mailman.jlab.org/mailman/listinfo/halld-offline<br class=""></blockquote><br class=""></blockquote><br class=""><br class=""></blockquote><br class=""><br class=""></div></div></blockquote></div><br class=""></div></body></html>