<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="">
Michael,
<div class=""><br class="">
</div>
<div class="">Just to briefly follow up on this discussion: You mentioned that the kinematic fitter assumes the photons come from the center of the target which is located at the center of CLAS. Does this also apply to the p pi(+)pi(-) final state that, among
others, you have been analyzing? </div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class=""><br class="">
</div>
<div class=""> Volker</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Nov 20, 2018, at 9:38 AM, Michael Kunkel <<a href="mailto:mkunkel@jlab.org" class="">mkunkel@jlab.org</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="auto" class="">Greetings,
<div class=""><br class="">
</div>
<div class="">“<span style="background-color: rgba(255, 255, 255, 0);" class="">Keep in mind that the Kinematic fitter and EC photons assumes the photons come from the center of the target which is located at the center of CLAS”</span></div>
<div class="">Not the reconstruction code.</div>
<div class="">See KinLine.C, I think, for how photons are corrected.</div>
<div class=""><br class="">
</div>
<div class="">Also my runIndex, with run = 10, should be used with photons. It should be good for all particles, as it’s just run 56855 information except EC and CC constants, which are set to run 10.<br class="">
<br class="">
<div class="">
<div style="direction: inherit;" class=""><span style="background-color: rgba(255, 255, 255, 0);" class="">BR<br class="">
MK <br class="">
----------------------------------------<br class="">
Michael C. Kunkel, USMC, PhD<br class="">
Forschungszentrum Jülich<br class="">
Nuclear Physics Institute and Juelich Center for Hadron Physics<br class="">
Experimental Hadron Structure (IKP-1)</span></div>
<div style="direction: inherit;" class=""><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.fz-2Djuelich.de_ikp&d=DwMGaQ&c=HPMtquzZjKY31rtkyGRFnQ&r=JNNs9K-auSa6rh56y2esuQ&m=s2hUdBUQl0aqeNaLGlKm_jtnYmV22fvK1SyT9KuO8cc&s=ysUlUJoNbpnHFwNUey1TLOSnu6p9TK1sWPcAA82Nrf0&e=" class="">www.fz-juelich.de/ikp</a> </div>
</div>
<div class=""><br class="">
On Nov 20, 2018, at 15:30, Andrea Celentano <<a href="mailto:andrea.celentano@ge.infn.it" class="">andrea.celentano@ge.infn.it</a>> wrote:<br class="">
<br class="">
</div>
<blockquote type="cite" class="">
<div class="">HI Michael,
<div class="">thanks for your detailed email.</div>
<div class=""><br class="">
</div>
<div class="">Can you explain me why the EC assumes the photons come from the target center at z=0? If I look at the
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__makePART.cc&d=DwMGaQ&c=HPMtquzZjKY31rtkyGRFnQ&r=JNNs9K-auSa6rh56y2esuQ&m=s2hUdBUQl0aqeNaLGlKm_jtnYmV22fvK1SyT9KuO8cc&s=MXu6m2_nNxAMOmhAOiJ9y3kwlQn9F4KQOAP84eBDhEs&e=" class="">
makePART.cc</a> file in the PID package, function GammaID(), I see:</div>
<div class=""><br class="">
</div>
<div class="">
<pre style="word-wrap: break-word; white-space: pre-wrap;" class="">int GammaID(bid_t * tbid, int tbid_index, part_t * part){
int IFail = 0;
vector3_t zero3 = {0.0,0.0,0.0};
float gamma=0;
float p, mag;
clasECHB_t *ECHB = (clasECHB_t *)getBank(&bcs_, "ECHB");
echb_t *echb;
clasEC1R_t *EC1R = getGroup(&bcs_,"EC1R",tbid->sec);
clasMVRT_t *MVRT = NULL;
/*look for MVRT bank. liji*/
if(MVRT = (clasMVRT_t *)getBank(&bcs_, "MVRT")){
zero3 = MVRT->mvrt[0].vert;
}else{
/*get target center position form TGEO bank. liji*/
zero3.z = target_z();
}</pre>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">It seems that, in case an MVRT bank exists, then the “zero3” is taken from it - aka the reaction vertex. In my analysis, I take events with 3 charged particles as seen by PART (p, pip, pim), hence MVRT is defined.</div>
<div class=""><br class="">
</div>
<div class="">Anyhow, we can discuss about this at the next g12 meeting. </div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class="">Andrea</div>
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On Nov 20, 2018, at 15:08, Michael C. Kunkel <<a href="mailto:mkunkel@jlab.org" class="">mkunkel@jlab.org</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">Greetings Andrea,<br class="">
<br class="">
I am cc'ing g12 this message as it pertains to g12 data analysis.<br class="">
<br class="">
Back in 2010, g12 investigated the photons detected in the EC. Diane Schott and I independently worked out a set of corrections for the MC that was with .1% of each other. However, these corrections never made it into the g12 analysis procedures because at
the time, no one was using EC photons. At the time, the kinematic fitter was not yet mature enough nor understood on the g12 dataset, so we hadn't used it for this study.<br class="">
Later, after the fitter was understood, g12 wanted to perform the EC photon corrections, but never got around to it.<br class="">
Keep in mind that the Kinematic fitter and EC photons assumes the photons come from the center of the target which is located at the center of CLAS. In g12 this is not that case at the target as not at CLAS center, so if you pursue this avenue, you should consider
this aspect.<br class="">
<br class="">
However, for the rudimentary corrections Diane and I derived, I remember the procedure and it was simply plotting the ratio of the reconstructed photon energy with the generated energy (for that event), vs the generated energy. You should see an exponential
decay, very much like the tail of a Poisson Distribution.<br class="">
<br class="">
The data correction was different, it was derived from gamma p -> p pi0 where pi0->e+e-gamma.<br class="">
Here, I cut the MxP -> pi0 and by using several 4-vector conservations, I was able to get a correction factor for gamma to create e+e-gamma = pi0.<br class="">
I can try to hunt the results and study, but I think it's easier to just reproduce them if you wanted to.<br class="">
<br class="">
Whichever you decide, I can help in either way.<br class="">
I know most of this information is somewhere on the g12 wiki, either on my page, which also links to the g7 page, or maybe Dianes page has it.<br class="">
<br class="">
BR<br class="">
MK<br class="">
----------------------------------------<br class="">
Michael C. Kunkel, PhD<br class="">
Forschungszentrum Jülich<br class="">
Nuclear Physics Institute and Juelich Center for Hadron Physics<br class="">
Experimental Hadron Structure (IKP-1)<br class="">
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.fz-2Djuelich.de_ikp&d=DwMGaQ&c=HPMtquzZjKY31rtkyGRFnQ&r=JNNs9K-auSa6rh56y2esuQ&m=s2hUdBUQl0aqeNaLGlKm_jtnYmV22fvK1SyT9KuO8cc&s=ysUlUJoNbpnHFwNUey1TLOSnu6p9TK1sWPcAA82Nrf0&e=" class="">www.fz-juelich.de/ikp</a><br class="">
<br class="">
On 11/20/18 12:02 PM, Andrea Celentano wrote:<br class="">
<blockquote type="cite" class="">Dear Michael,<br class="">
I write you to ask if in g12 anyone has already looked at the corrections to apply to photons measured in EC.<br class="">
<br class="">
I see that, for g11, a robust procedure has been set up by Dustin (see: <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__misportal.jlab.org_ul_Physics_Hall-2DB_clas_viewFile.cfm_2011-2D004.pdf-3FdocumentId-3D632&d=DwMGaQ&c=HPMtquzZjKY31rtkyGRFnQ&r=JNNs9K-auSa6rh56y2esuQ&m=s2hUdBUQl0aqeNaLGlKm_jtnYmV22fvK1SyT9KuO8cc&s=DlK6gnFclpelLAThM0B2acNbs5R7pIHEh2iwGabyshE&e=" class="">
https://misportal.jlab.org/ul/Physics/Hall-B/clas/viewFile.cfm/2011-004.pdf?documentId=632</a>). This procedure is based on kinematic fit of the gamma p -> p pi+ pi- gamma reaction, correlating the measured gamma 4-momentum with the one obtained to the kinfit
of p pi+ pi- with missing gamma hypothesis.<br class="">
<br class="">
I'd like to repeat the same exercise for g12, benefiting from the kinematic fit code you sent around as part of your efficiency studies. However, before doing so, I'd like to check that this has not yet been done. I'd like to join next g12 meeting (next week)
to discuss this, but to gain some time I write you in advance to know this.<br class="">
<br class="">
Thanks,<br class="">
<br class="">
Bests,<br class="">
Andrea<br class="">
<br class="">
</blockquote>
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
</div>
_______________________________________________<br class="">
G12 mailing list<br class="">
<a href="mailto:G12@jlab.org" class="">G12@jlab.org</a><br class="">
https://urldefense.proofpoint.com/v2/url?u=https-3A__mailman.jlab.org_mailman_listinfo_g12&d=DwICAg&c=HPMtquzZjKY31rtkyGRFnQ&r=JNNs9K-auSa6rh56y2esuQ&m=s2hUdBUQl0aqeNaLGlKm_jtnYmV22fvK1SyT9KuO8cc&s=9IhTt9kTHMPkA0lWdyHRKczsjzj4rKw2FIHlWf9geWQ&e=<br class="">
</div>
</blockquote>
</div>
<div class="">
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;">
<br class="">
</div>
</div>
</div>
</body>
</html>