Setting the correct parameters for Madgraph 5
Dear Experts,
I am trying to run the following processes with Madgraph 5 v2.2.2 with the following settings:
*******
import model Higgs_Effective
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
generate g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > l+ vl), ( t~ > b~ w-, w- > j j)
add process g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > j j), ( t~ > b~ w-, w- > l- vl~)
add process g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > j j), ( t~ > b~ w-, w- > j j)
add process g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > l+ vl), ( t~ > b~ w-, w- > l- vl~)
output /export/
set cut_decays F
set pdlabel 'lhapdf'
set lhaid 10800
set clusinfo F
set pdfwgt F
set ptj 0
set pta 0
set ptl 0
set etaj -1
set etaa -1
set etal -1
set drjj 0
set drll 0
set draa 0
set draj 0
set drjl 0
set dral 0
set MT 172.5
set nevents 1
set ebeam1 6500
set ebeam2 6500
***************
Even after setting the parameters for QED/QCD by hand I see strange output as below:
load MG5 configuration from input/mg5_
set lhapdf to /cvmfs/
fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
Please set the 'fastjet'variable to the full (absolute) /PATH/TO/
MG5_aMC> set fastjet /PATH/TO/
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/
import /export/
The import format was not given, so we guess it as command
import model Higgs_Effective
INFO: Change particles name to pass to MG5 convention
Defined multiparticle p = g u d s u~ d~ s~ c c~
Defined multiparticle j = g u d s u~ d~ s~ c c~
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 d s u~ d~ s~ a ve vm vt e- ve~ vm~ vt~ e+ c t b c~ t~ b~ z w+ h h1 w- mu- ta- mu+ ta+
define l+ = e+ mu+ ta+
Defined multiparticle l+ = e+ mu+ ta+
define l- = e- mu- ta-
Defined multiparticle l- = e- mu- ta-
generate g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > l+ vl), ( t~ > b~ w-, w- > j j)
INFO: Trying process: g g > t t~ HIW=0 QCD=2 QED=4 HIG=0
INFO: Process has 3 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t > b w+ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w+ > e+ ve WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > e+ vm WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > e+ vt WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > mu+ ve WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > mu+ vm WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > mu+ vt WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > ta+ ve WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > ta+ vm WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > ta+ vt WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t~ > b~ w- WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w- > d u~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > d c~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > s u~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > s c~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
1 processes with 12 diagrams generated in 0.128 s
Total: 1 processes with 12 diagrams
add process g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > j j), ( t~ > b~ w-, w- > l- vl~)
INFO: Trying process: g g > t t~ HIW=0 QCD=2 QED=4 HIG=0
INFO: Process has 3 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t > b w+ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w+ > u d~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > u s~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > d~ c WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > s~ c WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t~ > b~ w- WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w- > e- ve~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > e- vm~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > e- vt~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > mu- ve~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > mu- vm~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > mu- vt~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > ta- ve~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > ta- vm~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > ta- vt~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
1 processes with 12 diagrams generated in 0.108 s
Total: 2 processes with 24 diagrams
add process g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > j j), ( t~ > b~ w-, w- > j j)
INFO: Trying process: g g > t t~ HIW=0 QCD=2 QED=4 HIG=0
INFO: Process has 3 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t > b w+ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w+ > u d~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > u s~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > d~ c WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > s~ c WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t~ > b~ w- WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w- > d u~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > d c~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > s u~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > s c~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
1 processes with 13 diagrams generated in 0.068 s
Total: 3 processes with 37 diagrams
add process g g > t t~ HIG=0 HIW=0 QED=4 QCD=2, ( t > b w+, w+ > l+ vl), ( t~ > b~ w-, w- > l- vl~)
INFO: Trying process: g g > t t~ HIW=0 QCD=2 QED=4 HIG=0
INFO: Process has 3 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t > b w+ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w+ > e+ ve WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > e+ vm WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > e+ vt WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > mu+ ve WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > mu+ vm WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w+ > mu+ vt WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > ta+ ve WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > ta+ vm WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w+ > ta+ vt WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: t~ > b~ w- WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: w- > e- ve~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > e- vm~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > e- vt~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > mu- ve~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > mu- vm~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
INFO: Trying process: w- > mu- vt~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > ta- ve~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > ta- vm~ WEIGHTED=2 HIW=1 HIG=1
INFO: Trying process: w- > ta- vt~ WEIGHTED=2 HIW=1 HIG=1
INFO: Process has 1 diagrams
1 processes with 11 diagrams generated in 0.057 s
Total: 4 processes with 48 diagrams
I am not sure, from where it picks up the parameters HIG/HIW as 1 even though I set it to 0.
Is their any other way where I can fix this?
Much Thanks for any help
Ruchika
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask ruchika for more information if necessary.