<html><head></head><body>Dear Maurik,<div><br></div><div>One issue we have is that we turned off all multiple scattering in gsim, but it appears that some interactions remain. Do we need a collection of gsim cards to turn off all interactions? If so, what would they all be? I would hate to invest much time into this just to find that "XXX" was our problem.</div><div><br></div><div>Paul</div><div><br></div><div><br><div>
<div><div><div>--</div><div>Prof. Paul Eugenio</div><div>Florida State University</div><div>Department of Physics</div><div>Tallahassee, Florida, USA 32306</div><div><br></div><div>(850) 644-2585</div><div><a href="mailto:eugenio@fsu.edu">eugenio@fsu.edu</a></div><div><br></div></div></div><br class="Apple-interchange-newline">
</div>
<br><div><div>On Sep 2, 2010, at 11:57 AM, Maurik Holtrop wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear Craig,<div><br></div><div>You may need to take a closer look at the processes that are part of the GSIM run, which may tell you why that 1% of tracks is lost. It could very simply be one of the normal processes, such as multiple scattering in the target or any other part of the detector, that causes you to loose these tracks. Depending on the setting that you use for the simulation, I do quite agree with the statement: "since these tracks have already been detected by CLAS in real life, after passing through GSIM and reconstruction I should recover close to 100% of them. "<div><div>This statement would only be correct for a geometrical acceptance, and not for one that includes multiple scattering, nuclear processes, and less than 100% efficient detector elements. Add to that a less than 100% efficient tracking algorithm ( due to smearing, timing fluctuations etc) and I can easily see a loss of 1%.</div><div><br></div><div>How does the 99% efficiency of your GSIM run compare with the expected detector efficiency in the phase space you are studying?</div><div><br></div><div>Best regards,</div><div><br></div><div> Maurik </div><div><br></div><div><br></div><div>On Sep 2, 2010, at 11:19 AM, Craig Bookwalter wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>GPP knocking out a track due to holes in the DC is a possibility, except <br>for two things:<br>* I run gpp with a process flag of 0x20 which only processes the <br>tagger--there should be no DC efficiency stuff done.<br>* The events that I am feeding in were already detected by CLAS in real <br>life, which means that if the track went through a DC hole in real life, <br>it wouldn't be reconstructed and thus wouldn't be in my event sample. <br>And if the real-life track has been reconstructed, then it must not have <br>gone through an area of low efficiency in the DC, and thus GPP should <br>not knock it out. Of course the efficiency map could be messed up in <br>this case as well. In any case, when I look at the GSIM output before <br>going into GPP, I don't see hits, so this isn't the problem...thank you <br>though!!<br><br><br><br>Zhiwen Zhao wrote:<br><blockquote type="cite">If you use gpp with DC efficiency ON, would some inefficient DC wires <br></blockquote><blockquote type="cite">lose tracks?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Zhiwen<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On 09/01/2010 04:53 PM, Craig Bookwalter wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><blockquote type="cite">No problem man, this is a hard problem to explain. g12 had the target<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">centered at -90 cm, that's why the vertices are so far back. A better<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">way to understand it is with the pictures I attached--"input.png" is an<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">event that goes in, and "output.png" is what it looks like after going<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">through GSIM and a1c. As you can see, one of the tracks just disappears,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">GSIM generates no hits for it. Not sure how that happens.<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Kijun Park wrote:<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Craig Bookwalter wrote:<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">I'm not clear on what you mean by trigger time--for MC the event<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">start time is always 0 minus the photon propogation time to the event<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">vertex, and it's useful to keep it that way for debugging stuff like<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">this. Along with times from the ST and and TOF, that is the only time<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">I need for reconstruction. I think, at least.<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Kijun Park wrote:<br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Craig Bookwalter wrote:<br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Hi offliners,<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">I am in the middle of a study where I take real events from<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">reconstructed g12 data and feed those four-vectors back through<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">GSIM, with the hypothesis that since these tracks have already been<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">detected by CLAS in real life, after passing through GSIM and<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">reconstruction I should recover close to 100% of them. Currently I<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">am having a problem with GSIM not generating DC hits for some<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">tracks which pass through the DC volumes. When I look at these<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">events in gsim_int, I see that GSIM propogates them along strange<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">trajectories--strange because I have decays and multiple scattering<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">turned off. I have a PART file with a few events exhibiting this<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">problem located here:<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">/work/clas/clasg12/craigb/tmp/gsim_feedback/clas_offline/dc_no_hits.part<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Only the first seven events in this file are examples of what I am<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">having trouble with; I had to add more events to the end to get<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">some BOS buffer to flush and actually write the events to file.<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Also, I created this file to try and get a concentrated sample of<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">my troublesome events, but now when I run my gsim script on this<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">file, some of them have successfully-generated hits, even though I<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">have the RNDM card set. I am guessing this is because RNDM just<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">sets the random seed, which means changing an event's position<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">within a file changes your position within the random number<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">sequence. If there's a more intelligent/convenient way to share<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">this information with everyone, please let me know.<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">So in short, I would appreciate it if you all might take a look at<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">this with me and offer any suggestions as to why GSIM would just<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">not produce DC hits for what looks like 1% of tracks. Also, if I<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">use the "NOSEC 'ALL'" card this no-DC-hits behavior is enhanced by<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">a factor of 2 or so. Thanks in advance.<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">--cb<br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Hi Craig.<br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">I saw you didn't use the -T in your gpp. waht is the "trigbit:<br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">0xffff" at HEAD bank in your dc_no_hits.gpp ?<br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">How do you know trigger time ?<br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Kijun<br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Well, I thought your problem is out of scale of z-vertex (since I am<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">not g12 guy) from you first seven event in that file.<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">If not, could you point out your trouble from your file. Sorry for<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">slow understanding.<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">Kijun<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">_______________________________________________<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Clas_offline mailing list<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite">_______________________________________________<br></blockquote><blockquote type="cite">Clas_offline mailing list<br></blockquote><blockquote type="cite"><a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br></blockquote><blockquote type="cite"><a href="https://mailman.jlab.org/mailman/listinfo/clas_offline">https://mailman.jlab.org/mailman/listinfo/clas_offline</a><br></blockquote><blockquote type="cite"><br></blockquote><br><br>-- <br>-------------------------------------------------------------------------<br> Craig Bookwalter FSU Office: (850) 644 3808<br> Department of Physics JLab Office: (757) 269 5465<br> Florida State University <a href="mailto:craigb@hadron.physics.fsu.edu">craigb@hadron.physics.fsu.edu</a><br> Tallahasse, FL 32306 <a href="mailto:craigb@jlab.org">craigb@jlab.org</a> <br><br><br>"One toke? You poor fool. Just wait till you see those (expletive) bats."<br>-------------------------------------------------------------------------<br><br>_______________________________________________<br>Clas_offline mailing list<br><a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br>https://mailman.jlab.org/mailman/listinfo/clas_offline<br></div></blockquote></div><br></div></div>_______________________________________________<br>Clas_offline mailing list<br><a href="mailto:Clas_offline@jlab.org">Clas_offline@jlab.org</a><br>https://mailman.jlab.org/mailman/listinfo/clas_offline</blockquote></div><br></div></body></html>