Time remaining information

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

Hi Sherpas,

Trying out your nice new release I noticed the following (It is rather cosmetic, but I thought I'd report it anyway):

The time remaining information on the one or two last lines before completion of the cross section calculation seems to be nonsense. See log snippets below.

This is when generating QCD multi-jet processes.

Cheers,

Karl-Johan Grahn

===
Process_Group::CalculateTotalXSec(): Calculate xs for '2_2__j__j__j__j' (Comix)
Starting the calculation. Lean back and enjoy ... .
4.01755e+09 pb +- ( 5.47945e+08 pb = 13.6388 % ) 5000 ( 17466 -> 28.6 % )
full optimization:  ( 2s(2s) elapsed / 1m 53s(1m 51s) left )   
3.92822e+09 pb +- ( 4.06256e+08 pb = 10.342 % ) 10000 ( 27014 -> 52.3 % )
full optimization:  ( 4s(3s) elapsed / 1m 45s(1m 45s) left )   
...
3.35124e+09 pb +- ( 6.72926e+07 pb = 2.00799 % ) 540000 ( 583891 -> 98.2 % )
integration time:  ( 3m 19s(3m 18s) elapsed / 1s(1s) left )   
3.38557e+09 pb +- ( 6.73724e+07 pb = 1.98999 % ) 560000 ( 604251 -> 98.2 % )
integration time:  ( 3m 27s(3m 26s) elapsed / 1158050441d 7h 12s(1158050441d 7h 13s) left )   
2_2__j__j__j__j : 3.38557e+09 pb +- ( 6.73724e+07 pb = 1.98999 % )  exp. eff: 0.0272938 %
Process_Group::CalculateTotalXSec(): Calculate xs for '2_3__j__j__j__j__j' (Comix)
Starting the calculation. Lean back and enjoy ... .
1.21434e+08 pb +- ( 2.16819e+07 pb = 17.8549 % ) 5000 ( 18228 -> 27.4 % )
full optimization:  ( 7s(7s) elapsed / 6m 47s(6m 45s) left )   
1.90025e+08 pb +- ( 2.89626e+07 pb = 15.2415 % ) 10000 ( 28639 -> 48 % )
full optimization:  ( 13s(13s) elapsed / 6m 16s(6m 15s) left )   
...
1.99037e+08 pb +- ( 6.04239e+06 pb = 3.03581 % ) 1960000 ( 2045361 -> 97.3 % )
integration time:  ( 44m 16s(44m 12s) elapsed / 1m 14s(1m 13s) left )   
1.98593e+08 pb +- ( 5.99923e+06 pb = 3.02087 % ) 1980000 ( 2065907 -> 97.3 % )
integration time:  ( 44m 43s(44m 38s) elapsed / 40s(41s) left )   
1.98984e+08 pb +- ( 5.97083e+06 pb = 3.00066 % ) 2000000 ( 2086440 -> 97.3 % )
integration time:  ( 45m 10s(45m 5s) elapsed / 1158050441d 7h 12s(1158050441d 7h 13s) left )   
1.98865e+08 pb +- ( 5.93208e+06 pb = 2.98297 % ) 2020000 ( 2106978 -> 97.3 % )
integration time:  ( 45m 39s(45m 35s) elapsed / 1158050441d 6h 59m 33s(1158050441d 6h 59m 32s) left )   
===

[Frank Siegert]

Hi Karl-Johan,

thanks for the report. I can't actually reproduce it, but I can guess where it comes from. Can you try the fairly simple patch which I'm attaching? As you say it's just a cosmetic problem, so we won't provide this as an official patch that has to go into GENSER's installation.

Cheers, Frank

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

Hi Frank,

Thanks for the patch! I have tested it, but unfortunately the problem persists. Below is the contents of Run.dat for reference.

Naturally I understand if you don't give this highest priority!

Karl-Johan

(run){
  #EVENTS = 100000
  EVENTS = 10000
  SHERPA_OUTPUT = output.sherpa
  #FILE_SIZE = 10000
  FILE_SIZE = 1000000

  RANDOM_SEED = 1 1
}(run)

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

(processes){
  Process 93 93 -> 93 93 93{4};
  Order_EW 0;
  Max_N_Quarks 4;

  CKKW sqr(15/E_CMS);

  Enhance_Factor 15 {3}
  Enhance_Factor 150 {4};
  Enhance_Factor 1000 {5};
  Enhance_Factor 1000 {6};

  Enhance_Function exp(18.2*(1-exp(-1.7e-2*PPerp(p[2])))) {2};
  Enhance_Function exp(18.2*(1-exp(-1.7e-2*max(PPerp(p[2]),PPerp(p[3]),PPerp(p[4]))))) {3};
  Enhance_Function exp(18.2*(1-exp(-1.7e-2*max(PPerp(p[2]),PPerp(p[3]),PPerp(p[4]),PPerp(p[5]))))) {4};
  Enhance_Function exp(18.2*(1-exp(-1.7e-2*max(PPerp(p[2]),PPerp(p[3]),PPerp(p[4]),PPerp(p[5]),PPerp(p[6]))))) {5};
  Enhance_Function exp(18.2*(1-exp(-1.7e-2*max(PPerp(p[2]),PPerp(p[3]),PPerp(p[4]),PPerp(p[5]),PPerp(p[6]),PPerp(p[7]))))) {6};

  Integration_Error 0.02 {2};
  Integration_Error 0.03 {3};
  Integration_Error 0.04 {4};
  Integration_Error 0.05 {5};
  Integration_Error 0.06 {6};

  Print_Graphs Graphs;

  End process;
}(processes)

(selector){
  NJetFinder  2  10.0  0.0  0.4  -1
}(selector)

(me){
  ME_SIGNAL_GENERATOR = Comix
  RESULT_DIRECTORY = Results
  EVENT_GENERATION_MODE = Weighted
  EVENT_NORM = One
}(me)

(mi){
  MI_HANDLER = Amisic  # None or Amisic
}(mi)

#(isr){
#  PDF_LIBRARY = CTEQ6Sherpa
#  PDF_SET = cteq6l
#}(isr)

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

Apologies, had the path set up wrong so that I was re-running the same binary. Will test again.

Karl-Johan

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

I have now tried again with the right patched binary and the problem is indeed solved.

Thanks again!

Karl-Johan

[Frank Siegert]

Ok, thanks for the feedback. The fix is also committed to our SVN trunk, so it will definitely be in the next version.

[Stefan Hoeche]

Milestone old deleted