C
C  EXTERNAL CARDS FOR THE MONTE CARLO GENERATOR RIDI 
C  ----------------------------------------------------
C
C
C ------------------  Cards valid from version 2.0 --------------------
C
C ** Radiative corrections:  TRUE = yes
C    default = FALSE
C
RADCOR            TRUE
C
C ** Minimum energy of the radiated photon in case of radiative events
C    (it should be >0)
C    default = 0.03 GeV
C
EPHMIN            0.03
C
C ** Minimum angle of the radiated photon w.r.t. the incoming lepton 
C    (FSR events only; it should be >0)
C    default = 0.0001 radiant
C
THPHMN            0.0001
C
C ** Charm flavour only:  TRUE = yes
C    default = FALSE
C
ONLYCH            FALSE
C
C
C ------------------  Cards valid from version 1.0 --------------------
C
C ** Generation with weitghs:  TRUE = yes
C    default = FALSE
C
WEIGHT            TRUE
C
C ** q-qbar and/or q-qbar-gluon final state:  TRUE = yes
C    default = TRUE
C
QQBAR             TRUE
QQBARG            TRUE
C
C ** Transverse and/or longitudinal cross section:  TRUE = yes
C    default = TRUE
C
TRANSV            TRUE
LONGIT            TRUE
C
C ** Light flavours only:  TRUE = yes
C    default = FALSE
C
LIGHTF            FALSE
C
C ** High k_t only (q-qbar final state):  TRUE = yes
C    default = FALSE
C    (change it only if you really know what you are doing!)
C 
THLIM             FALSE
C
C ** Fragmentation:  TRUE = yes
C    default = TRUE
C
FRAGM             TRUE
C
C ** Enlarged output with check variables:  TRUE = yes
C    default = FALSE
C
CHECK             FALSE
C
C ** User output:  TRUE = yes
C    default = FALSE
C
USER              FALSE
C
C ** Number of final events 
C    default = 10000
C
NEVENT            10000
C
C ** Number of events to be used for the cross section determination 
C    (for generation of unweighted events only)
C    default = 1000
C
NEVWEI            1000
C
C ** Number of events to be printed
C    default = 0
C
NEVPRI            1
C
C ** K-factor choice: 0 = constant k-factors
C                     1 = exp[alphas(MX^2)*pi*C_F(A)]
C                     2 = exp[alphas(MX^2/4)*pi*C_F]       for qqbar
C                     2 = exp[alphas(MX^2/4)*pi*C_A]/3 +   for qqbarg
C                         exp[alphas(MX^2/4)*pi*C_F]*2/3
C                     3 = (1.+alphas(MX^2/4)*pi*C_F(A)/2.)**2
C    default = 2 (change it only if you really know what you are doing!)
C
IKFACT            2
C
C ** If IKFACT = 0 choose values of the constant K-factors:
C    first value for the q-qbar cross section  
C    second value for the q-qbar-gluon cross section
C    default = 2.1
C              4.
C    (change it only if you really know what you are doing!)
C
KFS             2.1            4.
C KFS             4            1.
C
C ** w_g=(1-z')^3*(2z'+1)^2 with z'=x_Bj/((ex_Bj/x_gluon)*xpom)
C    takes into account a more precise form of the triple pomeron vertex
C    at not too small values of z' (Levin&Wusthoff 1994)
C    0 = w_g NOT in the cross section evaluation 
C    1 = w_g in the cross section evaluation 
C    default = 1 (change it only if you really know what you are doing!)
C
IWG               1
C
C ** Type of parton shower for the q-qbar-gluon final state:
C    0 = parton shower only for the q-qbar system :qmax=sqrt(m_f^2+p_t^2)
C    1 = parton shower separately for the gluon jet           |
C        (qmax=sqrt(m_g^2+k_t^2)) and the q-qbar system     <-|
C    default = 1 (change it only if you really know what you are doing!)
C
IPSTYP            1
C
C ** Lepton choice:  1=positron; -1=electron
C    default = 1
C
LFLAG             1
C
C ** EBEAM: Lepton beam momentum (positive!)
C    default = 27.5
C
EBEAM             27.5
C
C ** PBEAM: Proton beam momentum (positive!)
C    default = 820.
C
PBEAM             820.
C
C ** Kinematical limits on variables: default values set to the limits.
C
C ..  Maximum x_Bjorken for diffractive events 
C       xbjddc < 0.1 in HERA kinematics:  xbjddc < xpom_max
C
XBJDDC            0.05
C
C ..  Q^2 , -(invariant mass of the virtual exchanged photon)
C       allowed range:    0.1 < qsqlim(1) < qsqlim(2) < 10000. GeV^2
C       (note that for values of Q^2<2 GeV^2 results may be effected
C        by large correction factors)
C
Q2LIM                         4.          1000.
C
C ..  t=(p-p')^2 , NEGATIVE invariant mass of the pomeron
C       allowed range:    -1. < tlim(1) < tlim(2) < -0.00001 GeV^2
C
TLIM                         -1.     -0.00001
C
C ..  Mx^2 , invariant mass squared of the gamma-pomeron system
C       allowed range:    3. < mx2lim(1) < mx2lim(2) < 10000. GeV^2
C
MX2LIM                        3.         10000.
C
C ..  x_pom , fraction of the pomeron momentum w.r.t. the proton one
C       allowed range:    0.00001 < xpomlim(1) < xpomlim(2) < 0.1
C
XPOLIM                           0.00001        0.05
C
C ..  k'_t^2 , transverse momentum of the hard gluon inside the pomeron
C              (for q-qbar-gluon final states only)
C       it should be  k'_t^2 > 0.001
C       default range:    0.001 < kt2lim(1) < kt2lim(2) < 100. GeV^2
C       (change it only if you really know what you are doing!)
C
KT2LIM                          0.001       100.
C
C ** Slope of the p_t^2 distribution of the scattered proton
C    default = 6. GeV^-2
C
BSLOPE            7.
C
C ** Charm mass 
C    default = 1.35 GeV (it should be chmass<mass of the D^+-)
C
CHMASS            1.35
C
C ** Gluon effective mass 
C    default = 0.75 GeV 
C    (change it only if you really know what you are doing!)
C
GMASS             0.75
C
C ** Cutoff for alpha_s: if(alpha_s.gt.alpscut)alpha_s=alpscut
C    (reasonable range [0.1-0.8])
C    default = 0.7
C    (change it only if you really know what you are doing!)
C
ALPSC           0.7
C
C ** Infrared cutoff for k'_t^2 (reasonable value in the range [0.04,0.2]) 
C    default = 0.06 GeV^2
C    (change it only if you really know what you are doing!)
C
KT2IC             0.06
C
C ** Gluon structure function type:  0 = user simple xg(x)
C                                    1 = use PDFLIB
C    default = 1
C
IGLUON            1
C
C ** If IGLUON = 0 choose parameters of the function 
C    xg(x)=gpar(1)*(1-x)**gpar(2)*(x/gpar(3))**gpar(4)
C    defaults: gpar(1)=3.; gpar(2)=5.; gpar(2)=0.01; gpar(2)=-0.4
C
GPAR              3.            5.          0.01          -0.4  
C
C ** If IGLUON = 1 choose gluon distribution according to PDFLIB 
C    (please consult the manual; special attention to Q^2_min):
C    gvalue(1) defines Ngroup; gvalue(2) defines Nset
C    defaults: gvalue(1)=5.; gvalue(2)=6. (GRV HO MSbar scheme)
C
GVAL                 3.            43.
C GVAL                 5.            6.
C
C ** Cutoff for the gluon structure function:
C    linear extrapolation for k'_t^2 < k_0^2
C    (reasonable value in the range [0.5-1.5])
C    default = 1. GeV^2
C    (change it only if you really know what you are doing!)
C
C With MRS-A k_0^2=0.625 is ok, since it's the lowest valid q^2 
GSFK02            0.625
C
C ** To optimaze variable generation in case of unweighted events
C    change here constants and powers of the generation functions 
C    (change them only if you really know what you are doing!)
C    
C ..  q-qbar diagram:
C
C     1/[x_bj^(PP+1)],    default = 0.5
C
PP                0.5
C
C     1/[Q^2^(RR+1)],    default = 1.
C
RR                1.
C
C     1/[(Mx^2+Q^2)^(SS+1)],    default = 1.
C
SS                1.
C
C ..  q-qbar-gluon diagram:
C
C
C     1/[(1-x_F)^(AA+1)],    default = 0.05
C     (x_F = p_z'/p_z)
C
AA                0.05
C
C     1/[(k'_t^2+m_g^2/BBB)^(BB+1)],    defaults BB = 1.5 
C                                                BBB = 0.025
BB                1.5
BBB               0.025
C
C     1/[(Q^2+CCC)^(CC+1)],    defaults CC = 1.5
C                                       CCC = 5.
CC                1.5
CCC               5.
C
C     1/[exbj*(exbj+EEE)],    default = 0.03
C     (exbj = x_bj in the electron-pomeron cms)
C
EEE               0.03
C
C     (1-z)^FF,    default = 0.
C     (z = exbj/x_gluon)
C
FF                0.
C
C **** JETSET data cards 
C
C MSTU     12=1
C MSTJ
C PARU
C PARJ
C
C
STOP