Weights in event generation Delphes
Dear all,
I have some doubts regarding the general workflow using MadGraph. Probably this question would better be asked at the Delphes team, but I do not get an answer since many months so I try here as well.
After running a simulation with MadGraph, using Pythia8 for the parton shower and Delphes for the detector simulation, I have a .root file as an output which contains all the information about the measurements.
At this point, I impose some cuts on the kinematics of the events and obtain some histograms. My question is the following:
1. I know the weight of each event corresponds to the value of the cross section of the specific event generated (running of parameters and weight of phase space) and therefore I have to fill the histograms weighting each event by its corresponding weight. However, I would like to re-scale the diagrams generated by theoretical computations (production cross section and branching ratios). Is there a way to extrapolate only the relative weight of one event w.r.t. to the other?
Also, I am using MadSpin to decay some final states of the generation. I was wondering whether the weight contains the branching ratios for the additional forced decays by MadSpin or if it corresponds to the value of the cross section of the event before using MadSpin.
I thank you very much for any help and wish you a nice day!
Best regards,
Guglielmo Coloretti
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- Last query:
- Last reply:
|
Revision history for this message
|
#1 |
Actually, your question sounds more related to us than to Delphes.
"Is there a way to extrapolate only the relative weight of one event w.r.t. to the other?"
Here the question is too vague to really answer. I would see that reweighting would be simpler to implement before running Delphes and also likely before running pythia8. At that stage, you have for sure the full information needed in order to perform the proper re-weighting. Obviously if your re-weighting is detector dependent then, this would be done after Delphes and at that stage i have no clue how you can do it, but again you should have the needed information at the Delphes level to do it (I guess).
"I was wondering whether the weight contains the branching ratios for the additional forced decays by MadSpin "
MadSpin does update the weight to add the associated branching ratio.
Cheers,
Olivier
|
Revision history for this message
|
#2 |
Dear Olivier,
Thank you very much for your reply. Now things are much more clear.
As far as when to weight the events, I am not sure how to have access to this information. Let's say I am running MadGraph at NLO with switcher on for Pythia and Delphes. In this case, is MadGraph re-weighting the events before feeding them to Pythia for the shower? And is there a way to check this at the end?
Thank you very much!
Best regards,
Guglielmo Coloretti
|
Revision history for this message
|
#3 |
You can use the re-weighting option of MG5aMC and/or the systematics option of MG5aMC to add additional weight to your events before running the parton-shower. The first option is entended for BSM physics while the second is intended to add weight for the pdf and scale variation.
You can also create your own reweighting algorithm via the plugin option.
If this is not enough, then will likely need to stop your generation after the step where you want to add your weight and then run on your new lhe (or hepmc) file the following tools of the simulation chain.
Cheers,
Olivier
|
Revision history for this message
|
#4 |
Dear Olivier,
Thank you very much for your detailed reply. I would be happy to use both option. However I am not too familiar with writing such algorithms. Anyways, while I am running a generation, at the switcher stage, the option re-weighting says 'not available'.
I am using version MG5_aMC_v3_4_1.
Thank you again!
Best regards,
Guglielmo Coloretti
|
Revision history for this message
|
#5 |
This means that you likely miss the dependency to make it work.
For reweighting, we need to have f2py install in your machine (the easiest is to do pip install numy, and then put the f2py executable within your path or indicate the path to that executable within MG5aMC)
Cheers,
Olivier
> On 13 Jul 2023, at 15:10, Guglielmo Coloretti <email address hidden> wrote:
>
> Question #707261 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Guglielmo Coloretti is still having a problem:
> Dear Olivier,
>
> Thank you very much for your detailed reply. I would be happy to use
> both option. However I am not too familiar with writing such algorithms.
> Anyways, while I am running a generation, at the switcher stage, the
> option re-weighting says 'not available'.
>
> I am using version MG5_aMC_v3_4_1.
>
> Thank you again!
>
> Best regards,
> Guglielmo Coloretti
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#6 |
This means that you likely miss the dependency to make it work.
For reweighting, we need to have f2py install in your machine (the easiest is to do pip install numy, and then put the f2py executable within your path or indicate the path to that executable within MG5aMC)
Cheers,
Olivier
> On 13 Jul 2023, at 15:10, Guglielmo Coloretti <email address hidden> wrote:
>
> Question #707261 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Guglielmo Coloretti is still having a problem:
> Dear Olivier,
>
> Thank you very much for your detailed reply. I would be happy to use
> both option. However I am not too familiar with writing such algorithms.
> Anyways, while I am running a generation, at the switcher stage, the
> option re-weighting says 'not available'.
>
> I am using version MG5_aMC_v3_4_1.
>
> Thank you again!
>
> Best regards,
> Guglielmo Coloretti
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#7 |
Dear Olivier,
Thank you very much. I already have installed f2py3 and f2py3.10. I changed the mg5_configurati
f2py_compiler_py3 = /usr/bin/f2py3.10
or
f2py_compiler_py3 = /usr/bin/f2py3
Now the switcher can be turned ON but I get the following error (I paste the .log file):
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.4.2 20xx-xx-xx *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
original_
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
reweight_
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
out = f(self, *args, **opt)
File "/home/
reweight_name = self.options[
AttributeError: 'NoneType' object has no attribute 'rsplit'
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.4.2 2023-01-20 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set group_subprocesses Auto
set ignore_
set low_mem_
set complex_mass_scheme False
set include_
set gauge unitary
set loop_optimized_
set loop_color_flows False
set max_npoint_
set default_
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion 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~
generate p p > g g
output scratch/
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
#######
#######
## 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
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(
#######
## INFORMATION FOR SMINPUTS
#######
BLOCK SMINPUTS #
1 1.325070e+02 # aewm1
2 1.166390e-05 # gf
3 1.300000e-01 # as (note that parameter not used if you use a pdf set)
#######
## INFORMATION FOR YUKAWA
#######
BLOCK YUKAWA #
5 4.700000e+00 # ymb
6 1.730000e+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 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
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# 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=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=elastic photon of proton,*
# +/-3=PDF of electron/positron beam *
# +/-4=PDF of muon/antimuon beam *
#******
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 alpha_s and its evol. *
# pdlabel: lhapdf=LHAPDF (installation needed) [1412.7420] *
# iww=Improved Weizsaecker-
# eva=Effective W/Z/A Approx. [21yy.zzzzz] *
# none=No PDF, same as lhapdf with lppx=0 *
#******
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_
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.
2 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
! 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.
#******
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#******
# Parton level cuts definition *
#******
0.0 = dsqrt_shat ! minimal shat for full process
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# 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 *
#******
0.4 = drjj ! min distance between jets
-1.0 = drjjmax ! max distance between jets
#******
# Minimum and maximum invariant mass for pairs *
#******
0.0 = mmjj ! min invariant mass of a jet pair
-1.0 = mmjjmax ! max invariant mass of a jet pair
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_
#******
# Inclusive cuts *
#******
0.0 = xptj ! minimum pt for at least one jet
#*****
# Control the pt's of the jets sorted by pt *
#*****
0.0 = ptj1min ! minimum pt for the leading jet in pt
0.0 = ptj2min ! minimum pt for the second jet in pt
-1.0 = ptj1max ! maximum pt for the leading jet in pt
-1.0 = ptj2max ! maximum pt for the second jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#*****
# Control the Ht(k)=Sum of k leading jets *
#*****
0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)
0.0 = ihtmin !inclusive Ht for all partons (including b)
-1.0 = ihtmax !inclusive Ht for all partons (including b)
#*****
# WBF cuts *
#*****
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# 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_
Thank you!!!
Best regards,
Guglielmo Coloretti
|
Revision history for this message
|
#8 |
Hi,
This has been fixed in 3.5.1 version release, can you try that version?
Cheers,
Olivier
> On 13 Jul 2023, at 15:45, Guglielmo Coloretti <email address hidden> wrote:
>
> Question #707261 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Guglielmo Coloretti is still having a problem:
> Dear Olivier,
>
> Thank you very much. I already have installed f2py3 and f2py3.10. I
> changed the mg5_configurati
>
> f2py_compiler_py3 = /usr/bin/f2py3.10
>
> or
>
> f2py_compiler_py3 = /usr/bin/f2py3
>
>
> Now the switcher can be turned ON but I get the following error (I paste the .log file):
>
> #******
> #* MadGraph5_
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.4.2 20xx-xx-xx *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https:/
> #* *
> #******
> #* *
> #* Command File for MadEvent *
> #* *
> #* run as ./bin/madevent.py filename *
> #* *
> #******
> generate_events
> Traceback (most recent call last):
> File "/home/
> return self.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> self.run_
> File "/home/
> original_fct(obj, *args, **opts)
> File "/home/
> self.exec_
> File "/home/
> stop = Cmd.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> reweight_
> File "/home/
> self.exec_cmd(line, precmd=True)
> File "/home/
> stop = Cmd.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> out = f(self, *args, **opt)
> File "/home/
> reweight_name = self.options[
> AttributeError: 'NoneType' object has no attribute 'rsplit'
> Run Options
> -----------
> stdout_level : None
>
> MadEvent Options
> ----------------
> automatic_
> notification_center : True
> run_mode : 2
> cluster_queue : None
> cluster_time : None
> cluster_size : 100
> cluster_memory : None
> nb_core : 16 (user set)
> cluster_temp_path : None
>
> Configuration Options
> -------
> pythia8_path : /home/gcolor/
> hwpp_path : None (user set)
> thepeg_path : None (user set)
> hepmc_path : None (user set)
> madanalysis_path : None (user set)
> madanalysis5_path : /home/gcolor/
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : /home/gcolor/
> exrootanalysis_path : /home/gcolor/
> syscalc_path : None (user set)
> rivet_path : None
> yoda_path : None
> lhapdf : /home/gcolor/
> lhapdf_py2 : None
> lhapdf_py3 : /home/gcolor/
> timeout : 60
> f2py_compiler : /usr/bin/f2py3.10 (user set)
> f2py_compiler_py2 : None
> f2py_compiler_py3 : /usr/bin/f2py3.10 (user set)
> web_browser : None
> eps_viewer : None
> text_editor : None
> fortran_compiler : None
> cpp_compiler : None
> auto_update : 7 (user set)
> cluster_type : condor
> cluster_
> cluster_nb_retry : 1
> cluster_local_path : None
> cluster_retry_wait : 300
> #******
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.4.2 2023-01-20 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https:/
> #* *
> #******
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #******
> set group_subprocesses Auto
> set ignore_
> set low_mem_
> set complex_mass_scheme False
> set include_
> set gauge unitary
> set loop_optimized_
> set loop_color_flows False
> set max_npoint_
> set default_
> set max_t_for_channel 99
> set zerowidth_tchannel True
> set nlo_mixed_expansion 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~
> generate p p > g g
> output scratch/
> #######
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> #######
> #######
> ## 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
> 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(
> #######
> ## INFORMATION FOR SMINPUTS
> #######
> BLOCK SMINPUTS #
> 1 1.325070e+02 # aewm1
> 2 1.166390e-05 # gf
> 3 1.300000e-01 # as (note that parameter not used if you use a pdf set)
> #######
> ## INFORMATION FOR YUKAWA
> #######
> BLOCK YUKAWA #
> 5 4.700000e+00 # ymb
> 6 1.730000e+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 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
>
> #******
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat MadEvent *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/
> # *
> # 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=
> #******
> # Collider type and energy *
> # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=elastic photon of proton,*
> # +/-3=PDF of electron/positron beam *
> # +/-4=PDF of muon/antimuon beam *
> #******
> 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 alpha_s and its evol. *
> # pdlabel: lhapdf=LHAPDF (installation needed) [1412.7420] *
> # iww=Improved Weizsaecker-
> # eva=Effective W/Z/A Approx. [21yy.zzzzz] *
> # none=No PDF, same as lhapdf with lppx=0 *
> #******
> 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_
> 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.
> #******
> None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
> {} = bias_parameters ! Specifies the parameters of the module.
> #
> #******
> # Parton level cuts definition *
> #******
> 0.0 = dsqrt_shat ! minimal shat for full process
> #
> #
> #******
> # BW cutoff (M+/-bwcutoff*
> #******
> 15.0 = bwcutoff ! (M+/-bwcutoff*
> #******
> # 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 *
> #******
> 0.4 = drjj ! min distance between jets
> -1.0 = drjjmax ! max distance between jets
> #******
> # Minimum and maximum invariant mass for pairs *
> #******
> 0.0 = mmjj ! min invariant mass of a jet pair
> -1.0 = mmjjmax ! max invariant mass of a jet pair
> {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
> {'default': False} = mxx_only_
> ! to pairs of particle/
> #******
> # Inclusive cuts *
> #******
> 0.0 = xptj ! minimum pt for at least one jet
> #******
> # Control the pt's of the jets sorted by pt *
> #******
> 0.0 = ptj1min ! minimum pt for the leading jet in pt
> 0.0 = ptj2min ! minimum pt for the second jet in pt
> -1.0 = ptj1max ! maximum pt for the leading jet in pt
> -1.0 = ptj2max ! maximum pt for the second jet in pt
> 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
> #******
> # Control the Ht(k)=Sum of k leading jets *
> #******
> 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
> 0.0 = ihtmin !inclusive Ht for all partons (including b)
> -1.0 = ihtmax !inclusive Ht for all partons (including b)
> #******
> # WBF cuts *
> #******
> 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
> 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
> #******
> # 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_
>
>
>
>
> Thank you!!!
>
> Best regards,
> Guglielmo Coloretti
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#9 |
Dear Olivier,
Sorry it took me a bit long to install it since I had to install LHAPDF independently. Anyways, not it does not break but it tells me the reweighting is failing:
ERROR: Reweighting failed
This is the full terminal output:
gcolor@
*******
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 3.5.1 2023-07-11 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* and *
* http://
* *
* Type 'help' for in-line help. *
* Type 'tutorial' to learn how MG5 works *
* Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
* Type 'tutorial MadLoop' to learn how MadLoop works *
* *
*******
load MG5 configuration from input/mg5_
set fastjet to /home/gcolor/
eMELA-config does not seem to correspond to a valid eMELA-config executable.
Please set the 'fastjet'variable to the full (absolute) /PATH/TO/
MG5_aMC> set eMELA /PATH/TO/
set ninja to /home/gcolor/
set collier to /home/gcolor/
set lhapdf to /home/gcolor/
set lhapdf to /home/gcolor/
Using default text editor "vi". Set another one in ./input/
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
Loading default model: sm
INFO: Restrict model sm with file models/
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = e+ mu+
Defined multiparticle l- = e- mu-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h w- ta- ta+
MG5_aMC>generate p p > z [QCD]
The current model sm does not allow to generate loop corrections of type ['QCD'].
MG5_aMC now loads 'loop_sm'.
import model loop_sm
INFO: Restrict model loop_sm with file models/
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ b t b~ t~ z w+ h w- ta- ta+
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS QCD+2*QED
Setting the born squared orders automatically in the process definition to QED^2<=2 QCD=0.
If this is not what you need, please regenerate with the correct orders.
WARNING: Use of multiparticles is non-trivial for NLO process generation and depends on the orders included, the process considered, as well as the PDF set chosen. See appendix D of arXiv:1804.10017 [hep-ph] for some guidance.
INFO: Generating FKS-subtracted matrix elements for born process: u u~ > z QCD^2=2 QED^2=2 (1 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: c c~ > z QCD^2=2 QED^2=2 (2 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: d d~ > z QCD^2=2 QED^2=2 (3 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: s s~ > z QCD^2=2 QED^2=2 (4 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: u~ u > z QCD^2=2 QED^2=2 (5 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: c~ c > z QCD^2=2 QED^2=2 (6 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: d~ d > z QCD^2=2 QED^2=2 (7 / 8)
INFO: Generating FKS-subtracted matrix elements for born process: s~ s > z QCD^2=2 QED^2=2 (8 / 8)
INFO: Generating virtual matrix element with MadLoop for process: u u~ > z QCD^2=2 QED^2=2 (1 / 8)
INFO: Generating virtual matrix element with MadLoop for process: c c~ > z QCD^2=2 QED^2=2 (2 / 8)
INFO: Generating virtual matrix element with MadLoop for process: d d~ > z QCD^2=2 QED^2=2 (3 / 8)
INFO: Generating virtual matrix element with MadLoop for process: s s~ > z QCD^2=2 QED^2=2 (4 / 8)
INFO: Generating virtual matrix element with MadLoop for process: u~ u > z QCD^2=2 QED^2=2 (5 / 8)
INFO: Generating virtual matrix element with MadLoop for process: c~ c > z QCD^2=2 QED^2=2 (6 / 8)
INFO: Generating virtual matrix element with MadLoop for process: d~ d > z QCD^2=2 QED^2=2 (7 / 8)
INFO: Generating virtual matrix element with MadLoop for process: s~ s > z QCD^2=2 QED^2=2 (8 / 8)
INFO: Generated 8 subprocesses with 48 real emission diagrams, 8 born diagrams and 8 virtual diagrams
MG5_aMC>output scratch
INFO: Writing out the aMC@NLO code, using optimized Loops
INFO: directory /home/gcolor/
If you continue this directory will be deleted and replaced.
Do you want to continue? [y, n][60s to answer]
>
INFO: initialize a new directory: scratch
INFO: remove old information in scratch
INFO: Generating real emission matrix-elements...
INFO: Generating Helas calls for FKS process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (1 / 8)
INFO: Processing color information for process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: c c~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (2 / 8)
INFO: Generating Helas calls for FKS process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (3 / 8)
INFO: Reusing existing color information for process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: s s~ > z [ all = QCD QED ] QCD^2=2 QED^2=2 (4 / 8)
INFO: Generating Helas calls for FKS process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2 (5 / 8)
INFO: Processing color information for process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: c~ c > z [ all = QCD QED ] QCD^2=2 QED^2=2 (6 / 8)
INFO: Generating Helas calls for FKS process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2 (7 / 8)
INFO: Reusing existing color information for process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating Helas calls for FKS process: s~ s > z [ all = QCD QED ] QCD^2=2 QED^2=2 (8 / 8)
INFO: ... Done
INFO: Writing files in P0_uux_z (1 / 4)
INFO: Creating files in directory V0_uux_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: u u~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Writing files in P0_ddx_z (2 / 4)
INFO: Creating files in directory V0_ddx_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: d d~ > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Writing files in P0_uxu_z (3 / 4)
INFO: Creating files in directory V0_uxu_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: u~ u > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Writing files in P0_dxd_z (4 / 4)
INFO: Creating files in directory V0_dxd_z
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2
INFO: Generating born Feynman diagrams for Process: d~ d > z [ all = QCD QED ] QCD^2=2 QED^2=2
History written to /home/gcolor/
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 16 routines in 0.312 s
save configuration file to /home/gcolor/
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
^[[AINFO: Generate web pages
^[[AType "launch" to generate events from this process, or see
/home/gcolor/
Run "open index.html" to see more information about this process.
MG5_aMC>launch
INFO: *******
* *
* W E L C O M E to M A D G R A P H 5 *
* a M C @ N L O *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 3.5.1 2023-07-11 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /home/gcolor/
INFO: load configuration from /home/gcolor/
INFO: load configuration from /home/gcolor/
INFO: set f2py compiler to /usr/bin/f2py3
INFO: set f2py compiler to /usr/bin/f2py3
Using default text editor "vi". Set another one in ./input/
launch auto
The following switches determine which programs are run:
/================== Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = HERWIG6 | OFF|PYTHIA6Q|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON|NLO|NLO_TREE|LO |
| 6. Run MadAnalysis5 on the events generated | madanalysis = Not Avail | Please install module |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>3
The following switches determine which programs are run:
/================== Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = OFF | PYTHIA6Q|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON|NLO|NLO_TREE|LO |
| 6. Run MadAnalysis5 on the events generated | madanalysis = Not Avail | Please install module |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>3
The following switches determine which programs are run:
/================== Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA6Q | PYTHIA6PT|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON|NLO|NLO_TREE|LO |
| 6. Run MadAnalysis5 on the events generated | madanalysis = Not Avail | Please install module |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>3
The following switches determine which programs are run:
/================== Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA6PT | PYTHIA8|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON|NLO|NLO_TREE|LO |
| 6. Run MadAnalysis5 on the events generated | madanalysis = Not Avail | Please install module |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>3
The following switches determine which programs are run:
/================== Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA8 | HERWIG6|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON|NLO|NLO_TREE|LO |
| 6. Run MadAnalysis5 on the events generated | madanalysis = Not Avail | Please install module |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>4
The following switches determine which programs are run:
/================== Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA8 | HERWIG6|
| 4. Decay onshell particles | madspin = ON | onshell|OFF |
| 5. Add weights to events for new hypp. | reweight = OFF | ON|NLO|NLO_TREE|LO |
| 6. Run MadAnalysis5 on the events generated | madanalysis = Not Avail | Please install module |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>5
The following switches determine which programs are run:
/================== Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA8 | HERWIG6|
| 4. Decay onshell particles | madspin = ON | onshell|OFF |
| 5. Add weights to events for new hypp. | reweight = ON | NLO|NLO_TREE|LO|OFF |
| 6. Run MadAnalysis5 on the events generated | madanalysis = Not Avail | Please install module |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>
INFO: will run in mode: aMC@NLO
INFO: modify parameter parton_shower of the run_card.dat to PYTHIA8
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_card.dat |
| 4. reweight : reweight_card.dat |
| 5. shower : shower_card.dat |
\------
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ][90s to answer]
>set nevents 10
INFO: modify parameter nevents of the run_card.dat to 10
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_card.dat |
| 4. reweight : reweight_card.dat |
| 5. shower : shower_card.dat |
\------
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ]
>
INFO: NLO reweighting is on ON. Automatically set store_rwgt_info to True
INFO: modify parameter store_rwgt_info of the run_card.dat to True
INFO: modify parameter extralibs of the shower_card.dat to pythia8 z dl stdc++
INFO: modify parameter extrapaths of the shower_card.dat to ../lib /home/gcolor/
INFO: Update the dependent parameter of the param_card.dat
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.119
INFO: Starting run
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
INFO: For gauge cancellation, the width of 'Z' has been set to zero.
INFO: Compiling source...
Gtk-Message: 18:18:34.443: Not loading module "atk-bridge": The functionality is provided by GTK natively. Please try to not load it.
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 16 cores
INFO: Compiling P0_uux_z...
INFO: Compiling P0_ddx_z...
INFO: Compiling P0_uxu_z...
INFO: Compiling P0_dxd_z...
INFO: P0_uxu_z done.
INFO: P0_dxd_z done.
INFO: P0_uux_z done.
INFO: P0_ddx_z done.
INFO: Checking test output:
INFO: P0_uux_z
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_ddx_z
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uxu_z
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxd_z
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 16 cores
INFO: Cleaning previous results
INFO: Doing NLO matched to parton shower
INFO: Setting up grids
INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 18h19 ]
sum of cpu time of last step: 0 second
INFO: Determining the number of unweighted events per channel
Intermediate results:
Random seed: 33
Total cross section: 2.359e+06 +- 1.2e+04 pb
Total abs(cross section): 2.662e+06 +- 1.2e+04 pb
INFO: Computing upper envelope
INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 18h19 ]
INFO: Idle: 0, Running: 3, Completed: 1 [ current time: 18h19 ]
INFO: Idle: 0, Running: 2, Completed: 2 [ 0.021s ]
INFO: Idle: 0, Running: 1, Completed: 3 [ 0.027s ]
INFO: Idle: 0, Running: 0, Completed: 4 [ 0.037s ]
sum of cpu time of last step: 0 second
INFO: Updating the number of unweighted events per channel
Intermediate results:
Random seed: 33
Total cross section: 2.358e+06 +- 6.9e+03 pb
Total abs(cross section): 2.668e+06 +- 6.6e+03 pb
INFO: Generating events
INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 18h19 ]
sum of cpu time of last step: 0 second
INFO: Doing reweight
INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 18h19 ]
INFO: Collecting events
INFO:
----
Summary:
Process p p > z [QCD]
Run at p-p collider (6500.0 + 6500.0 GeV)
Number of events generated: 10
Total cross section: 2.358e+06 +- 6.9e+03 pb
----
INFO: The /home/gcolor/
INFO: Events generated
reweight -from_cards
ERROR: Reweighting failed
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 > z
INFO: options:
INFO: detected model: loop_sm. Loading...
set f2py compiler to /usr/bin/f2py3
set f2py compiler to /usr/bin/f2py3
set fastjet to /home/gcolor/
set ninja to /home/gcolor/
set collier to /home/gcolor/
set lhapdf to /home/gcolor/
set lhapdf to /home/gcolor/
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
launch
INFO: Will use seed 540515021
WARNING: set the mass of the b-quark to its value in the param_card.dat: 4.7 GeV
INFO: We need to recalculate the branching fractions for t~,z,t,w+,w-
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /home/gcolor/
INFO: Run "set stdout_level DEBUG" before import for more information.
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 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 > z --no_warning=
INFO: Done 0.8882
INFO: generating the full matrix element squared (with decay)
INFO: generate p p > z, z > all all QCD=99 --no_warning=
INFO: Done 1.972
INFO: generate matrix element for decay only (1 - > N).
INFO: output standalone_msF /home/gcolor/
INFO: Done 0.9608
INFO: Compiling code
INFO: detect independant decays
INFO: Done in 0.0305597782135
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: 18h20
INFO: Event 6/75 : 0.12s
INFO: Not enough events for at least one production mode.
INFO: This is ok as long as you don't reuse the max weight for other generations.
INFO:
INFO: Decaying the events...
INFO: Total number of events written: 10/10
INFO: Average number of trial points per production event: 6.3
INFO: Branching ratio to allowed decays: 1
INFO: Number of events with weights larger than max_weight: 0
INFO: Number of subprocesses 5
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/gcolor/
INFO: MadSpin Done
INFO: Preparing MCatNLO run
INFO: Using LHAPDF v6.5.4 interface for PDFs
INFO: Compiling MCatNLO for PYTHIA8...
INFO: ... done
INFO: Showering events...
INFO: (Running in /home/gcolor/
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 18h20 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 0.25s ]
INFO: The file /home/gcolor/
It contains showered and hadronized events in the HEPMC format obtained by showering the parton-level event file /home/gcolor/
INFO: Run complete
INFO:
quit
INFO:
MG5_aMC>
Thank you very much!
Best regards,
Guglielmo Coloretti
|
Revision history for this message
|
#10 |
So this seems to crash quite early. Not that clear why...
Maybe an f2py issue?
can you run in debug mode (adding --debug)
This might add some printout to help find what the issue is.
Cheers,
Olivier
|
Revision history for this message
|
#11 |
Thanks Olivier Mattelaer, that solved my question.
