CRITICAL: The following initial/final state ((1, 21), (1, 25)) can not be found in the new model/process.

Asked by SHUBHAM YADAV

Hi ,
I am getting the following error and i am also sharing the run summary. Please suggest me what should i do now??

CRITICAL: The following initial/final state ((1, 21), (1, 25)) can not be found in the new model/process. If you want to set the weights of such events to zero use "change allow_missing_finalstate False"
Command "generate_events run_13" interrupted with error:

=== Results Summary for run: run_13 tag: tag_1 ===

     Cross-section : 13.88 +- 0.0793 pb
     Nb of events : 100

INFO: No version of lhapdf. Can not run systematics computation
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: Running Reweighting
change mode LO # Define type of Reweighting. For LO sample this command
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 1.254
WARNING: REWEIGHTING: original card and new card are identical.
INFO: starts to compute weight for events with the following modification to the param_card:
INFO:
INFO: Event nb 0 current time: 22h31
INFO: Event nb 10 current time: 22h31
INFO: Event nb 20 current time: 22h31
INFO: Event nb 30 current time: 22h31
INFO: Event nb 40 current time: 22h31
INFO: Event nb 50 current time: 22h31
INFO: Event nb 60 0.022s
INFO: Event nb 70 0.025s
INFO: Event nb 80 0.027s
INFO: Event nb 90 0.03s
INFO: All event done (nb_event: 100) 0.033s
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/Events/run_13/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 13.8823 pb (indicative error: 0.079297 pb)
change model loop_sm
change process g g > H g [sqrvirt=QCD]
change process g u > H u [sqrvirt=QCD] --add
change process g u~ > H u~ [sqrvirt=QCD] --add
change process g d~ > H d~ [sqrvirt=QCD] --add
change process g c > H c [sqrvirt=QCD] --add
change process g c~ > H c~ [sqrvirt=QCD] --add
change process g s > H s [sqrvirt=QCD] --add
change process g s~ > H s~ [sqrvirt=QCD] --add
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 1.198
INFO: detected model: loop_sm. Loading...
INFO: generating the square matrix element for reweighting (second model and/or processes)
INFO: generate g g > H g --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g g > H g pert_QCD --no_warning=duplicate;add process g u > H u --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g u > H u pert_QCD --no_warning=duplicate;add process g u~ > H u~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g u~ > H u~ pert_QCD --no_warning=duplicate;add process g d~ > H d~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g d~ > H d~ pert_QCD --no_warning=duplicate;add process g c > H c --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g c > H c pert_QCD --no_warning=duplicate;add process g c~ > H c~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g c~ > H c~ pert_QCD --no_warning=duplicate;add process g s > H s --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g s > H s pert_QCD --no_warning=duplicate;add process g s~ > H s~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g s~ > H s~ pert_QCD --no_warning=duplicate;
INFO: RETRY with generate g g > H g [sqrvirt=QCD] ;add process g u > H u [sqrvirt=QCD] ;add process g u~ > H u~ [sqrvirt=QCD] ;add process g d~ > H d~ [sqrvirt=QCD] ;add process g c > H c [sqrvirt=QCD] ;add process g c~ > H c~ [sqrvirt=QCD] ;add process g s > H s [sqrvirt=QCD] ;add process g s~ > H s~ [sqrvirt=QCD] ;
INFO: Done 4.695
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: change model loop_sm
KEEP: change process g g > H g [sqrvirt=QCD]
KEEP: change process g u > H u [sqrvirt=QCD]
KEEP: change process g u~ > H u~ [sqrvirt=QCD]
KEEP: change process g d~ > H d~ [sqrvirt=QCD]
KEEP: change process g c > H c [sqrvirt=QCD]
KEEP: change process g c~ > H c~ [sqrvirt=QCD]
KEEP: change process g s > H s [sqrvirt=QCD]
KEEP: change process g s~ > H s~ [sqrvirt=QCD]
INFO: Event nb 0 current time: 22h32
INFO: Event nb 10 0.057s
INFO: Event nb 20 0.081s
INFO: Event nb 30 0.1s
CRITICAL: The following initial/final state ((1, 21), (1, 25)) can not be found in the new model/process. If you want to set the weights of such events to zero use "change allow_missing_finalstate False"
Command "generate_events run_13" interrupted with error:
Exception :
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj/run_13_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
gzipping output file: unweighted_events.lhe
INFO: Done
INFO:
INFO: Original cross-section: 13.88234 +- 0.07929705 pb (cross-section from sum of weights: 13.88233500000001)
INFO: Computed cross-section:
INFO: rwgt_1 : 13.88233500000001 +- 0.07929705 pb
quit
INFO:
MG5_aMC>

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
Olivier Mattelaer (olivier-mattelaer) said :
#1

The problem is that even if you have the line:
change process g u > H u [sqrvirt=QCD] --add

you do not have the following one:
change process u g > H u [sqrvirt=QCD] --add

Cheers,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#2

What is the solution??

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

To add (all) the symmetric process for the initial state.

Cheers,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#4

please provide me a command

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

So in your reweight_card add

change process u g > H u [sqrvirt=QCD] --add
change process u~ g > H u~ [sqrvirt=QCD] --add
change process d~ g > H d~ [sqrvirt=QCD] --add
change process c g > H c [sqrvirt=QCD] --add
change process c~ g > H c~ [sqrvirt=QCD] --add
change process s g > H s [sqrvirt=QCD] --add
change process s~ g > H s~ [sqrvirt=QCD] --add

Cheers,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#6
Revision history for this message
SHUBHAM YADAV (sy388405) said :
#7

Thankyou for your fast reply.
as ypou suggested i have modified and added that in reweight card but i'm still getting the following error

Command "generate_events run_33" interrupted with error:
TypeError : int() argument must be a string, a bytes-like object or a real number, not 'NoneType'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj/run_33_decayed_1_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run

 === Results Summary for run: run_33 tag: tag_1 ===

     Cross-section : 13.95 +- 0.07945 pb
     Nb of events : 10

INFO: No version of lhapdf. Can not run systematics computation
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: Running Reweighting
change mode LO # Define type of Reweighting. For LO sample this command
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 1.282
WARNING: REWEIGHTING: original card and new card are identical.
INFO: starts to compute weight for events with the following modification to the param_card:
INFO:
INFO: Event nb 0 current time: 02h36
INFO: All event done (nb_event: 10) current time: 02h36
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/Events/run_33/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 13.9502 pb (indicative error: 0.0794473 pb)
change model loop_sm
change process g g > H g [sqrvirt=QCD]
change process g u > H u [sqrvirt=QCD] --add
change process g u~ > H u~ [sqrvirt=QCD] --add
change process g d~ > H d~ [sqrvirt=QCD] --add
change process g c > H c [sqrvirt=QCD] --add
change process g c~ > H c~ [sqrvirt=QCD] --add
change process g s > H s [sqrvirt=QCD] --add
change process g s~ > H s~ [sqrvirt=QCD] --add
change process u g > H u [sqrvirt=QCD] --add
change process u~ g > H u~ [sqrvirt=QCD] --add
change process d~ g > H d~ [sqrvirt=QCD] --add
change process c g > H c [sqrvirt=QCD] --add
change process c~ g > H c~ [sqrvirt=QCD] --add
change process s g > H s [sqrvirt=QCD] --add
change process s~ g > H s~ [sqrvirt=QCD] --add
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 1.171
INFO: detected model: loop_sm. Loading...
INFO: generating the square matrix element for reweighting (second model and/or processes)
INFO: generate g g > H g --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g g > H g pert_QCD --no_warning=duplicate;add process g u > H u --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g u > H u pert_QCD --no_warning=duplicate;add process g u~ > H u~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g u~ > H u~ pert_QCD --no_warning=duplicate;add process g d~ > H d~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g d~ > H d~ pert_QCD --no_warning=duplicate;add process g c > H c --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g c > H c pert_QCD --no_warning=duplicate;add process g c~ > H c~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g c~ > H c~ pert_QCD --no_warning=duplicate;add process g s > H s --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g s > H s pert_QCD --no_warning=duplicate;add process g s~ > H s~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g s~ > H s~ pert_QCD --no_warning=duplicate;add process u g > H u --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process u g > H u pert_QCD --no_warning=duplicate;add process u~ g > H u~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process u~ g > H u~ pert_QCD --no_warning=duplicate;add process d~ g > H d~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process d~ g > H d~ pert_QCD --no_warning=duplicate;add process c g > H c --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process c g > H c pert_QCD --no_warning=duplicate;add process c~ g > H c~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process c~ g > H c~ pert_QCD --no_warning=duplicate;add process s g > H s --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process s g > H s pert_QCD --no_warning=duplicate;add process s~ g > H s~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process s~ g > H s~ pert_QCD --no_warning=duplicate;
INFO: RETRY with generate g g > H g [sqrvirt=QCD] ;add process g u > H u [sqrvirt=QCD] ;add process g u~ > H u~ [sqrvirt=QCD] ;add process g d~ > H d~ [sqrvirt=QCD] ;add process g c > H c [sqrvirt=QCD] ;add process g c~ > H c~ [sqrvirt=QCD] ;add process g s > H s [sqrvirt=QCD] ;add process g s~ > H s~ [sqrvirt=QCD] ;add process u g > H u [sqrvirt=QCD] ;add process u~ g > H u~ [sqrvirt=QCD] ;add process d~ g > H d~ [sqrvirt=QCD] ;add process c g > H c [sqrvirt=QCD] ;add process c~ g > H c~ [sqrvirt=QCD] ;add process s g > H s [sqrvirt=QCD] ;add process s~ g > H s~ [sqrvirt=QCD] ;
INFO: Done 6.222
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: change model loop_sm
KEEP: change process g g > H g [sqrvirt=QCD]
KEEP: change process g u > H u [sqrvirt=QCD]
KEEP: change process g u~ > H u~ [sqrvirt=QCD]
KEEP: change process g d~ > H d~ [sqrvirt=QCD]
KEEP: change process g c > H c [sqrvirt=QCD]
KEEP: change process g c~ > H c~ [sqrvirt=QCD]
KEEP: change process g s > H s [sqrvirt=QCD]
KEEP: change process g s~ > H s~ [sqrvirt=QCD]
KEEP: change process u g > H u [sqrvirt=QCD]
KEEP: change process u~ g > H u~ [sqrvirt=QCD]
KEEP: change process d~ g > H d~ [sqrvirt=QCD]
KEEP: change process c g > H c [sqrvirt=QCD]
KEEP: change process c~ g > H c~ [sqrvirt=QCD]
KEEP: change process s g > H s [sqrvirt=QCD]
KEEP: change process s~ g > H s~ [sqrvirt=QCD]
INFO: Event nb 0 current time: 02h36
INFO: All event done (nb_event: 10) 0.028s
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/Events/run_33/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 4.51908 pb (indicative error: 2.90525 pb)
/home/shubham/Cards/param_card_loop_sm.dat
Command "/home/shubham/Cards/param_card_loop_sm.dat" not recognized, please try again
INFO: Original cross-section: 13.95017 +- 0.07944727 pb (cross-section from sum of weights: 13.950166)
INFO: Computed cross-section:
INFO: rwgt_1 : 13.950166 +- 0.07944727 pb
INFO: rwgt_2 : 4.519084040351715 +- 2.9052451133005235 pb
INFO: quit rwgt
decay_events -from_cards
INFO: Running MadSpin
INFO: This functionality allows for the decay of resonances
INFO: in a .lhe file, keeping track of the spin correlation effets.
INFO: BE AWARE OF THE CURRENT LIMITATIONS:
INFO: (1) Only a succession of 2 body decay are currently allowed
************************************************************
* *
* W E L C O M E to M A D S P I N *
* *
************************************************************
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: detected model: heft. Loading...
No MadAnalysis5 version number could be read from the path supplied '/home/shubham/work/MG5_aMC_v2_9_20/HEPTools/madanalysis5/madanalysis5'.
The specified version of MadAnalysis5 will not be active in your session.
set ninja to /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/lib
set collier to /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/lib
lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

Set group_subprocesses to Auto
Note that you need to regenerate all processes
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
set spinmode onshell
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
launch
INFO: Will use seed 730526214
INFO: generating the production square matrix element for onshell
INFO: generating the full matrix element squared (with decay)
INFO: generate matrix element for decay only (1 - > N).
INFO: output standalone /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/madspin_me
INFO: Done 2.866
INFO: Compiling code
INFO: Compiling code
INFO: Estimating the maximum weight
INFO: *****************************
INFO: Probing the first 75 events with 400 phase space points
INFO: Event 1/75 : 0.062552s
INFO: Event 6/75 : 0.279485s
INFO: Event 11/75 : 0.495521s
INFO: Event 16/75 : 0.712336s
INFO: Event 21/75 : 0.928856s
INFO: Event 26/75 : 1.145483s
INFO: Event 31/75 : 1.362548s
INFO: Event 36/75 : 1.579722s
INFO: Event 41/75 : 1.797284s
INFO: Event 46/75 : 2.014822s
INFO: Event 51/75 : 2.232498s
INFO: Event 56/75 : 2.450166s
INFO: Event 61/75 : 2.668479s
INFO: Event 66/75 : 2.886761s
INFO: Event 71/75 : 3.104699s
INFO: The decayed event file has been moved to the following location:
INFO: /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/Events/run_33_decayed_1/unweighted_events.lhe
INFO: MadSpin Done
INFO: Creating root files
Warning in <TTree::Bronch>: Using split mode on a class: TRootWeight with a custom Streamer
** Reading /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/Events/run_33_decayed_1/unweighted_events.lhe
** [################################################################] (100.00%)
** Exiting...
INFO: Running Pythia8 [arXiv:1410.3012]
Splitting .lhe event file for PY8 parallelization...
Submitting Pythia8 jobs...
Pythia8 shower jobs: 1 Idle, 0 Running, 0 Done [0 second]
Pythia8 shower jobs: 0 Idle, 0 Running, 1 Done [1 seconds]
Merging results from the split PY8 runs...
Command "generate_events run_33" interrupted with error:
TypeError : int() argument must be a string, a bytes-like object or a real number, not 'NoneType'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj/run_33_decayed_1_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
INFO: Done
INFO:
quit
INFO:
MG5_aMC>

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

So now you have the re-weighting working fine.
And face an issue within running pythia8 (which might just be because you do not have enough events?).

But here I really need to have the content of the debug file in order to be able to comment (if the error is not an internal pythia8 one, which does not seem to be the case)

A work around here is to run pythia8 manually outside of MG5aMC

Cheers,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#9

please explain this --
A work around here is to run pythia8 manually outside of MG5aMC

content of debug file---

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.9.20 2024-06-17 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_33
Traceback (most recent call last):
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1541, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/madevent_interface.py", line 2474, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/common_run_interface.py", line 7126, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/madevent_interface.py", line 2584, in run_generate_events
    self.exec_cmd('shower --no_default', postcmd=False, printcmd=False)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/madevent_interface.py", line 3927, in do_shower
    self.exec_cmd('%s %s'%(shower,' '.join(args)),
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/madevent_interface.py", line 4669, in do_pythia8
    new_djr_HwU = histograms.HwUList(djr_file,run_id=0)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/various/histograms.py", line 1659, in __init__
    self.parse_histos_from_PY8_XML_stream(stream, run_id,
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/various/histograms.py", line 1830, in parse_histos_from_PY8_XML_stream
    wgt_label[mandatory_attribute] = int(wgt_label[mandatory_attribute])
TypeError: int() argument must be a string, a bytes-like object or a real number, not 'NoneType'
                              Run Options
                              -----------
               stdout_level : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
                   run_mode : 2
              cluster_queue : None (user set)
               cluster_time : None (user set)
               cluster_size : 100
             cluster_memory : 100 (user set)
                    nb_core : 20 (user set)
          cluster_temp_path : None

                      Configuration Options
                      ---------------------
               pythia8_path : /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
                thepeg_path : None (user set)
                 hepmc_path : /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/hepmc (user set)
           madanalysis_path : None (user set)
          madanalysis5_path : /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/madanalysis5/madanalysis5 (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /home/shubham/work/MG5_aMC_v2_9_20/Delphes (user set)
        exrootanalysis_path : /home/shubham/work/MG5_aMC_v2_9_20/ExRootAnalysis (user set)
               syscalc_path : None (user set)
                     lhapdf : lhapdf-config
                 lhapdf_py2 : None
                 lhapdf_py3 : None
                    timeout : 60
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                web_browser : None
                 eps_viewer : None
                text_editor : None
           fortran_compiler : None
               cpp_compiler : None
                auto_update : 7 (user set)
               cluster_type : condor
      cluster_status_update : (600, 30)
           cluster_nb_retry : 1
         cluster_local_path : None
         cluster_retry_wait : 300
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.9.20 2024-06-17 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set group_subprocesses Auto
set ignore_six_quark_processes False
set low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model heft
generate p p > h j
output pp_hj
######################################################################
## PARAM_CARD AUTOMATICALLY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118800e+01 # MZ
   25 1.250000e+02 # MH
  9000006 1.250001e+02 # MP
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000e+00 # d : 0.0
  2 0.000000e+00 # u : 0.0
  3 0.000000e+00 # s : 0.0
  4 0.000000e+00 # c : 0.0
  11 0.000000e+00 # e- : 0.0
  12 0.000000e+00 # ve : 0.0
  13 0.000000e+00 # mu- : 0.0
  14 0.000000e+00 # vm : 0.0
  16 0.000000e+00 # vt : 0.0
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  24 8.041900e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.325070e+02 # aEWM1
    2 1.166390e-05 # Gf
    3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set)

 ###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.200000e+00 # ymb
    6 1.645000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
DECAY 9000006 6.382339e-03 # WH1
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 9000006 # h1
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0 # Particle/Antiparticle distinction (0=own anti)
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update to_full *
#*********************************************************************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#*********************************************************************
  10 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron, 4=photon from muon *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     6500.0 = ebeam1 ! beam 1 total energy in GeV
     6500.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     nn23lo1 = pdlabel ! PDF set
     230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
 False = fixed_fac_scale ! if .true. use fixed fac scale
 91.188 = scale ! fixed ren scale
 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
 1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# Phase-Space Optimization strategy (basic options)
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
                             ! 0: sum over helicity, 1: importance sampling
   2 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
                             ! 1 is old strategy (using amp square)
        ! 2 is new strategy (using only the denominator)
# To see advanced option for Phase-Space optimization: type "update psoptim"
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
 None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
 {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 20.0 = ptj ! minimum pt for the jets
 -1.0 = ptjmax ! maximum pt for the jets
 {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
 5.0 = etaj ! max rap for the jets
 {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
 {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
 {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0.0 = ptheavy ! minimum pt for at least one heavy final state
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   True = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#10

i have run for 1000 events then it's started to show the same reweight error. with this i am qlso sharing its debug file's content

CRITICAL: The following initial/final state ((1, 21), (1, 25)) can not be found in the new model/process. If you want to set the weights of such events to zero use "change allow_missing_finalstate False"
Command "generate_events run_35" interrupted with error:
Exception :

 === Results Summary for run: run_35 tag: tag_1 ===

     Cross-section : 13.93 +- 0.08011 pb
     Nb of events : 1000

INFO: No version of lhapdf. Can not run systematics computation
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: Running Reweighting
change mode LO # Define type of Reweighting. For LO sample this command
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 1.171
WARNING: REWEIGHTING: original card and new card are identical.
INFO: starts to compute weight for events with the following modification to the param_card:
INFO:
INFO: Event nb 0 current time: 03h37
INFO: Event nb 10 current time: 03h37
INFO: Event nb 20 current time: 03h37
INFO: Event nb 30 current time: 03h37
INFO: Event nb 40 0.02s
INFO: Event nb 50 0.025s
INFO: Event nb 60 0.028s
INFO: Event nb 70 0.031s
INFO: Event nb 80 0.034s
INFO: Event nb 90 0.037s
INFO: Event nb 100 0.04s
INFO: Event nb 200 0.068s
INFO: Event nb 300 0.096s
INFO: Event nb 400 0.12s
INFO: Event nb 500 0.15s
INFO: Event nb 600 0.18s
INFO: Event nb 700 0.21s
INFO: Event nb 800 0.24s
INFO: Event nb 900 0.26s
INFO: All event done (nb_event: 1000) 0.29s
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/Events/run_35/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 13.928 pb (indicative error: 0.0801058 pb)
change model loop_sm
change process g g > H g [sqrvirt=QCD]
change process g u > H u [sqrvirt=QCD] --add
change process g u~ > H u~ [sqrvirt=QCD] --add
change process g d~ > H d~ [sqrvirt=QCD] --add
change process g c > H c [sqrvirt=QCD] --add
change process g c~ > H c~ [sqrvirt=QCD] --add
change process g s > H s [sqrvirt=QCD] --add
change process g s~ > H s~ [sqrvirt=QCD] --add
change process u g > H u [sqrvirt=QCD] --add
change process u~ g > H u~ [sqrvirt=QCD] --add
change process d~ g > H d~ [sqrvirt=QCD] --add
change process c g > H c [sqrvirt=QCD] --add
change process c~ g > H c~ [sqrvirt=QCD] --add
change process s g > H s [sqrvirt=QCD] --add
change process s~ g > H s~ [sqrvirt=QCD] --add
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 1.206
INFO: detected model: loop_sm. Loading...
INFO: generating the square matrix element for reweighting (second model and/or processes)
INFO: generate g g > H g --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g g > H g pert_QCD --no_warning=duplicate;add process g u > H u --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g u > H u pert_QCD --no_warning=duplicate;add process g u~ > H u~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g u~ > H u~ pert_QCD --no_warning=duplicate;add process g d~ > H d~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g d~ > H d~ pert_QCD --no_warning=duplicate;add process g c > H c --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g c > H c pert_QCD --no_warning=duplicate;add process g c~ > H c~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g c~ > H c~ pert_QCD --no_warning=duplicate;add process g s > H s --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g s > H s pert_QCD --no_warning=duplicate;add process g s~ > H s~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process g s~ > H s~ pert_QCD --no_warning=duplicate;add process u g > H u --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process u g > H u pert_QCD --no_warning=duplicate;add process u~ g > H u~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process u~ g > H u~ pert_QCD --no_warning=duplicate;add process d~ g > H d~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process d~ g > H d~ pert_QCD --no_warning=duplicate;add process c g > H c --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process c g > H c pert_QCD --no_warning=duplicate;add process c~ g > H c~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process c~ g > H c~ pert_QCD --no_warning=duplicate;add process s g > H s --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process s g > H s pert_QCD --no_warning=duplicate;add process s~ g > H s~ --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process s~ g > H s~ pert_QCD --no_warning=duplicate;
INFO: RETRY with generate g g > H g [sqrvirt=QCD] ;add process g u > H u [sqrvirt=QCD] ;add process g u~ > H u~ [sqrvirt=QCD] ;add process g d~ > H d~ [sqrvirt=QCD] ;add process g c > H c [sqrvirt=QCD] ;add process g c~ > H c~ [sqrvirt=QCD] ;add process g s > H s [sqrvirt=QCD] ;add process g s~ > H s~ [sqrvirt=QCD] ;add process u g > H u [sqrvirt=QCD] ;add process u~ g > H u~ [sqrvirt=QCD] ;add process d~ g > H d~ [sqrvirt=QCD] ;add process c g > H c [sqrvirt=QCD] ;add process c~ g > H c~ [sqrvirt=QCD] ;add process s g > H s [sqrvirt=QCD] ;add process s~ g > H s~ [sqrvirt=QCD] ;
INFO: Done 6.1
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: change model loop_sm
KEEP: change process g g > H g [sqrvirt=QCD]
KEEP: change process g u > H u [sqrvirt=QCD]
KEEP: change process g u~ > H u~ [sqrvirt=QCD]
KEEP: change process g d~ > H d~ [sqrvirt=QCD]
KEEP: change process g c > H c [sqrvirt=QCD]
KEEP: change process g c~ > H c~ [sqrvirt=QCD]
KEEP: change process g s > H s [sqrvirt=QCD]
KEEP: change process g s~ > H s~ [sqrvirt=QCD]
KEEP: change process u g > H u [sqrvirt=QCD]
KEEP: change process u~ g > H u~ [sqrvirt=QCD]
KEEP: change process d~ g > H d~ [sqrvirt=QCD]
KEEP: change process c g > H c [sqrvirt=QCD]
KEEP: change process c~ g > H c~ [sqrvirt=QCD]
KEEP: change process s g > H s [sqrvirt=QCD]
KEEP: change process s~ g > H s~ [sqrvirt=QCD]
INFO: Event nb 0 current time: 03h38
INFO: Event nb 10 0.024s
CRITICAL: The following initial/final state ((1, 21), (1, 25)) can not be found in the new model/process. If you want to set the weights of such events to zero use "change allow_missing_finalstate False"
Command "generate_events run_35" interrupted with error:
Exception :
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj/run_35_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
gzipping output file: unweighted_events.lhe
INFO: Done
INFO:
INFO: Original cross-section: 13.92799 +- 0.08010582 pb (cross-section from sum of weights: 13.927987999999868)
INFO: Computed cross-section:
INFO: rwgt_1 : 13.927987999999868 +- 0.08010582 pb
quit
INFO:

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

Content of its debug file

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.9.20 2024-06-17 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_35
Traceback (most recent call last):
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/reweight_interface.py", line 1236, in calculate_matrix_element
    orig_order, Pdir, hel_dict = self.id_to_path_second[tag]
KeyError: ((1, 21), (1, 25))

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1541, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/madevent_interface.py", line 2474, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/common_run_interface.py", line 7126, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/madevent_interface.py", line 2572, in run_generate_events
    self.exec_cmd('reweight -from_cards', postcmd=False)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/common_run_interface.py", line 2178, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1689, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/various/misc.py", line 109, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/reweight_interface.py", line 565, in do_launch
    weight = self.calculate_weight(event)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/reweight_interface.py", line 1023, in calculate_weight
    w_new = self.calculate_matrix_element(new_event, 1)
  File "/home/shubham/work/MG5_aMC_v2_9_20/madgraph/interface/reweight_interface.py", line 1242, in calculate_matrix_element
    raise Exception
Exception
                              Run Options
                              -----------
               stdout_level : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
                   run_mode : 2
              cluster_queue : None (user set)
               cluster_time : None (user set)
               cluster_size : 100
             cluster_memory : 100 (user set)
                    nb_core : 20 (user set)
          cluster_temp_path : None

                      Configuration Options
                      ---------------------
               pythia8_path : /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
                thepeg_path : None (user set)
                 hepmc_path : /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/hepmc (user set)
           madanalysis_path : None (user set)
          madanalysis5_path : /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/madanalysis5/madanalysis5 (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /home/shubham/work/MG5_aMC_v2_9_20/Delphes (user set)
        exrootanalysis_path : /home/shubham/work/MG5_aMC_v2_9_20/ExRootAnalysis (user set)
               syscalc_path : None (user set)
                     lhapdf : lhapdf-config
                 lhapdf_py2 : None
                 lhapdf_py3 : None
                    timeout : 60
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                web_browser : None
                 eps_viewer : None
                text_editor : None
           fortran_compiler : None
               cpp_compiler : None
                auto_update : 7 (user set)
               cluster_type : condor
      cluster_status_update : (600, 30)
           cluster_nb_retry : 1
         cluster_local_path : None
         cluster_retry_wait : 300
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.9.20 2024-06-17 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set group_subprocesses Auto
set ignore_six_quark_processes False
set low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model heft
generate p p > h j
output pp_hj
######################################################################
## PARAM_CARD AUTOMATICALLY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118800e+01 # MZ
   25 1.250000e+02 # MH
  9000006 1.250001e+02 # MP
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000e+00 # d : 0.0
  2 0.000000e+00 # u : 0.0
  3 0.000000e+00 # s : 0.0
  4 0.000000e+00 # c : 0.0
  11 0.000000e+00 # e- : 0.0
  12 0.000000e+00 # ve : 0.0
  13 0.000000e+00 # mu- : 0.0
  14 0.000000e+00 # vm : 0.0
  16 0.000000e+00 # vt : 0.0
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  24 8.041900e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.325070e+02 # aEWM1
    2 1.166390e-05 # Gf
    3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set)

 ###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.200000e+00 # ymb
    6 1.645000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
DECAY 9000006 6.382339e-03 # WH1
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 9000006 # h1
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0 # Particle/Antiparticle distinction (0=own anti)
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update to_full *
#*********************************************************************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#*********************************************************************
  1000 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron, 4=photon from muon *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     6500.0 = ebeam1 ! beam 1 total energy in GeV
     6500.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     nn23lo1 = pdlabel ! PDF set
     230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
 False = fixed_fac_scale ! if .true. use fixed fac scale
 91.188 = scale ! fixed ren scale
 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
 1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# Phase-Space Optimization strategy (basic options)
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
                             ! 0: sum over helicity, 1: importance sampling
   2 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
                             ! 1 is old strategy (using amp square)
        ! 2 is new strategy (using only the denominator)
# To see advanced option for Phase-Space optimization: type "update psoptim"
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
 None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
 {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 20.0 = ptj ! minimum pt for the jets
 -1.0 = ptjmax ! maximum pt for the jets
 {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
 5.0 = etaj ! max rap for the jets
 {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
 {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
 {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0.0 = ptheavy ! minimum pt for at least one heavy final state
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   True = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#11
Revision history for this message
SHUBHAM YADAV (sy388405) said :
#12
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#13
Revision history for this message
SHUBHAM YADAV (sy388405) said :
#14

Hi, I am facing this issue now Command "generate_events run_43" interrupted with error:
KeyError : ((21, 21), (25, 21))
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj/run_43_tag_1_debug.log'.
Please attach this file to your report.

  === Results Summary for run: run_43 tag: tag_1 ===

     Cross-section : 13.88 +- 0.07881 pb
     Nb of events : 10

INFO: No version of lhapdf. Can not run systematics computation
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: Running Reweighting
change mode LO # Define type of Reweighting. For LO sample this command
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 1.187
WARNING: REWEIGHTING: original card and new card are identical.
INFO: starts to compute weight for events with the following modification to the param_card:
INFO:
INFO: Event nb 0 current time: 23h46
INFO: All event done (nb_event: 10) current time: 23h46
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj/Events/run_43/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 13.8777 pb (indicative error: 0.0788103 pb)
change keep_ordering True
launch
WARNING: REWEIGHTING: original card and new card are identical.
INFO: starts to compute weight for events with the following modification to the param_card:
INFO:
INFO: Event nb 0 current time: 23h46
Command "generate_events run_43" interrupted with error:
KeyError : ((21, 21), (25, 21))
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj/run_43_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
gzipping output file: unweighted_events.lhe
INFO: Done
INFO:
INFO: Original cross-section: 13.87767 +- 0.07881034 pb (cross-section from sum of weights: 13.877672)
INFO: Computed cross-section:
INFO: rwgt_1 : 13.877672 +- 0.07881034 pb
quit
INFO:
MG5_aMC>

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

Can you provide the various file that you are using, such that I will better understand what you are doing, and maybe even try to reproduce your issue.

Cheers,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#16

various file like what ???

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

The script/command that you use to generate the process
the run_card.dat
the param_card.dat
the reweight_card.dat
and potential other files/information, that would be useful for me to try to understand/reproduce your issue.

Cheers,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#18

import model heft

generate p p > h j

output pp pp_hj

launch

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

######################################################################
## PARAM_CARD AUTOMATICALLY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118800e+01 # MZ
   25 1.250000e+02 # MH
  9000006 1.250001e+02 # MP
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000e+00 # d : 0.0
  2 0.000000e+00 # u : 0.0
  3 0.000000e+00 # s : 0.0
  4 0.000000e+00 # c : 0.0
  11 0.000000e+00 # e- : 0.0
  12 0.000000e+00 # ve : 0.0
  13 0.000000e+00 # mu- : 0.0
  14 0.000000e+00 # vm : 0.0
  16 0.000000e+00 # vt : 0.0
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  24 8.041900e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.325070e+02 # aEWM1
    2 1.166390e-05 # Gf
    3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set)

 ###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.200000e+00 # ymb
    6 1.645000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
DECAY 9000006 6.382339e-03 # WH1
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 9000006 # h1
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0 # Particle/Antiparticle distinction (0=own anti)

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update to_full *
#*********************************************************************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#*********************************************************************
  10 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron, 4=photon from muon *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     6500.0 = ebeam1 ! beam 1 total energy in GeV
     6500.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     nn23lo1 = pdlabel ! PDF set
     230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
 False = fixed_fac_scale ! if .true. use fixed fac scale
 91.188 = scale ! fixed ren scale
 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
 1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# Phase-Space Optimization strategy (basic options)
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
                             ! 0: sum over helicity, 1: importance sampling
   2 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
                             ! 1 is old strategy (using amp square)
        ! 2 is new strategy (using only the denominator)
# To see advanced option for Phase-Space optimization: type "update psoptim"
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
 None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
 {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 20.0 = ptj ! minimum pt for the jets
 -1.0 = ptjmax ! maximum pt for the jets
 {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
 5.0 = etaj ! max rap for the jets
 {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
 {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
 {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0.0 = ptheavy ! minimum pt for at least one heavy final state
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   True = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

#*************************************************************************
# Reweight Module *
# Matrix-Element reweighting at LO/NLO *
# Mattelaer Olivier arxiv:1607.00763 *
#*************************************************************************
#
# Note:
# 1) the value of alpha_s will be used from the event so the value in
# the param_card is not taken into account.
# 2) It is (in general) dangerous/wrong to change parameters by a large
# amount, if this changes the shape of the matrix elements a lot.
# (For example, changing a particle's mass by much more than its
# width leads to very inaccurate result). In such a case, separate
# event generation runs are needed.
#
#************************************************************************
# ENTER YOUR COMMANDS BELOW.
#************************************************************************

change mode LO # Define type of Reweighting. For LO sample this command
                   # has no effect since only "LO" mode is allowed.
launch

change keep_ordering True
launch

change model loop_sm
change process g g > H g [sqrvirt=QCD]
change process g u > H u [sqrvirt=QCD] --add
change process g u~ > H u~ [sqrvirt=QCD] --add
change process g d~ > H d~ [sqrvirt=QCD] --add
change process g c > H c [sqrvirt=QCD] --add
change process g c~ > H c~ [sqrvirt=QCD] --add
change process g s > H s [sqrvirt=QCD] --add
change process g s~ > H s~ [sqrvirt=QCD] --add
change process u g > H u [sqrvirt=QCD] --add
change process u~ g > H u~ [sqrvirt=QCD] --add
change process d~ g > H d~ [sqrvirt=QCD] --add
change process c g > H c [sqrvirt=QCD] --add
change process c~ g > H c~ [sqrvirt=QCD] --add
change process s g > H s [sqrvirt=QCD] --add
change process s~ g > H s~ [sqrvirt=QCD] --add
launch
~/Cards/param_card_loop_sm.dat
# SPECIFY A PATH OR USE THE SET COMMAND LIKE THIS:
# set sminputs 1 130 # modify 1/alpha_EW

#************************************************************************
# Manual: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Reweight
#************************************************************************
#
# Example of (standard) code for the computation of two weights:
#
# launch ! tag to start the computation of the first weight
# set BLOCKNAME ID VALUE ! rule to modify the current param_card
# set BLOCKNAME ID VALUE ! rule to modify the current param_card
# launch ! start to compute a second weight
# /home/Cards/param_card_2.dat ! you can also enter a path to a valid card
#
# Note: The command to specify the parameter are AFTER the associated "launch"
#
# Possible options:
# You can enter one of the following lines to customize the reweighting
# procedure. These need to be given before the 'launch' command.
#
# change model NAME : use another model for the matrix-elements to reweight
# with. In this case you need to provide the path to a correct
# param_card for the new model; you cannot modify the original one
# with the 'set' command.
# change process DEF [--add]: change the process by which you reweight.
# The initial and final state particles of the new process should
# be exactly identical to the ones in the original process.
# change helicity False: perform the reweighting by helicity summed
# matrix-elements even if the events have been written with a
# single helicity state.
# change mode XXX: change the type of reweighting performed.
# allowed values: LO, NLO, LO+NLO
# - This command has no effect for reweighting an .lhe event file with LO accuracy.
# In that case LO mode is always used (whatever entry is set).
# - When the .lhe file reweighted is at NLO accuracy, then all modes are allowed.
# * "LO" is an approximate leading order method
# * "NLO" is the NLO accurate method
# * "LO+NLO" runs both
# - "NLO" and "LO+NLO" modes requires 'store_rwgt_info' equals True (run_card.dat)
# If the reweighting is done at generation level this parameter will
# automatically be set on True.
#************************************************************************

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#* Manual: *
#* cp3.irmp.ucl.ac.be/projects/madgraph/wiki/MadSpin *
#* *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
# set spinmode onshell # Use one of the madspin special mode
 set max_weight_ps_point 400 # number of PS to estimate the maximum for each event

# specify the decay for the final state particles
 set spinmode onshell
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
# running the actual code
launch

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
pythia8_card
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
!...Parton showering on or off
      MSTP(61)=1
      MSTP(71)=1

!...Fragmentation/hadronization on or off
      MSTJ(1)=1

!...Multiple interactions on or off
      MSTP(81)=20

!...Don't stop execution after 10 errors
      MSTU(21)=1

!...PDFset if MG set not supported by pythia-pgs package (set in lhapdf5 or higher)
      LHAID= 10041

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

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

I have tested the following reweight_card.dat

 change mode LO # Define type of Reweighting. For LO sample this command
                    # has no effect since only "LO" mode is allowed.

 #change keep_ordering True
 change model loop_sm
 change process p p > H j [sqrvirt=QCD]
 launch
 ~/Cards/param_card_loop_sm.dat

and it worked out of the box.

Now I did run without madspin first, because your madspin card does not specify a higgs decay mode into that file.
Then I tested with di photon decay and no issue.
And finally your pythia8 card is not a pythia8 card but a pythia6 card (and pythia6 is really a legacy code that should not be used anymore)

Cheers,

Olivier

So the exact script that I did run:

 import model heft
 generate p p > h j
 output
 launch
 madspin=ON
 reweight=ON
 ./rwgt.dat
 set use_syst F
 decay h > a a

with the file rwgt containing the lines mentioned above

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#20

Dear Madgraph expert,
As you have mentioned above-
 "set use_syst F" Where should i add this and what it means....

but i have installed pythia8 from Madgraph interface, and it is by default as it is...
I did'nt make any changes in it. then how it becomes card of pythia6...

As you suggest i have do the changes. I did run without Hadronization. Its working. But When I on the Hadronization. It show the following Error --

Merging results from the split PY8 runs...
Command "generate_events run_13" interrupted with error:
TypeError : int() argument must be a string, a bytes-like object or a real number, not 'NoneType'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/run_13_decayed_1_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run

  === Results Summary for run: run_13 tag: tag_1 ===

     Cross-section : 13.88 +- 0.0793 pb
     Nb of events : 1000

INFO: No version of lhapdf. Can not run systematics computation
INFO: Creating Plots for parton level
INFO: Plots for parton level generated, see /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/HTML/run_13/plots_parton.html
INFO: End Plots for parton level
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: Running Reweighting
change mode LO # Define type of Reweighting. For LO sample this command
change model loop_sm
change process p p > H j [sqrvirt=QCD]
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 0.6355
INFO: detected model: loop_sm. Loading...
INFO: generating the square matrix element for reweighting (second model and/or processes)
INFO: generate p p > H j --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > H j pert_QCD --no_warning=duplicate;
INFO: RETRY with generate p p > H j [sqrvirt=QCD] ;
INFO: Done 3.728
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: change model loop_sm
KEEP: change process p p > H j [sqrvirt=QCD]
INFO: Event nb 0 current time: 16h08
INFO: Event nb 10 current time: 16h08
INFO: Event nb 20 current time: 16h08
INFO: Event nb 30 0.026s
INFO: Event nb 40 0.033s
INFO: Event nb 50 0.042s
INFO: Event nb 60 0.053s
INFO: Event nb 70 0.062s
INFO: Event nb 80 0.071s
INFO: Event nb 90 0.08s
INFO: Event nb 100 0.09s
INFO: Event nb 200 0.18s
INFO: Event nb 300 0.28s
INFO: Event nb 400 0.37s
INFO: Event nb 500 0.46s
INFO: Event nb 600 0.55s
INFO: Event nb 700 0.64s
INFO: Event nb 800 0.72s
INFO: Event nb 900 0.81s
INFO: All event done (nb_event: 1000) 0.9s
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/Events/run_13/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 12.0403 pb (indicative error: 0.225895 pb)
/home/shubham/Cards/param_card_loop_sm.data
Command "/home/shubham/Cards/param_card_loop_sm.data" not recognized, please try again
set use_syst F
WARNING: Invalid Syntax. The command 'set' should be placed after the 'launch' one. Continuing by adding automatically 'launch'
launch
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: change model loop_sm
KEEP: change process p p > H j [sqrvirt=QCD]
INFO: Event nb 0 current time: 16h08
INFO: Event nb 10 current time: 16h08
INFO: Event nb 20 current time: 16h08
INFO: Event nb 30 0.028s
INFO: Event nb 40 0.035s
INFO: Event nb 50 0.044s
INFO: Event nb 60 0.055s
INFO: Event nb 70 0.065s
INFO: Event nb 80 0.074s
INFO: Event nb 90 0.084s
INFO: Event nb 100 0.094s
INFO: Event nb 200 0.18s
INFO: Event nb 300 0.28s
INFO: Event nb 400 0.39s
INFO: Event nb 500 0.48s
INFO: Event nb 600 0.58s
INFO: Event nb 700 0.67s
INFO: Event nb 800 0.76s
INFO: Event nb 900 0.85s
INFO: All event done (nb_event: 1000) 0.94s
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/Events/run_13/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 12.0403 pb (indicative error: 0.225895 pb)
INFO: Original cross-section: 13.88234 +- 0.07929705 pb (cross-section from sum of weights: 13.882335000000138)
INFO: Computed cross-section:
INFO: rwgt_1 : 12.04033301184774 +- 0.2258947930255331 pb
INFO: rwgt_2 : 12.04033301184774 +- 0.2258947930255331 pb
INFO: quit rwgt
decay_events -from_cards
INFO: Running MadSpin
INFO: This functionality allows for the decay of resonances
INFO: in a .lhe file, keeping track of the spin correlation effets.
INFO: BE AWARE OF THE CURRENT LIMITATIONS:
INFO: (1) Only a succession of 2 body decay are currently allowed
************************************************************
* *
* W E L C O M E to M A D S P I N *
* *
************************************************************
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: detected model: heft. Loading...
No MadAnalysis5 version number could be read from the path supplied '/home/shubham/work/MG5_aMC_v2_9_20/HEPTools/madanalysis5/madanalysis5'.
The specified version of MadAnalysis5 will not be active in your session.
set ninja to /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/lib
set collier to /home/shubham/work/MG5_aMC_v2_9_20/HEPTools/lib
lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

Set group_subprocesses to Auto
Note that you need to regenerate all processes
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a
launch
INFO: Will use seed 410470592
INFO: We need to recalculate the branching fractions for t,h,w+,w-,t~,z
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /home/shubham/work/MG5_aMC_v2_9_20/models/heft with file models/heft/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
WARNING: The LO estimate for the width of particle 25
WARNING: differs from the one in the banner by 33 percent
INFO:
INFO: decay channels for t : ( width = 1.4915 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b w+
INFO:
INFO:
INFO: decay channels for z : ( width = 2.441756 GeV )
INFO: BR d1 d2
INFO: 1.523658e-01 d~ d
INFO: 1.523658e-01 s~ s
INFO: 1.507440e-01 b~ b
INFO: 1.188141e-01 u~ u
INFO: 1.188141e-01 c~ c
INFO: 6.793744e-02 ve~ ve
INFO: 6.793744e-02 vm~ vm
INFO: 6.793744e-02 vt~ vt
INFO: 3.438713e-02 e+ e-
INFO: 3.438713e-02 mu+ mu-
INFO: 3.430972e-02 ta+ ta-
INFO:
INFO:
INFO: decay channels for w+ : ( width = 2.04793 GeV )
INFO: BR d1 d2
INFO: 3.333610e-01 d~ u
INFO: 3.333610e-01 s~ c
INFO: 1.111195e-01 e+ ve
INFO: 1.111195e-01 mu+ vm
INFO: 1.110390e-01 ta+ vt
INFO:
INFO:
INFO: decay channels for h : ( width = 0.004768364 GeV )
INFO: BR d1 d2
INFO: 9.027834e-01 b~ b
INFO: 5.426389e-02 ta+ ta-
INFO: 4.091760e-02 g g
INFO: 2.035121e-03 a a
INFO:
INFO:
INFO: decay channels for t~ : ( width = 1.4915 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b~ w-
INFO:
INFO:
INFO: decay channels for w- : ( width = 2.04793 GeV )
INFO: BR d1 d2
INFO: 3.333610e-01 d u~
INFO: 3.333610e-01 s c~
INFO: 1.111195e-01 e- ve~
INFO: 1.111195e-01 mu- vm~
INFO: 1.110390e-01 ta- vt~
INFO:
INFO: generating the production square matrix element
INFO: generate p p > h j;
INFO: Done 0.7758
INFO: generating the full matrix element squared (with decay)
INFO: generate p p > h j, h > a a QCD=99 --no_warning=duplicate;
INFO: Done 0.7781
INFO: generate matrix element for decay only (1 - > N).
INFO: output standalone_msF /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/decay_me
INFO: Done 0.4685
INFO: Compiling code
INFO: detect independant decays
INFO:
INFO: Estimating the maximum weight
INFO: *****************************
INFO: Probing the first 75 events
INFO: with 400 phase space points
INFO:
INFO: Event 1/75 : current time: 16h09
INFO: Event 6/75 : 0.059s
INFO: Event 11/75 : 0.11s
INFO: Event 16/75 : 0.16s
INFO: Event 21/75 : 0.21s
INFO: Event 26/75 : 0.26s
INFO: Event 31/75 : 0.31s
INFO: Event 36/75 : 0.36s
INFO: Event 41/75 : 0.42s
INFO: Event 46/75 : 0.48s
INFO: Event 51/75 : 0.54s
INFO: Event 56/75 : 0.6s
INFO: Event 61/75 : 0.66s
INFO: Event 66/75 : 0.71s
INFO: Event 71/75 : 0.77s
INFO:
INFO: Decaying the events...
INFO: Total number of events written: 1000/1000
INFO: Average number of trial points per production event: 1.04
INFO: Branching ratio to allowed decays: 0.00203512
INFO: Number of events with weights larger than max_weight: 0
INFO: Number of subprocesses 3
INFO: Number of failures when restoring the Monte Carlo masses: 0
INFO: The decayed event file has been moved to the following location:
INFO: /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/Events/run_13_decayed_1/unweighted_events.lhe.gz
INFO: MadSpin Done
INFO: Creating Plots for parton level
INFO: Plots for parton level generated, see /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/HTML/run_13_decayed_1/plots_parton.html
INFO: End Plots for parton level
INFO: Running Pythia8 [arXiv:1410.3012]
Splitting .lhe event file for PY8 parallelization...
Submitting Pythia8 jobs...
Pythia8 shower jobs: 1 Idle, 9 Running, 0 Done [0 second]
Pythia8 shower jobs: 0 Idle, 9 Running, 1 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 8 Running, 2 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 7 Running, 3 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 6 Running, 4 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 5 Running, 5 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 4 Running, 6 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 3 Running, 7 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 2 Running, 8 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 1 Running, 9 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 0 Running, 10 Done [1 seconds]
Merging results from the split PY8 runs...
Command "generate_events run_13" interrupted with error:
TypeError : int() argument must be a string, a bytes-like object or a real number, not 'NoneType'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/run_13_decayed_1_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
INFO: Done
INFO:
quit
INFO:
MG5_aMC>

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

So form the help command:

INFO: *** HELP MESSAGE ***
## Information on parameter use_syst from class RunCardLO
## current value: True (parameter should be of type <class 'bool'>)
## Add in the lhef file information needed for the computation of systematic uncertainty (scale variation and pdf)
## This value is considered as been set by the user

So it does not matter if it is True or False, I just set it to False to avoid issue related to lhapdf which is compiled on my machine for python3.12 and here I need to run with python3.10 making the code to crash if the systematic was on (Technically I could have kept the writting of the information and just turn the running of the systematics)

For the where, I have put that command in the main script, but this is equivalent to edit the run_card by hand.

For your question about pythia8, I do not know what you did to use a pythia6 card, nobody report a similar issue.
In the cards directory you should have a _default card for pythia8 does that file also contains pythia6 information?
For sure pythia8 is going to crash if you pass to him a pythia6 card...

In case of problem with the integration of pythia8, the best is typically to run pythia8 outside madgraph and if the problem still persist then contact pythia8 authors.

Cheers,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#22

hi, have you run with on state of hadronization?
i make run with on state of hadronization via pythia8, it shows the following error--

Command "generate_events run_20" interrupted with error:
TypeError : int() argument must be a string, a bytes-like object or a real number, not 'NoneType'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/run_20_tag_1_debug.log'.
Please attach this file to your report.

 === Results Summary for run: run_20 tag: tag_1 ===

     Cross-section : 13.97 +- 0.08095 pb
     Nb of events : 1000

INFO: No version of lhapdf. Can not run systematics computation
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
INFO: Extracting the banner ...
INFO: process: p p > h j
INFO: options:
INFO: Running Reweighting
change mode LO # Define type of Reweighting. For LO sample this command
change model loop_sm
change process p p > H j [sqrvirt=QCD]
launch
INFO: detected model: heft. Loading...
INFO: generating the square matrix element for reweighting
INFO: generate p p > h j ;
INFO: Done 0.5929
INFO: detected model: loop_sm. Loading...
INFO: generating the square matrix element for reweighting (second model and/or processes)
INFO: generate p p > H j --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > H j pert_QCD --no_warning=duplicate;
INFO: RETRY with generate p p > H j [sqrvirt=QCD] ;
INFO: Done 3.77
INFO: starts to compute weight for events with the following modification to the param_card:
INFO: change model loop_sm
KEEP: change process p p > H j [sqrvirt=QCD]
INFO: Event nb 0 current time: 23h21
INFO: Event nb 10 current time: 23h21
INFO: Event nb 20 0.024s
INFO: Event nb 30 0.03s
INFO: Event nb 40 0.04s
INFO: Event nb 50 0.049s
INFO: Event nb 60 0.06s
INFO: Event nb 70 0.07s
INFO: Event nb 80 0.08s
INFO: Event nb 90 0.089s
INFO: Event nb 100 0.1s
INFO: Event nb 200 0.19s
INFO: Event nb 300 0.28s
INFO: Event nb 400 0.36s
INFO: Event nb 500 0.45s
INFO: Event nb 600 0.54s
INFO: Event nb 700 0.63s
INFO: Event nb 800 0.71s
INFO: Event nb 900 0.81s
INFO: All event done (nb_event: 1000) 0.89s
INFO: Event /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/Events/run_20/unweighted_events.lhe have now the additional weight
INFO: new cross-section is : 12.1382 pb (indicative error: 0.231714 pb)
/home/shubham/Cards/param_card_loop_sm.data
Command "/home/shubham/Cards/param_card_loop_sm.data" not recognized, please try again
INFO: Original cross-section: 13.97449 +- 0.08094698 pb (cross-section from sum of weights: 13.974489000000144)
INFO: Computed cross-section:
INFO: rwgt_1 : 12.138221767985428 +- 0.23171421756984312 pb
INFO: quit rwgt
decay_events -from_cards
INFO: Creating root files
Warning in <TTree::Bronch>: Using split mode on a class: TRootWeight with a custom Streamer
** Reading /home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/Events/run_20/unweighted_events.lhe
** [################################################################] (100.00%)
** Exiting...
INFO: Running Pythia8 [arXiv:1410.3012]
Splitting .lhe event file for PY8 parallelization...
Submitting Pythia8 jobs...
Pythia8 shower jobs: 1 Idle, 9 Running, 0 Done [0 second]
Pythia8 shower jobs: 0 Idle, 9 Running, 1 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 8 Running, 2 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 7 Running, 3 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 6 Running, 4 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 5 Running, 5 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 4 Running, 6 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 3 Running, 7 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 2 Running, 8 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 1 Running, 9 Done [1 seconds]
Pythia8 shower jobs: 0 Idle, 0 Running, 10 Done [1 seconds]
Merging results from the split PY8 runs...
Command "generate_events run_20" interrupted with error:
TypeError : int() argument must be a string, a bytes-like object or a real number, not 'NoneType'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/run_20_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
INFO: Done
INFO:
quit
INFO:
MG5_aMC>

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

Hi,

Indeed it did not work for me this time,
I think that you were lucky since for sure it was working last week, I guess the reason is that I did install gnuplot last night for studying a bug report related to that tool, and since here the issue that I face was plot related it likely face such issue only if gnuplot is installed on the machine. Anyway this is the patch that I have include to fix the issue (and that will be part of 2.9.21 release: https://github.com/mg5amcnlo/mg5amcnlo/commit/f85c4de987a895d9c0b999025d35ddd84a3ec51a):

diff --git a/madgraph/various/histograms.py b/madgraph/various/histograms.py
index 094c1c674..593582fcd 100755
--- a/madgraph/various/histograms.py
+++ b/madgraph/various/histograms.py
@@ -1824,9 +1824,9 @@ class HwUList(histograms_PhysicsObjectList):
         # Also cast them in the proper type
         for wgt_label in all_weights:
             for mandatory_attribute in ['PDF','MUR','MUF','MERGING','ALPSFACT']:
- if mandatory_attribute not in wgt_label:
+ if mandatory_attribute not in wgt_label or wgt_label[mandatory_attribute] is None:
                     wgt_label[mandatory_attribute] = '-1'
- if mandatory_attribute=='PDF':
+ elif mandatory_attribute=='PDF':
                     wgt_label[mandatory_attribute] = int(wgt_label[mandatory_attribute])
                 elif mandatory_attribute in ['MUR','MUF','MERGING','ALPSFACT']:
                     wgt_label[mandatory_attribute] = float(wgt_label[mandatory_attribute])

With that patch, i did not had any issue anymore.

Thanks for the report,

Olivier

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#24

Hi, now what should be my next step ?

Revision history for this message
Sihyun Jeon (shjeon) said :
#25

fix the lines in "madgraph/various/histograms.py" in your madgraph repository as Olivier did above and rerun it
ones with - are deleted lines and + are added lines

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#26

Hi, I have do the suggested modifications and did rerun, now it shows the following error--

Command "generate_events run_26" interrupted with error:
IndentationError : unindent does not match any outer indentation level (histograms.py, line 1825)
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/shubham/work/MG5_aMC_v2_9_20/pp_hj2/run_26_decayed_1_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
INFO: Done
INFO:
quit
INFO:
MG5_aMC>

and I'm also adding the line no. 1825--

   for wgt_label in all_weights:

Revision history for this message
SHUBHAM YADAV (sy388405) said :
#27

I am sharing the content of line no. 1825 to 1832

   for wgt_label in all_weights:
          for mandatory_attribute in ['PDF','MUR','MUF','MERGING','ALPSFACT']:
          if mandatory_attribute not in wgt_label or wgt_label[mandatory_attribute] is None:
             wgt_label[mandatory_attribute] = '-1'
             elif mandatory_attribute=='PDF':
         wgt_label[mandatory_attribute] = int(wgt_label[mandatory_attribute])
        elif mandatory_attribute in ['MUR','MUF','MERGING','ALPSFACT']:
         wgt_label[mandatory_attribute] = float(wgt_label[mandatory_attribute])

Revision history for this message
Sihyun Jeon (shjeon) said :
#28

this is a python error that is easy to google with keyword "IndentationError : unindent does not match any outer indentation level "

https://stackoverflow.com/questions/26720841/python-indentationerror-unindent-does-not-match-any-outer-indentation-level

your 2nd line for and 3rd line if are aligned which shouldn't be aligned
your 3rd line if and 5th line elif are not aligned which should be aligned
and so on
just make the left blanks to be properly formatted

Revision history for this message
Sihyun Jeon (shjeon) said :
#29
Revision history for this message
SHUBHAM YADAV (sy388405) said :
#30

Thankyou, finally it works. I hope I will disturb you soon .......

Can you help with this problem?

Provide an answer of your own, or ask SHUBHAM YADAV for more information if necessary.

To post a message you must log in.