Total inclusive jet cross section

[Generic User (don't modify these fields)]

Dear Sherpa Team,

I run the Sherpa for inclusive jet production @7TeV using the following Run.dat

(fragmentation){

  FRAGMENTATION = Ahadic
  DECAYMODEL    = Hadrons

}(fragmentation)

(me){
  ME_SIGNAL_GENERATOR   = Comix
  EVENT_GENERATION_MODE = Weighted
}(me)

(beam){
  BEAM_1          = 2212
  BEAM_ENERGY_1   = 3500.
  BEAM_2          = 2212
  BEAM_ENERGY_2   = 3500.
}(beam)
(mi){
  MI_HANDLER    = Amisic
  SCALE_MIN     = 2.0745
}(mi)

(processes){
Process 93 93 -> 93 93 93{2};
   Order_EW 0; CKKW sqr(20/E_CMS);
   Max_N_Quarks 4;

   Integration_Error 0.02 {2};
   Enhance_Function pow(PPerp2(p[2]),2) {2};

   Integration_Error 0.03 {3};
   Enhance_Function pow(max(PPerp2(p[2]),PPerp2(p[3]),PPerp2(p[4])),2) {3};
   Enhance_Factor 4 {3};

   Integration_Error 0.04 {4};
   Enhance_Function pow(max(PPerp2(p[2]),PPerp2(p[3]),PPerp2(p[4]),PPerp2(p[5])),2) {4};
   Enhance_Factor 16 {4};

End process;
}(processes);

(selector){
  NJetFinder 2 25.0 0 1
}(selector)

I get slightly different results for total cross section, if I run it with MI_HANDLER=Amisic and MI_HANDLER=None options. The difference is ~2.2%, but the accuracy of calculation is ~0.7%. So, it look like a systematic effect, rather than statistical fluctuation.

So, my question is whether such a difference is expected or there is a bug in my steering?

I have attached the logs from the two runs.

Thank you in advance, Pavel

[Generic User (don't modify these fields)]

log files

[Frank Siegert]

I'm trying whether I can reproduce this with Sherpa trunk over night.