Hi,

In order to use MadSpin you need a model which can handle both the production and the decay which is clearly not your case.
On the other hand you can create such model easily with the following command:
import model YOURMODEL
add model plugin hgg_plugin

This will create a model named YOURMODEL__hgg_plugin.
Model that you can use with MadSpin and/or directly to generate your events (at LO only).

Cheers,

Olivier






On 28 Apr 2015, at 08:43, Dario <email address hidden> wrote:

> Here is the production model param card
> 
> ** Attachment added: "param_card_TopFCNC.dat"
>   https://bugs.launchpad.net/mg5amcnlo/+bug/1449521/+attachment/4386572/+files/param_card_TopFCNC.dat
> 
> -- 
> You received this bug notification because you are subscribed to
> MadGraph5_aMC@NLO.
> https://bugs.launchpad.net/bugs/1449521
> 
> Title:
>  Problems bypassing param Card in MadSpin
> 
> Status in MadGraph5_aMC@NLO Generator:
>  New
> 
> Bug description:
>  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
> 
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/mg5amcnlo/+bug/1449521/+subscriptions