Madgraph through batch
Hello,
I want to run Madgraph through a batch system (and have done previously). The submission script settings I used are as follows:
import model heft
define l = l+ l-
define nu = vl vl~
generate p p > h j, h > j j l nu
output higgs_lvjj_200
launch higgs_lvjj_200
pythia=ON
pgs=OFF
done
set nevents 5000
set ptj 200
done
My problem is that this does not set the number of events or transverse jet momentum to the required values - it just ignores them. Note, this worked as desired when I had 'pythia=OFF'. Any thoughts as to how to solve this? Also note: the output from the batch of what happened does not show the option where you change the cards (it did previously).
Many thanks,
James
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- JamesB
- Solved:
- Last query:
- Last reply:
Revision history for this message
|
#1 |
Hi James,
Thanks for the information. This is indeed a bug of the new version.
One quick workaround is to remove your first "done" in the script. i.e. do
import model heft
define l = l+ l-
define nu = vl vl~
generate p p > h j, h > j j l nu
output higgs_lvjj_200
launch higgs_lvjj_200
pythia=ON
pgs=OFF
set nevents 5000
set ptj 200
done
I will fix this for the next version.
Thanks a lot for the information,
Olivier
Revision history for this message
|
#2 |
Hi,
Just to be complete, here is the patch that fix the issue:
=== modified file 'madgraph/
--- madgraph/
+++ madgraph/
@@ -2365,7 +2365,8 @@
if 'allow_arg' in opts:
del opts['allow_arg']
-
+
+ allowed_args +=["0", "done"]
Thanks,
Olivier
Revision history for this message
|
#3 |
Hi Oliver,
I am repeating this process, with the inputs exactly as above:
import model heft
define l = l+ l-
define nu = vl vl~
generate p p > h j, h > j j l nu
output higgs_lvjj_200
launch higgs_lvjj_200
pythia=ON
pgs=OFF
set nevents 5000
set ptj 200
done
But this no longer runs pythia. I only get the unweighted_
Thanks,
James
Revision history for this message
|
#4 |
Hi,
This runs for me nicely (but with 2.6.1 and the development version which has the above patch).
What is your screen output?
Cheers,
Olivier
On 29 Jan 2018, at 17:57, JamesB <<email address hidden>
import model heft
define l = l+ l-
define nu = vl vl~
generate p p > h j, h > j j l nu
output higgs_lvjj_200
launch higgs_lvjj_200
pythia=ON
pgs=OFF
set nevents 5000
set ptj 200
done
Revision history for this message
|
#5 |
Hello Oliver,
The output from the batch job has the following line:
'pythia-pgs not available. Ignore commmand'
I re-downloaded Madgraph (MG5_aMC_v2_6_1) and installed Pythia8 through the interface as well but I still can't get Pythia to run.
Thoughts,
James
Revision history for this message
|
#6 |
Hi James,
Pythia=ON is an old syntax that is still supported but which is not recommended.
Due to that reason, when you use that version, we assume that you still want to have the same stack of program that were
activated when such command were mainstream. (i.e. this command ask to run Pythia6 and not Pythia8)
This is why if you want to keep such script, you first have to install the pythia-pgs package.
If you want to use PY8, you should replace that line by shower=Pythia8
Cheers,
Olivier
> On 1 Feb 2018, at 15:27, JamesB <email address hidden> wrote:
>
> Question #663012 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> JamesB is still having a problem:
> Hello Oliver,
>
> The output from the batch job has the following line:
>
> 'pythia-pgs not available. Ignore commmand'
>
> I re-downloaded Madgraph (MG5_aMC_v2_6_1) and installed Pythia8 through
> the interface as well but I still can't get Pythia to run.
>
> Thoughts,
> James
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#7 |
Thanks Olivier Mattelaer, that solved my question.
Revision history for this message
|
#8 |
In fact, this still does not seem to work.
The output is:
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program ESC[1mdetectorE
| 3. Choose an analysis package (plot/convert) ESC[1manalysisE
| 4. Decay onshell particles ESC[1m madspinESC[0m = ESC[31mOFFESC[0m |
| 5. Add weights to events for new hypp. ESC[1mreweightE
\======
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=OFF' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
set nevents 20000
ESC[1;60mINFO: modify parameter nevents of the run_card.dat to 20000 ESC[0m
set ptj 200
ESC[1;60mINFO: modify parameter ptj of the run_card.dat to 200.0 ESC[0m
INFO: Update the dependent parameter of the param_card.dat
Generating 20000 events with run name run_01
What do you reckon?
Regards,
James
Revision history for this message
|
#9 |
Hi,
Sorry I should miss a point, but now you have correctly
> showerESC[0m = ESC[32mPythia8E
so I do not see anything kind of wrong in that output.
Cheers,
Olivier
> On 1 Feb 2018, at 18:02, JamesB <email address hidden> wrote:
>
> Question #663012 on MadGraph5_aMC@NLO changed:
> https:/
>
> JamesB posted a new comment:
> In fact, this still does not seem to work.
>
> The output is:
>
>
> The following switches determine which programs are run:
> /======
> | 1. Choose the shower/
> | 2. Choose the detector simulation program ESC[1mdetectorE
> | 3. Choose an analysis package (plot/convert) ESC[1manalysisE
> | 4. Decay onshell particles ESC[1m madspinESC[0m = ESC[31mOFFESC[0m |
> | 5. Add weights to events for new hypp. ESC[1mreweightE
> \======
> Either type the switch number (1 to 5) to change its setting,
> Set any switch explicitly (e.g. type 'shower=OFF' at the prompt)
> Type 'help' for the list of all valid option
> Type '0', 'auto', 'done' or just press enter when you are done.
> set nevents 20000
> ESC[1;60mINFO: modify parameter nevents of the run_card.dat to 20000 ESC[0m
> set ptj 200
> ESC[1;60mINFO: modify parameter ptj of the run_card.dat to 200.0 ESC[0m
> INFO: Update the dependent parameter of the param_card.dat
> Generating 20000 events with run name run_01
>
> What do you reckon?
> Regards,
> James
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#10 |
Hi Oliver,
I'm still having the same problem. Madgraph runs fine through the batch but Pythia is not used, so only unweighted_
Many thanks,
James
My run settings are as follows:
import model heft
define l = l+ l-
define nu = ve ve~ vm vm~
generate p p > h j, h > j j l nu
output lvjj
launch lvjj
shower=Pythia8
set nevents 100
set ptj1min 200
done
Revision history for this message
|
#11 |
Hi,
Here is after a run on my university cluster (this slurm scheduler).
Everything seems fine here.
run_01_
run_shower.sh tag_1_pts.dat tag_1_pythia8_
What scheduler are you using?
Cheers,
Olivier
Revision history for this message
|
#12 |
Hello,
Thanks for looking into this.
The problem is that Pythia does not run:
<location>
I believe this error is due to lhapdf.
In MG5_aMC_
set lhapdf to <location>
I have googled this error but I don't really understand it. Any ideas on how to fix it?
Regards,
James
Revision history for this message
|
#13 |
Do you have some log file for PY8?
This might help to understand why PY8 is crashing.
Cheers,
Olivier
> On 23 Feb 2018, at 15:47, JamesB <email address hidden> wrote:
>
> Question #663012 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> JamesB is still having a problem:
> Hello,
>
> Thanks for looking into this.
>
> The problem is that Pythia does not run:
> <location>
>
> I believe this error is due to lhapdf.
>
> In MG5_aMC_
> set lhapdf to <location>
>
> I have googled this error but I don't really understand it. Any ideas on
> how to fix it?
>
> Regards,
> James
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#14 |
This is the debug.log:
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * i* * * *
#* * * *
#* *
#* *
#* VERSION 2.6.1 2017-12-12 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events run_01
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
postcmd=False, printcmd=False)
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
raise Exception, self.fail_msg
Exception: program /data/atlas/
automatic_
cluster_
Revision history for this message
|
#15 |
That's the MG log,
If you look in Events/run_01
you might see some file (and/or directory) related to py8.
Which should contain the py8 log.
Cheers,
Olivier
> On 23 Feb 2018, at 16:07, JamesB <email address hidden> wrote:
>
> Question #663012 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> JamesB is still having a problem:
> This is the debug.log:
>
>
> #******
> #* MadGraph5_
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * i* * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.1 2017-12-12 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https:/
> #* *
> #******
> #* *
> #* Command File for MadEvent *
> #* *
> #* run as ./bin/madevent.py filename *
> #* *
> #******
> generate_events run_01
> Traceback (most recent call last):
> File "/home/
> return self.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> self.exec_
> File "/home/
> stop = Cmd.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> postcmd=False, printcmd=False)
> File "/home/
> stop = Cmd.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> self.cluster.
> File "/home/
> raise Exception, self.fail_msg
> Exception: program /data/atlas/
> Run Options
> -----------
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_
> notification_center : True
> cluster_temp_path : None
> cluster_memory : None (user set)
> cluster_size : 100
> cluster_queue : None (user set)
> nb_core : 24 (user set)
> cluster_time : 24 (user set)
> run_mode : 2
>
> Configuration Options
> -------
> text_editor : None
> cluster_local_path : None
> cluster_
> pythia8_path : /home/spet3829/
> hwpp_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : /home/spet3829/
> thepeg_path : None (user set)
> cluster_type : condor
> madanalysis5_path : /home/spet3829/
> cluster_nb_retry : 1
> eps_viewer : None
> web_browser : None
> syscalc_path : None (user set)
> madanalysis_path : None (user set)
> lhapdf : /home/spet3829/
> f2py_compiler : None
> hepmc_path : None (user set)
> cluster_retry_wait : 300
> fortran_compiler : None
> auto_update : 7 (user set)
> exrootanalysis_path : /home/spet3829/
> timeout : 60
> cpp_compiler : None
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Revision history for this message
|
#16 |
Apologies.
Here is PY8_log.txt
./MG5aMC_
./MG5aMC_
./MG5aMC_
./MG5aMC_
./MG5aMC_
Revision history for this message
|
#17 |
Hi,
This seems to be a path/ld_librarypath isssue related to C++.
I would suggest to reinstall pythia8, if that does not help it might indicate some incompatibility between mounting point where you install the code and where you run it. In that case you should contact the IT responsible of your cluster for instructions.
Cheers,
Olivier
Revision history for this message
|
#18 |
Hi Oliver,
Problem solved. There was an environment issue with my ROOT setup.
Thanks for your help,
James