SMEFT NLO Reweighting Error

Asked by Binish Batool on 2019-08-06

Hello,

I have one question about the NLO Reweighting using SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO-NLO and I am producing

generate p p > t t~ a QED=1 QCD=2 NP=0 [QCD]

using reweight card (please find below) and I get the following error

INFO:
   --------------------------------------------------------------
      Summary:
      Process p p > t t~ a QED=1 QCD=2 NP=0 [QCD]
      Run at p-p collider (6500.0 + 6500.0 GeV)
      Number of events generated: 1000
      Total cross section: 2.357e+00 +- 1.5e-02 pb
   --------------------------------------------------------------

INFO: The /home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/ttgamma_NLO_EFT/Events/run_02/events.lhe.gz file has been generated.

INFO: Events generated
reweight -from_cards
INFO: split the event file in bunch of 2500 events
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 14h04 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 1m 4s ]
INFO: Collect and combine the various output file.
Command "launch auto " interrupted with error:
Exception : Some of the run failed: Please read /home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/ttgamma_NLO_EFT/Events/run_02/events.lhe.gz_0.lhe_tag_1_debug.log
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/ttgamma_NLO_EFT/run_02_tag_1_debug.log'.
Please attach this file to your report.

######## Reweight CARD #####################
change model NLO

launch

path-to-param_card
  set dim62f 22 4.0
#############################

the log information from

Log1
#########home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/ttgamma_NLO_EFT/Events/run_02/events.lhe.gz_0.lhe_tag_1_debug.log####

  File "/home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 1606, in create_standalone_directory
    mymod = __import__('%s.Source.rwgt2py' % data['paths'][1], globals(), locals(), [],-1)
ImportError: libLHAPDF.so: cannot open shared object file: No such file or directory

and

the other log information from

#####Log2
/home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/ttgamma_NLO_EFT/run_02_tag_1_debug.log

  File "/home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 2101, in do_reweight
    raise Exception, "Some of the run failed: Please read %s_%s_debug.log" % (f, self.run_tag)
Exception: Some of the run failed: Please read /home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/ttgamma_NLO_EFT/Events/run_02/events.lhe.gz_0.lhe_tag_1_debug.log
Value of current Options:

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

I have checked and tested if I switched back the NLO to LO in reweight card, it can successfully produce the event with reweighted id. So why I am not able to implement the reweighted event at NLO ?

Cheers
Binish

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2019-08-06
Last reply:
2019-08-16

Looks like you have an issue with lhapdf library (which one are you using?)
LO re-weighting does not required the lhapdf library while the NLO one needs it.

I guess that you found the reason of the issue with your lhapdf library due to your other question on launchpad.

Cheers,

Olivier

Binish Batool (binishbatool) said : #2

Yes I fixed the issue as well but after fixing this bug I could not see NLO reweighting either. So what happened MadGraph detected the reweight card but it could not generate the weighted events. Is this problem known to you? How can I solve it? I used MadGraph 2.5.6 and SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO-NLO model for the process under study.

generate p p > t t~ a NP=0 QED=1 QCD=2 [QCD]

What is the log information now?

2.5.6? your log shows 2.6.5 (2.5.6 is very very old)

Cheers,

Olivier

Binish Batool (binishbatool) said : #4

sorry I swapped the number in version. yes it's 2.6.5

>> 03:56:32 INFO: Running Reweighting
>> 03:56:32 INFO: generating the square matrix element for reweighting
>> 03:56:32 INFO: generate p p > t t~ a QED=1 QCD=2
>> --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add
>> process p p > t t~ a pert_QCD EFT=1 QED=1 QCD=3 --no_warning=duplicate;
>> 03:56:32 INFO: Done 21.14
>> 03:56:32 INFO: generate p p > t t~ a EFT=1 QED=1 QCD=2 [ virt=QCD] ;
>> 03:56:32 INFO: Done 212.2
>> 03:56:32 INFO: Original cross-section: 2.8497059 +- 0.017886814 pb
>> 03:56:32 INFO: Computed cross-section:

"Computed cross-section: it should not be empty "

and then it finished with usual way of generation.

Can you try to run in debugging mode? This should give additional information.

just add the flag "--debug" when you run your executable 9mg5.mg5_aMC or any othe executable of MG5aMC in principle)

Cheers,

Olivier

> On 7 Aug 2019, at 22:13, Binish Batool <email address hidden> wrote:
>
> Question #682641 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/682641
>
> Binish Batool posted a new comment:
> sorry I swapped the number in version. yes it's 2.6.5
>
>
>>> 03:56:32 INFO: Running Reweighting
>>> 03:56:32 INFO: generating the square matrix element for reweighting
>>> 03:56:32 INFO: generate p p > t t~ a QED=1 QCD=2
>>> --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add
>>> process p p > t t~ a pert_QCD EFT=1 QED=1 QCD=3 --no_warning=duplicate;
>>> 03:56:32 INFO: Done 21.14
>>> 03:56:32 INFO: generate p p > t t~ a EFT=1 QED=1 QCD=2 [ virt=QCD] ;
>>> 03:56:32 INFO: Done 212.2
>>> 03:56:32 INFO: Original cross-section: 2.8497059 +- 0.017886814 pb
>>> 03:56:32 INFO: Computed cross-section:
>
>
> "Computed cross-section: it should not be empty "
>
> and then it finished with usual way of generation.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Binish Batool (binishbatool) said : #6

I produced the events in Debug mode. But I did not see any additional information than earlier.

###########################################3

INFO: Events generated
reweight -from_cards
REWEIGHT: load configuration from /home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/NLO_RW/Cards/amcatnlo_configuration.txt
REWEIGHT: load configuration from /home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/input/mg5_configuration.txt
REWEIGHT: load configuration from /home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/NLO_RW/Cards/amcatnlo_configuration.txt
REWEIGHT: Extracting the banner ...
REWEIGHT: process: p p > t t~ a QED=1 QCD=2
REWEIGHT: options: NP=0
REWEIGHT: Running Reweighting
REWEIGHT: detected model: SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO-NLO. Loading...
REWEIGHT: Restrict model SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO with file models/SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO/restrict_NLO.dat .
REWEIGHT: Change particles name to pass to MG5 convention
REWEIGHT: Restrict model /home/binish/MG5_aMC_v2.6.5/MG5_aMC_v2_6_5/models/SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO with file models/SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO/restrict_NLO.dat .
REWEIGHT: Change particles name to pass to MG5 convention
REWEIGHT: generating the square matrix element for reweighting
REWEIGHT: generate p p > t t~ a QED=1 QCD=2 NP=0 --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > t t~ a pert_QCD QED=1 QCD=3 NP=0 --no_warning=duplicate;
REWEIGHT: Done 5.453
REWEIGHT: generate p p > t t~ a QED=1 QCD=2 [ virt=QCD] NP=0 ;
REWEIGHT: Done 45.18
REWEIGHT: Original cross-section: 2.3565432 +- 0.014708088 pb
REWEIGHT: Computed cross-section:
REWEIGHT: gzipping output file: events.lhe
REWEIGHT:
decay_events -from_cards

#########################################

Again "Computer cross-section" is empty.

Binish

Binish Batool (binishbatool) said : #7

  </mgrwgt>
  <rwgt>
   <wgt id='1001'> 0.45297E+01 </wgt>
   <wgt id='1002'> 0.35215E+01 </wgt>
   <wgt id='1003'> 0.59681E+01 </wgt>
   <wgt id='1004'> 0.41672E+01 </wgt>
   <wgt id='1005'> 0.32396E+01 </wgt>
   <wgt id='1006'> 0.54904E+01 </wgt>
   <wgt id='1007'> 0.49352E+01 </wgt>
   <wgt id='1008'> 0.38368E+01 </wgt>
   <wgt id='1009'> 0.65024E+01 </wgt>
  </rwgt>

you see nothing has happened. Copying weight id's from LHE file.

Binish Batool (binishbatool) said : #8

Sorry I did not hear anything about it?

Hi,

I would need to take a deep look, and that is taking too much time for my one hour of internet per day during my hollidays.

Cheers,

Olivier

> On 13 Aug 2019, at 11:47, Binish Batool <email address hidden> wrote:
>
> Question #682641 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/682641
>
> Binish Batool posted a new comment:
> Sorry I did not hear anything about it?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Do you have any file like this one?
Events/run_01/events.lhe.gz_0.lhe_tag_1_debug.log

Thanks,

Olivier

> On 9 Aug 2019, at 18:17, Binish Batool <email address hidden> wrote:
>
> Question #682641 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/682641
>
> Binish Batool posted a new comment:
> </mgrwgt>
> <rwgt>
> <wgt id='1001'> 0.45297E+01 </wgt>
> <wgt id='1002'> 0.35215E+01 </wgt>
> <wgt id='1003'> 0.59681E+01 </wgt>
> <wgt id='1004'> 0.41672E+01 </wgt>
> <wgt id='1005'> 0.32396E+01 </wgt>
> <wgt id='1006'> 0.54904E+01 </wgt>
> <wgt id='1007'> 0.49352E+01 </wgt>
> <wgt id='1008'> 0.38368E+01 </wgt>
> <wgt id='1009'> 0.65024E+01 </wgt>
> </rwgt>
>
> you see nothing has happened. Copying weight id's from LHE file.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,

When trying to reproduce your issue, I face a bug linked to your model which as nothing to do with the re-weighting.
So I create a bug report for it:
https://bugs.launchpad.net/mg5amcnlo/+bug/1840504

Did you apply any special trick to avoid this? Or might it be the same problem but that for some reason you never realize that issue with standard run before?

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Binish Batool for more information if necessary.

To post a message you must log in.