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

#ifndef __DTEVT1__
#define __DTEVT1__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     This module is required for interfacing DPMJET (and rQMD-2.4)    *
//     to FLUKA                                                         *
//                                                                      *
//----------------------------------------------------------------------*
//

//c --------- /HKKEVT/ ------------------------------------------------
const int    NMXHKK = 250000;
extern struct {
	int	nhkk;
	int	nevhkk;
	int	isthkk[NMXHKK];
	int	idhkk[NMXHKK];
	int	jmohkk[NMXHKK][2];
	int	jdahkk[NMXHKK][2];
	double	phkk[NMXHKK][5];
	double	vhkk[NMXHKK][4];
	double	whkk[NMXHKK][4];
} dtevt1_;

//
//                       WHKK(4,NMXHKK) GIVES POSITIONS AND TIMES IN
//                       PROJECTILE FRAME, THE CHAINS ARE CREATED ON
//                       THE POSITIONS OF THE PROJECTILE NUCLEONS
//                       IN THE PROJECTILE FRAME (TARGET NUCLEONS IN
//                       TARGET FRAME) BOTH POSITIONS ARE THREFORE NOT
//                       COMPLETELY CONSISTENT. THE TIMES IN THE
//                       PROJECTILE FRAME HOWEVER ARE OBTAINED BY
//                       LORENTZ TRANSFORMING FROM THE LAB SYSTEM.
//
//  Based on the proposed standard COMMON block (Sjostrand Memo 17.3,89)
//
// NMXHKK: maximum numbers of entries (partons/particles) that can be
//    stored in the commonblock.
//
// NHKK: the actual number of entries stored in current event. These are
//    found in the first NHKK positions of the respective arrays below.
//    Index IHKK, 1 <= IHKK <= NHKK, is below used to denote a given
//    entry.
//
// ISTHKK(IHKK): status code for entry IHKK, with following meanings:
//    = 0 : null entry.
//    = 1 : an existing entry, which has not decayed or fragmented.
//        This is the main class of entries which represents the
//        "final state" given by the generator.
//    = 2 : an entry which has decayed or fragmented and therefore
//        is not appearing in the final state, but is retained for
//        event history information.
//    = 3 : a documentation line, defined separately from the event
//        history. (incoming reacting
//        particles, etc.)
//    = 4 - 10 : undefined, but reserved for future standards.
//    = 11 - 20 : at the disposal of each model builder for constructs
//        specific to his program, but equivalent to a null line in the
//        context of any other program. One example is the cone defining
//        vector of HERWIG, another cluster or event axes of the JETSET
//        analysis routines.
//    = 21 - : at the disposal of users, in particular for event tracking
//        in the detector.
//
// IDHKK(IHKK) : particle identity, according to the Particle Data Group
//    standard.
//
// JMOHKK(1,IHKK) : pointer to the position where the mother is stored.
//    The value is 0 for initial entries.
//
// JMOHKK(2,IHKK) : pointer to position of second mother. Normally only
//    one mother exist, in which case the value 0 is used. In cluster
//    fragmentation models, the two mothers would correspond to the q
//    and qbar which join to form a cluster. In string fragmentation,
//    the two mothers of a particle produced in the fragmentation would
//    be the two endpoints of the string (with the range in between
//    implied).
//
// JDAHKK(1,IHKK) : pointer to the position of the first daughter. If an
//    entry has not decayed, this is 0.
//
// JDAHKK(2,IHKK) : pointer to the position of the last daughter. If an
//    entry has not decayed, this is 0. It is assumed that the daughters
//    of a particle (or cluster or string) are stored sequentially, so
//    that the whole range JDAHKK(1,IHKK) - JDAHKK(2,IHKK) contains
//    daughters. Even in cases where only one daughter is defined (e.g.
//    K0 -> K0S) both values should be defined, to make for a uniform
//    approach in terms of loop constructions.
//
// PHKK(1,IHKK) : momentum in the x direction, in GeV/c.
//
// PHKK(2,IHKK) : momentum in the y direction, in GeV/c.
//
// PHKK(3,IHKK) : momentum in the z direction, in GeV/c.
//
// PHKK(4,IHKK) : energy, in GeV.
//
// PHKK(5,IHKK) : mass, in GeV/c**2. For spacelike partons, it is allowed
//    to use a negative mass, according to PHKK(5,IHKK) = -sqrt(-m**2).
//
// VHKK(1,IHKK) : production vertex x position, in mm.
//
// VHKK(2,IHKK) : production vertex y position, in mm.
//
// VHKK(3,IHKK) : production vertex z position, in mm.
//
// VHKK(4,IHKK) : production time, in mm/c (= 3.33*10**(-12) s).
//********************************************************************

#ifdef __cplusplus
}
#endif
#endif