Dear experts,
I am trying to study top plus higgs production from pp collisions in MG5. The idea is to generate p p > h t at LO ( latter on at NLO) and then decay h > a a and the top via MadSpin. The model I am using for production (Top FCNC effective interactions) does not include the higgs gamma gamma effective interaction, so I was wondering if I could use it just for production and then switch to heft for the Higgs decay in MadSpin. If I naively insert in the MadSpin card the lines:
import model heft /Users/dario/MG5/MadSpin/param_card.dat
decay h > a a
launch
(where param_card.dat is the generic heft one with some minor modifications in the masses -look at attached file-) and run pp > h t in MG5, I get the (expected) param card incompatibility error:
Command "generate_events run_01" interrupted with error:
InvalidCmd : The two param_card seems very different.
So we prefer not to proceed. If you are sure about what you are doing,
you can use the command 'import model MODELNAME PARAM_CARD_PATH --bypass_check'
I added the piece "--bypass_check" in the MadSpin card in order but now I get the following error in MadSpin:
Command "generate p p > t h --no_warning=duplicate" interrupted with error:
InvalidCmd : No amplitudes generated from process Process: g/u/c/d/s/u~/c~/d~/s~ g/u/c/d/s/u~/c~/d~/s~ > t h HIG=1 HIW=1 WEIGHTED=4. Please enter a valid process
Command "generate_events run_01" interrupted with error:
InvalidCmd : No processes generated. Please generate a process first.
I have never used MadSpin before, so I'm not sure what is happening. I would appreciate your help.
Thank you
Dario
MG5_aMC>generate p p > t h
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED=6
INFO: Trying coupling order WEIGHTED=7
INFO: Trying coupling order WEIGHTED=8
INFO: Trying process: g u > t h NP=2 WEIGHTED=8
INFO: Process has 4 diagrams
INFO: Trying process: g c > t h NP=2 WEIGHTED=8
INFO: Process has 4 diagrams
INFO: Process u g > t h added to mirror process g u > t h
INFO: Process c g > t h added to mirror process g c > t h
2 processes with 8 diagrams generated in 0.167 s
Total: 2 processes with 8 diagrams
MG5_aMC>output
INFO: initialize a new directory: PROC_TopFCNC_1
INFO: remove old information in PROC_TopFCNC_1
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g u > t h NP=2 WEIGHTED=8
INFO: Processing color information for process: g u > t h NP=2
INFO: Generating Helas calls for process: g c > t h NP=2 WEIGHTED=8
INFO: Reusing existing color information for process: g c > t h NP=2
INFO: Creating files in directory P0_gq_th
INFO: Generating Feynman diagrams for Process: g u > t h NP=2 WEIGHTED=8
INFO: Generating Feynman diagrams for Process: g c > t h NP=2 WEIGHTED=8
INFO: Finding symmetric diagrams for subprocess group gq_th
History written to /Users/dario/MG5/PROC_TopFCNC_1/Cards/proc_card_mg5.dat
Generated helas calls for 2 subprocesses (8 diagrams) in 0.016 s
Wrote files for 22 helas calls in 0.057 s
Export UFO model to MG4 format
ALOHA: aloha creates FFS12 routines
ALOHA: aloha creates FFVS3 routines
ALOHA: aloha creates FFV1 routines
ALOHA: aloha creates FFV4 routines
ALOHA: aloha creates FFVS15 routines
ALOHA: aloha creates FFV9 routines
ALOHA: aloha creates FFVS10 routines
ALOHA: aloha creates FFVS22 routines
ALOHA: aloha creates FFS3 routines
ALOHA: aloha creates FFV12 routines
ALOHA: aloha creates FFV17 routines
save configuration file to /Users/dario/MG5/PROC_TopFCNC_1/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /Users/dario/MG5/PROC_TopFCNC_1 done.
Type "launch" to generate events from this process, or see
/Users/dario/MG5/PROC_TopFCNC_1/README
Run "open index.html" to see more information about this process.
MG5_aMC>launch
************************************************************
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.2.3 2015-02-10 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /Users/dario/MG5/PROC_TopFCNC_1/Cards/me5_configuration.txt
INFO: load configuration from /Users/dario/MG5/input/mg5_configuration.txt
INFO: load configuration from /Users/dario/MG5/PROC_TopFCNC_1/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
1 Run the pythia shower/hadronization: pythia=OFF
2 Run PGS as detector simulator: pgs=OFF
3 Run Delphes as detector simulator: delphes=OFF
4 Decay particles with the MadSpin module: madspin=ON
5 Add weight to events based on coupling parameters: reweight=OFF
Either type the switch number (1 to 5) to change its default setting,
or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
Type '0', 'auto', 'done' or just press enter when you are done.
[0, 1, 2, 3, 4, 5, auto, done, pythia=ON, ... ][60s to answer]
>done
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
9 / plot : plot_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/MadWidth/...).
Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madspin, 9, enter path, ... ][60s to answer]
>done
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
P0_gq_th
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 11h47 ]
INFO: End survey
refine 10000
Creating Jobs
INFO: Refine results to 10000
P0_gq_th
INFO: Idle: 2, Running: 8, Completed: 2 [ 2.5s ]
INFO: Idle: 0, Running: 0, Completed: 12 [ 7.5s ]
INFO: Combining runs
INFO: finish refine
refine 10000
Creating Jobs
INFO: Refine results to 10000
P0_gq_th
INFO: Combining runs
INFO: finish refine
combine_events
INFO: Combining Events
INFO: Creating root files
** Reading unweighted_events.lhe
** [################################################################] (100.00%)
** Exiting...
=== Results Summary for run: run_01 tag: tag_1 ===
Cross-section : 115.4 +- 0.2205 pb
Nb of events : 10000
running syscalc on mode parton
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
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 > t h
INFO: options:
INFO: detected model: TopFCNC. Loading...
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
import model heft /Users/dario/MG5/MadSpin/param_card.dat --bypass_check
INFO: detected model: heft. Loading...
INFO: Restrict model heft with file models/heft/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
Removed obsolete multiparticles all
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- m- ve~ vm~ vt~ e+ m+ t b t~ b~ z w+ h h1 w- tt- tt+
decay h > a a
launch
INFO: Will use seed 638177697
INFO: generating Madspin matrix element
INFO: We need to recalculate the branching fractions for h
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model heft with file models/heft/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
** Reading unweighted_events.lhe
** [################################################################] (100.00%)
** Exiting...
WARNING: The LO estimate for the width of particle 25
WARNING: differs from the one in the banner by 14 percent
INFO:
INFO: decay channels for h : ( width = 0.004770471 GeV )
INFO: BR d1 d2
INFO: 9.023637e-01 b~ b
INFO: 5.423783e-02 tt+ tt-
INFO: 4.119509e-02 g g
INFO: 2.203346e-03 a a
INFO:
INFO: generating the production square matrix element
INFO: generate p p > t h --no_warning=duplicate;
Command "generate p p > t h --no_warning=duplicate" interrupted with error:
InvalidCmd : No amplitudes generated from process Process: g/u/c/d/s/u~/c~/d~/s~ g/u/c/d/s/u~/c~/d~/s~ > t h HIG=1 HIW=1 WEIGHTED=4. Please enter a valid process
Command "generate_events run_01" interrupted with error:
InvalidCmd : No processes generated. Please generate a process first.
quit
INFO:
INFO:
more information in /Users/dario/MG5/PROC_TopFCNC_1/index.html
Here is the production model param card