[Eg2] RC code instructions

Hayk Hakobyan hayk at jlab.org
Tue Apr 12 18:02:55 EDT 2011


Dear Eg2,

Here I send you some instructions to use the HAPRAD_CPP code.
https://github.com/usm-data-analysis/HAPRAD_cpp

For some reason in the production version of HAPRAD_CPP there was the
code with old HAPRAD2 parameters. Now the production version is
updated so it has the write version in it.

To process pi+'s one should create in the test directory text file
called "parameters" and put the following inside:

In the case of the deuterium:
par0 = 5.54600e-02;
par1 = 1.01711e+00;
par2 = -8.47055e-01;
par3 = 7.76589e-02;
par4 = 5.86747e-11;
A1 = 2.51386e-01;
B1 = 1.04605e+00;
C1 = -1.16972e+00;
D1 = -1.40120e+00;
E1 = 1.97121e+00;

In the case of the carbon:
par0 = 1.05430e-01;
par1 = 1.33184e+00;
par2 = -1.25147e+00;
par3 = -2.61048e-01;
par4 = 3.67470e-01;
A1 = 9.79762e-02;
B1 = 5.62038e-01;
C1 = -4.64066e-01;
D1 = -2.27958e-01;
E1 = 2.89126e-01;

In the case of the iron:
par0 = 1.07049e-01
par1 = 1.44052e+00
par2 = -1.37651e+00
par3 = -3.70841e-01
par4 = 6.24209e-01
A1 = 8.66764e-02
B1 = 6.02301e-01
C1 = -3.76230e-01
D1 = -2.42141e-01
E1 = 3.98611e-01


In the case of the lead:
par0 = 1.07485e-01;
par1 = 1.46133e+00;
par2 = -1.41410e+00;
par3 = -3.36003e-01;
par4 = 5.16617e-01;
A1 = 1.00957e-01;
B1 = 6.29042e-01;
C1 = -5.31384e-01;
D1 = -2.48889e-01;
E1 = 3.21765e-01;

Besides that one should take care of the fraction of protons in the nuclear
for exclusive tail by changing it from 1/2 which is for case of deuterium to
the corresponding number for the solid target in the function called
GetFactor3(). This requires the program to be recompiled. I'll let you
know when this part also will be simplified for use.

To use the code for pi- and pi0 there are some modifications which
still should be done. I'll do it pretty soon and let you know too.

Please feel free to ask any question you may have.

Thanks,
Hayk.




More information about the Eg2 mailing list