Generating double higgs with specific final states

Asked by Ana Sc

Dear Madgraph experts,

I am sorry if my questions has already been answered somewhere, I went through the answers and didn't find a full answer to my problem.
I am trying to generate dihiggs decaying in wwzz, in final states: for w all, for z leptonic + hadronic
What I am currently doing is:

generate p p > h h [QCD]
launch
madspin=none
decay h > w+ w- > l+ vl l- vl~
decay h > w+ w- > l+ vl j j
decay h > w+ w- > j j l- vl~
decay h > w+ w- > j j j j
decay h > z z > l+ l- l+ l-
decay h > z z > l+ l- j j
decay h > z z > j j j j

1)To my understanding in this way I can have all the possible combinations of hh decays (wwww.wwzz,zzzz) is this correct? Is there a way to ask only for hh>wwzz (w decaying to any final state, z only to jj or l+l-)?

2) When I look into lhe file (decayed) and take an example of event:
14 0 +1.6871474e-04 6.65252500e+02 7.54677100e-03 9.77857300e-02
       21 -1 0 0 502 501 +0.0000000000e+00 +0.0000000000e+00 +2.2687198516e+03 2.2687198516e+03 0.0000000000e+00 0.0000e+00 -1.0000e+00
       21 -1 0 0 501 502 -0.0000000000e+00 -0.0000000000e+00 -4.8767683136e+01 4.8767683136e+01 0.0000000000e+00 0.0000e+00 1.0000e+00
       25 2 1 2 0 0 +2.7887079943e+02 +1.2325786866e+02 +1.2678493149e+03 1.3099726373e+03 1.2500000000e+02 0.0000e+00 0.0000e+00
       25 2 1 2 0 0 -2.7887079943e+02 -1.2325786866e+02 +9.5210285351e+02 1.0075148975e+03 1.2500000000e+02 0.0000e+00 0.0000e+00
        4 1 4 4 503 0 -4.1108205784e+01 -5.3524322331e+00 +1.9480574219e+02 1.9916779434e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
      -24 2 4 4 0 0 -2.2022680381e+02 -9.5121201751e+01 +6.5509098151e+02 7.0217404512e+02 7.9726453834e+01 0.0000e+00 0.0000e+00
        1 1 6 6 504 0 -9.8434697131e+01 -6.6440595383e+01 +4.1195111631e+02 4.2872772775e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
       -2 1 6 6 0 504 -1.2179210668e+02 -2.8680606367e+01 +2.4313986520e+02 2.7344631736e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       -3 1 4 4 0 503 -1.7535789289e+01 -2.2784234435e+01 +1.0220612795e+02 1.0617305607e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       24 2 3 3 0 0 +1.8547353971e+02 +6.6781720478e+01 +8.3746259722e+02 8.6407905988e+02 8.0179722273e+01 0.0000e+00 0.0000e+00
        2 1 10 10 505 0 +2.1453510778e+01 +2.2228207394e+01 +5.2396799295e+01 6.0825742124e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       -1 1 10 10 0 505 +1.6402002893e+02 +4.4553513082e+01 +7.8506579792e+02 8.0325331776e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 +9.3113660177e+01 +4.3045236028e+01 +3.6500543450e+02 3.7914642722e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
      -14 1 3 3 0 0 +2.8359757468e-01 +1.3430911283e+01 +6.5381274214e+01 6.6747140937e+01 0.0000000000e+00 0.0000e+00 1.0000e+00
<mgrwt>
Here I see gg(21) and hh(25) (which is okay), but only one w+(24) and one w-(-24), while I was expecting to see either 2w+ 2w-, or w+w- 2z (etc..). But, I see 8 final state particles, 6 quarks and mu vmu, which would mean that 3 w decayed hadronicaly and 1 leptonicaly. I am a bit conufsed by this and I am not able to find an event where I have w+w-z z in lhe file.

Thank you in advance for your help!

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Ana Sc
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Concerning the "missing" W/Z in the lhef, this is normal (and it would be WRONG to have both here). I have attached the associated FAQ which should explains that point.
FAQ #2173: “Intermediate Boson missing from events in LHE output file”.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#2

For the second question,

Doing this,

decay h > w+ w- > l+ vl l- vl~
decay h > w+ w- > l+ vl j j
decay h > w+ w- > j j l- vl~
decay h > w+ w- > j j j j
decay h > z z > l+ l- l+ l-
decay h > z z > l+ l- j j
decay h > z z > j j j j

assumes no interference between ww and zz which is likely not true for the four jet final state (But I have not clue how important it is)

1)To my understanding in this way I can have all the possible combinations of hh decays (wwww.wwzz,zzzz) is this correct? Is there a way to ask only for hh>wwzz (w decaying to any final state, z only to jj or l+l-)?

If you have exactly two decay line, this should normally trigger the fact that each higgs will pick one specific line.
so using should in principle be what you want. (please double check in case spinmode none does not follow that convention)
define zdec = l+ l- j
define wdec = l+ vl l- vl~ j
decay h > w+ w- > 4wdec
decay h > z > 4zdec

Cheers,

Olivier

Revision history for this message
Ana Sc (anasc) said :
#3

Thank you for your very fast and accurate answer.
Indeed, having exactly two decay lines produces only wwzz case.