[Eg2] RC code instructions

[Hayk Hakobyan]

Hi Raphaël,

These parameters are used in TSemiInclusiveModel in HAPRAD_CPP.
The function you can find there too.

To extract them you can use my program which does it using RooFit.
The program is located here:
https://github.com/usm-data-analysis/HAPRAD_fit
it has instructions to use.

The idea is to find the description of the Pt,Z,Xb distribution with
some analytical function. (Pt-transverse momentum, Z=Eh/nu, Xb-Bjorken
X).

Let me know if you have any questions.

Best regards,
Hayk.





> Hi Hayk,
>
> Could you explain the parameters? What formula did you use to extract
> those
> from data?
>
> Thanks
>
> On Tue, Apr 12, 2011 at 5:02 PM, Hayk Hakobyan <hayk at jlab.org> wrote:
>
>>
>> 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.
>>
>>
>> _______________________________________________
>> Eg2 mailing list
>> Eg2 at jlab.org
>> https://mailman.jlab.org/mailman/listinfo/eg2
>>
>
>
>
> --
> Raphaël Dupré
>
> Web Page <http://www.jlab.org/%7Edupre/>
>
> Argonne National Laboratory
>
> Université Claude Bernard Lyon 1
>