/*
 * DO NOT MANUALLY UPDATE THIS FILE.
 * If needed: make cpp_headers
 *
 * fortran2c.py - automatic conversion
 *       atffac.inc -> atffac.h
 */

#ifndef __ATFFAC__
#define __ATFFAC__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== ATFFAC ===========================================================*
//
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     Include file atffac:  ATomic Form FACtors                        *
//                                                                      *
//     Created on  18  March 1992   by    Alfredo Ferrari & Paola Sala  *
//                                                   Infn - Milan       *
//                                                                      *
//     Gmoliz(iz) = Z^1/3 / 121                                         *
//     Algmlz(iz) = Log (Gmoliz(iz))                                    *
//     Xsielz(iz) = asymptotic contribution of atomic electrons to pair *
//                  and bremsstrahlung                                  *
//     Fclmbz(iz) = Coulomb correction                                  *
//     Aagelz(iz) = a for the Tsai fit to the atomic elastic form factor*
//                  [1-F^2(q)]=(aq)^4/[1+(aq)^2]^2, [a] = [MeV/c]^-1    *
//                  For Z>=5 is given by a = 111.7 / (Z^1/3 me)         *
//     Apginz(iz) = a' for the Tsai fit to the atomic inelastic form    *
//                  factor S(q)=(a'q)^4/[1+(a'q)^2]^2, [a'] = [MeV/c]^-1*
//                  For Z>=5 is given by a' = 724.2 / (Z^2/3 me)        *
//                                                                      *
//     Actually a and a' are stored already squared !!                  *
//                                                                      *
//     Ph10o4(iz) = Phi1(0)/4, full screening value of the Phi1 elas-   *
//                  tic screening function, Phi1(0)/4 = Log(184.15) or  *
//                  equivalent for low Z                                *
//     Ps10o4(iz) = Psi1(0)/4, full screening value of the Psi1 inelas- *
//                  tic screening function, Psi1(0)/4 = Log(1194) or    *
//                  equivalent for low Z                                *
//     Alradz(iz) = Log(184.15/Z^1/3) or equivalent for low Z           *
//     Asqzft(iz) = a parameter for the fit to S(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     Bsqzft(iz) = b parameter for the fit to S(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     Csqzft(iz) = c parameter for the fit to S(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     Dsqzft(iz) = d parameter for the fit to S(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     Esqzft(iz) = e parameter for the fit to S(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     Afqzft(iz) = a parameter for the fit to F(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     Bfqzft(iz) = b parameter for the fit to F(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     Cfqzft(iz) = c parameter for the fit to F(q,Z) computed with the *
//                  Hartree-Fock method                                 *
//     X0fqzt(iz) = Log(X0) parameter for the F(q,Z) tabulations        *
//                                                                      *
//     Xrfqzt(iz) = Log of Xr ratio for the fit to F(q,Z) tabulations   *
//                 (from 1st to N_th value = X0 x Xr^(N-1))             *
//   Xifqzt(i,iz) = Log of X extra (N+i_th) points for the F(q,Z) tab-  *
//                  ulations                                            *
//     Nfqfzt(iz) = number of "ordinary" tabulation points for the      *
//                  F(q,Z) tabulations                                  *
//                                                                      *
//----------------------------------------------------------------------*
//
//  111.7 = 184.15 / sqrt(e)
const double A121 = 121.E+00;
//  724.2 = 1194   / sqrt(e)
const double A111P7 = 111.7E+00;
const double A724P2 = 724.2E+00;
const double A184 = 184.15E+00;
const double A1194 = 1194.E+00;
//
const int    NFQLST = 5;
extern struct {
	double	gmoliz[100];
	double	algmlz[100];
	double	xsielz[100];
	double	fclmbz[100];
	double	aagelz[100];
	double	apginz[100];
	double	alradz[100];
	double	asqzft[100];
	double	bsqzft[100];
	double	csqzft[100];
	double	dsqzft[100];
	double	esqzft[100];
	double	afqzft[100];
	double	bfqzft[100];
	double	cfqzft[100];
	double	x0fqzt[100];
	double	xrfqzt[100];
	double	a0fqzt[100];
	double	a1fqzt[100];
	double	ph10o4[100];
	double	ps10o4[100];
	double	xifqzt[100][NFQLST+1];
	int	kfqfzt[100];
	int	nfqfzt[100];
} atffac_;
const int atffac_xifqzt_base = 0;


#ifdef __cplusplus
}
#endif
#endif