Add class description
[u/mrichter/AliRoot.git] / PDF / npdf / eksarp.F
CommitLineData
82cedcfe 1#include "pdf/pilot.h"
2 function eksarp(x,v2,x1,x116,ptm,indpi)
3#include "pdf/impdp.inc"
4 dimension ptm(10)
5
6 z=v2
7 xx=x-x1
8 if (x.le.x116) then
9 eksarp=ptm(1)+ptm(2)*z+ptm(3)*z**2+ptm(4)*xx+ptm(5)*xx**2
10 else
11 z1= log(x116/x1)
12 xx1=x116-x1
13 ff0=ptm(1)+ptm(2)*z1+ptm(3)*z1**2+ptm(4)*xx1+ptm(5)*xx1**2
14 xx16=x-x116
15 z16= log(x/x116)
16 qexp=19.d0+(indpi-1)*(indpi-2)*8.d0/2.d0
17 eksarp=ff0+ptm(6)*xx16+ptm(7)*xx16**2+ptm(8)*xx16**3
18 # +ptm(9)*xx16**qexp+ptm(10)*z16
19 endif
2011 return
21 end
22
23