!
!   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
!
MSEL=0     ! Your fate is in the following cards 

!  Upsilon (1S)

MSUB(86) = 1 ! gg->J/Psi g
KFPR(86,1) = 553 ! request Upsilon(1S)
MSUB(106) = 1 ! gg->J/Psi gamma
KFPR(106,1) = 553 ! request Upsilon(1S)
MSUB(107) = 1 ! g gamma->J/Psi g
KFPR(107,1) = 553 ! request Upsilon(1S)
PARP(38) = 15.73 ! GeV^3; from Upsilon(1S) Width

!  Upsilon (2S)

!MSUB(86) = 1 ! gg->J/Psi g
!KFPR(86,1) = 100553 ! request Upsilon(2S)
!MSUB(106) = 1 ! gg->J/Psi gamma
!KFPR(106,1) = 100553 ! request Upsilon(2S)
!MSUB(107) = 1 ! g gamma->J/Psi g
!KFPR(107,1) = 100553 ! request Upsilon(2S)
!PARP(38) = 14.79 ! GeV^3; from Upsilon(2S) Width

!  Upsilon (3S)

!MSUB(86) = 1 ! gg->J/Psi g
!KFPR(86,1) = 200553 ! request Upsilon(3S)
!MSUB(106) = 1 ! gg->J/Psi gamma
!KFPR(106,1) = 200553 ! request Upsilon(3S)
!MSUB(107) = 1 ! g gamma->J/Psi g
!KFPR(107,1) = 200553 ! request Upsilon(3S)
!PARP(38) = 9.44 ! GeV^3; from Upsilon(3S) Width

! Chi_0,1,2b

!MSUB(87) = 1 ! gg->Chi0c g
!KFPR(87,1) = 10551 ! request chi_0b
!MSUB(104) = 1 ! gg->Chi0c
!KFPR(104,1) = 10551 ! request chi_0b
!PARP(39) = 0.135 ! GeV^3; guessed from Upsilon width (???)

!MSUB(88) = 1 ! gg->Chi1c g
!KFPR(88,1) = 20553 ! request chi_1b

!MSUB(89) = 1 ! gg->Chi2c g
!KFPR(89,1) = 555 ! request chi_2b
!MSUB(105) = 1 ! gg->Chi2c 
!KFPR(105,1) = 555 ! request chi_2b


! Other Kinematic Inputs

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.
!
MSTP(81)=0     !TURN OFF MULTIPLE INTERACTIONS

!
! Force Upsilon(1S) decays to ee
!
MDME(1114,1) = 1 !   0    0.025200    e-              e+     
MDME(1115,1) = 0 !   0    0.024800    mu-             mu+    
MDME(1116,1) = 0 !   0    0.026700    tau-            tau+   
MDME(1117,1) = 0 !  32    0.015000    d               dbar   
MDME(1118,1) = 0 !  32    0.045000    u               ubar   
MDME(1119,1) = 0 !  32    0.015000    s               sbar   
MDME(1120,1) = 0 !  32    0.045000    c               cbar   
MDME(1121,1) = 0 !   4    0.774300    g               g      
MDME(1122,1) = 0 !   4    0.029000    gamma           g      
!
! Force Upsilon(2S) decays to ee
!
MDME(1658,1) = 1 !   0    0.014000    e-              e+                     
MDME(1659,1) = 0 !   0    0.014000    mu-             mu+                    
MDME(1660,1) = 0 !   0    0.014000    tau-            tau+                   
MDME(1661,1) = 0 !  32    0.008000    d               dbar                   
MDME(1662,1) = 0 !  32    0.024000    u               ubar                   
MDME(1663,1) = 0 !  32    0.008000    s               sbar                   
MDME(1664,1) = 0 !  32    0.024000    c               cbar                   
MDME(1665,1) = 0 !   4    0.425000    g               g               g      
MDME(1666,1) = 0 !   4    0.020000    gamma           g               g      
MDME(1667,1) = 0 !   0    0.185000    Upsilon         pi+             pi-    
MDME(1668,1) = 0 !   0    0.088000    Upsilon         pi0             pi0    
MDME(1669,1) = 0 !   0    0.043000    chi_0b          gamma                  
MDME(1670,1) = 0 !   0    0.067000    chi_1b          gamma                  
MDME(1671,1) = 0 !   0    0.066000    chi_2b          gamma                  
!
! Force Upsilon(3S) decays to ee
!
!   If I could, I would !


!
! Force Chi_0b decays to Upsilon(1S)
!
MDME(1600,1) = 1 ! 0.02  chi_b0 -> Upsilon gamma
MDME(1601,1) = 0 ! 0.98  chi_b0 -> g g
!
! Force Chi_1b decays to Upsilon(1S)
!
MDME(1123,1) = 1 ! 0.35  chi_b0 -> Upsilon gamma
MDME(1124,1) = 0 ! 0.65  chi_b0 -> g g
!
! Force Chi_2b decays to Upsilon(1S)
!
MDME(1600,1) = 1 ! 0.22  chi_b0 -> Upsilon gamma
MDME(1601,1) = 0 ! 0.78  chi_b0 -> g g
!
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).