MadGraph example
[1]:
import hepi
print(hepi.__version__)
import smpl
import numpy as np
import hepi.madgraph as mg
import hepi.util as util
import matplotlib.pyplot as plt
mg.set_path("/opt/MG5_aMC_v2_7_0/")
print (mg.get_path())
0.1.4.42+dirty
/opt/MG5_aMC_v2_7_0/
[ ]:
No on-shell subtraction through madstr
[6]:
params = [
"mastercode_with_gm2.in",
]
pss = [
(1000011,-1000011),
]
for pa,pb in pss:
for param in params:
i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"cteq6l1","cteq66",1., 1.)
li = [i]
li = hepi.mass_scan([i],pa, np.linspace(100,1000,7+8))
mg_dl = mg.run(li,skip=False,madstr=False)
hepi.mass_plot(mg_dl,"LO",pa,logy=True)
hepi.mass_plot(mg_dl,"NLO",pa,logy=True)
hepi.title(plt.gca(),li[0],scenario="mastercode")
Running: 15 jobs
./output/4bec207f684440baa1701f222feffd61ef9ed63a626da77e10441653c9a42084.out
No module named madgraph
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 5.2.7.3 20xx-xx-xx *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/apn/data/de.neuwirthinformatik.Alexander/Development/git/hepi/docs/source/examples/output/4bec207f684440baa1701f222feffd61ef9ed63a626da77e10441653c9a42084.bdir/Cards/amcatnlo_configuration.txt
INFO: load configuration from /opt/MG5_aMC_v2_7_0/input/mg5_configuration.txt
INFO: load configuration from /home/apn/data/de.neuwirthinformatik.Alexander/Development/git/hepi/docs/source/examples/output/4bec207f684440baa1701f222feffd61ef9ed63a626da77e10441653c9a42084.bdir/Cards/amcatnlo_configuration.txt
Using default eps viewer "gv". Set another one in ./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
calculate_xsect -f
INFO: will run in mode: NLO
INFO: Starting run
INFO: Compiling the code
write ./param_card.dat
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
INFO: Using LHAPDF v6.3.0 interface for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
WARNING: Could not compile StdHEP because its source directory could not be found in the SOURCE folder.
Check the MG5_aMC option 'output_dependencies'.
This will prevent the use of HERWIG6/Pythia6 shower.
INFO: Compiling directories...
INFO: Compiling on 16 cores
INFO: Compiling P0_uux_elmelp...
INFO: Compiling P0_ddx_elmelp...
INFO: Compiling P0_uxu_elmelp...
INFO: Compiling P0_dxd_elmelp...
INFO: P0_dxd_elmelp done.
INFO: P0_uux_elmelp done.
INFO: P0_ddx_elmelp done.
INFO: P0_uxu_elmelp done.
INFO: Checking test output:
INFO: P0_uux_elmelp
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_ddx_elmelp
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uxu_elmelp
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxd_elmelp
INFO: Result for test_ME:
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 fixed order NLO
INFO: Setting up grids
INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 02h24 ]
INFO:
Results after grid setup:
Total cross section: 2.708e-01 +- 1.6e-03 pb
INFO: Refining results, step 1
INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 02h24 ]
INFO:
--------------------------------------------------------------
Final results and run summary:
Process p p > 1000011 -1000011 [QCD]
Run at p-p collider (6500.0 + 6500.0 GeV)
Total cross section: 2.692e-01 +- 9.7e-04 pb
--------------------------------------------------------------
INFO: The results of this run and the HwU and GnuPlot files with the plots have been saved in /home/apn/data/de.neuwirthinformatik.Alexander/Development/git/hepi/docs/source/examples/output/4bec207f684440baa1701f222feffd61ef9ed63a626da77e10441653c9a42084.bdir/Events/run_01
INFO: Run complete
INFO:
quit
INFO:
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
On-shell subtraction through madstr
[7]:
params = [
"mastercode_with_gm2.in",
]
pss = [
(1000022,1000022),
]
for pa,pb in pss:
for param in params:
i = hepi.Input(hepi.Order.NLO,13000,pa,pb,param,"cteq6l1","cteq66",1., 1.,model_path="/opt/MG5_aMC_v2_7_0/models/EWKino_NLO_UFO_py3")
li = [i]
li = hepi.mass_scan([i],pa, np.linspace(100,1000,7+8))
mg_dl = mg.run(li,skip=False,madstr=True)
hepi.mass_plot(mg_dl,"LO",pa,logy=True)
hepi.mass_plot(mg_dl,"NLO",pa,logy=True)
hepi.title(plt.gca(),li[0],scenario="mastercode")
Running: 15 jobs
./output/bc70cb3af5e83ab35b937bca29b2536487b5662454482b4f2315ada2001d6774.out
No module named madgraph
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 5.2.7.3 20xx-xx-xx *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/apn/data/de.neuwirthinformatik.Alexander/Development/git/hepi/docs/source/examples/output/bc70cb3af5e83ab35b937bca29b2536487b5662454482b4f2315ada2001d6774.bdir/Cards/amcatnlo_configuration.txt
INFO: load configuration from /opt/MG5_aMC_v2_7_0/input/mg5_configuration.txt
INFO: load configuration from /home/apn/data/de.neuwirthinformatik.Alexander/Development/git/hepi/docs/source/examples/output/bc70cb3af5e83ab35b937bca29b2536487b5662454482b4f2315ada2001d6774.bdir/Cards/amcatnlo_configuration.txt
Using default eps viewer "gv". Set another one in ./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
calculate_xsect -f
INFO: will run in mode: NLO
INFO: Starting run
INFO: Compiling the code
write ./param_card.dat
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
INFO: MadSTR: Forcing width MDL_WSUL to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSUL to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSCL to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSCL to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSUR to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSUR to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSCR to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSCR to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSDL to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSDL to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSSL to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSSL to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSBL to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSBL to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSDR to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSDR to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSSR to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSSR to zero inside param_card.inc
INFO: MadSTR: Forcing width MDL_WSBR to zero inside param_card.inc
INFO: MadSTR: Forcing width MP__MDL_WSBR to zero inside param_card.inc
WARNING: The replacements above ensure poles cancelation, and affect all widths
EXCEPT those which enter the resonance-treatment counterterms, which
are taken from the param_card.
Do NOT set these widths to zero in the param_card.
INFO: Using LHAPDF v6.3.0 interface for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
WARNING: Could not compile StdHEP because its source directory could not be found in the SOURCE folder.
Check the MG5_aMC option 'output_dependencies'.
This will prevent the use of HERWIG6/Pythia6 shower.
INFO: Compiling directories...
INFO: Compiling on 16 cores
INFO: Compiling P0_uux_n1n1...
INFO: Compiling P0_ccx_n1n1...
INFO: Compiling P0_ddx_n1n1...
INFO: Compiling P0_uxu_n1n1...
INFO: Compiling P0_cxc_n1n1...
INFO: Compiling P0_ssx_n1n1...
INFO: Compiling P0_dxd_n1n1...
INFO: Compiling P0_sxs_n1n1...
INFO: Compiling P0_bbx_n1n1...
INFO: Compiling P0_bxb_n1n1...
INFO: P0_bxb_n1n1 done.
INFO: P0_uxu_n1n1 done.
INFO: P0_ccx_n1n1 done.
INFO: P0_dxd_n1n1 done.
INFO: P0_bbx_n1n1 done.
INFO: P0_cxc_n1n1 done.
INFO: P0_sxs_n1n1 done.
INFO: P0_ddx_n1n1 done.
INFO: P0_uux_n1n1 done.
INFO: P0_ssx_n1n1 done.
INFO: Checking test output:
INFO: P0_uux_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_ccx_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_ddx_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_ssx_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uxu_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_cxc_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxd_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_sxs_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_bbx_n1n1
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_bxb_n1n1
INFO: Result for test_ME:
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 fixed order NLO
INFO: Setting up grids
INFO: Idle: 0, Running: 10, Completed: 0 [ current time: 02h31 ]
INFO: Idle: 0, Running: 9, Completed: 1 [ 5.1s ]
INFO: Idle: 0, Running: 8, Completed: 2 [ 5.5s ]
INFO: Idle: 0, Running: 7, Completed: 3 [ 6s ]
INFO: Idle: 0, Running: 6, Completed: 4 [ 6.2s ]
INFO: Idle: 0, Running: 5, Completed: 5 [ 7.6s ]
INFO: Idle: 0, Running: 4, Completed: 6 [ 7.6s ]
INFO: Idle: 0, Running: 3, Completed: 7 [ 7.8s ]
INFO: Idle: 0, Running: 2, Completed: 8 [ 9.5s ]
INFO: Idle: 0, Running: 1, Completed: 9 [ 10.8s ]
INFO: Idle: 0, Running: 0, Completed: 10 [ 11.2s ]
INFO:
Results after grid setup:
Total cross section: 7.307e-06 +- 2.4e-07 pb
INFO: Refining results, step 1
INFO: Idle: 0, Running: 10, Completed: 0 [ current time: 02h31 ]
INFO: Idle: 0, Running: 9, Completed: 1 [ 4.7s ]
INFO: Idle: 0, Running: 8, Completed: 2 [ 4.9s ]
INFO: Idle: 0, Running: 7, Completed: 3 [ 5s ]
INFO: Idle: 0, Running: 6, Completed: 4 [ 5.9s ]
INFO: Idle: 0, Running: 5, Completed: 5 [ 6.4s ]
INFO: Idle: 0, Running: 4, Completed: 6 [ 6.5s ]
INFO: Idle: 0, Running: 3, Completed: 7 [ 7.6s ]
INFO: Idle: 0, Running: 2, Completed: 8 [ 8s ]
INFO: Idle: 0, Running: 1, Completed: 9 [ 9.1s ]
INFO: Idle: 0, Running: 0, Completed: 10 [ 9.9s ]
INFO:
--------------------------------------------------------------
Final results and run summary:
Process p p > 1000022 1000022 [QCD]
Run at p-p collider (6500.0 + 6500.0 GeV)
Total cross section: 7.402e-06 +- 1.2e-07 pb
--------------------------------------------------------------
INFO: The results of this run and the HwU and GnuPlot files with the plots have been saved in /home/apn/data/de.neuwirthinformatik.Alexander/Development/git/hepi/docs/source/examples/output/bc70cb3af5e83ab35b937bca29b2536487b5662454482b4f2315ada2001d6774.bdir/Events/run_01
INFO: Run complete
INFO:
quit
INFO:
stty: 'standard input'stty: 'standard input': Inappropriate ioctl for device
: Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: stty: 'standard input': Inappropriate ioctl for device
'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
stty: 'standard input': Inappropriate ioctl for device
[ ]:
[ ]: