!
!   PROTON STRUCTURE FUNCTION SELECTION FROM PDFLIB
!
!   SET MSTP(52) = 2 FOR PDFLIB CHOICE
!   SET MSTP(51) = 1000 * NGROUP + NSET
!  EXAMPLES:
!  MRSD-=3031 ; MRSG=3041
!  CTEQ2M=4018 ; CTEQ2PM=4024 ; CTEQ3M=4030
!  GRVH0=5002 ; GRV94H0=5006
!
MSTP(51)=4030
MSTP(52)=2
!
! CHANGING RANDOM NUMBER GENERATOR SEED
!
MRPY(1) = RANSEED2      ! RANDOM SEED FIELD

! This set of cards roughly corresponds to the models explored by
! J. Lykken & K. Matchev, hep-ph/9903238


MSEL=0      ! choose your own fate

MSUB(271)=1 !    squark - squark processes 
MSUB(272)=1 ! 
MSUB(273)=1 ! 
MSUB(274)=1 ! 
MSUB(275)=1 ! 
MSUB(276)=1 ! 
MSUB(277)=1 ! 
MSUB(278)=1 ! 
MSUB(279)=1 ! 
MSUB(280)=1 ! 

! set SUSY model:

IMSS(1)=1   ! SUGRA [2] or general spectrum [1] 
IMSS(2)=0   ! Gaugino masses are input by hand
IMSS(3)=1   ! M3 treated as the gluino pole mass
IMSS(4)=0   ! Higgs sector determined by PYTHIA from ma and tanb
IMSS(7)=0   ! No extra D-term contributions 
IMSS(8)=0   ! [0]-Staus are treated separately
IMSS(9)=1   ! Right-handed squarks are input by hand
IMSS(10)=0   ! [0]-N2 decays correctly
IMSS(11)=0   ! [0]-SUGRA; [1]-GMed


!RMSS array: (note that all of this will be recomputed by Pythia
!             taking into account mixings & stuff - so to get the 
!             physical masses & BR you need to print out the whole
!             thing with pylist)

RMSS(1)=75.0              ! U(1) gaugino mass
RMSS(2)=141.473648        ! SU(2) gaugino mass
RMSS(3)=454.29306         ! SU(3) gaugino mass
RMSS(4)=245.527969        ! mu, higgsino mass parameter
RMSS(5)=10.               ! tan beta
RMSS(6)=200.0             ! left slepton mass
RMSS(7)=200.0             ! right slepton mass
RMSS(8)=150.0             ! left squark mass
RMSS(9)=150.0             ! right squark mass
RMSS(10)=400.0            ! left squark mass - 3rd generation
RMSS(11)=400.0            ! right sbottom mass
RMSS(12)=400.0            ! right stop mass
RMSS(13)=450.0            ! left stau mass
RMSS(14)=450.0            ! right stau mass
RMSS(15)=-438.160217      ! bottom trilinear coupling A_b
RMSS(16)=-335.252197      ! top trilinear coupling A_t
RMSS(17)=-53.6412468      ! tau trilinear coupling A_tau
RMSS(18)=-0.0268513039    ! Higgs mixing angle alpha - not used in our case
RMSS(19)=500              ! M_A
RMSS(20)=  0.0417031534   ! GUT scale coupling constant
RMSS(21)=1.               ! Gravitino mass in eV
RMSS(22)=436.741058       ! Right up squark mass
RMSS(23)=0.               ! D_x term contribution (not used here)
RMSS(24)=0.               ! D_y term contribution (not used here)
RMSS(25)=0.               ! D_s term contribution (not used here)


!
CKIN(3)=PTLOW           ! Minimum Pt for hard 2->2 process.
CKIN(4)=PTHIGH             ! Maximum Pt for hard 2->2 process.
CKIN(9)=-4.2              ! Rapidity range for product with
CKIN(10)=4.2              ! largest rapidity in hard 2->2 process.

!
! SELECTING PYTHIA KINEMATIC CUTS
!
MSTP(81)=0     !TURN OFF MULTIPLE INTERACTIONS
!

!MDCY(C111,1)=0    !SETS PI ZERO STABLE (LIMITS SIZE OF OUTPUT)
MDCY(15,1)=0 ! sets tau stable - this is important (for tauola)!
! The next line is _mandatory_ for RunII production
MSTJ(22)=2 !A particle is decayed only if its proper
!lifetime is smaller than PARJ(71)
! default should be ok. (strange particles stable, charm & bottom
! decay).