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

#ifndef __PARNUC__
#define __PARNUC__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== parnuc ===========================================================*
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     Created on  2 August 1991    by    Alfredo Ferrari & Paola Sala  *
//                                                   Infn - Milan       *
//                                                                      *
//     included in :                                                    *
//                                                                      *
//                   Bdpree                                             *
//                   Bimsel                                             *
//                   Cmanmo                                             *
//                   Couset                                             *
//                   Muoabs                                             *
//                   Nucdcy                                             *
//                   Nucnuc                                             *
//                   Nwisel                                             *
//                   Peanut                                             *
//                   Phdset                                             *
//                   Phdwll                                             *
//                   Pioabs                                             *
//                   Pionuc                                             *
//                   Pmprab                                             *
//                   Pncset                                             *
//                   Pncset                                             *
//                   Prepre                                             *
//                   Sbcomp                                             *
//                   Umofin                                             *
//                                                                      *
//          Kionpr(i) = i0 + i1 * 1000 + i2 * 1000                      *
//                      i0 = particle type                              *
//                         0 = normal "target-frame"     secondary      *
//                         1 = normal "projectile-frame" secondary      *
//                         2 = original "target-frame"     ion projec-  *
//                             tile nucleon                             *
//                         3 = original "projectile-frame" ion projec-  *
//                             nucleon                                  *
//                      i1 = id of the originating nucleon if from an   *
//                           ion projectile                             *
//                      i2 = not yet used                               *
//          Krsnce(i) = id of (possible) resonance for i_th particle    *
//        Krsdcy(j,k) = j_th property of the k_th resonance (j=1 is the *
//                      id)                                             *
//                                                                      *
//----------------------------------------------------------------------*
//
//  Standard parameter for central density: it is exactly r0nucl of
//  Parevt
const double ALEVEL = 8.E-03;
//  "Naive" r0 for sigmas
const double RCNUCL = 1.12E+00;
//  r0 for sigmas according to Kikuchi & Kawai
const double R0SIG = 1.3E+00;
//  r0 for Coulomb barrier according to Kikuchi & Kawai
const double R0SIGK = 1.5E+00;
//  Parameter for Coulomb barrier according to Dostrovski (r0c = 1.63 fm)
//  The original paper quotes 1.7 fm rather than 1.63
const double RCOULB = 1.5E+00;
//  Central density in nuclei
const double COULBH = 0.88235E-03;
//  For Taufor, Aefrav & Rhonuc the combinations 1 and 3 are ok. If the
//  Hannes expression for Taufor is used we must use Taufor from
//  combination 3 and Aefrav and Rhonuc from combination 1
//  Combination 3:
//3    PARAMETER ( TAUFO0 = 3.0D+00 )
//  Combination 1:
//     PARAMETER ( TAUFO0 = 5.0 D+00 )
const double RHONU0 = 0.75E+00/PIPIPI/RCNUCL/RCNUCL/RCNUCL;
//  Thresholds for explicit interactions:
//    Primary projectile:
const double TAUFO0 = 10.0E+00;
//    Secondary nucleons:
const double EKEEXP = 0.03E+00;
//    Just for Pauli checks:
const double EKREXP = 0.05E+00;
//
const double EKEMNM = 0.01E+00;
const int    NCPMX = 1000;
const int    NRSMX = 10;
const int    NDCMX = 5;
extern struct {
	double	ekori;
	double	pxori;
	double	pyori;
	double	pzori;
	double	ptori;
	double	taufor;
	double	tfrabs;
	double	tfrdcy;
	double	textau;
	double	texqel;
	double	texres;
	double	texnut;
	double	ennuc[NCPMX];
	double	pnucl[NCPMX];
	double	ekfnuc[NCPMX];
	double	xstnuc[NCPMX];
	double	ystnuc[NCPMX];
	double	zstnuc[NCPMX];
	double	xlcnuc[NCPMX];
	double	ylcnuc[NCPMX];
	double	zlcnuc[NCPMX];
	double	pxnucl[NCPMX];
	double	pynucl[NCPMX];
	double	pznucl[NCPMX];
	double	rstnuc[NCPMX];
	double	freepa[NCPMX];
	double	crrpan[NCPMX];
	double	crrpap[NCPMX];
	double	bstnuc[NCPMX];
	double	agenuc[NCPMX];
	double	taufpa[NCPMX];
	double	slcnuc[NCPMX];
	double	rhnucl[2][NCPMX];
	double	amssrs[NRSMX];
	double	gammrs[NRSMX];
	double	enrsdc[NRSMX][NDCMX];
	double	pxrsdc[NRSMX][NDCMX];
	double	pyrsdc[NRSMX][NDCMX];
	double	pzrsdc[NRSMX][NDCMX];
	double	bndgav;
	double	defmin;
	double	frpmby;
	double	ftfrdc;
	double	amrsft;
	double	angxin;
	double	angyin;
	double	angzin;
	int	kpnucl[NCPMX];
	int	krfnuc[NCPMX];
	int	isfnuc[NCPMX];
	int	ilinuc[2][NCPMX];
	int	inucts[2][NCPMX];
	int	nnucts[2];
	int	kionpr[NCPMX];
	int	krsnce[NCPMX];
	logical	latauf[NCPMX];
	int	krsdcy[NRSMX][NDCMX];
	int	kpori;
	int	ibori;
	int	ichori;
	int	ibnucl;
	int	npnuc;
	int	npanmo;
} parnuc_;


#ifdef __cplusplus
}
#endif
#endif