External LHE file not working with Pythia

Hi Sarah,

I never try that.
The problem is that you don't have the information in the lhe file that MG use in order to run pythia.
You need to add the run_card (with the correct tag) in the banner in order to have the translation of some MG parameter
correctly passed to pythia (like if we have to run matching,…)

Cheers,

Olivier

On Dec 28, 2013, at 1:35 PM, Sarah Malik <email address hidden> wrote:

> New question #241369 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/madgraph5/+question/241369
>
> Hi,
>
> I'm trying to run an external LHE file produced with POWHEG through the pythia+delphes framework in Madgraph (v 5.1.5.15). I copy the LHE file to a run_01 directory and rename it unweighted_events.lhe and then gzip it.
> Then I run pythia on it but get the following error:
> At line 148 of file ME2pythia.f
> Fortran runtime error: Bad integer for item 1 in list input
> 1
>
> I think its not able to read the LHE file. When I run over an LHE file produced by Madgraph then everything works fine.
> The LHE file from POWHEG that I'm trying to run over is below. Any suggestions on what the problem might be and how I can modify the LHE file so its readable?
>
> Thanks
> Sarah
>
> -----------------------------------------
> <LesHouchesEvents version="1.0">
> <!--
> file generated with POWHEG-BOX version 1.0
> Input file powheg.input contained:
> # minimum suggested values for ncall and itmx:
> # (ncall1,itmx1,ncall2,itmx2) = (250K,5,1M,5)
>
>
> decaymode 1 ! 1 for S-mediated fermionic DM
> ! 2 for PS-mediated fermionic DM
> ! -1, -2 for full theory (mediator mass and width needed)
>
> DMmass 20 ! DM-candidate mass
> DMLambda 500 ! UV scale
>
> #DMphimass 500
> #DMphiwidth 100
>
> testplots 1 ! (default 0, do not) do NLO and PWHG distributions
> bornonly 0 ! (default 0) if 1 do Born only
>
> runningscale 3 ! (default 3) 0, 1, 2, 3 correspond to 2mX (fixed), ptj, XXbar invmass, Ht/2
> ! Details in Born_phsp.f
>
> renscfact 1 ! (default 1d0) ren scale factor: muren = muref * renscfact
> facscfact 1 ! (default 1d0) fac scale factor: mufact = muref * facscfact
>
> mass_low 20 ! minimum value for XXbar invariant mass
> mass_high 8000 ! maximum value for XXbar invariant mass
>
> bornktmin 1 ! (default 0d0) generation cut. Minimum kt in underlying Born
> ###
> bornsuppfact 80 ! (default 0d0) mass param for Born suppression factor. If < 0 suppfact = 1
> ###
> withnegweights 1 ! (default 0) if 1, use negative weights
> smartsig 1
>
> numevts 5 ! number of events to be generated
> ih1 1 ! hadron 1 (1 for protons, -1 for antiprotons)
> ih2 1 ! hadron 2 (1 for protons, -1 for antiprotons)
> ebeam1 4000d0 ! energy of beam 1
> ebeam2 4000d0 ! energy of beam 2
>
> ! To be set only if using internal (mlm) pdfs
> ndns1 131 ! pdf set for hadron 1 (mlm numbering)
> ndns2 131 ! pdf set for hadron 2 (mlm numbering)
>
> ! To be set only if using LHA pdfs
> ! 21000 MSTW2008 (LO central)
> ! 21100 MSTW2008 (NLO central)
>
> lhans1 21100 ! pdf set for hadron 1 (LHA numbering)
> lhans2 21100 ! pdf set for hadron 2 (LHA numbering)
>
> ! Parameters to allow or not the use of stored data
> use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
> use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <>
>
> ncall1 50000 ! number of calls for initializing the integration grid
> itmx1 4 ! number of iterations for initializing the integration grid
> ncall2 50000 ! number of calls for computing the integral and finding upper bounding envelope
> itmx2 4 ! number of iterations for computing the integral and finding upper bouding envelope
> foldcsi 1 ! number of folds on csi integration
> foldy 1 ! number of folds on y integration
> foldphi 1 ! number of folds on phi integration
> nubound 100000 ! number of calls to set up the upper bounding norms for radiation
>
> ! OPTIONAL PARAMETERS
> xupbound 2d0 ! increase upper bound for radiation generation by this factor
>
> #iseed 5437 ! Start the random number generator with seed iseed
> #rand1 0 ! skipping rand2*100000000+rand1 numbers.
> #rand2 0 ! (see RM48 short writeup in CERNLIB)
>
> #manyseeds 1 ! Used to perform multiple runs with different random
> ! seeds in the same directory.
> ! If set to 1, the program asks for an integer j.
> ! The file pwgseeds.dat at line j is read, and the
> ! integer at line j is used to initialize the random
> ! sequence for the generation of the events.
> ! The event file is called pwgevents-'j'.lhe
>
> End of powheg.input content
> PDF package: lha
> Random number generator initialized with: 54217137 18503643 0
> -->
> <init>
> 2212 2212 4.00000E+03 4.00000E+03 -1 -1 -1 -1 -4 1
> -1.00000E+00 -1.00000E+00 1.00000E+00 10001
> </init>
> <event>
> 7 10001 3.42641E-05 3.94731E+01 -1.00000E+00 1.48520E-01
> 4 -1 0 0 511 0 0.000000000E+00 0.000000000E+00 6.502557531E+01 6.502557531E+01 0.000000000E+00 0.00000E+00 9.000E+00
> -4 -1 0 0 0 502 0.000000000E+00 0.000000000E+00 -1.395013273E+03 1.395013273E+03 0.000000000E+00 0.00000E+00 9.000E+00
> 23 2 1 2 0 0 2.589726567E+01 5.611445066E+01 -6.206042143E+02 7.229823010E+02 3.656970287E+02 0.00000E+00 9.000E+00
> 19 1 3 3 0 0 -1.105067427E+02 1.421821233E+02 -4.604244326E+02 4.947910210E+02 2.000000000E+01 0.00000E+00 9.000E+00
> -19 1 3 3 0 0 1.364040084E+02 -8.606767262E+01 -1.601797817E+02 2.281912799E+02 2.000000000E+01 0.00000E+00 9.000E+00
> 21 1 1 2 501 502 1.176247515E+01 -6.794116046E+01 -6.893515132E+02 6.927913581E+02 7.629394531E-06 0.00000E+00 9.000E+00
> 21 1 1 2 511 501 -3.765974083E+01 1.182670980E+01 -2.003197011E+01 4.426518915E+01 2.384185791E-07 0.00000E+00 9.000E+00
> #rwgt 1 27 4.458279237580892E-007 54217137 18503663 0
> </event>
> <event>
> 7 10001 1.95542E-05 9.79019E+00 -1.00000E+00 2.03118E-01
> -5 -1 0 0 0 502 0.000000000E+00 0.000000000E+00 1.320324351E+02 1.320324351E+02 0.000000000E+00 0.00000E+00 9.000E+00
> 21 -1 0 0 502 511 0.000000000E+00 0.000000000E+00 -1.437605840E+03 1.437605840E+03 0.000000000E+00 0.00000E+00 9.000E+00
> 23 2 1 2 0 0 1.190509490E+02 1.913476912E+01 -7.998704366E+01 2.469110107E+02 2.000696184E+02 0.00000E+00 9.000E+00
> 19 1 3 3 0 0 1.221357406E+02 3.435915059E+00 -1.368507279E+02 1.845455672E+02 2.000000000E+01 0.00000E+00 9.000E+00
> -19 1 3 3 0 0 -3.084791616E+00 1.569885406E+01 5.686368424E+01 6.236544351E+01 2.000000000E+01 0.00000E+00 9.000E+00
> -5 1 1 2 0 501 -1.171095870E+02 -2.873054508E+01 -2.632058127E+01 1.234215240E+02 0.000000000E+00 0.00000E+00 9.000E+00
> 21 1 1 2 501 511 -1.941361983E+00 9.595775956E+00 -1.199265780E+03 1.199305740E+03 2.157918644E-05 0.00000E+00 9.000E+00
> #rwgt 1 1 1.310511802431645E-006 54217137 18504554 0
> </event>
> <event>
> 7 10001 1.40949E-05 1.71134E+02 -1.00000E+00 1.16189E-01
> 5 -1 0 0 501 0 0.000000000E+00 0.000000000E+00 1.193437829E+03 1.193437829E+03 0.000000000E+00 0.00000E+00 9.000E+00
> 2 -1 0 0 502 0 0.000000000E+00 0.000000000E+00 -1.059155188E+03 1.059155188E+03 0.000000000E+00 0.00000E+00 9.000E+00
> 23 2 1 2 0 0 3.649319228E+02 -2.616708000E+02 9.057325966E+02 1.292992199E+03 8.061205701E+02 0.00000E+00 9.000E+00
> 19 1 3 3 0 0 1.869818869E+02 -5.506289544E+02 5.998478490E+02 8.356864923E+02 2.000000000E+01 0.00000E+00 9.000E+00
> -19 1 3 3 0 0 1.779500359E+02 2.889581544E+02 3.058847476E+02 4.573057064E+02 2.000000000E+01 0.00000E+00 9.000E+00
> 5 1 1 2 502 0 -3.801460611E+02 9.121474866E+01 -6.087516405E+02 7.234706061E+02 0.000000000E+00 0.00000E+00 9.000E+00
> 2 1 1 2 501 0 1.521413825E+01 1.704560514E+02 -1.626983151E+02 2.361302123E+02 0.000000000E+00 0.00000E+00 9.000E+00
> #rwgt 1 30 8.916683405805559E-006 54217137 18513050 0
> </event>
> <event>
> 7 10001 2.01168E-05 2.43846E+01 -1.00000E+00 1.63610E-01
> 21 -1 0 0 501 511 0.000000000E+00 0.000000000E+00 4.817909493E+02 4.817909493E+02 0.000000000E+00 0.00000E+00 9.000E+00
> 21 -1 0 0 502 501 0.000000000E+00 0.000000000E+00 -2.073755117E+02 2.073755117E+02 0.000000000E+00 0.00000E+00 9.000E+00
> 23 2 1 2 0 0 5.775404542E+01 -8.357506263E+01 1.536440670E+02 2.868107193E+02 2.198489675E+02 0.00000E+00 9.000E+00
> 19 1 3 3 0 0 1.638091468E+01 -1.387525077E+02 1.616069443E+02 2.145632709E+02 2.000000000E+01 0.00000E+00 9.000E+00
> -19 1 3 3 0 0 4.137313074E+01 5.517744509E+01 -7.962877241E+00 7.224744845E+01 2.000000000E+01 0.00000E+00 9.000E+00
> 5 1 1 2 502 0 -7.898850597E+01 9.556292377E+01 -1.127390907E+02 1.675755324E+02 0.000000000E+00 0.00000E+00 9.000E+00
> -5 1 1 2 0 511 2.123446054E+01 -1.198786115E+01 2.335104613E+02 2.347802093E+02 0.000000000E+00 0.00000E+00 9.000E+00
> #rwgt 1 30 6.080127828600106E-006 54217137 18515932 0
> </event>
> <event>
> 7 10001 1.93103E-05 9.12593E+01 -1.00000E+00 1.28102E-01
> 21 -1 0 0 501 511 0.000000000E+00 0.000000000E+00 3.082500010E+02 3.082500010E+02 0.000000000E+00 0.00000E+00 9.000E+00
> -5 -1 0 0 0 502 0.000000000E+00 0.000000000E+00 -6.980161432E+02 6.980161432E+02 0.000000000E+00 0.00000E+00 9.000E+00
> 23 2 1 2 0 0 -2.243618221E+01 -7.413103914E+01 -6.098931790E+02 6.488592217E+02 2.074844732E+02 0.00000E+00 9.000E+00
> 19 1 3 3 0 0 -1.248250686E+00 -9.827786529E+01 -1.153010892E+02 1.528215891E+02 2.000000000E+01 0.00000E+00 9.000E+00
> -19 1 3 3 0 0 -2.118793152E+01 2.414682615E+01 -4.945920898E+02 4.960376326E+02 2.000000000E+01 0.00000E+00 9.000E+00
> 21 1 1 2 501 502 8.023825240E+01 1.447510435E+02 2.072809853E+02 2.652479002E+02 1.618438983E-05 0.00000E+00 9.000E+00
> -5 1 1 2 0 511 -5.780207019E+01 -7.062000432E+01 1.284605154E+01 9.215902218E+01 0.000000000E+00 0.00000E+00 9.000E+00
> #rwgt 1 29 2.349376421846550E-006 54217137 18520911 0
> </event>
> </LesHouchesEvents>
> #Random number generator exit values: 54217137 18522497 0
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

Thanks very much for your response, it fixed that particular problem I was using.
However, now I seem to be having issues with LHAPDF. I get the following error when running on Pythia with a PDF set other than CTEQ6l (and some others). I had a look through previous questions regarding this issue and I think this may be due to Pythia having its own version of LHAPDF and its using the PDFsets that are in there (MSTW08 isn't one of them and this is what I'm trying to run with). I believe one of the answers suggested in the past was to link pythia with the correct LHAPDF, but any ideas how I do that? Do I have to modify the makefile in pythia-pgs/src/makefile?

    MSTP(98) changed from 0 to 1
 ==== PYTHIA WILL USE LHAPDF ====
 WRONG LHAPDF set number = 21100 given! STOP EXE!

Thanks!
Sarah

Hi,

> I believe one of the answers suggested in the past was to link pythia with the correct LHAPDF, but any ideas how I do that? Do I have to modify the makefile in pythia-pgs/src/makefile?

Sorry, I have no idea how to do that.

Cheers,

Olivier

On Jan 8, 2014, at 1:51 PM, Sarah Malik <email address hidden> wrote:

> Question #241369 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/madgraph5/+question/241369
>
> Sarah Malik posted a new comment:
> Hi Olivier,
>
> Thanks very much for your response, it fixed that particular problem I was using.
> However, now I seem to be having issues with LHAPDF. I get the following error when running on Pythia with a PDF set other than CTEQ6l (and some others). I had a look through previous questions regarding this issue and I think this may be due to Pythia having its own version of LHAPDF and its using the PDFsets that are in there (MSTW08 isn't one of them and this is what I'm trying to run with). I believe one of the answers suggested in the past was to link pythia with the correct LHAPDF, but any ideas how I do that? Do I have to modify the makefile in pythia-pgs/src/makefile?
>
> MSTP(98) changed from 0 to 1
> ==== PYTHIA WILL USE LHAPDF ====
> WRONG LHAPDF set number = 21100 given! STOP EXE!
>
> Thanks!
> Sarah
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

Ok, one of the threads I found which I think refers to the same issue I'm having is:
https://bugs.launchpad.net/madgraph5/+bug/1169027
I'll try the suggestion you made (https://bugs.launchpad.net/madgraph5/+bug/1169027/comments/4) and see if it works.

Thanks
Sarah