Single Z boson production using e+e- collision
Hi
I want to generate this process e+ e- --> gamma*/Z, Z --> b b~ in Madgraph. What I did is
generate e+ e- > Z, Z > b b~
But it returned an error resulting to zero cross section after executing launch:
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program | detector = OFF | PGS|Delphes |
| 3. Choose an analysis package (plot/convert) | analysis = MadAnalysis5 | OFF|ExRoot |
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON |
\======
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia6' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. pythia8 : pythia8_card.dat |
| 4. madanalysis5_parton : madanalysis5_
| 5. madanalysis5_hadron : madanalysis5_
\------
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, pythia8, 4, enter path, ... ][90s to answer]
>
INFO: Update the dependent parameter of the param_card.dat
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: P1_ll_z_z_bbx
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 12h13 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 3.4s ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 3.4s ]
INFO: End survey
refine 10000
Creating Jobs
INFO: Refine results to 10000
INFO: Generating 10000.0 unweighted events.
INFO: Effective Luminosity 1.2e+103 pb^-1
INFO: need to improve 0 channels
Survey return zero cross section.
Typical reasons are the following:
1) A massive s-channel particle has a width set to zero.
2) The pdf are zero for at least one of the initial state particles
or you are using maxjetflavor=4 for initial state b:s.
3) The cuts are too strong.
Please check/correct your param_card and/or your run_card.
Zero result detected: See https:/
I clearly made some mistakes at this point. How do I correctly generate this process in Madgraph?
Sincerely,
Dennis
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- 2020-04-28
- Last query:
- 2020-04-28
- Last reply:
- 2020-04-28
Hi Dennis,
Your syntax
e+ e- > Z, Z > b b~
requires a ONSHELL Z.
So you have two issues:
1. You do not have the gamma* contribution
2. If your beam energy is not close of the Z peak then the onshell condition related to your syntax will not be possible leading to a 0 result. The onshell condition is that the invariant mass of the Z is on the pole mass up to a given number of times the width. (that number can be set in the run_card: bwcutoff parameter.
The syntax that you are looking for is likely
e+ e- > b b~
Cheers,
Olivier
> On 27 Apr 2020, at 06:18, Dennis Arogancia <email address hidden> wrote:
>
> New question #690230 on MadGraph5_aMC@NLO:
> https:/
>
> Hi
>
> I want to generate this process e+ e- --> gamma*/Z, Z --> b b~ in Madgraph. What I did is
>
> generate e+ e- > Z, Z > b b~
>
> But it returned an error resulting to zero cross section after executing launch:
>
> The following switches determine which programs are run:
> /======
> | 1. Choose the shower/
> | 2. Choose the detector simulation program | detector = OFF | PGS|Delphes |
> | 3. Choose an analysis package (plot/convert) | analysis = MadAnalysis5 | OFF|ExRoot |
> | 4. Decay onshell particles | madspin = OFF | ON|onshell |
> | 5. Add weights to events for new hypp. | reweight = OFF | ON |
> \======
> Either type the switch number (1 to 5) to change its setting,
> Set any switch explicitly (e.g. type 'shower=Pythia6' at the prompt)
> Type 'help' for the list of all valid option
> Type '0', 'auto', 'done' or just press enter when you are done.
>>
> Do you want to edit a card (press enter to bypass editing)?
> /------
> | 1. param : param_card.dat |
> | 2. run : run_card.dat |
> | 3. pythia8 : pythia8_card.dat |
> | 4. madanalysis5_parton : madanalysis5_
> | 5. madanalysis5_hadron : madanalysis5_
> \------
> 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, pythia8, 4, enter path, ... ][90s to answer]
>>
> INFO: Update the dependent parameter of the param_card.dat
> Generating 10000 events with run name run_01
> survey run_01
> INFO: compile directory
> compile Source Directory
> Using random number seed offset = 21
> INFO: Running Survey
> Creating Jobs
> Working on SubProcesses
> INFO: P1_ll_z_z_bbx
> INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 12h13 ]
> INFO: Idle: 0, Running: 0, Completed: 1 [ 3.4s ]
> INFO: Idle: 0, Running: 0, Completed: 1 [ 3.4s ]
> INFO: End survey
> refine 10000
> Creating Jobs
> INFO: Refine results to 10000
> INFO: Generating 10000.0 unweighted events.
> INFO: Effective Luminosity 1.2e+103 pb^-1
> INFO: need to improve 0 channels
> Survey return zero cross section.
> Typical reasons are the following:
> 1) A massive s-channel particle has a width set to zero.
> 2) The pdf are zero for at least one of the initial state particles
> or you are using maxjetflavor=4 for initial state b:s.
> 3) The cuts are too strong.
> Please check/correct your param_card and/or your run_card.
> Zero result detected: See https:/
>
>
> I clearly made some mistakes at this point. How do I correctly generate this process in Madgraph?
>
>
> Sincerely,
>
> Dennis
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Dennis Arogancia (dcarogancia) said : | #2 |
Hi Olivier,
Thank you for your reply. You were right I was looking e+ e- > b b~. I was thinking of the process syntax in Pythia for single boson production: WeakSingleBoson
generate_events run_02
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/
return self.run_
File "/home/
"reconstructed event file for reconstruction '%s'."%
MadGraph5Error: MadAnalysis5 failed to produce the reconstructed event file for reconstruction 'BasicReco'.
I encountered no error during the installation of MadAnalysis package.
Sincerely,
Dennis
That is something that you should report to the MadAnalysis5 team.
Cheers,
Olivier
Dennis Arogancia (dcarogancia) said : | #4 |
Thanks Olivier Mattelaer, that solved my question.