Issue while decaying the double higgs production

Asked by Ramkrishna on 2020-03-06

Dear experts,

I am trying to generate the double Higgs production along with its decay chain (p p > h h > a a w+ w- > a a j j j j ).

Since I was getting an error while producing this. So, I tried to first generate only p p > h h > a a w+ w-. I did this like:

1. Proc card:

define p = g u c d s u~ c~ d~ s~
define j = g u c d s b u~ c~ d~ s~ b~
define h1 = h
define h2 = h
generate g g > h1 h2 [QCD]
output sm_hh_decay_madspin

2. Madspin card:

set spinmode none
# set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
decay h1 > a a
decay h2 > w+ w-
launch

If I do like this then it works without any issue. But, the final LHE file does not have the decay information. You can see my log here[2] and also I copied the one event information from decayed LHE file [1].

I went through all the links suggested by the launchpad before posting this question, but I am unable to find the solution of my problem. Sorry, if I missed any link where there was a solution.

Please let me know if you need more information.

with regards,
Ram

-------------------------------------------

[1]
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 82 # gh
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
</slha>
<MGGenerationInfo>
# Number of Events : 10000
# Integrated weight (pb) : 0.014504
</MGGenerationInfo>
<madspin>
set seed 706301291
set spinmode none
decay h1 > a a
decay h2 > w+ w-
launch
</madspin>
</header>
<init>
2212 2212 6.500000e+03 6.500000e+03 0 0 247000 247000 -4 1
   +1.4504000e-02 +3.7773000e-05 +1.4504000e-02 0
<generator name='MadGraph5_aMC@NLO' version='2.7.0'>please cite 1405.0301 </generator>
</init>
<event>
 4 0 +1.4504000e-02 4.38307500e+02 7.54677100e-03 1.03126700e-01
       21 -1 0 0 502 501 +0.0000000000e+00 +0.0000000000e+00 +5.8594273064e+02 5.8594273064e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       21 -1 0 0 501 502 -0.0000000000e+00 -0.0000000000e+00 -8.1967693431e+01 8.1967693431e+01 0.0000000000e+00 0.0000e+00 1.0000e+00
       25 1 1 2 0 0 -1.4918265873e+02 +5.4888242495e+01 +3.8070633355e+02 4.3108061571e+02 1.2500000000e+02 0.0000e+00 0.0000e+00
       25 1 1 2 0 0 +1.4918265873e+02 -5.4888242495e+01 +1.2326870365e+02 2.3682980836e+02 1.2500000000e+02 0.0000e+00 0.0000e+00
<mgrwt>
<rscale> 2 0.43830754E+03</rscale>
<asrwt>0</asrwt>
<pdfrwt beam="1"> 1 21 0.90145033E-01 0.43830754E+03</pdfrwt>
<pdfrwt beam="2"> 1 21 0.12610415E-01 0.43830754E+03</pdfrwt>
<totfact> 0.51756454E+04</totfact>
</mgrwt>
</event>

-------------------------------------------

[2]

(base) rams-mbp:sm_hh_decay_madspin ram$ ./bin/generate_events
************************************************************
* *
* 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.7.0 20xx-xx-xx *
* *
* 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/ramkrishna/cernbox/post_doc_ihep/MG5_aMC_v2_7_0/sm_hh_decay_madspin/Cards/me5_configuration.txt
INFO: load configuration from /Users/ramkrishna/cernbox/post_doc_ihep/MG5_aMC_v2_7_0/input/mg5_configuration.txt
INFO: load configuration from /Users/ramkrishna/cernbox/post_doc_ihep/MG5_aMC_v2_7_0/sm_hh_decay_madspin/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
generate_events
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|======== other options ========\
| 1. Choose the shower/hadronization program | shower = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = ExRoot | MadAnalysis5|OFF |
| 4. Decay onshell particles | madspin = ON | onshell|OFF |
| 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=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>3
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|======== other options ========\
| 1. Choose the shower/hadronization program | shower = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = MadAnalysis5 | OFF|ExRoot |
| 4. Decay onshell particles | madspin = ON | onshell|OFF |
| 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=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>3
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|======== other options ========\
| 1. Choose the shower/hadronization program | shower = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = OFF | ExRoot|MadAnalysis5 |
| 4. Decay onshell particles | madspin = ON | onshell|OFF |
| 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=Pythia8' 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. madspin : madspin_card.dat |
| 4. MadLoopParams : MadLoopParams.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, 4, enter path, ... ][90s to answer]
>
INFO: Update the dependent parameter of the param_card.dat
Generating 10000 events with run name run_03
survey run_03
INFO: compile directory
compile Source Directory
Using random number seed offset = 39
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: P0_gg_hh
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 19h41 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 7s ]
INFO: P0_gg_hh/G1 is at 0.015 +- 0.000522 pb. Now submitting iteration #2.
INFO: Idle: 1, Running: 0, Completed: 5 [ 17.3s ]
INFO: P0_gg_hh/G1 is at 0.01462 +- 0.000135 pb. Now submitting iteration #3.
INFO: Idle: 0, Running: 1, Completed: 7 [ 36.8s ]
INFO: Survey finished for P0_gg_hh/G1 at 0.01458 +- 4.53e-05 pb
INFO: Idle: 0, Running: 0, Completed: 9 [ 37.3s ]
INFO: End survey
refine 10000
Creating Jobs
INFO: Refine results to 10000
INFO: Generating 10000.0 unweighted events.
INFO: Effective Luminosity 788697.620191 pb^-1
INFO: need to improve 1 channels
Current estimate of cross-section: 0.014581 +- 4.5306e-05
INFO: Idle: 2, Running: 2, Completed: 0 [ current time: 19h42 ]
INFO: Idle: 0, Running: 1, Completed: 3 [ 9.4s ]
INFO: P0_gg_hh/G1 is at 1670/11367 (15%) event. Resubmit 8 job at iteration 1.
INFO: Idle: 3, Running: 4, Completed: 6 [ 19.5s ]
INFO: Idle: 0, Running: 3, Completed: 10 [ 29.5s ]
INFO: P0_gg_hh/G1 is at 5092/11456 (44%) event. Resubmit 15 job at iteration 2.
INFO: Idle: 10, Running: 4, Completed: 15 [ 39.3s ]
INFO: Idle: 6, Running: 4, Completed: 19 [ 48.9s ]
INFO: Idle: 1, Running: 4, Completed: 24 [ 58.6s ]
INFO: Idle: 0, Running: 2, Completed: 27 [ 1m 6s ]
INFO: found enough event for P0_gg_hh/G1
INFO: Idle: 0, Running: 0, Completed: 30 [ 1m 6s ]
INFO: Combining runs
INFO: finish refine
refine 10000
Creating Jobs
INFO: Refine results to 10000
INFO: Generating 10000.0 unweighted events.
INFO: Effective Luminosity 792884.721456 pb^-1
INFO: need to improve 0 channels
Current estimate of cross-section: 0.014504 +- 3.7773e-05
INFO: Idle: 0, Running: 0, Completed: 0 [ current time: 19h43 ]
INFO: Combining runs
INFO: finish refine
INFO: Combining Events

  === Results Summary for run: run_03 tag: hh_decay_max_weight_ps_point_off ===

     Cumulative sequential time for this run: 5m03s
     Cross-section : 0.0145 +- 3.777e-05 pb
     Nb of events : 10000

INFO: No version of lhapdf. Can not run systematics computation
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: g g > h1 h2
INFO: options:
INFO: detected model: loop_sm. Loading...
/Users/ramkrishna/cernbox/post_doc_ihep/MG5_aMC_v2_7_0/HEPTools/lhapdf6/bin/lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

set ninja to /Users/ramkrishna/cernbox/post_doc_ihep/MG5_aMC_v2_7_0/HEPTools/lib
set collier to /Users/ramkrishna/cernbox/post_doc_ihep/MG5_aMC_v2_7_0/HEPTools/lib
set spinmode none
decay h1 > a a
decay h2 > w+ w-
launch
INFO: Will use seed 706301291
decaying event number 100 [0.027773141861 s]
decaying event number 200 [0.0612690448761 s]
decaying event number 300 [0.0885779857635 s]
decaying event number 400 [0.118600130081 s]
decaying event number 500 [0.147078990936 s]
decaying event number 600 [0.178290128708 s]
decaying event number 700 [0.20571398735 s]
decaying event number 800 [0.233953952789 s]
decaying event number 900 [0.264091014862 s]
decaying event number 1000 [0.295778989792 s]
decaying event number 2000 [0.604248046875 s]
decaying event number 3000 [0.893077135086 s]
decaying event number 4000 [1.1714220047 s]
decaying event number 5000 [1.48783111572 s]
decaying event number 6000 [1.77781009674 s]
decaying event number 7000 [2.06083607674 s]
decaying event number 8000 [2.34033894539 s]
decaying event number 9000 [2.61759305 s]
INFO: The decayed event file has been moved to the following location:
INFO: /Users/ramkrishna/cernbox/post_doc_ihep/MG5_aMC_v2_7_0/sm_hh_decay_madspin/Events/run_03_decayed_1/unweighted_events.lhe.gz
INFO: MadSpin Done
INFO: storing files of previous run
INFO: Done
quit
INFO:
(base) rams-mbp:sm_hh_decay_madspin ram$

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2020-03-06
Last reply:
2020-03-06

Dear Ram,

Looks like MadSpin is confused about the particles h1 h2.
Seems that the definition of such label are not passed to madspin correctly.
Actually they are no point to use a multi-particle name for a single particle.
So you should use everywhere "h" instead of h1/h2

Secondly, the decay
decay h > w+ w-
can not be handle by MadSpin since this decay is IMPOSSIBLE
(your syntax requires that the three particles are onshell wich can not be the case in the standard model).

Third, the decay
decay h > a a
is not possible in the standard model
you should either use the heft model (which use an EFT vertex)

So in order to test I have use the following madspin_card
set spinmode none
import model heft ParamCardPath
decay h > a a
decay h > b b~

Where ParamCardPath is the path to the param_card related to the heft model that I need to use
to compute the decay matrix-element.

and then I get the expected result:

8 0 +2.9581156e-05 3.61096100e+02 7.54677100e-03 1.05815100e-01
       21 -1 0 0 502 501 +0.0000000000e+00 +0.0000000000e+00 +3.6195806723e+01 3.6195806723e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       21 -1 0 0 501 502 -0.0000000000e+00 -0.0000000000e+00 -9.0059035447e+02 9.0059035447e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
       25 2 1 2 0 0 +1.5705905215e+01 +3.8079372625e+00 -7.6756738389e+02 7.7784694166e+02 1.2500000000e+02 0.0000e+00 0.0000e+00
       25 2 1 2 0 0 -1.5705905215e+01 -3.8079372625e+00 -9.6827163865e+01 1.5893921954e+02 1.2500000000e+02 0.0000e+00 0.0000e+00
       22 1 4 4 0 0 +9.5074282576e+00 -4.5391088567e+01 +4.4527402125e+00 4.6589365834e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       22 1 4 4 0 0 -2.5213333472e+01 +4.1583151305e+01 -1.0127990407e+02 1.1234985370e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
        5 1 3 3 503 0 +3.3579489441e+01 +1.4201333886e+01 -7.4828393576e+02 7.4918635767e+02 4.7000000000e+00 0.0000e+00 1.0000e+00
       -5 1 3 3 0 503 -1.7873584223e+01 -1.0393396623e+01 -1.9283448294e+01 2.8660584161e+01 4.7000000000e+00 0.0000e+00 1.0000e+00

Cheers,

Olivier

> On 6 Mar 2020, at 20:08, Ramkrishna <email address hidden> wrote:
>
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s b u~ c~ d~ s~ b~
> define h1 = h
> define h2 = h
> generate g g > h1 h2 [QCD]
> output sm_hh_decay_madspin

Actually the webpage for madspin says that you can use multi-particle label for final state.
Not for initial states. So the line
decay h1 > a a
decay h2 > w+ w-
were just ignored.
It will not be trivial to support multiparticle at that level so I will not do it.

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Ramkrishna for more information if necessary.

To post a message you must log in.