&CASCADE_input
NrEvents = 1000      ! Nr of events to process
PBE1 =-27.5          ! Beam1 energy    
KBE1 =-11            ! Flavor beam1 -11: positron, 22: photon  2212: proton 
IRE1 = 0             ! 0: beam 1 has no structure 
PBE2 = 920.          ! Beam2 energy
KBE2 = 2212          ! Flavor beam2 11: electron, 22: photon 2212: proton
IRE2 = 1             ! 0:  beam 2 has no structure 
QMIN = 2.            ! (D=5.0) minimum Q2 to be generated
! QMAX = 100.        ! (D=10.D8) maximum Q2 to be generated
! YMIN = 0.01        ! (D=0.) YMIN of photon
! YMAX = 0.8         ! (D=1.) YMAX of photon
! THEI = 0.          ! (D=0.) minimum scattering angle for e
! THEM = 180.        ! (D=180.) maximum scattering angle for e
NFLA = 5             ! (D=5) nr of flavours used in str.fct
MULT = 0             ! (D=0) multiple scattering
Process_Id = 10      ! gamma g -> qq_bar
IRPA = 1             !  gamma g --> qq_bar
IRPB = 0             !  g + g --> g + g
IRPC = 0             !  g + q --> g + q
IHFL = 5             ! (D=4) produced flavour 
! PTCU = 9           ! (D=0) p_t **2 cut for process
Hadronisation = 1    ! Hadronisation on (=1)
SpaceShower = 1      ! Space-like Parton Shower
TimeShower = 1       ! Time-like Parton Shower
ScaleTimeShower = 4  ! Scale choice for Time-like Shower
!                    1: 2(m^2_1t+m^2_2t)    
!                    2: shat     
!                    3: 2(m^2_1+m^2_2)     
!                    4: 2*scalup (from lhe file)     
!ScaleFactorFinalShower = 1.               ! scale factor for 
!                                            Final State Parton Shower  
PartonEvolution = 1                        ! type of parton evolution 
!                                            in Space-like Shower
!                                            1: CCFM
!                                            2: full all flavor TMD evolution
! EnergyShareRemnant = 4                   ! energy sharing in proton remnant
!                                            1: (a+1)(1-z)**a <z>=1/(a+2)=1/3
!                                            2: (a+1)(1-z)**a
!                                                <z>=1/(a+2)=mq/(mq+mQ
!                                            3: N/(z(1-1/z-c/(1-z))**2)
!                                               c=(mq/mQ)**2
!                                            4: PYZDIS: KFL1=1                   IRAM = 0             ! (D=0) Running of alpha_em(Q2)
IRAS = 1             ! (D=1) Running of alpha_s(MU2)
IQ2S = 3             ! (D=1) Scale MU2 of alpha_s
SCAL = 1.            !  scale factor for renormalisation scale 
SCAF = 1.            !  scale factor for factorisation scale
PartonDensity = 1201 ! use ONLY internal CCFM set
! TMDDensityPath= '../UPDF'                ! Path to TMD density 
!                                            for internal files
! IINT = 1           ! (D=0) select MC integration procedure
NCAL = 20000         ! (D=20000) Nr of calls per iteration for bases
ACC1 = 1.            ! (D=1) relative prec.(%) for grid optimisation
ACC2 = 0.5           ! (0.5) relative prec.(%) for integration
! INTE = 0           ! Interaction type (D=0)
! KT1 = 0.44         ! (D=0.0) intrinsic kt for beam 1
! KT2 = 0.44         ! (D=0.0) intrinsic kt for beam 2
! KTRE = 0.35        ! (D=0.35) primordial kt when non-trivial target remnant is split into two particles
!HepMCOUTPUT = TRUE
!HepMCFormat = 2
!HepMCFile = 'output.hepmc'
!RivetOUT = TRUE
!RivetRun = 'somerun'
!RivetFile = 'steering_ep_output.yoda'
!RivetAnalyses='MC_XS','MC_PDFS','MC_DIS_Check'
&End 

&PYTHIA6_input
! P6_MSTJ(41) = 2     !  (D = 2) type of branchings allowed in shower. =1 only QCD
P6_PMAS(4,1)= 1.6     ! charm mass
P6_PMAS(5,1)= 4.75    ! bottom mass
P6_PMAS(25,1)= 125.   ! higgs mass
P6_PARU(112) = 0.2    ! lambda QCD
P6_MSTU(111) =  2     ! = 0 : alpha_s is fixed, =1 first order; =2 2nd order;
P6_MSTU(112) = 4      ! nr of flavours wrt lambda_QCD
P6_MSTU(113) = 3      ! min nr of flavours for alphas
P6_MSTU(114) = 5      ! max nr of flavours for alphas
P6_MSTJ(41) = 2       !  (D = 2) type of branchings allowed in shower. =1 only P6_MSTJ(48) = 1       !   1=max angle def. in PARJ(85)
!  from Professor tune 28/3/2010
P6_MSTJ(11)  = 5      !  Bowler
P6_PARJ(47)  = 0.8    ! fixed parameter
P6_PARJ(1) = 0.073    ! P(qq)/P(q) (D=0.1)
P6_PARJ(2) = 0.20     ! P(s)/P(u) (D=0.3)
P6_PARJ(3) = 0.94     ! P(us)/P(ud)/P(s)/P(d) (D=0.4)
P6_PARJ(4) = 0.032    ! (1/3)P(ud_1)/P(ud_0) (D=0.050)
P6_PARJ(11) = 0.31    ! P(S=1)d,u (D=0.5)
P6_PARJ(12) = 0.4     ! P(S=1)s   (D=0.6)
P6_PARJ(13) = 0.54    ! P(S=1)c,b (D=0.75)
P6_PARJ(21) = 0.325   ! sigma_q (D=0.360)
P6_PARJ(25) = 0.63    ! Extra eta suppression (D=1.)
P6_PARJ(26) = 0.12    ! Extra eta' suppression (D=0.4)
P6_PARJ(41) = 0.5     ! Lund symm. fragm.: a (D=0.3)
P6_PARJ(42) = 0.6     ! Lund symm. fragm.: b (D=0.58)
P6_PARJ(47) = 0.67    ! rb (D=1.0)
P6_PARJ(81) = 0.29    ! LAMBDA for alpha_s in ps (D=0.29)
P6_PARJ(82) = 1.65    ! Parton shower cut-off (D=1.0)
&End
