SubProcess doesn't have available phase space. Please check mass spectrum

Asked by Guglielmo Coloretti

Dear all,

I tried looking in all other similar questions but did not find any solutions there.

I am using a BSM model where I am introducing three additional scalars with charges (SU3 x SU2 x UY1):

six1: (6bar, 1, 4/3 )
six2: (6bar, 1, -2/3 )
six3: (1, 1, -2 )

Anyways, while using
generate p p > six2, (six2 > six1 six3, (six1 > u u), (six3 > e- e-))

I got as error:

******************************************************************************************
P1_dd_six2_six2_six1six3_six1_cc_six3_emem,P1_ss_six2_six2_six1six3_six1_cc_six3_emem SubProcesses doesn't have available phase-space.
Please check mass spectrum.

Error when reading /home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/PROC_6b6b1_UFO_0/SubProcesses/P1_dd_six2_six2_six1six3_six1_uu_six3_emem/G1/results.dat
Command "generate_events run_01" interrupted with error:
FileNotFoundError : [Errno 2] No such file or directory: '/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/PROC_6b6b1_UFO_0/SubProcesses/P1_dd_six2_six2_six1six3_six1_uu_six3_emem/G1/results.dat'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/PROC_6b6b1_UFO_0/run_01_tag_1_debug.log'.
Please attach this file to your report.

And the log file is:

generate_events run_01
Traceback (most recent call last):
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/madevent_interface.py", line 2404, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/common_run_interface.py", line 7634, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/madevent_interface.py", line 2602, in run_generate_events
    self.exec_cmd('refine %s' % nb_event, postcmd=False)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/interface/madevent_interface.py", line 3601, in do_refine
    x_improve.launch() # create the ajob for the refinment.
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/madevent/gen_ximprove.py", line 1050, in launch
    self.results = sum_html.collect_result(self.cmd,
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/madevent/sum_html.py", line 753, in collect_result
    P_comb.add_results(os.path.basename(G), path, mfactors[G])
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/madevent/sum_html.py", line 434, in add_results
    oneresult.read_results(filepath)
  File "/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/madgraph/madevent/sum_html.py", line 281, in read_results
    finput = open(filepath)
FileNotFoundError: [Errno 2] No such file or directory: '/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/PROC_6b6b1_UFO_0/SubProcesses/P1_dd_six2_six2_six1six3_six1_uu_six3_emem/G1/results.dat'
Related File: /home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/PROC_6b6b1_UFO_0/SubProcesses/P1_dd_six2_six2_six1six3_six1_uu_six3_emem/G1/results.dat
******************************************************************************************

However, masses for scalars are:
six1: 600 GeV
six2: 2200 GeV
six3: 600 GeV

I tried to isolate the problem, but all parts of the process spit out a cross section independently. For instance, the following processes work:
generate p p > six2, (six2 > six1 six3)
generate six1 > u u
generate six3 > e- e-

I am using version 3.4.1. I tried to use 2.9.13 as well, but got the same error.

I was wondering whether someone has any idea why my generation is not working. I thank you all in advance!

Best regards,
Guglielmo Coloretti

Question information

Language:
English Edit question
Status:
Expired
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:

This question was reopened

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

Hi,

Do you have a log file within the following directory:
/home/gcolor/Desktop/Softwares/MG5_aMC_v3_4_1/PROC_6b6b1_UFO_0/SubProcesses
/P1_dd_six2_six2_six1six3_six1_uu_six3_emem/G1/

If yes what is the content of that file?

Cheers,

Olivier

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#2

Dear Olivier,

Thank you for your swift reply. The content is the following:

 ******************************************************************************************
 Process in group number 1
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.13000000000000003
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23_lo_as_0130_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23lo1: 0.13000000000000000
 Define smin to 3.7478884000000008E-005
 Warning: pt or E min of a jet should in general be >0
 Warning: pt or E min of a jet should in general be >0
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

  External Params
  ---------------------------------

 mdl_qs1 = 0.10000000000000001
 mdl_qs2 = 0.10000000000000001
 mdl_qs0 = 0.10000000000000001
 mdl_lams = 0.10000000000000001
 mdl_lams1 = 0.10000000000000001
 mdl_lams2 = 0.10000000000000001
 mdl_lams3 = 0.10000000000000001
 mdl_lams4 = 0.10000000000000001
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663700000000000E-005
 aS = 0.13000000000000003
 mdl_ymdo = 5.0400000000000002E-003
 mdl_ymup = 2.5500000000000002E-003
 mdl_yms = 0.10100000000000001
 mdl_ymc = 1.2700000000000000
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 172.00000000000000
 mdl_yme = 5.1099999999999995E-004
 mdl_ymm = 0.10566000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MZ = 91.187600000000003
 mdl_Me = 5.1099999999999995E-004
 mdl_MMU = 0.10566000000000000
 mdl_MTA = 1.7769999999999999
 mdl_MU = 2.5500000000000002E-003
 mdl_MC = 1.2700000000000000
 mdl_MT = 172.00000000000000
 mdl_MD = 5.0400000000000002E-003
 mdl_MS = 0.10100000000000001
 mdl_MB = 4.7000000000000002
 mdl_MH = 125.00000000000000
 mdl_MSIX1 = 600.00000000000000
 mdl_MSIX2 = 2200.0000000000000
 mdl_MSIX3 = 600.00000000000000
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_WT = 1.5083359999999999
 mdl_WH = 4.0699999999999998E-003
 mdl_WSIX1 = 0.40000000000000002
 mdl_WSIX2 = 0.40000000000000002
 mdl_WSIX3 = 0.40000000000000002
  Internal Params
  ---------------------------------

 mdl_MZ__exp__2 = 8315.1783937600012
 mdl_MZ__exp__4 = 69142191.720053151
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15625.000000000000
 mdl_complexi = (0.0000000000000000,1.0000000000000000)
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824359746197842
 mdl_sqrt__aEW = 8.8422894590285753E-002
 mdl_ee = 0.31345100004952897
 mdl_MW__exp__2 = 6371.9284088904105
 mdl_sw2 = 0.23369913342182447
 mdl_cw = 0.87538612427783857
 mdl_sqrt__sw2 = 0.48342438232036300
 mdl_sw = 0.48342438232036300
 mdl_g1 = 0.35807170271074895
 mdl_gw = 0.64839716719502682
 mdl_vev = 246.22056907348590
 mdl_vev__exp__2 = 60624.568634871241
 mdl_lam = 0.12886689630821144
 mdl_yb = 2.6995322804122722E-002
 mdl_yc = 7.2944808428161397E-003
 mdl_ydo = 2.8948175943144366E-005
 mdl_ye = 2.9350233942354699E-006
 mdl_ym = 6.0687783137949077E-004
 mdl_ys = 5.8011225600348830E-004
 mdl_yt = 0.98791394091683138
 mdl_ytau = 1.0206529494239589E-002
 mdl_yup = 1.4646398542662327E-005
 mdl_muH = 88.388347648318444
 mdl_I1a11 = (2.89481759431443664E-005,0.0000000000000000)
 mdl_I1a22 = (5.80112256003488297E-004,0.0000000000000000)
 mdl_I1a33 = (2.69953228041227219E-002,0.0000000000000000)
 mdl_I2a11 = (1.46463985426623273E-005,0.0000000000000000)
 mdl_I2a22 = (7.29448084281613968E-003,0.0000000000000000)
 mdl_I2a33 = (0.98791394091683138,0.0000000000000000)
 mdl_I3a11 = (1.46463985426623273E-005,0.0000000000000000)
 mdl_I3a22 = (7.29448084281613968E-003,0.0000000000000000)
 mdl_I3a33 = (0.98791394091683138,0.0000000000000000)
 mdl_I4a11 = (2.89481759431443664E-005,0.0000000000000000)
 mdl_I4a22 = (5.80112256003488297E-004,0.0000000000000000)
 mdl_I4a33 = (2.69953228041227219E-002,0.0000000000000000)
 mdl_ee__exp__2 = 9.8251529432049817E-002
 mdl_sw__exp__2 = 0.23369913342182450
 mdl_cw__exp__2 = 0.76630086657817542
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.36055512754639896
 mdl_G__exp__2 = 1.6336281798666927
  Couplings of 6b6b1_UFO
  ---------------------------------

        GC_36 = 0.00000E+00 0.10000E+00
        GC_41 = 0.00000E+00 0.20000E+00
        GC_42 = 0.00000E+00 0.56569E+00
        GC_43 = 0.00000E+00 0.56569E+00

 Collider parameters:
 --------------------

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = nn23lo1
 alpha_s(Mz)= 0.1300 running at 2 loops.
 alpha_s(Mz)= 0.1300 running at 2 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations 1000 5 3
Enter desired fractional accuracy:
 Desired fractional accuracy: 0.10000000000000001
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Explicitly summing over helicities
Enter Configuration Number:
Running Configuration Number: 1
 Not subdividing B.W.
 Attempting mappinvarients 1 6
 Determine nb_t
 T-channel found: 0
 Completed mapping 6
 about to integrate 10 1000 5 3 10 1
 Using non-zero grid deformation.
 10 dimensions 1000 events 10 invarients 5 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Error opening grid
 Using Uniform Grid! 24
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Masses: 0.504E-02 0.504E-02 0.255E-02 0.255E-02 0.511E-03 0.511E-03
 Set CM energy to 13000.00
 Mapping Graph 1 to config 1
 Determine nb_t
 T-channel found: 0
 Setting BW -1 1 600.00000000000000
 Setting BW -2 2 600.00000000000000
 Setting PDF BW 9 3 2200.0000000000000
 Transforming s_hat BW 2.8639053254437875E-002 5.2071005917159773E-006
Using random seed offsets 1 : 1
  with seed 24
 Ranmar initialization seeds 31059 9398

 ********************************************
 * You are using the DiscreteSampler module *
 * part of the MG5_aMC framework *
 * Author: Valentin Hirschi *
 ********************************************

  Particle 3 4 5 6
      Et > 0.0 0.0 0.0 0.0
       E > 0.0 0.0 0.0 0.0
     Eta < -1.0 -1.0 -1.0 -1.0
   xqcut: 0.0 0.0 0.0 0.0
d R # 3 > -0.0 0.0 0.0 0.0
d R # 4 > -0.0 -0.0 0.0 0.0
d R # 5 > -0.0 -0.0 -0.0 0.0
s min # 3> 0.0 0.0 0.0 0.0
s min # 4> 0.0 0.0 0.0 0.0
s min # 5> 0.0 0.0 0.0 0.0
xqcutij # 3> 0.0 0.0 0.0 0.0
xqcutij # 4> 0.0 0.0 0.0 0.0
xqcutij # 5> 0.0 0.0 0.0 0.0
 No cut BW -1 1
 No cut BW -2 2
 No cut BW -3 3
 alpha_s for scale 2200.1879308311254 is 8.5183497487988619E-002
 Added good helicity 49 63.999999997931965 in event 1 local: 1
 RESET CUMULATIVE VARIABLE
 RESET CUMULATIVE VARIABLE
 Iteration 1 Mean: 0.4783E-06 Abs mean: 0.4783E-06 Fluctuation: 0.148E-07 0.120E-05 93.5%
  1 0.4783E-06 0.4783E-06 +- 0.1480E-07 0.98
 Writing out events 1.0225683281406867E-010 0.97873103998917310
 Relative summed weights:
  0.1000E+01 0.0000E+00
 Relative number of events:
  0.1000E+01 0.0000E+00
 Events:
        1015 0
 Error: failed to reduce to color indices: 0 0

ls status:
input_app.txt
run1_app.log

 ******************************************************************************************

I can see that the issue is with the color flow, although for separate processes and vertices no error is raised.

Best regards,
Guglielmo Coloretti

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#3
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#4

Hi,

> I can see that the issue is with the color flow, although for separate processes and vertices no error is raised.

Those color computation are done in different way and therefore are sensitive to different type of information coming from the UFO. So this does not surprise me too much.

What is surprising is that it fails to combine two uncolor particle together.
So it seems that either six1 or six3 are providing wrong color information.
I would say that the issue is more likely to be six3.
Can you check within the UFO or within the param_card which color is assigned to those particles?

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said (last edit ):
#5

Dear Olivier,

Thank you very much for your help. So the param_card tells me the following information (which are the expected ones):

###################################
## INFORMATION FOR QNUMBERS 9000005
###################################
BLOCK QNUMBERS 9000005 # six1
      1 4 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 6 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000006
###################################
BLOCK QNUMBERS 9000006 # six2
      1 -2 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 6 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000007
###################################
BLOCK QNUMBERS 9000007 # six3
      1 -6 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)

They seem to have the right quantum numbers to me.

Best regards,
Guglielmo Coloretti

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#6

Sorry I am keeping closing and re-opening the question....

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said (last edit ):
#7

I found that at least the processes:

generate p p > six2
generate six2 > six3 six1, (six3 > u u), (six1 > e- e-)

are independently working. Though Pythia8 cannot simulate the showering in neither of them.

Best regards,
Guglielmo Coloretti

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

Everything seems fine.
Can I ask you to put the model online such that I can try to reproduce the error?
Or if you prefer you can send me the model by email: <email address hidden> (please mention the question number 705553)

Cheers,

Olivier

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

Hi,

This seems to fix the issue.
I'm just not sure about what is the LHEF convention for sextet within that file and this seems to not be documented in the reference file. So if you have issue with parton-shower do not hesitate to report it here.

diff --git a/madgraph/iolibs/template_files/addmothers.f b/madgraph/iolibs/template_files/addmothers.f
index 03d5ae0eb..f1151dbdf 100644
--- a/madgraph/iolibs/template_files/addmothers.f
+++ b/madgraph/iolibs/template_files/addmothers.f
@@ -1025,7 +1025,7 @@ c enddo
          do k = 1,2
             found=.false.
             do j=1,nexternal
- if(icol(k,i).eq.icol(1,j).or.icol(k,i).eq.icol(2,j))then
+ if(abs(icol(k,i)).eq.abs(icol(1,j)).or.abs(icol(k,i)).eq.abs(icol(2,j)))then
                   found=.true.
                   goto 10 ! break
                endif

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#10

Hi Olivier,

Thank you very much for your help!! I was wondering where exactly should I add these lines.

Best regards,
Guglielmo Coloretti

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#11
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#12

Either in madgraph/iolibs/template_files/addmothers.f (in that case, you have the exact line number to edit in the above patch)
and then regenerate the code.

Or if you want to edit directly one of the generated code, then you need to edit the file
SubProcesses/addmothers.f (but in that case the line number is likely wrong)

Cheers,

Olivier

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#13

Hi Olivier,

Again, thank you very much for your help. Unfortunately I am not really sure what to here. As far as I understood, the changes are to be made at line 1025 (7 positions from the beginning of it).

If I add the above content there, however, I obtain several errors while launching the generations.

My /madgraph/iolibs/template_files/addmothers.f looks like this at line 1025:

      do k = 1,2
            found=.false.
            do j=1,nexternal
               if(icol(k,i).eq.icol(1,j).or.icol(k,i).eq.icol(2,j))then
                  found=.true.
                  goto 10 ! break
               endif
            enddo
 10 continue
            if (.not.found)then
               call correct_external_flow_epsilon(icol, jpart, maxcolor,
     & icol(k,i))
            endif
         enddo
 20 continue
      enddo
      return
      end

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

then use git with the following commit:

Cheers,

Olivier

> On 6 Mar 2023, at 11:30, Guglielmo Coloretti <email address hidden> wrote:

>
> Question #705553 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/705553
>
> Status: Answered => Open
>
> Guglielmo Coloretti is still having a problem:
> Hi Olivier,
>
> Again, thank you very much for your help. Unfortunately I am not really
> sure what to here. As far as I understood, the changes are to be made at
> line 1025 (7 positions from the beginning of it).
>
> If I add the above content there, however, I obtain several errors while
> launching the generations.
>
> My /madgraph/iolibs/template_files/addmothers.f looks like this at line
> 1025:
>
> do k = 1,2
> found=.false.
> do j=1,nexternal
> if(icol(k,i).eq.icol(1,j).or.icol(k,i).eq.icol(2,j))then
> found=.true.
> goto 10 ! break
> endif
> enddo
> 10 continue
> if (.not.found)then
> call correct_external_flow_epsilon(icol, jpart, maxcolor,
> & icol(k,i))
> endif
> enddo
> 20 continue
> enddo
> return
> end
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#15

Hi Olivier,

Sorry, for the very last time :) I believe your forgot to paste the git commit.

Best regards,
Guglielmo Coloretti

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

Hi,

I actually did not forget it but my mail convert it automatically to an image that was discarded by launchpad.
Thanks for letting me know.
So here it is:

https://github.com/mg5amcnlo/mg5amcnlo/commit/e868becdffcbb752b8f01bcea2513c8ec5e52aad

Cheers,

Olivier

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#17

Dear Olivier,

Thank you very very much!

Now the generation is working and it computes the cross section. However, although Pythia8 and Delphes run without errors, the output files (e.g. the .root file) are empty. Do you have any suggestions where to look at in the code?

Best regards,
Guglielmo Coloretti

Revision history for this message
Launchpad Janitor (janitor) said :
#18

This question was expired because it remained in the 'Open' state without activity for the last 15 days.