Hi Hayk,<br><br>Could you explain the parameters? What formula did you use to extract those from data?<br><br>Thanks<br><br><div class="gmail_quote">On Tue, Apr 12, 2011 at 5:02 PM, Hayk Hakobyan <<a href="mailto:hayk@jlab.org">hayk@jlab.org</a>> wrote:<br>
<blockquote class="gmail_quote"><br>
Dear Eg2,<br>
<br>
Here I send you some instructions to use the HAPRAD_CPP code.<br>
<a href="https://github.com/usm-data-analysis/HAPRAD_cpp">https://github.com/usm-data-analysis/HAPRAD_cpp</a><br>
<br>
For some reason in the production version of HAPRAD_CPP there was the<br>
code with old HAPRAD2 parameters. Now the production version is<br>
updated so it has the write version in it.<br>
<br>
To process pi+'s one should create in the test directory text file<br>
called "parameters" and put the following inside:<br>
<br>
In the case of the deuterium:<br>
par0 = 5.54600e-02;<br>
par1 = 1.01711e+00;<br>
par2 = -8.47055e-01;<br>
par3 = 7.76589e-02;<br>
par4 = 5.86747e-11;<br>
A1 = 2.51386e-01;<br>
B1 = 1.04605e+00;<br>
C1 = -1.16972e+00;<br>
D1 = -1.40120e+00;<br>
E1 = 1.97121e+00;<br>
<br>
In the case of the carbon:<br>
par0 = 1.05430e-01;<br>
par1 = 1.33184e+00;<br>
par2 = -1.25147e+00;<br>
par3 = -2.61048e-01;<br>
par4 = 3.67470e-01;<br>
A1 = 9.79762e-02;<br>
B1 = 5.62038e-01;<br>
C1 = -4.64066e-01;<br>
D1 = -2.27958e-01;<br>
E1 = 2.89126e-01;<br>
<br>
In the case of the iron:<br>
par0 = 1.07049e-01<br>
par1 = 1.44052e+00<br>
par2 = -1.37651e+00<br>
par3 = -3.70841e-01<br>
par4 = 6.24209e-01<br>
A1 = 8.66764e-02<br>
B1 = 6.02301e-01<br>
C1 = -3.76230e-01<br>
D1 = -2.42141e-01<br>
E1 = 3.98611e-01<br>
<br>
<br>
In the case of the lead:<br>
par0 = 1.07485e-01;<br>
par1 = 1.46133e+00;<br>
par2 = -1.41410e+00;<br>
par3 = -3.36003e-01;<br>
par4 = 5.16617e-01;<br>
A1 = 1.00957e-01;<br>
B1 = 6.29042e-01;<br>
C1 = -5.31384e-01;<br>
D1 = -2.48889e-01;<br>
E1 = 3.21765e-01;<br>
<br>
Besides that one should take care of the fraction of protons in the nuclear<br>
for exclusive tail by changing it from 1/2 which is for case of deuterium to<br>
the corresponding number for the solid target in the function called<br>
GetFactor3(). This requires the program to be recompiled. I'll let you<br>
know when this part also will be simplified for use.<br>
<br>
To use the code for pi- and pi0 there are some modifications which<br>
still should be done. I'll do it pretty soon and let you know too.<br>
<br>
Please feel free to ask any question you may have.<br>
<br>
Thanks,<br>
Hayk.<br>
<br>
<br>
_______________________________________________<br>
Eg2 mailing list<br>
<a href="mailto:Eg2@jlab.org">Eg2@jlab.org</a><br>
<a href="https://mailman.jlab.org/mailman/listinfo/eg2">https://mailman.jlab.org/mailman/listinfo/eg2</a><br>
</blockquote></div><br><br clear="all"><br>-- <br>Raphaël Dupré<br><br><a href="http://www.jlab.org/%7Edupre/">Web Page</a><br><br>Argonne National Laboratory<br><br>Université Claude Bernard Lyon 1<br>