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

#ifndef __PHOTEL__
#define __PHOTEL__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== photel ===========================================================*
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     Include file photel:                                             *
//             created on 24-October-1990 by   A. Ferrari & P. Sala     *
//                                                   INFN - Milan       *
//                                                                      *
//     Included in the following routines:                              *
//                                                                      *
//                Photo   (new version)                                 *
//                Phoflu  (version with fluorescence of photo)          *
//                Bdphfl  (block data for photoelectrons and            *
//                           fluorescence data)                         *
//     Variable:                                                        *
//                Umuphe (i,j) = tabulated emission cosines of emitted  *
//                               photoelectrons at given cumulative pro-*
//                               bability ( prob = j / nphumu ) and at  *
//                               given electron beta**2 ( beta2 =       *
//                               (i-1) x dphbt2 ), computed according   *
//                               to the relativistic theory of F. Von   *
//                               Sauter                                 *
//                                                                      *
//                Dumphe (i,j) = tabulated derivatives of the inverse   *
//                               function computed at umuphe            *
//                                                                      *
//----------------------------------------------------------------------*
//
const double PHB2MX = 9.800000000000000E-01;
const double DPHBT2 = 7.000000000000000E-02;
const double BETA00 = 2.645751311064591E-01;
const int    NPHUMU = 30;
//
const int    NPHBET = 15;
extern struct {
	double	umuphe[NPHBET][NPHUMU+1];
	double	dumphe[NPHBET][NPHUMU+1];
} photel_;
const int photel_umuphe_base = 0;
const int photel_dumphe_base = 0;


#ifdef __cplusplus
}
#endif
#endif