&CASCADE_input
NrEvents = -1                              ! Nr of events to process
!NrEvents = 5                              ! Nr of events to process
!NrEvents = 1000                            ! Nr of events to process
Process_Id = -1                            ! Read LHE file
!Hadronisation = 1                         ! Hadronisation on (=1)
Hadronisation = 0                          ! Hadronisation on (=1)
pythia6 = 1                                ! Hadronisation with PYTHIA6 (=0 PYTHIA8)
SpaceShower = 1                            ! Space-like Parton Shower
SpaceShowerOrderAlphas=2                   ! Order alphas in Space Shower
IPSForDis = 1                              ! switch needed for initial state PS for DIS with LHE file
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 = 2                        ! 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
PartonDensity = 102200                     ! use TMDlib: PB-TMDNLO-set2
! PartonDensity = 102100                   ! use TMDlib: PB-TMDNLO-set1
! TMDDensityPath= './share'                ! Path to TMD density 
!                                            for internal files
lheInput = 'KaTie-DIS-qpm-offshell-q2gt100.lhe.gz'             ! LHE input file
!lheInput = 'KaTie-DIS-qpm-offshell-q2gt3.5.lhe.gz'             ! LHE input file
! lheInput = 'KaTie-PBset1-DIS-qpm-offshell-q2gt100.lhe'
lheHasOnShellPartons = 0                   ! = 0 LHE file has 
!                                              off-shell parton configuration
lheReweightTMD = 0                         ! Reweight with new TMD 
!                                            given in PartonDensity
lheScale = 2                               ! Scale defintion for TMD
!                                            0: use scalup
!                                            1: use shat
!                                            2: use 1/2 Sum mt of final parton/particles
!                                            3: use shat for Born and 1/2 Sum mt of final parton(particle)
!                                            4: use shat for Born and max pt of most forward/backward parton(particle)
!ScaleFactorMatchingScale = 1.              ! Scale factor for matching scale when including TMDs
lheNBornpart = 2                           ! Nr of hard partons (particles) (Born process)
! lheWeightId = 0                          ! use weight Id = ... 
!                                            as weight for LHE file
Uncertainty_TMD = 0                        ! calculate and store uncertainty TMD pdfs
&End 

&CASCADE_PS
IPS_Start_Scale = 1 	! =0(default): use LHEscale (=1 use shat) as starting scale for initial stata PS
!HepMCOUTPUT = TRUE
!HepMCFormat = 2
!HepMCFile = 'output.hepmc'
!RivetOUT = TRUE
!RivetRun = 'somerun'
!RivetFile = 'output.yoda'
!RivetAnalyses='MC_XS','MC_PDFS','MC_DIS_Check'
&End 


&PYTHIA6_input
P6_Itune = 370          ! Retune of Perugia 2011 w CTEQ6L1          (Oct 2012)
P6_MSTJ(45) = 4         ! Nr of flavors in final state shower: g->qqbar 
P6_MSTJ(41) = 1         !  (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) =  1       ! = 0 : alpha_s is fixed, =1 first order; =2 2nd order;
P6_MSTJ(48) =   1       ! (D=0), 0=no max. angle, 1=max angle def. in PARJ(85)
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
&End
