bin/HARDCOL/parameter.f

0001 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
0002 c...this subroutine is used to set the necessary parameters for      c
0003 c...the initialization for hard color singlet exchange.              c
0004 c...to use the program youd need to make a directory: (data) to      c
0005 c...save all the obtained data-files.                                c
0006 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
0007 c   to have a better understanding of setting the parameters         c
0008 c   you may see the README file to get more detailed information.    c
0009 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
0010 c copyright (c) Rikard Enberg, Gunnar Ingelman, Leszek Motyka        c
0011 c reference: hep-ph/0111090                                          c
0012 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
0013
0014         SUBROUTINE SETPARAMETERS
0015 c...preamble: declarations.
0016         IMPLICIT DOUBLE PRECISION(A-H, O-Z)
0017         IMPLICIT INTEGER(I-N)
0018
0019       include "hardcol_set_par.inc"
0020
0021 c...user process event common block.
0022
0023       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0024       COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
0025       COMMON/PYDATR/MRPY(6),RRPY(100)
0026       COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
0027       COMMON/HCLPAR/ECM,NEV
0028       SAVE /PYPARS/,/PYDAT1/,/PYDATR/,/PYSUBS/,/HCLPAR/
0029
0030
0031 c       logical wronginput
0032 c
0033 c... Read parameter setting from file hardcol_set_par.nam
0034 c
0035         Call Read_parameter_settings
0036 c
0037 c... change the intitial state of the random number
0038 c        mrpy(1) = 19780503             ! default value
0039          mrpy(1) = irandom
0040          write( 6, * ) ' '
0041          write( 6, * ) ' Change default value of random, mrpy(1), to ',
0042      +   mrpy(1)
0043          write( 6, * ) ' '
0044 c... end random change
0045
0046 C       pi = dacos(-1.0d0)
0047
0048       ecm = ENERGYOFLHC
0049       nev = NUMOFEVENTS
0050       MSTP(2) = MSTP2
0051       CKIN(3) = CKIN3
0052       MSEL = MSEL0
0053       MSUB(406) = MSUB406
0054       MSUB(407) = MSUB407
0055       MSUB(408) = MSUB408
0056       MSTP(198) = MSTP198
0057
0058
0059 c...error message.
0060 c      wronginput=.false.
0061 c       CALL uperror(wronginput)
0062 c       if(wronginput) stop '-----input error! stop the program !!!'
0063
0064 c       CALL parameters()
0065 c      CALL dparameters()
0066 c      CALL coupling()
0067
0068         return
0069         end
0070
0071
0072 c***************************************
0073 c***************************************
0074         SUBROUTINE Read_parameter_settings
0075 c
0076 c... Get parameters from namelist
0077         implicit double precision(a-h, o-z)
0078         implicit integer(i-n)
0079
0080         Namelist / hardcol_set_par / ENERGYOFLHC, NUMOFEVENTS,
0081      +  MSTP2, CKIN3, MSEL0, MSUB406, MSUB407, MSUB408, MSTP198,
0082      +  irandom
0083
0084       include "hardcol_set_par.inc"
0085 c
0086 c-------------------------------------------------------------------------------
0087 c
0088         open( unit=1, file='hardcol_set_par.nam',Status='Old',Err=99)
0089         read( 1, nml=hardcol_set_par, err=90)
0090         write( 6, * ) ' '
0091         write( 6, * ) ' Contents of namelist *hardcol_set_par*: '
0092         write( 6, nml=hardcol_set_par)
0093         write( 6, * ) ' '
0094         Close( 1 )
0095         Return
0096 c
0097   90    Write( 6, * ) ' !!!!! Unable to read namelist hardcol_set_par '
0098         Call Exit
0099   99    Write( 6, * ) ' !!!!! Unable to open hardcol_set_par.nam'
0100         Call Exit
0101         End
0102