<div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div>Dear Mauri and dear all,</div><div><br></div><div>in the mail below you can find a summary of the procedure I developed to implement the cc matching cuts on e1f.</div><div><br></div><div>Let me know if you have any question.</div><div><br></div><div>Best regards,</div><div><br></div><div>silvia</div><div><div class="h5"><div><br></div><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Silvia Pisano</b> <span dir="ltr"><<a href="mailto:silvia.pisano@lnf.infn.it" target="_blank">silvia.pisano@lnf.infn.it</a>></span><br>Date: 8 April 2016 at 17:00<br>Subject: cc cuts<br>To: Harut Avagyan <<a href="mailto:avakian@jlab.org" target="_blank">avakian@jlab.org</a>>, Nathan Harrison <<a href="mailto:nathanh@jlab.org" target="_blank">nathanh@jlab.org</a>>, Aram Movsisyan <<a href="mailto:maram@mail.desy.de" target="_blank">maram@mail.desy.de</a>><br><br><br><div dir="ltr">Dear all,
<div><br></div><div>after some weeks I managed to implement and test the Osipenko's cuts on phi (from Nathan), theta and time. For theta and time I took definitions (and inspiration!) from what Aram did (thanks a lot for your help with it, Aram!), and developed an automatic fit procedure to define the cuts.</div><div><br></div><div>I wrote a set of function that I will put on github, so that if someone has to implement those cuts can use it (Nathan, if you will be asked to do it for the analysis review just let me know and I will help you with it).</div><div><br></div><div>For the theta fit I used a polynomial background plus a gaussian signal, while for time I used a polynomial background plus N gaussians in order to account for multiple reflection peaks.</div><div><br></div><div>This is an example for theta fits:</div><div><br></div><div><a href="http://www.lnf.infn.it/~pisanos/dh/cc_theta_cuts_Sect3.pdf" target="_blank">http://www.lnf.infn.it/~pisano<wbr>s/dh/cc_theta_cuts_Sect3.pdf</a><br></div><div><br></div><div>(other sectors and in the same page, changing the corresponding number)</div><div><br></div><div>and for time:</div><div><br></div><div><a href="http://www.lnf.infn.it/~pisanos/dh/cc_time_cuts_Sect3.pdf" target="_blank">http://www.lnf.infn.it/~pisano<wbr>s/dh/cc_time_cuts_Sect3.pdf</a><br></div><div><br></div><div>The number of peaks is determined through the root class TSpectrum, and it seems to work in a pretty good way (it failures just in few cases, if any).</div><div><br></div><div>I followed the procedure delined by Wes, cutting around three sigma from the mean for theta and applying just a lower cut on time, so to account for the multiple relfection. However, the combined action of the three cuts - theta, time and phi - seems to be quite effective in selecting just the gaussian signal, as you can see from, e.g.,:</div><div><br></div><div><a href="http://www.lnf.infn.it/~pisanos/dh/cc_time_Sect_4_check.pdf" target="_blank">http://www.lnf.infn.it/~pisano<wbr>s/dh/cc_time_Sect_4_check.pdf</a><br></div><div><br></div><div>for time and in</div><div><br></div><div><a href="http://www.lnf.infn.it/~pisanos/dh/cc_theta_Sect_4_check.pdf" target="_blank">http://www.lnf.infn.it/~pisano<wbr>s/dh/cc_theta_Sect_4_check.pdf</a><br></div><div><br></div><div>for theta, where the empty histogram shows the variable distribution before the cuts, while the green, solid filled histogram represents the distribution after the cuts.</div><div><br></div><div>The global effects on nphe is</div><div><br></div><div><a href="http://www.lnf.infn.it/~pisanos/dh/cc_nphe_check.pdf" target="_blank">http://www.lnf.infn.it/~pisano<wbr>s/dh/cc_nphe_check.pdf</a><br></div><div><br></div><div>I am now going to redo the higher-level analysis to see what happens to the observables, and I will keep you updated on it.</div><div><br></div><div>Thanks again to all of you for your help,</div><div><br></div><div>silvia</div><div><br></div><div><br></div></div></div></div></div></div></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"> "Life was just chemistry, <br>  and chemistry was just physics." - Jonah Lehrer</div><div dir="ltr"><br></div></div>
</div>