Unable to feed lhe file to pythia8

Asked by PARTHA KUMAR PAUL

I am trying to feed pythia8 one lhe file which I have created from calchep. The process is DM DM -> H H. The DM is a singlet scalar. I am writing the following code in pythia8,

#include"Pythia8/Pythia.h"
#include<iostream>
#include<sstream>
#include<cmath>
#include<fstream>
#include<numeric>

#include "Pythia8Plugins/FastJet3.h"
#include "fastjet/PseudoJet.hh"
#include "fastjet/ClusterSequence.hh"

using namespace Pythia8;
using namespace std;

using namespace fastjet;

int main()
{

int nEvent = 100000; //number of events we want to generate

Pythia pythia;

Event& event = pythia.event;
pythia.readString("Beams:frameType=4");
pythia.readString("Beams:LHEF=sstohh.lhe");
pythia.readString("LesHouches:setLifetime=2");
pythia.readString("SLHA:minMassSM=0.");

pythia.readString("PartonLevel:FSR=off");
pythia.readString("PartonLevel:ISR=off");
pythia.readString("PartonLevel:MPI=off");

pythia.readString("HadronLevel:Hadronize=on");
pythia.readString("HadronLevel:Decay=on");
pythia.readString("SigmaProcess:renormScale2 = 4");
pythia.readString("SigmaProcess:factorScale2 = 4");

//initialise pythia
pythia.init();

I am getting the following error,

PYTHIA Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off
 PYTHIA Info from SLHAinterface::initSLHA: importing MASS entries for id = {50800000}
 PYTHIA Warning in SLHAinterface::initSLHA: adding off shell DECAY mode 25 -> 24 -24
 PYTHIA Warning in SLHAinterface::initSLHA: adding off shell DECAY mode 25 -> 23 23
 PYTHIA Info from SLHAinterface::initSLHA: importing DECAY tables for id = {25}
 PYTHIA Error in Pythia::init: could not set up PDF for beam A
 PYTHIA Abort from Pythia::init: PDF initialization failed
 PYTHIA Abort from Pythia::next: not properly initialized so cannot generate events

*------- PYTHIA Event and Cross Section Statistics -------------------------------------------------------------*
 | |
 | Subprocess Code | Number of events | sigma +- delta |
 | | Tried Selected Accepted | (estimated) (mb) |
 | | | |
 |-----------------------------------------------------------------------------------------------------------------|
 | | | |
 | | | |
 | sum | 0 0 0 | 0.000e+00 0.000e+00 |
 | |
 *------- End PYTHIA Event and Cross Section Statistics ----------------------------------------------------------*

 *------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------*
 | |
 | times message |
 | |
 | 1 Abort from Pythia::init: PDF initialization failed |
 | 10000 Abort from Pythia::next: not properly initialized so cannot generate events |
 | 1 Error in Pythia::init: could not set up PDF for beam A |
 | 1 Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off |
 | 1 Info from SLHAinterface::initSLHA: importing DECAY tables |
 | 1 Info from SLHAinterface::initSLHA: importing MASS entries |
 | 2 Warning in SLHAinterface::initSLHA: adding off shell DECAY mode |
 | |
 *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------*

Please let me know what am I missing.

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
PoolsAPK (poolsapk123) said (last edit ):
#1

It seems that you are encountering an error while trying to feed an LHE file to Pythia8. According to my search results, the error message “could not set up PDF for beam A” usually occurs when Pythia8 is unable to initialize the parton distribution function (PDF) for the beam. This could be due to a number of reasons, such as an incorrect path to the PDF data files or an outdated version of Pythia8.
To resolve this issue, you can try the following steps:
Check the path to the PDF data files: Make sure that the path to the PDF data files is correct and that Pythia8 has read access to these files.
Update Pythia8: Make sure that you are using the latest version of Pythia8. You can download the latest version from the official Pythia website.
Check your code: Make sure that your code is correctly configured and that you are passing the correct arguments to Pythia8.
I hope this helps. Let me know if there’s anything else I can do for you.

Revision history for this message
PARTHA KUMAR PAUL (parthakumarpaul) said :
#2

Since I am using a BSM particle, a singlet scalar which PDF do I have to use?
In the same code when I am feeding the pp->h process, it is working fine.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Should not you ask such question on the gitlab of pythia?
Sorry but here we are not expert in that code.

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask PARTHA KUMAR PAUL for more information if necessary.

To post a message you must log in.