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

#ifndef __THRSCM__
#define __THRSCM__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== Thrscm ===========================================================*
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     THReShold (for FLUKA) CoMmon: for FLUKA9x/.../200x               *
//                                                                      *
//     Created on    15 May 1990    by    Alfredo Ferrari & Paola Sala  *
//                                                   Infn - Milan       *
//                                                                      *
//        Glbthr = global  threshold kinetic energy in GeV for particle *
//                 production and transport                             *
//        Ethrin = threshold kinetic energy in GeV for inelastic inter. *
//        Ethrel = threshold kinetic energy in GeV for elastic inter.   *
//        Ethstr = threshold kinetic energy in GeV for star scoring     *
//        Ethrij = threshold kinetic energy in GeV for each kind of     *
//                 particle                                             *
//        Pthrij = threshold momentum in GeV/c for each kind of         *
//                 particle                                             *
//        Pthneu = neutron momentum at the kinetic energy threshold     *
//        Ennmax = switch energy in GeV for the low energy neutron      *
//                 group treatment                                      *
//        Khvtrn = max. index for full heavy particle transport         *
//        Llownt = .true. if low energy neutron production and transport*
//                 has to be activated                                  *
//        Llmncp = .true. if low energy neutron production and transport*
//                 with MCNP has to be activated, .false. for default   *
//                 multigroup transport (never implemented!)            *
//        Lhvtrn = flag for full heavy particle transport               *
//        Liontr = flag for full transport of ions heavier than alphas  *
//                                                                      *
//----------------------------------------------------------------------*
//
//   Minimum absolute energy for transport/ranging out:
//   Minimum sensible energy/nucleon for transport/ranging out:
const double ETHRMN = 2.0E-06;
//
const double EOAMNM = 1.0E-06;
extern struct {
	double	glbthr;
	double	ethrin;
	double	ethrel;
	double	ethstr;
	double	ethrij[NALLWP+6+1];
	double	pthrij[NALLWP+6+1];
	double	pthneu;
	double	ennmax;
	int	khvtrn;
	logical	llownt;
	logical	llmcnp;
	logical	lhvtrn;
	logical	liontr;
} thrscm_;
const int thrscm_ethrij_base = -6;
const int thrscm_pthrij_base = -6;


#ifdef __cplusplus
}
#endif
#endif