CalcHEP

Not able to use MSTW2008LO PDF in CalcHep

Asked by Harinder

Hello,

I had a batch file which was used 1 year back with older version of CalcHep where I could successfully used MSTW2008LO PDF as

#########################################
# PDF's #
#########################################
pdf1 : LHA:MSTW2008lo68cl.LHgrid:0:1
pdf2 : LHA:MSTW2008lo68cl.LHgrid:0:1

But now I am trying to run same batch file using CalcHep 3.6.18 and LHAPDF-6.1.4 but Its giving an error
I am getting error: (Just when it started calculating sigma)

&&&&&&&&&&&&&&&&&&&&&&&&&&
.......RemTV........
Kr: {(1*2^0) N^0 (N-1)^0 (N+1)^0
[0/t33 (0,2,1)][1/t33 (0,1,2)]}ok
Kr: {(1*2^0) N^0 (N-1)^0 (N+1)^0
[1/t33 (0,1,1)]}ok
Kr: {(1*2^0) N^1 (N-1)^0 (N+1)^0
}ok
ERROR in strfun
ERROR in strfun
Error in calculation of cross section for u,u->u,u
 for run MWX200.
Please have a look at cs.o and cs.out in Processes/m4/unitary/p1/MWX200.
Killing all jobs and quitting.

&&&&&&&&&&&&&&&&&&&&

If I change my PDF from above to

pdf1 : cteq6l(proton)
pdf2 : cteq6l(proton)

it runs perfectly fine.

Can anyone tell me what could be the mistake? Any suggestion is appreciated.

Question information

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

Which version LHAPDF is using by you?
Best
    Alexander Pukhov

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

I see , you use version 6.
As I remember version 6 has another interface.
So, nothing is changed in CalcHEP. It still supports version 5 LHAPDF.
I'll try to implement version 6.

Revision history for this message
Harinder (harinder-singh-bawa) said : 		#3

Hello,

Thanks for the info. Yes, I am using LHAPDF6, In the meantime, I can try to use LHAPDF 5

Harinder

Revision history for this message
Harinder (harinder-singh-bawa) said : 		#4

Hello, I tried using 5.9.1 LHAPDF with CalcHEP 3.6.18 instead of version 6 as I was getting errors above:

But here too in LHAPDF version 5. I am getting somewhat different error:

File $Home/Work/lhapdf-5.9.1/../local/share/lhapdf/PDFsets/MSTW2008lo68cl.LHgrid cannot be opened !
 If you have not already done so:
 Use the bin/lhapdf-getdata script to download the file.
 Set the environmental variable LHAPATH to specify the directory if not as default (above).

I can see
 21000 1 0 0 MSTW2008lo68cl.LHgrid 0 1.0 1000000000. 1.E-06 1. MSTW2008-lo(central)
in
$Home/Work/lhapdf-5.9.1/../local/share/lhapdf/PDFsets.index

Why it doesnt recognise then?

\

FYI
In my Batch script: I am using
#########################################
# PDF's #
#########################################
pdf1 : LHA:MSTW2008lo68cl.LHgrid:0:1
pdf2 : LHA:MSTW2008lo68cl.LHgrid:0:1

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

I have tested 5.951 and confirms that it works fine.
I recommend to try GUI interface before launching batch.

Record in share/lhapdf/PDFsets.index means nothing. One has to go to LHAPDF WEB site and download corresponding file.
I guess file
 $Home/Work/lhapdf-5.9.1/../local/share/lhapdf/PDFsets/MSTW2008lo68cl.LHgrid
indeed is absent on your computer.

Best
   Alexander Pukhov

Revision history for this message
Best Alexander Belyaev (alexander.belyaev) said : 		#6

Dear Harinder,

I confirm also, that CalcHEP works with the LHAPDF5 versions,
( so please download the respective PDFset),
and I also confirm that CalcHEP will not work at the moment with the
LHAPDF6 versions since the interface to LHAPDF6 has been qualitatively changed
and has been written on C++ versus fortran version for LHAPDF5.

So, our advice is to use LHAPDF5 version with which CalcHEP works fine as we can see on our side.
Please, confirm this.

Alexander Belyaev

Revision history for this message
Harinder (harinder-singh-bawa) said : 		#7

Hello,

Thanks! Yes, Its indeed works for LHAPDF 5 after downloading PDF. Thanks for the answers

Harinder

Revision history for this message
Harinder (harinder-singh-bawa) said : 		#8

Thanks Alexander Belyaev, that solved my question.