From 5ebdc6076f7031dd3bbde98768e50d931092ce8a Mon Sep 17 00:00:00 2001
From: Chris Green <greenc@fnal.gov>
Date: Tue, 9 Oct 2018 13:56:43 -0500
Subject: [PATCH] New package: pythia6 (#9340)
Includes optional extra code for compatibility with Root.
---
.../builtin/packages/pythia6/CMakeLists.txt | 115 ++++++++
.../repos/builtin/packages/pythia6/package.py | 172 +++++++++++
.../builtin/packages/pythia6/pythia6.patch | 278 ++++++++++++++++++
3 files changed, 565 insertions(+)
create mode 100644 var/spack/repos/builtin/packages/pythia6/CMakeLists.txt
create mode 100644 var/spack/repos/builtin/packages/pythia6/package.py
create mode 100644 var/spack/repos/builtin/packages/pythia6/pythia6.patch
diff --git a/var/spack/repos/builtin/packages/pythia6/CMakeLists.txt b/var/spack/repos/builtin/packages/pythia6/CMakeLists.txt
new file mode 100644
index 0000000000..9362f03668
--- /dev/null
+++ b/var/spack/repos/builtin/packages/pythia6/CMakeLists.txt
@@ -0,0 +1,115 @@
+# ======================================================================
+# pythia6 main build file
+#
+# setup cmake
+# cd .../path/to/build/directory
+# cmake [-DCMAKE_INSTALL_PREFIX=/install/path]
+# [-DCMAKE_BUILD_TYPE=<RelWithDebInfo|Debug|Release|MinSizeRel> ]
+# [-Drun_long_tests:BOOL=ON]
+# .../path/to/pythia/source
+# make
+# make test
+# make install
+# ======================================================================
+
+# use cmake 2.6 or later
+cmake_minimum_required (VERSION 2.6)
+
+project(pythia6 C Fortran)
+message(STATUS "pythia version is ${PYTHIA6_VERSION}")
+
+enable_testing()
+
+#build all libraries in a single directory to enable testing
+set(LIBRARY_OUTPUT_PATH ${PROJECT_BINARY_DIR}/lib)
+
+# set compiler flags
+# default GNU compiler flags:
+# DEBUG -g
+# RELEASE -O3 -DNDEBUG
+# MINSIZEREL -Os -DNDEBUG
+# RELWITHDEBINFO -O2 -g
+set( CMAKE_C_FLAGS_DEBUG "-g -O0" )
+set( CMAKE_C_FLAGS_RELWITHDEBINFO "-O3 -g -DNDEBUG -fno-omit-frame-pointer" )
+set( CMAKE_C_FLAGS_RELEASE "-O3 -DNDEBUG" )
+set( CMAKE_Fortran_FLAGS_DEBUG "-g -O0 -fno-second-underscore" )
+set( CMAKE_Fortran_FLAGS_RELWITHDEBINFO "-O3 -g -DNDEBUG -fno-omit-frame-pointer -fno-second-underscore" )
+set( CMAKE_Fortran_FLAGS_RELEASE "-O3 -DNDEBUG -fno-second-underscore" )
+
+message(STATUS "CMAKE_Fortran_COMPILER_INIT = ${CMAKE_Fortran_COMPILER_INIT}")
+message(STATUS "CMAKE_Fortran_COMPILER_FULLPATH = ${CMAKE_Fortran_COMPILER_FULLPATH}")
+message(STATUS "CMAKE_Fortran_COMPILER = ${CMAKE_Fortran_COMPILER}")
+
+if(NOT CMAKE_BUILD_TYPE)
+ set(CMAKE_BUILD_TYPE RelWithDebInfo CACHE STRING "" FORCE)
+endif()
+message(STATUS "cmake build type set to ${CMAKE_BUILD_TYPE}")
+
+message("ENV_FLAGS = $ENV{FFLAGS}")
+string(TOUPPER ${CMAKE_BUILD_TYPE} BTYPE_UC )
+if( ${BTYPE_UC} MATCHES "DEBUG")
+ set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS_DEBUG}")
+ set(CMAKE_Fortran_FLAGS "${CMAKE_Fortran_FLAGS_DEBUG}")
+elseif( ${BTYPE_UC} MATCHES "RELEASE")
+ set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS_RELEASE}")
+ set(CMAKE_Fortran_FLAGS "${CMAKE_Fortran_FLAGS_RELEASE}")
+elseif( ${BTYPE_UC} MATCHES "RELWITHDEBINFO")
+ set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS_RELWITHDEBINFO}")
+ set(CMAKE_Fortran_FLAGS "${CMAKE_Fortran_FLAGS_RELWITHDEBINFO}")
+endif()
+set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS} $ENV{CFLAGS}")
+set(CMAKE_Fortran_FLAGS "${CMAKE_Fortran_FLAGS} $ENV{FFLAGS}")
+message("CMAKE_C_FLAGS = ${CMAKE_C_FLAGS}")
+message("CMAKE_Fortran_FLAGS = ${CMAKE_Fortran_FLAGS}")
+
+# source
+# shared library
+file(GLOB src_files *.c *.F *.f)
+add_library(Pythia6 SHARED ${src_files})
+
+# Installation:
+# Library.
+install(TARGETS Pythia6 DESTINATION lib)
+
+# Include-able file.
+install( FILES example/main60.f
+ DESTINATION include )
+
+# Documentation.
+install(DIRECTORY doc DESTINATION .)
+
+# Examples
+install(DIRECTORY example DESTINATION .)
+
+# tests
+macro( pythia_test testname )
+ set ( package_library_list ${PROJECT_BINARY_DIR}/lib/${CMAKE_STATIC_LIBRARY_PREFIX}CLHEP-${PACKAGE}-${VERSION}${CMAKE_STATIC_LIBRARY_SUFFIX} )
+ link_libraries( Pythia6 )
+ message( STATUS "building ${testname} from ${CMAKE_CURRENT_SOURCE_DIR}/example in ${CMAKE_CURRENT_BINARY_DIR}" )
+ add_executable(${testname} example/${testname}.f)
+ add_test( ${testname} ${CMAKE_CURRENT_BINARY_DIR}/${testname} )
+endmacro( pythia_test )
+
+pythia_test( main61 )
+pythia_test( main63 )
+pythia_test( main66 )
+pythia_test( main67 )
+pythia_test( main68 )
+pythia_test( main69 )
+pythia_test( main71 )
+pythia_test( main72 )
+pythia_test( main73 )
+pythia_test( main75 )
+pythia_test( main77 )
+pythia_test( main78 )
+pythia_test( main81 )
+
+# these examples take a while to run
+if( run_long_tests )
+ pythia_test( main62 )
+ pythia_test( main64 )
+ pythia_test( main65 )
+ pythia_test( main70 )
+ pythia_test( main74 )
+ pythia_test( main79 )
+endif( run_long_tests )
diff --git a/var/spack/repos/builtin/packages/pythia6/package.py b/var/spack/repos/builtin/packages/pythia6/package.py
new file mode 100644
index 0000000000..7967092934
--- /dev/null
+++ b/var/spack/repos/builtin/packages/pythia6/package.py
@@ -0,0 +1,172 @@
+##############################################################################
+# Copyright (c) 2013-2016, Lawrence Livermore National Security, LLC.
+# Produced at the Lawrence Livermore National Laboratory.
+#
+# This file is part of Spack.
+# Created by Todd Gamblin, tgamblin@llnl.gov, All rights reserved.
+# LLNL-CODE-647188
+#
+# For details, see https://github.com/llnl/spack
+# Please also see the LICENSE file for our notice and the LGPL.
+#
+# This program is free software; you can redistribute it and/or modify
+# it under the terms of the GNU Lesser General Public License (as
+# published by the Free Software Foundation) version 2.1, February 1999.
+#
+# This program is distributed in the hope that it will be useful, but
+# WITHOUT ANY WARRANTY; without even the IMPLIED WARRANTY OF
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the terms and
+# conditions of the GNU Lesser General Public License for more details.
+#
+# You should have received a copy of the GNU Lesser General Public
+# License along with this program; if not, write to the Free Software
+# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+from spack import *
+import numbers
+import os
+from six import iteritems
+from six.moves.urllib.parse import urlparse
+
+
+def _is_integral(x):
+ """Accepts only integral values."""
+ try:
+ return isinstance(int(x), numbers.Integral) and \
+ (not isinstance(x, bool)) and \
+ int(x) == x
+ except ValueError:
+ return False
+
+
+class Pythia6(CMakePackage):
+ """PYTHIA is a program for the generation of high-energy physics events,
+ i.e. for the description of collisions at high energies between elementary
+ particles such as e+, e-, p and pbar in various combinations.
+
+ PYTHIA6 is a Fortran package which is no longer maintained: new
+ prospective users should use Pythia8 instead.
+
+ This recipe includes patches required to interoperate with Root.
+ """
+
+ homepage = 'https://pythiasix.hepforge.org/'
+ url = 'http://www.hepforge.org/archive/pythiasix/pythia-6.4.28.tgz'
+
+ version('6.4.28',
+ sha256='01cbff47e99365b5e46f6d62c1735d3cae1932c4710604850d59f538cb758020')
+
+ # Root's TPythia6 interface requires extra sources to be built into
+ # the Pythia6 library.
+ variant('root', default=False,
+ description='Build extra (non OEM) code to allow use by Root.')
+
+ # The maximum number of particles (NMXHEP) supported by the arrays
+ # in the /HEPEVT/ COMMON block may need tweaking if pythia6 is
+ # intended to be used with other code with different requirements.
+ variant('nmxhep', default=4000, values=_is_integral, description='Extent of particle arrays in the /HEPEVT/ COMMON block.')
+
+ # In the unlikely event of new versions >6.4.28,
+ # pythia6_common_address.c should be checked for accuracy against
+ # the definitions of the relevant COMMON blocks in the Pythia6
+ # Fortran source, and patched if necessaary.
+ resource(
+ name='root-pythia6-shim',
+ url='https://root.cern.ch/download/pythia6.tar.gz',
+ sha256='d613dcb27c905710e2f13a934913cc5545e3e5d0e477e580107385d9ef260056',
+ when='+root',
+ destination='.',
+ placement={'pythia6_common_address.c': 'pythia6_common_address.c',
+ 'tpythia6_called_from_cc.F': 'tpythia6_called_from_cc.F'}
+ )
+
+ # Download examples separately.
+ examples \
+ = {'main60.f':
+ 'd713b8b267c4405cc9d31c58bba267ae3378902a26fa52393003bf35fd56902c',
+ 'main61.f':
+ 'e2a3d5524e43d16f60d9edc6e7198d41006d1ba127fb7b0e265aa509e13128b4',
+ 'main62.f':
+ 'dce822a72fe2d6cfb6d43c479ba98928fb0a39290a6ee26fdcacc66229313045',
+ 'main63.f':
+ 'b2dd343b3cd7969979b80c564d82b92e0d776d66bb19d346b52f2af27adeb62d',
+ 'main64.f':
+ 'a35f2f232e6e0d68d67fd350d4d46b0a353f5c7811de0c2db47ae16d17ed1843',
+ 'main65.f':
+ '03c81e0bbd77710b0461e18265e80e3bd51360b9f416c86013401f882ac39a5e',
+ 'main66.f':
+ '50dd9221a7e84ee7c5005db6758e5880d190eab8cce8a52e7c7b29e9fee8d3da',
+ 'main67.f':
+ '1984aa90fe4e3d628c3bcceaa6fca1b08231d835158d975fa171337d55ca4a2f',
+ 'main68.f':
+ 'c8d6def1298477ffec6a1d98c7e02dcee0debe6badc7c63f752f9194b82f212d',
+ 'main69.f':
+ 'd14399d43f8c4b670907558849d3e5a4d7625d027de3c10002185c58b20b061a',
+ 'main71.f':
+ '2e47af778003b0596e8999f0914033c6eda7335211b9e96ac3075d45a3cde12e',
+ 'main72.f':
+ 'e27ce2af68b40436c51c65767ebb5ff0955ab8dfdfc5fc5c217ae73cd53070da',
+ 'main73.f':
+ '567db2d1a66896ce5103ffa7e10742442b0e934088883e91339536e0249772c4',
+ 'main75.f':
+ 'b850986c43a5af1e7d13b66d22b01584e3c68bb338be32eac39e31f971b80be4',
+ 'main77.f':
+ '0679852c4f35719531ad38dc1dbb374b884181eb5e483c36d8867ccb449177a4',
+ 'main78.f':
+ '5babc59fe6a0bd57d97ec398cf01745bc9b72ce6ce0711e934d53c7821e21912',
+ 'main79.f':
+ '27ca84d6d0877f3605cbc1b865c3e1f571e7d2c9301094a4122e726a903dbead',
+ 'main81.f':
+ 'b02fecd1cd0f9ba16eaae53e9da0ba602569fdf0e46856cccdfb4c5b7ba33e8b',
+ 'ttbar.lhe':
+ 'fb0d43175cc392b19c2b6633dcf673d0b56229b60bec92df4aa782c7196b149c'}
+
+ for example, checksum in iteritems(examples):
+ resource(name=example,
+ url='http://pythiasix.hepforge.org/examples/{0}'.
+ format(example),
+ sha256=checksum,
+ expand=False,
+ destination='example',
+ placement={example: example}
+ )
+
+ # Docs.
+ docs \
+ = {'http://www.hepforge.org/archive/pythiasix/update_notes-6.4.28.txt':
+ 'a229be4ba9a4eb65a9d53600a5f388b620038d56694c6cb4671c2be224b67751',
+ 'http://home.thep.lu.se/~torbjorn/pythia6/lutp0613man2.pdf':
+ '03d637310ea80f0d7aea761492bd38452c602890d8cf913a1ec9edacd79fa43d',
+ 'https://pythiasix.hepforge.org/pythia6-announcement.txt':
+ '2a52def41f0c93e32e0db58dbcf072b987ebfbd32e42ccfc1f9382fcf65f1271'}
+
+ for docurl, checksum in iteritems(docs):
+ doc = os.path.basename(urlparse(docurl).path)
+ resource(name=doc,
+ url=docurl,
+ sha256=checksum,
+ expand=False,
+ destination='doc',
+ placement={doc: doc}
+ )
+
+ # The included patch customizes some routines provided in dummy form
+ # by the original source to be useful out of the box in the vast
+ # majority of cases. If your case is different, platform- or
+ # variant-based adjustments should be made.
+ patch('pythia6.patch', level=0)
+
+ def patch(self):
+ # Use our provided CMakeLists.txt. The Makefile provided with
+ # the source is GCC (gfortran) specific, and would have required
+ # additional patching for the +root variant.
+ llnl.util.filesystem.copy(os.path.join(os.path.dirname(__file__),
+ 'CMakeLists.txt'),
+ self.stage.source_path)
+ # Apply the variant value at the relevant place in the source.
+ filter_file(r'^(\s+PARAMETER\s*\(\s*NMXHEP\s*=\s*)\d+',
+ r'\1{0}'.format(self.spec.variants['nmxhep'].value),
+ 'pyhepc.f')
+
+ def cmake_args(self):
+ args = ['-DPYTHIA6_VERSION={0}'.format(self.version.dotted)]
+ return args
diff --git a/var/spack/repos/builtin/packages/pythia6/pythia6.patch b/var/spack/repos/builtin/packages/pythia6/pythia6.patch
new file mode 100644
index 0000000000..5d45541e07
--- /dev/null
+++ b/var/spack/repos/builtin/packages/pythia6/pythia6.patch
@@ -0,0 +1,278 @@
+diff -Naur pytime.f pytime.f
+--- pytime.f 2018-08-23 11:24:27.000000000 -0500
++++ pytime.f 2018-08-23 11:46:08.189999826 -0500
+@@ -1,4 +1,3 @@
+-
+ C*********************************************************************
+
+ C...PYTIME
+@@ -59,17 +58,18 @@
+ C IDATI(6)=ISEC
+
+ C...Example 4: GNU LINUX libU77, SunOS.
+-C CALL IDATE(IDTEMP)
+-C IDATI(1)=IDTEMP(3)
+-C IDATI(2)=IDTEMP(2)
+-C IDATI(3)=IDTEMP(1)
+-C CALL ITIME(IDTEMP)
+-C IDATI(4)=IDTEMP(1)
+-C IDATI(5)=IDTEMP(2)
+-C IDATI(6)=IDTEMP(3)
++ CALL IDATE(IDTEMP)
++ IDATI(1)=IDTEMP(3)
++ IDATI(2)=IDTEMP(2)
++ IDATI(3)=IDTEMP(1)
++ CALL ITIME(IDTEMP)
++ IDATI(4)=IDTEMP(1)
++ IDATI(5)=IDTEMP(2)
++ IDATI(6)=IDTEMP(3)
+
+ C...Common code to ensure right century.
+ IDATI(1)=2000+MOD(IDATI(1),100)
+
+ RETURN
+ END
++
+diff -Naur upevnt.f upevnt.f
+--- upevnt.f 2018-08-23 11:24:27.000000000 -0500
++++ upevnt.f 2018-08-23 11:46:08.189999826 -0500
+@@ -1,56 +1,3 @@
+-
+-C...Old example: handles a simple Pythia 6.4 initialization file.
+-
+-c SUBROUTINE UPINIT
+-
+-C...Double precision and integer declarations.
+-c IMPLICIT DOUBLE PRECISION(A-H, O-Z)
+-c IMPLICIT INTEGER(I-N)
+-
+-C...Commonblocks.
+-c COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
+-c COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
+-c SAVE /PYDAT1/,/PYPARS/
+-
+-C...User process initialization commonblock.
+-c INTEGER MAXPUP
+-c PARAMETER (MAXPUP=100)
+-c INTEGER IDBMUP,PDFGUP,PDFSUP,IDWTUP,NPRUP,LPRUP
+-c DOUBLE PRECISION EBMUP,XSECUP,XERRUP,XMAXUP
+-c COMMON/HEPRUP/IDBMUP(2),EBMUP(2),PDFGUP(2),PDFSUP(2),
+-c &IDWTUP,NPRUP,XSECUP(MAXPUP),XERRUP(MAXPUP),XMAXUP(MAXPUP),
+-c &LPRUP(MAXPUP)
+-c SAVE /HEPRUP/
+-
+-C...Read info from file.
+-c IF(MSTP(161).GT.0) THEN
+-c READ(MSTP(161),*,END=110,ERR=110) IDBMUP(1),IDBMUP(2),EBMUP(1),
+-c & EBMUP(2),PDFGUP(1),PDFGUP(2),PDFSUP(1),PDFSUP(2),IDWTUP,NPRUP
+-c DO 100 IPR=1,NPRUP
+-c READ(MSTP(161),*,END=110,ERR=110) XSECUP(IPR),XERRUP(IPR),
+-c & XMAXUP(IPR),LPRUP(IPR)
+-c 100 CONTINUE
+-c RETURN
+-C...Error or prematurely reached end of file.
+-c 110 WRITE(MSTU(11),5000)
+-c STOP
+-
+-C...Else not implemented.
+-c ELSE
+-c WRITE(MSTU(11),5100)
+-c STOP
+-c ENDIF
+-
+-C...Format for error printout.
+-c 5000 FORMAT(1X,'Error: UPINIT routine failed to read information'/
+-c &1X,'Execution stopped!')
+-c 5100 FORMAT(1X,'Error: You have not implemented UPINIT routine'/
+-c &1X,'Dummy routine in PYTHIA file called instead.'/
+-c &1X,'Execution stopped!')
+-
+-c RETURN
+-c END
+-
+ C*********************************************************************
+
+ C...UPEVNT
+@@ -120,3 +67,54 @@
+
+ RETURN
+ END
++
++C...Old example: handles a simple Pythia 6.4 event file.
++
++c SUBROUTINE UPEVNT
++
++C...Double precision and integer declarations.
++c IMPLICIT DOUBLE PRECISION(A-H, O-Z)
++c IMPLICIT INTEGER(I-N)
++
++C...Commonblocks.
++c COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
++c COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
++c SAVE /PYDAT1/,/PYPARS/
++
++C...User process event common block.
++c INTEGER MAXNUP
++c PARAMETER (MAXNUP=500)
++c INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
++c DOUBLE PRECISION XWGTUP,SCALUP,AQEDUP,AQCDUP,PUP,VTIMUP,SPINUP
++c COMMON/HEPEUP/NUP,IDPRUP,XWGTUP,SCALUP,AQEDUP,AQCDUP,IDUP(MAXNUP),
++c &ISTUP(MAXNUP),MOTHUP(2,MAXNUP),ICOLUP(2,MAXNUP),PUP(5,MAXNUP),
++c &VTIMUP(MAXNUP),SPINUP(MAXNUP)
++c SAVE /HEPEUP/
++
++C...Read info from file.
++c IF(MSTP(162).GT.0) THEN
++c READ(MSTP(162),*,END=110,ERR=110) NUP,IDPRUP,XWGTUP,SCALUP,
++c & AQEDUP,AQCDUP
++c DO 100 I=1,NUP
++c READ(MSTP(162),*,END=110,ERR=110) IDUP(I),ISTUP(I),
++c & MOTHUP(1,I),MOTHUP(2,I),ICOLUP(1,I),ICOLUP(2,I),
++c & (PUP(J,I),J=1,5),VTIMUP(I),SPINUP(I)
++c 100 CONTINUE
++c RETURN
++C...Special when reached end of file or other error.
++c 110 NUP=0
++
++C...Else not implemented.
++c ELSE
++c WRITE(MSTU(11),5000)
++c STOP
++c ENDIF
++
++C...Format for error printout.
++c 5000 FORMAT(1X,'Error: You have not implemented UPEVNT routine'/
++c &1X,'Dummy routine in PYTHIA file called instead.'/
++c &1X,'Execution stopped!')
++
++c RETURN
++c END
++
+diff -Naur upinit.f upinit.f
+--- upinit.f 2018-08-23 11:24:27.000000000 -0500
++++ upinit.f 2018-08-23 11:46:08.190999819 -0500
+@@ -1,4 +1,3 @@
+-
+ C*********************************************************************
+
+ C...UPINIT
+@@ -64,3 +63,56 @@
+
+ RETURN
+ END
++
++C...Old example: handles a simple Pythia 6.4 initialization file.
++
++c SUBROUTINE UPINIT
++
++C...Double precision and integer declarations.
++c IMPLICIT DOUBLE PRECISION(A-H, O-Z)
++c IMPLICIT INTEGER(I-N)
++
++C...Commonblocks.
++c COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
++c COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
++c SAVE /PYDAT1/,/PYPARS/
++
++C...User process initialization commonblock.
++c INTEGER MAXPUP
++c PARAMETER (MAXPUP=100)
++c INTEGER IDBMUP,PDFGUP,PDFSUP,IDWTUP,NPRUP,LPRUP
++c DOUBLE PRECISION EBMUP,XSECUP,XERRUP,XMAXUP
++c COMMON/HEPRUP/IDBMUP(2),EBMUP(2),PDFGUP(2),PDFSUP(2),
++c &IDWTUP,NPRUP,XSECUP(MAXPUP),XERRUP(MAXPUP),XMAXUP(MAXPUP),
++c &LPRUP(MAXPUP)
++c SAVE /HEPRUP/
++
++C...Read info from file.
++c IF(MSTP(161).GT.0) THEN
++c READ(MSTP(161),*,END=110,ERR=110) IDBMUP(1),IDBMUP(2),EBMUP(1),
++c & EBMUP(2),PDFGUP(1),PDFGUP(2),PDFSUP(1),PDFSUP(2),IDWTUP,NPRUP
++c DO 100 IPR=1,NPRUP
++c READ(MSTP(161),*,END=110,ERR=110) XSECUP(IPR),XERRUP(IPR),
++c & XMAXUP(IPR),LPRUP(IPR)
++c 100 CONTINUE
++c RETURN
++C...Error or prematurely reached end of file.
++c 110 WRITE(MSTU(11),5000)
++c STOP
++
++C...Else not implemented.
++c ELSE
++c WRITE(MSTU(11),5100)
++c STOP
++c ENDIF
++
++C...Format for error printout.
++c 5000 FORMAT(1X,'Error: UPINIT routine failed to read information'/
++c &1X,'Execution stopped!')
++c 5100 FORMAT(1X,'Error: You have not implemented UPINIT routine'/
++c &1X,'Dummy routine in PYTHIA file called instead.'/
++c &1X,'Execution stopped!')
++
++c RETURN
++c END
++
+diff -Naur upveto.f upveto.f
+--- upveto.f 2018-08-23 11:24:27.000000000 -0500
++++ upveto.f 2018-08-23 11:46:08.190999819 -0500
+@@ -1,54 +1,3 @@
+-
+-C...Old example: handles a simple Pythia 6.4 event file.
+-
+-c SUBROUTINE UPEVNT
+-
+-C...Double precision and integer declarations.
+-c IMPLICIT DOUBLE PRECISION(A-H, O-Z)
+-c IMPLICIT INTEGER(I-N)
+-
+-C...Commonblocks.
+-c COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
+-c COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
+-c SAVE /PYDAT1/,/PYPARS/
+-
+-C...User process event common block.
+-c INTEGER MAXNUP
+-c PARAMETER (MAXNUP=500)
+-c INTEGER NUP,IDPRUP,IDUP,ISTUP,MOTHUP,ICOLUP
+-c DOUBLE PRECISION XWGTUP,SCALUP,AQEDUP,AQCDUP,PUP,VTIMUP,SPINUP
+-c COMMON/HEPEUP/NUP,IDPRUP,XWGTUP,SCALUP,AQEDUP,AQCDUP,IDUP(MAXNUP),
+-c &ISTUP(MAXNUP),MOTHUP(2,MAXNUP),ICOLUP(2,MAXNUP),PUP(5,MAXNUP),
+-c &VTIMUP(MAXNUP),SPINUP(MAXNUP)
+-c SAVE /HEPEUP/
+-
+-C...Read info from file.
+-c IF(MSTP(162).GT.0) THEN
+-c READ(MSTP(162),*,END=110,ERR=110) NUP,IDPRUP,XWGTUP,SCALUP,
+-c & AQEDUP,AQCDUP
+-c DO 100 I=1,NUP
+-c READ(MSTP(162),*,END=110,ERR=110) IDUP(I),ISTUP(I),
+-c & MOTHUP(1,I),MOTHUP(2,I),ICOLUP(1,I),ICOLUP(2,I),
+-c & (PUP(J,I),J=1,5),VTIMUP(I),SPINUP(I)
+-c 100 CONTINUE
+-c RETURN
+-C...Special when reached end of file or other error.
+-c 110 NUP=0
+-
+-C...Else not implemented.
+-c ELSE
+-c WRITE(MSTU(11),5000)
+-c STOP
+-c ENDIF
+-
+-C...Format for error printout.
+-c 5000 FORMAT(1X,'Error: You have not implemented UPEVNT routine'/
+-c &1X,'Dummy routine in PYTHIA file called instead.'/
+-c &1X,'Execution stopped!')
+-
+-c RETURN
+-c END
+-
+ C*********************************************************************
+
+ C...UPVETO
+@@ -98,3 +47,4 @@
+
+ RETURN
+ END
++
--
GitLab