Calchep3.6.23 and LHAPDF6.15: Problem with PDFs

Asked by Fernando Marroquim

I am using the Ubuntu 12.04, calchep3.6.23 , LHAPDF6.15 and MSTWL2008.
I got the following message:
LHAPDF 6.1.5 loading /home/fmarroquim/LHAPDF615/share/LHAPDF/MSTW2008lo90cl/MSTW2008lo90cl_0000.dat
MSTW2008lo90cl PDF set, member #0, version 3; LHAPDF ID = 21050
ERROR in strfun
LHAPDF 6.1.5 loading /home/fmarroquim/LHAPDF615/share/LHAPDF/MSTW2008lo90cl/MSTW2008lo90cl_0000.dat
MSTW2008lo90cl PDF set, member #0, version 3; LHAPDF ID = 21050

In the interactive numeric mode

The option
S.F.1: MSTW2008lo90cl:0:0
does not allow to change the last figure to 1 or -1 it is always set to 0
It allows to change the PDF set.

The option S.F.2: it is working fine.

Similar problem is occuring in the batch mode appearing the message: ERROR in strfun
when I put
pdf1: LHA:MSTW2008lo90cl:0:1
pdf2: LHA:MSTW2008lo90cl:0:1

Question information

Language:
English Edit question
Status:
Solved
For:
CalcHEP Edit question
Assignee:
No assignee Edit question
Solved by:
Fernando Marroquim
Solved:
Last query:
Last reply:
Revision history for this message
Alexander Pukhov (pukhov) said :
#1

At first screen of symbolic mode you should see
     Last correction March 14. ( line 22, file c_source/symb/lbl.c)
 If not, please, download the last improved version. We indeed had such bug.
If your version is the last one , please, let me know and I'll check it again.
Best
   Alexander Pukhov

Revision history for this message
Fernando Marroquim (fmarroquim) said :
#2

Dear Alexander
Thank you for your answer.

Below is the line 22 from file lbl.c
goto_xy(15,Y+1); print("Version 3.6.23: Last correction March 14, 2015");

Regards
Fernando

Revision history for this message
Alexander Pukhov (pukhov) said :
#3

Thank you for confirmation.
Unfortunately I don't see this bug on my CentOS Linux.
I'll re-read code carefully to see what could be a reason.

Revision history for this message
Alexander Pukhov (pukhov) said :
#4

If in results/session.dat you have

LHA:MSTW2008lo90cl:0:0

then indeed proton/antiProton switch will not work because the last number has to be 1 or -1.
The proton/antiProton switch changes its sign.

I guess that everything will work properly when you remove files results/session.dat and results/aux/session.dat and start new session.
   Indeed there is bug in the code which allows to set this zero. It appears when one press the Esc key in menu on LHA sets.
To remove this bug just delete the line 250 in
   c_source/num/strfun/sf_lha.c

sgn[i1]=0;

and recompile CaclHEP by "make". n_calchep should be recompiled as well.

For batch mode I still don't see a problem.

Thanks a lot for bug information.
    Alexander Pukhov

Revision history for this message
Fernando Marroquim (fmarroquim) said :
#5

Dear Alexander,
When setting
sgn[i1]=0; or sgn[i1]=1; in file c_source/num/strfun/sf_lha.c
it is working fine in interactive mode.

I still have problem in batch mode.
I am using the calchep example batch_file changing only the lines:
S.F.1: MSTW2008nlo90cl:0:1
S.F.1: MSTW2008nlo90cl:0:1

In batch mode the session. dat appears the following first lines:
#Subprocess 1 ( u, D -> W+, b, B )
#Session_number 1
#Initial_state inP1=4.000000E+03 inP2=4.000000E+03
 Polarizations= { 0.000000E+00 0.000000E+00 }
  StrFun1="LHA:MSTW2008nlo90cl:0:1"
  StrFun2="LHA:MSTW2008nlo90cl:0:1"

Then there are many
ERROR in strfun
messages
Thanks
Fernando

Revision history for this message
Alexander Pukhov (pukhov) said :
#6

I have launched it several times and in my computer it works fine.
Have you

export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:$<path to your lhapdf>

instruction in calchep_batch ?

Best
   Alexander Pukhov

Revision history for this message
Fernando Marroquim (fmarroquim) said :
#7

CORRECTIONof my last message
Dear Alexander,
When deleting line 250 sgn[i1]=0;
or
setting sgn[i1]=1; in file c_source/num/strfun/sf_lha.c
it is working fine in interactive mode.

I still have problem in batch mode.
I am using the calchep example batch_file changing only the lines:
S.F.1: MSTW2008nlo90cl:0:1
S.F.1: MSTW2008nlo90cl:0:1

In batch mode the session. dat appears the following first lines:
#Subprocess 1 ( u, D -> W+, b, B )
#Session_number 1
#Initial_state inP1=4.000000E+03 inP2=4.000000E+03
 Polarizations= { 0.000000E+00 0.000000E+00 }
  StrFun1="LHA:MSTW2008nlo90cl:0:1"
  StrFun2="LHA:MSTW2008nlo90cl:0:1"

Then there are many
ERROR in strfun
messages
Thanks
Fernan

Revision history for this message
Fernando Marroquim (fmarroquim) said :
#8

Dear Alexander,
It is working fine in both ways:
Interactive and batch modes.
Thank you very much for your support.
Fernando