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

#ifndef __MULHD__
#define __MULHD__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== mulhd ============================================================*
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//   Common mulhd for the new multiple scattering in Fluka              *
//                  created by A. Ferrari & P. Sala on 23-Oct-1989      *
//                                                                      *
//          included in:                                                *
//                        bdnopt                                        *
//                        deflts                                        *
//                        enion                                         *
//                        epilog                                        *
//                        fiprou                                        *
//                        flukam                                        *
//                        formfh                                        *
//                        hmsnsc                                        *
//                        hmulrt                                        *
//                        kaskad                                        *
//                        mulhad                                        *
//                        mulmix                                        *
//                        omegah                                        *
//                        stepop                                        *
//                        thnszh                                        *
//                        zeroin                                        *
//                                                                      *
//                        Mxxmdf = maximum number of media              *
//                                                                      *
//----------------------------------------------------------------------*
//
//     PARAMETER ( TIMESS = 3.00D+00 )
const int    MXXPT1 = 2;
//  Tmsrlx: times "relaxed" if the nuclear finite size correction must
//  not be taken into account
//     PARAMETER ( TMSRLX = 2.00D+00 )
const double TIMESS = 2.00E+00;
//     PARAMETER ( EPSINS = 0.10D+00 )
const double TMSRLX = 1.50E+00;
//  Epsrlx: epsins "relaxed" if the nuclear finite size correction must
//  not be taken into account
//     PARAMETER ( EPSRLX = 0.15D+00 )
const double EPSINS = 0.15E+00;
//  Sqepsn = sqrt (epsins)
//     PARAMETER ( SQEPSN = 0.3162277660168379 D+00 )
const double EPSRLX = 0.50E+00;
//  Sqepsr = sqrt (epsrlx)
//     PARAMETER ( SQEPSR = 0.3872983346207417 D+00 )
const double SQEPSN = 0.3872983346207417E+00;
const double SQEPSR = 0.7071067811865475E+00;
const double PARNSI = 1.732050807568877E+00*SQEPSN;
//     PARAMETER ( R0NCMS = 1.12 D+00 )
const double PRNSR0 = 1.732050807568877E+00*SQEPSR;
const double R0NCMS = 1.20E+00;
//
extern struct {
	double	blcc[MXXMDF];
	double	blccra[MXXMDF];
	double	fanolc[MXXMDF];
	double	zdfano[MXXMDF];
	double	xcc[MXXMDF];
	double	ztilde[MXXPT1+1][MXXMDF];
	double	alpztl[MXXPT1+1][MXXMDF];
	double	rldflu[MXXMDF];
	double	alpzt2[MXXPT1+1][MXXMDF];
	double	tfflu0[MXXMDF];
	double	xr0flu[MXXMDF];
	double	ecutm[2][NALLWP+6+1][MXXMDF];
	double	estepf[MXXMDF];
	double	hthnsz[NALLWP+6+1][MXXMDF];
	double	ae1o3[MXXMDF];
	double	parnsr[MXXMDF];
	double	heesli[MXXMDF];
	double	elmcsf[2][MXXMDF];
	double	thmspr;
	double	thmssc;
	double	hmsamp;
	double	hmreje;
	double	hmsngl[5+1];
	logical	lsrcrh[MXXMDF];
	logical	lnscrh[MXXMDF];
	logical	lnsceh[MXXMDF];
	logical	ltopt[MXXMDF];
	logical	lhfano[MXXMDF];
	int	ihorel[NALLWP+6+1];
	int	nfscat;
	int	nfllda;
	int	ihmsng;
	logical	lmlrgh;
	logical	lmcspr;
	logical	lmcssc;
} mulhd_;
const int mulhd_ztilde_base = 0;
const int mulhd_alpztl_base = 0;
const int mulhd_alpzt2_base = 0;
const int mulhd_ecutm_base = -6;
const int mulhd_hthnsz_base = -6;
const int mulhd_hmsngl_base = 0;
const int mulhd_ihorel_base = -6;


#ifdef __cplusplus
}
#endif
#endif