Computation fails when using [LOonly=QCD]

Asked by Hesham El Faham

Hi,

I am trying to generate the following process:

generate p p > t t t~ t~ NP^2==2 QCD<=4 QED<=4 YT<=4

in smeftatnlo (I adapted the model such that the Yukawa is explicit, hence YT).

Running this process as shown above is perfectly fine and it generates ~1000 independent diagrams. However, running the exact same process with [LOonly=QCD] or [LOonly=all] results in the computation failing. I checked the end of the log file in the sub-process directories, and this is what it reads:

"Too many gridpacks to keep track off"

--> Is there a way around that?

Another thing I noted that with this syntax is that more diagrams are generated (~2000 independent) compared to the normal generation,

--> is that understood? (I was not able to compare the cross-sections since the computation never succeeded)

Even though I want to study this process at LO, the reason I use the [LOonly=QCD] syntax is because I want to do a Fixed Order HwU Analysis (and separate smeft orders, etc.) and this is the only way, as far as I know, by which I can do that.

Your help is much appreciated.

Best,
Hesham

Question information

Language:
English Edit question
Status:
Open
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said :
#1

Hi Hesham,
what version of the code are you using?
I have tried the SM-version of this process very recently, and everything was running fine...
Cheers,

Marco

Revision history for this message
Hesham El Faham (helfaham) said :
#2

Hi Marco,

I used 3.0.4 and the latest one 3_3_1, and I get the same issue with both.

Yeah, I think the error appears when using smeftatnlo.

Best,
Hesham

Revision history for this message
marco zaro (marco-zaro) said :
#3

Ok. Do you confirm that with the SM it works?

Can you share the model?

Marco

Revision history for this message
Hesham El Faham (helfaham) said :
#4

Hi marco,

So I removed the YT explicit coupling to make sure things are simple, so in what follows I am just using the default models.

I did a run with sm-no_b_mass and with the command:
generate p p > t t t~ t~ QCD<=4 QED<=4 [LOonly=all] and it worked.

I then used the default smeftatnlo at leading order and did this:
generate p p > t t t~ t~ NP^2==2 QCD<=4 QED<=4 [LOonly=all]. It didn't work and failed the same way as I had mentioned early on.

In both runs I have used a restriction card that takes the 4 fermion operators only, and in both I used the latest version of the code.

Best,
Hesham

Can you help with this problem?

Provide an answer of your own, or ask Hesham El Faham for more information if necessary.

To post a message you must log in.