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