(run){
  ### basic settings ###
  SCL:=0.5;
  J1CUT:=80.0;
  J2CUT:=60.0;
  ETAMAX:=2.8;
  RJET:=0.4;
  ######################

  EVENTS 10;
  BATCH_MODE 2;
  OUTPUT 2;
  EVT_FILE_PATH=Event/V_HTSCL_mstw08_jJ1CUT_jJ2CUT_eETAMAX/

  SHERPA_LDADD SherpaLHOLE njet
  # For grid generation only use dummy virtuals
  USE_DUMMY_VIRTUAL=1;
}(run)

(beam){
  BEAM_1 = 2212; BEAM_ENERGY_1 = 3500;
  BEAM_2 = 2212; BEAM_ENERGY_2 = 3500;
}(beam)

(isr){
  PDF_LIBRARY MSTW08Sherpa;
  PDF_SET mstw2008nlo;
}(isr)

(processes){
  Process 93 93 -> 93 93;
  ## add extra jets e.g. Process 93 93 -> 93 93 93 for 3 jets
  NLO_QCD_Mode Fixed_Order;
  NLO_QCD_Part V;
  ## For Grid Generation
  Loop_Generator Internal;
  ## For Event generation
  #Loop_Generator LHOLE;
  Order_EW 0;
  End process;
}(processes)

## top loops are not included the current version of NJet
## to ensure the top quark is not included in the running
## of alpha_s we set the mass to a large value
(model){
  MASS[6] = 1e10
}(model)

(selector){
  ## asymmetric cuts : set the second parameter to the number
  ## of final state jets corresponding to the process above
  FastjetFinder   antikt  2  J2CUT  0.0  RJET  0.75  ETAMAX
  FastjetFinder   antikt  1  J1CUT  0.0  RJET  0.75  ETAMAX
}(selector)

(me){
  ME_SIGNAL_GENERATOR Comix LHOLE;
  EVENT_GENERATION_MODE Weighted;
  SCALES VAR{SCL*SCL*H_T2}
  RESULT_DIRECTORY=Result/V_HTSCL_mstw08_jJ1CUT_jJ2CUT_eETAMAX/
}(me)

(integration){
  INTEGRATOR 6;
  ERROR 0.3;
  PSI_ITMIN=10000;
  Finish_Optimization=On;
}(integration)