<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Greetings,<br>
<br>
If I understood Eugene correctly, then what I did was plot the beam
energy in 6 MeV bins for all events from the pi0.<br>
I also did this for the Glux, by creating a histogram of the output
of Johann's glfux calculator.<br>
<br>
I am attaching the results, but saving judgment until my mistakes
are pointed out.<br>
<br>
BR<br>
MK<br>
<br>
<br>
<div class="moz-cite-prefix">On 3/28/14 4:27 PM, Eugene Pasyuk
wrote:<br>
</div>
<blockquote
cite="mid:1729877211.136211.1396038429654.JavaMail.root@jlab.org"
type="cite">
<style type="text/css">p { margin: 0; }</style>
<div style="font-family: arial,helvetica,sans-serif; font-size:
12pt; color: #000000">reaction does not matter, pi0 or anything
else<br>
<br>
<div><span name="x"></span>-Eugene<span name="x"></span><br>
</div>
<br>
<hr id="zwchr">
<blockquote style="border-left:2px solid
#1010FF;margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>From:
</b>"Michael C. Kunkel" <a class="moz-txt-link-rfc2396E" href="mailto:mkunkel@jlab.org"><mkunkel@jlab.org></a><br>
<b>To: </b>"Eugene Pasyuk" <a class="moz-txt-link-rfc2396E" href="mailto:pasyuk@jlab.org"><pasyuk@jlab.org></a><br>
<b>Cc: </b>"g12" <a class="moz-txt-link-rfc2396E" href="mailto:g12@jlab.org"><g12@jlab.org></a>, "Johann Goetz"
<a class="moz-txt-link-rfc2396E" href="mailto:theodore.goetz@gmail.com"><theodore.goetz@gmail.com></a><br>
<b>Sent: </b>Friday, March 28, 2014 4:22:52 PM<br>
<b>Subject: </b>Re: [G12] GFlux<br>
<br>
Greetings,<br>
<br>
I am wondering if Johann's trip program takes into account the
tagger hit status to be 7 or 15 as well as does the GFlux
calculator when using the "good" option take this into
account?<br>
<br>
Also, Eugene, when you wrote to look at the flux and yield and
normalized yield per E-bin, would the yield of pi0 suffice?<br>
<br>
BR<br>
MK<br>
<br>
<br>
<div class="moz-cite-prefix">On 3/28/14 4:13 PM, Eugene Pasyuk
wrote:<br>
</div>
<blockquote
cite="mid:720603509.134399.1396037608071.JavaMail.root@jlab.org">
<style>p { margin: 0; }</style>
<div style="font-family: arial,helvetica,sans-serif;
font-size: 12pt; color: #000000">The other thing is to
make sure when you select your events and match them to
the tagger you use tagger hits with status 7 or 15 only.<br>
<br>
<div><span></span>-Eugene<span></span><br>
</div>
<br>
<hr id="zwchr">
<blockquote style="border-left:2px solid
#1010FF;margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>From:
</b>"Eugene Pasyuk" <a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:pasyuk@jlab.org" target="_blank"><pasyuk@jlab.org></a><br>
<b>To: </b>"Michael C. Kunkel" <a
moz-do-not-send="true" class="moz-txt-link-rfc2396E"
href="mailto:mkunkel@jlab.org" target="_blank"><mkunkel@jlab.org></a><br>
<b>Cc: </b>"g12" <a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:g12@jlab.org" target="_blank"><g12@jlab.org></a>,
"Johann Goetz" <a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:theodore.goetz@gmail.com" target="_blank"><theodore.goetz@gmail.com></a><br>
<b>Sent: </b>Friday, March 28, 2014 4:04:41 PM<br>
<b>Subject: </b>Re: [G12] GFlux<br>
<br>
<style>p { margin: 0; }</style>
<div style="font-family: arial,helvetica,sans-serif;
font-size: 12pt; color: #000000">The way to figure out
is to look at the flux and yield and normalized yield
per E-bin, that is the finest binning. If they are
consistent, the the problem is in re-binning
procedure. If they are not, then we have to look
further. Technically speaking there is not unusual
that you may have holes on photon spectrum, that would
just indicate that some of the counters are dead.
<div><br>
<div><span></span>-Eugene<span></span><br>
</div>
<br>
<hr id="zwchr">
<blockquote style="border-left:2px solid
#1010FF;margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>From:
</b>"Michael C. Kunkel" <a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:mkunkel@jlab.org" target="_blank"><mkunkel@jlab.org></a><br>
<b>To: </b>"Eugene Pasyuk" <a
moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:pasyuk@jlab.org" target="_blank"><pasyuk@jlab.org></a>,
"Johann Goetz" <a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:theodore.goetz@gmail.com"
target="_blank"><theodore.goetz@gmail.com></a><br>
<b>Cc: </b>"g12" <a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:g12@jlab.org" target="_blank"><g12@jlab.org></a><br>
<b>Sent: </b>Friday, March 28, 2014 3:55:45 PM<br>
<b>Subject: </b>Re: [G12] GFlux<br>
<br>
Greetings,<br>
<br>
I am wondering where we are at with the flux
issue?<br>
Meaning, what is going on with the flux in the 3.5
GeV region?<br>
<br>
I feel that the flux is in error for a few
reasons. <br>
First when I bin flux in 15 MeV, I see 0 flux at
3.515 - 3.53.<br>
Second, when I bin flux in 50 MeV binning I see a
drop in flux at 3.5 MeV<br>
Here is the output<br>
//<br>
3.45 3.5 1.06825209003e+12<br>
3.5 3.55 470202243457.0<br>
3.55 3.6 1.69234003548e+12<br>
//<br>
<br>
Third reason is that if you look at my attached
pdf, I show the pi0 cross section for various
energies, if you were to go to compare the cross
sections around 3.5 GeV to the neighboring Egamma
bins, you will notice a sizable discrepancy. See
G12_XSection_ALL.pdf<br>
<br>
Fourth reason, if I were to compare the pi0 cross
section in unit of dsigma/dt and compare it other
experiments at 4 and 5 GeV, our cross section is
off by a factor of 2, see t_vs_dsigdt_Compare.pdf<br>
<br>
In order for any g12 papers to be published, we
need verification of cross sections, which as of
this moment I do not believe we have.<br>
<br>
BR<br>
MK<br>
<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">On 3/25/14 10:16 PM,
Eugene Pasyuk wrote:<br>
</div>
<blockquote
cite="mid:1908667858.873792.1395800206351.JavaMail.root@jlab.org">
<style>p { margin: 0; }</style>
<div style="font-family:
arial,helvetica,sans-serif; font-size: 12pt;
color: #000000">This is not quite correct,
especially when bins are small and comparable
with the width of E-bin.
<div>Just imagine, You energy bin is
3.000-3.015 GeV. E-bin let's say covers
3.0121-3.0181 GeV. The center of E-bin will
be 3.01505 GeV. When you enter this into the
histogram it will fall to 3.015-3.030 bin
then. However, this E-bin should be split
between two adjacent bins in new binning
scheme. There are two possible ways of doing
this. Either randomize energy within the
width of the E-bin before you enter it to th
ehistogram, or split it between adjacent
bins using weight which depends on how the
bin boundary splits E-bin and have two
entries to the histogram with sum of the
weights equal one. You can imagine, if the
bin width is chosen to be less than the
widths of E-bin, then you will get
completely empty bins in resulting
histogram.<br>
<br>
<div><span></span>-Eugene<span></span><br>
</div>
<br>
<hr id="zwchr">
<blockquote style="border-left:2px solid
#1010FF;margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>From:
</b>"Johann Goetz" <a
moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:theodore.goetz@gmail.com"
target="_blank"><theodore.goetz@gmail.com></a><br>
<b>To: </b>"Eugene Pasyuk" <a
moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:pasyuk@jlab.org"
target="_blank"><pasyuk@jlab.org></a><br>
<b>Sent: </b>Tuesday, March 25, 2014
9:12:49 PM<br>
<b>Subject: </b>Re: [G12] GFlux<br>
<br>
<p dir="ltr">I consider the beam energy to
be the center of the energy paddle and
fill the histogram at that energy with
the number if counts in the gflux file.
Is that wrong?</p>
<div class="gmail_quote">On Mar 25, 2014
5:59 PM, "Eugene Pasyuk" <<a
moz-do-not-send="true"
href="mailto:pasyuk@jlab.org"
target="_blank">pasyuk@jlab.org</a>>
wrote:<br>
<blockquote class="gmail_quote"
style="margin:0 0 0
.8ex;border-left:1px #ccc
solid;padding-left:1ex">
<div>
<div
style="font-size:12pt;font-family:arial,helvetica,sans-serif">I
see four numbers in each row in
this file. I figured the first one
is E-bin number. What are the
other three? I guess the second
and the fourth are some energies.
What is the third number which is
integer and varies a lot.<br>
When you say you "histogram" does
it mean you take one fixed number
for each E-bin (beam centroid?)
and put it in the histogram, or
something else? <br>
<br>
<div><span></span>-Eugene<span></span><br>
</div>
<br>
<hr>
<blockquote
style="padding-left:5px;font-size:12pt;font-style:normal;margin-left:5px;font-family:Helvetica,Arial,sans-serif;text-decoration:none;font-weight:normal;border-left:2px
solid #1010ff"><b>From: </b>"Johann
Goetz" <<a
moz-do-not-send="true"
href="mailto:theodore.goetz@gmail.com"
target="_blank">theodore.goetz@gmail.com</a>><br>
<b>To: </b>"g12" <<a
moz-do-not-send="true"
href="mailto:g12@jlab.org"
target="_blank">g12@jlab.org</a>><br>
<b>Sent: </b>Tuesday, March 25,
2014 5:40:42 PM<br>
<b>Subject: </b>Re: [G12] GFlux<br>
<br>
<div dir="ltr">The output of
gflux includes a file that has
the energies of the individual
energy paddles and the counts.
These files are then used to
fill a histogram in the
g12-gflux script where the
range and bins are specified
on the command line. The
energies are slightly
"corrected" to the ones found
in the following file:
<div>
/group/clas/parms/pcor/g12/tagger_energies.dat</div>
<div>But, the correction is
very small and averaged over
the whole run period (i.e.
not one a run-by-run basis).</div>
<div>
<div><br>
</div>
</div>
<div>And before you ask, the
the summing of the flux is
valid to at least 4 figures
- probably 5. </div>
</div>
<div class="gmail_extra"><br
clear="all">
<div>-- <a
moz-do-not-send="true"
href="https://sites.google.com/site/theodoregoetz/"
target="_blank">Johann</a></div>
<br>
<br>
<div class="gmail_quote">On
Tue, Mar 25, 2014 at 4:28
PM, Michael C. Kunkel <span
dir="ltr"><<a
moz-do-not-send="true"
href="mailto:mkunkel@jlab.org"
target="_blank">mkunkel@jlab.org</a>></span>
wrote:<br>
<blockquote
class="gmail_quote"
style="margin:0 0 0
.8ex;border-left:1px #ccc
solid;padding-left:1ex">
Greetings,<br>
<br>
If the data shows a hole
at beam energies, how does
the gflux for this<br>
energy have events?<br>
Below is a table of Gflux
for 15MeV binning, which
is to fine for<br>
analysis but great for
illustrating this
question.<br>
<br>
Also, notice the 0 entries
at 3.5 GeV, and low events
in neighbors.<br>
FORMAT: Emin:Emax:Flux<br>
<br>
<br>
<br>
_______________________________________________<br>
G12 mailing list<br>
<a moz-do-not-send="true"
href="mailto:G12@jlab.org" target="_blank">G12@jlab.org</a><br>
<a moz-do-not-send="true"
href="https://mailman.jlab.org/mailman/listinfo/g12" target="_blank">https://mailman.jlab.org/mailman/listinfo/g12</a><br>
</blockquote>
</div>
<br>
</div>
<br>
_______________________________________________<br>
G12 mailing list<br>
<a moz-do-not-send="true"
href="mailto:G12@jlab.org"
target="_blank">G12@jlab.org</a><br>
<a moz-do-not-send="true"
href="https://mailman.jlab.org/mailman/listinfo/g12"
target="_blank">https://mailman.jlab.org/mailman/listinfo/g12</a><br>
</blockquote>
<br>
</div>
</div>
</blockquote>
</div>
</blockquote>
<br>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
G12 mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:G12@jlab.org" target="_blank">G12@jlab.org</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="https://mailman.jlab.org/mailman/listinfo/g12" target="_blank">https://mailman.jlab.org/mailman/listinfo/g12</a>
</pre>
</blockquote>
<br>
</blockquote>
<br>
</div>
</div>
</blockquote>
<br>
</div>
</blockquote>
<br>
</blockquote>
<br>
</div>
</blockquote>
<br>
</body>
</html>