************************************************************************
************************************************************************
* Stefan Kretzer (kretzer@pa.msu.edu) Nov 2001:                        *
*                                                                      *
* Library of FF parametrizations updated to be posted on               *
* http://www.pv.infn.it/~radici/FFdatabase/                            *
* site maintained by Marco Radici (radici@pv.infn.it)                  *
* and Rainer Jakob (rainer@theorie.physik.uni-wuppertal.de)            *
* within the EU network:                                               *
* "Hadronic Physics with High Energy Electromagnetic Probes"           *
*                                                                      *
* Also posted on that site are the individual FF sets included in this *
* library as distributed by the corresponding authors.                 *
************************************************************************
* K, KKP and BFGW fragmentation functions for                          *
* iparton -> ihadron                                                   *
*                                                                      *
* Subroutine "dlib" returns D(z,Q2) [NOT z*D(z,Q2)] as an array        *
* "dff(iparton)" of dimension dff(-5:5)                                *
*                                                                      *
* iparton = 5,4,3,2,1,0,-1,...,-5 means b,c,s,d,u,g,ubar,...,bbar      *
*                                                                      *
* [Q2] = GeV^2                                                         *
*                                                                      *
* ffset = 1,2,3 means K, KKP, BFGW                                     *
* Details are to be found in the corresponding references:             *
* S.Kretzer, Phys.Rev.D62, 054001 (2000)                               *
* B.A.Kniehl, G.Kramer, B.Potter, Nucl.Phys.B582, 514 (2000)           *
* L.Bourhis, M.Fontannaz, J.P.Guillet, M.Werlen,Eur.Phys.J.C19,89(2001)*
*                                                                      *
* fforder=0,1 is LO, NLO(MSbar)                                        *
* no LO set for BFGW                                                   *
*                                                                      *
* ihadron=1,2,3,4,5 is pi,K,h,p,n                                      *
* no pi,K sets for BFGW                                                *
* no p,n sets for K, BFGW                                              *
*                                                                      *
* icharge=0,1,2,3 is 0,+,-,+&-                                         *
* note: 3 = +&- = charge sum (NOT average) in this library             *
* for the neutral (icharge=0) particles K^0, n:                        *
* icp = 1,2,3 chooses between particle, anti-particle or sum of both   *
* icp is inactive for pi,h,p and charged Kaons                         *
*                                                                      *
* Following KKP, FFs into {pi^0;(anti-)K^0} and into {pi^+/-;K^+/-}    *
* are related by isospin.                                              *
*                                                                      *
* ipi = 1,2,3 is a flag for BFGW; inactive for K, KKP:                 *
* ipi = 1: best fit     (formula (8))  of BFGW                         *
* ipi = 2: large Ng set (formula (9))  of BFGW                         *
* ipi = 3: low Ng set   (formula (10)) of BFGW                         *
************************************************************************
* If options (ihadron, icharge etc.) are chosen which do not exist for *
* a given parametrization (ffset) or if a flag is chosen outside its   *
* range as defined above then the code returns a corresponding warning *
* and stops.                                                           *
************************************************************************
************************************************************************