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

#ifndef __BREANG__
#define __BREANG__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== breang ===========================================================*
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     Created on  22  March 1991   by    Alfredo Ferrari & Paola Sala  *
//                                                   Infn - Milan       *
//                                                                      *
//     Included in the following routines:                              *
//                                                                      *
//                      BDANDI                                          *
//                      BREMNW                                          *
//                      BREMS                                           *
//                                                                      *
//         Mxphtb = number of intervals on which phi(y) is tabulated    *
//         Yphimn = minimum value of y for which phi is taublated       *
//         Yphirt = ratio of y between two tabulated points             *
//         Yphimx = maximum value of y for which phi is taublated       *
//                  (=Yphimn*Yphirt**(Mxphtb-1))                        *
//         Phiytb = tabulated values of phi                             *
//         Alymin = natural logarithm of Yphimn                         *
//         Alyrat = natural logarithm of Yphirt                         *
//         Alytra = natural logarithm of the max. y allowed for Xsitra  *
//         Phia00 = used for asymptotic behaviour                       *
//         Phib00 = used for asymptotic behaviour                       *
//         Phic0  = used for asymptotic behaviour                       *
//         Phid0  = used for asymptotic behaviour                       *
//        ( for y > yphimx, Phi = Phia00 + Phib00 / y^2                 *
//          for y < yphimn, Phi = log(y) + Phic0 + Phid0 x y^2 )        *
//         Accrit = parameter used for the Migdal polarization effect,  *
//                  given by Nel x Lambda_compt^2 x r0 / pi             *
//         Zbrem  = approximate "effective" Z for bremsstrahlung        *
//         Fcoul  = Coulomb factor fc(Z)                                *
//         Gmolie = factor entering Moliere's expansion of Thomas-Fermi *
//                  form factor ( = Z^1/3/121 )                         *
//         Algmol = natural logarithm of Gmolie                         *
//                                                                      *
//----------------------------------------------------------------------*
//
//  Tpifsc = 2 x pi x fine structure constant
const double TPIFSC = 4.5850621648360624E-02;
const int    MXPHTB = 100;
const double YPHIMN = 1.0E-01;
//
const double YPHIRT = 1.07E+00;
extern struct {
	double	phiytb[MXPHTB];
	double	yphimx;
	double	alymin;
	double	alyrat;
	double	alytra;
	double	phia00;
	double	phib00;
	double	phic0;
	double	phid0;
	double	accrit[MXXMDE];
	double	zbrem[MXXMDE];
	double	fcoul[MXXMDE];
	double	gmolie[MXXMDE];
	double	algmol[MXXMDE];
} breang_;


#ifdef __cplusplus
}
#endif
#endif