*  Namelist to control input data
*
&InFiles
  ! Number of intput files
    NInputFiles = 15

  ! Input files:

     InputFileNames =
     'datafiles/hera/h1zeusCombined/inclusiveDis/1506.06042/HERA1+2_NCep_920-thexp.dat',
     'datafiles/hera/h1zeusCombined/inclusiveDis/1506.06042/HERA1+2_NCep_820-thexp.dat',
     'datafiles/hera/h1zeusCombined/inclusiveDis/1506.06042/HERA1+2_NCep_575-thexp.dat',
     'datafiles/hera/h1zeusCombined/inclusiveDis/1506.06042/HERA1+2_NCep_460-thexp.dat',
     'datafiles/hera/h1zeusCombined/inclusiveDis/1506.06042/HERA1+2_NCem-thexp.dat',
     'datafiles/hera/h1zeusCombined/inclusiveDis/1506.06042/HERA1+2_CCep-thexp.dat',
     'datafiles/hera/h1zeusCombined/inclusiveDis/1506.06042/HERA1+2_CCem-thexp.dat',

     'datafiles/lhc/atlas/zptfull8tev/zpt_y00_bwi.dat',
     'datafiles/lhc/atlas/zptfull8tev/zpt_y01_bwi.dat',
     'datafiles/lhc/atlas/zptfull8tev/zpt_y02_bwi.dat',
     'datafiles/lhc/atlas/zptfull8tev/zpt_y03_bwi.dat',
     'datafiles/lhc/atlas/zptfull8tev/zpt_y04_bwi.dat',
     'datafiles/lhc/atlas/zptfull8tev/zpt_y05_bwi.dat',
     'datafiles/lhc/atlas/zptfull8tev/zpt_y06_bwi.dat',     
     'datafiles/lhc/atlas/zptfull8tev/zpt_y07_bwi.dat',
&End

&InCorr
  ! Number of correlation (statistical, systematical or full) files
    NCorrFiles = 36

  ! Correlation files:
    CorrFileNames(1) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y0.corr',
    CorrFileNames(2) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y1.corr',
    CorrFileNames(3) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y2.corr',
    CorrFileNames(4) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y3.corr',
    CorrFileNames(5) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y4.corr',
    CorrFileNames(6) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y5.corr',
    CorrFileNames(7) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y6.corr',
    CorrFileNames(8) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y0_y7.corr',
    CorrFileNames(9) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y1_y1.corr',
    CorrFileNames(10) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y1_y2.corr',
    CorrFileNames(11) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y1_y3.corr',
    CorrFileNames(12) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y1_y4.corr',
    CorrFileNames(13) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y1_y5.corr',
    CorrFileNames(14) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y1_y6.corr',
    CorrFileNames(15) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y1_y7.corr',
    CorrFileNames(16) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y2_y2.corr',
    CorrFileNames(17) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y2_y3.corr',
    CorrFileNames(18) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y2_y4.corr',
    CorrFileNames(19) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y2_y5.corr',
    CorrFileNames(20) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y2_y6.corr',
    CorrFileNames(21) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y2_y7.corr',
    CorrFileNames(22) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y3_y3.corr',
    CorrFileNames(23) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y3_y4.corr',
    CorrFileNames(24) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y3_y5.corr',
    CorrFileNames(25) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y3_y6.corr',
    CorrFileNames(26) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y3_y7.corr',
    CorrFileNames(27) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y4_y4.corr',
    CorrFileNames(28) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y4_y5.corr',
    CorrFileNames(29) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y4_y6.corr',
    CorrFileNames(30) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y4_y7.corr',
    CorrFileNames(31) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y5_y5.corr',
    CorrFileNames(32) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y5_y6.corr',
    CorrFileNames(33) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y5_y7.corr',
    CorrFileNames(34) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y6_y6.corr',
    CorrFileNames(35) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y6_y7.corr',
    CorrFileNames(36) = 'datafiles/lhc/atlas/zptfull8tev/corr/zpt_y7_y7.corr',
&End

&ReduceSyst
    ! even with tolerance =0 the following flag may speed up calculations
  do_reduce = .false.  ! turn-on to simplify/speedup chi2 calculation.
    ! tolerance = 0.0 for exact calculation, > 0.0 for improved speed.
  tolerance = 0.0
    ! depending on blas library, toggling the following flag may improve chi2 computation speed:
  useBlas = .false.
&End

&CovarToNuisance
  ! Global switch for using nuisance param representation for covariance mat.
  LConvertCovToNui = .true.

   ! Tolerance -- zero means exact transformation
  Tolerance = 0.0

   ! The following lines allow to adjust error scaling properties (default: :M)
  DataName     = 'CMS electon Asymmetry rapidity', 'CMS W muon asymmetry', 'D0 W asymmetry 2013'
  DataSystType = ':A', ':A', ':A'
&End

*
* Main steering cards
*
&xFitter 
 ! XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
 !
 ! Chi2 definition. Following options are supported:
 !  
 ! -- Bias corrections for uncertainties --
 ! 'StatScale'    :  'Poisson',  'NoRescale' ( see also 'ExtraSystRescale' below )
 ! 'UncorSysScale':  'Poisson',  'Linear',  'NoRescale'
 ! 'CorSysScale'  :  'Linear',   'NoRescale'
 ! 
 ! -- Treatment of systematics in chi2 ---
 ! 'UncorChi2Type':  'Diagonal'  
 ! 'CorChi2Type'  :  'Hessian', 'Matrix', 'Offset'
 !
 ! -- Extra corrections ---
 !   are given as comma separated list for Chi2ExtraParam, they are off by default.
 !  'PoissonCorr'            : extra log correction accounting for changing uncertainties 
 !  'FirstIterationRescale' : re-scale uncertainties at the first iteration only 
 !  'ExtraSystRescale'      : additional re-scaling of stat. uncertainty to account for syst. shifts.

   CHI2SettingsName = 'StatScale', 'UncorSysScale', 'CorSysScale', 'UncorChi2Type', 'CorChi2Type'
   Chi2Settings     = 'Poisson'  , 'Linear',        'Linear'     , 'Diagonal'     , 'Hessian'
   Chi2ExtraParam = 'PoissonCorr'



 ! Debug flag
  LDEBUG     = False


 ! Quadratic approximation for asymmetric uncertainties
 ! AsymErrorsIterations = 10
&End


*
* Output steering cards
*
&Output 
  ! -- Q2 values at which the pdfs & errors are done (up to 20)
  Q2VAL = 1.9, 3.0, 4.0, 5., 10., 100., 6464, 8317 
!  Q2VAL = 1.9, 4., 10., 100., 6464, 8317 

  ! How many x points to write (standard = 101)
  OUTNX = 101

  ! x-range of output (standard = 1E-4 1.0)
  OUTXRANGE = 1E-4, 0.9999

  ! Write out LHAPDF5 output
  ! WriteLHAPDF5 = false
&End



*
* Process dependent cuts
*
&Cuts

  !--------------------- NC ep  --------------------------

  ! Rule #1: Q2 cuts
   ProcessName(1)     = 'NC e+-p'
   Variable(1)        = 'Q2'
   CutValueMin(1)     = 10
   CutValueMax(1)     = 1000000.0

  ! Rule #2: x cuts
   ProcessName(2)     = 'NC e+-p'
   Variable(2)        = 'x'
   CutValueMin(2)     = 0.000001 
   CutValueMax(2)     = 1.0

  !---------------------  CC ep  ------------------

   ProcessName(3)     = 'CC e+-p'
   Variable(3)        = 'Q2'
   CutValueMin(3)     = 10
   CutValueMax(3)     = 1000000.0

   ProcessName(4)     = 'CC e+-p'
   Variable(4)        = 'x'
   CutValueMin(4)     = 0.000001 
   CutValueMax(4)     = 1.0

&End

