MadGraph5_aMC@NLO

Crash for VBF Higgs production in aMC@NLO 2.1.1

Bug #1307493 reported by Josh Bendavid
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
MadGraph5_aMC@NLO
Invalid
Undecided
marco zaro

Bug Description

I get a crash when trying to generate vbf higgs production at nlo.

(same problem whether I use "generate p p > h j j [QCD]" or "generate p p > h j j $$ w+ w- z [QCD]" for the pure vbf contribution

INFO: Checking test output:
INFO: P0_uux_hssx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxu_hcsx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uxc_hsdx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uux_huux
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
Error detected in "launch "
write debug file /data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
aMCatNLOError : Poles do not cancel, run cannot continue
quit
INFO:

bash-4.1$ cat /data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/run_01_tag_1_debug.log
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.1.1 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* and *
#* http://amcatnlo.cern.ch *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#************************************************************
launch
Traceback (most recent call last):
  File "/data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/bin/internal/extended_cmd.py", line 872, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/bin/internal/extended_cmd.py", line 865, in onecmd_orig
    return func(arg, **opt)
  File "/data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/bin/internal/amcatnlo_run_interface.py", line 1184, in do_launch
    self.compile(mode, options)
  File "/data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/bin/internal/amcatnlo_run_interface.py", line 3283, in compile
    self.check_tests(test, this_dir)
  File "/data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/bin/internal/amcatnlo_run_interface.py", line 3297, in check_tests
    return self.parse_check_poles_log(pjoin(dir, '%s.log' % test))
  File "/data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/bin/internal/amcatnlo_run_interface.py", line 3332, in parse_check_poles_log
    raise aMCatNLOError('Poles do not cancel, run cannot continue')
aMCatNLOError: Poles do not cancel, run cannot continue
Value of current Options:
              web_browser : None
              text_editor : None
        cluster_temp_path : None
             syscalc_path : None
            cluster_queue : madgraph
         madanalysis_path : None
                   lhapdf : /data/bendavid/mgstandalone/lhapdfinst/bin/lhapdf-config
                 mg5_path : /data/bendavid/mgstandalone/MG5_aMC_v2_1_1
           cluster_memory : None
    cluster_status_update : (600, 30)
             cluster_time : None
               hepmc_path : None
             pythia8_path : None
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : None
          pythia-pgs_path : None
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
      exrootanalysis_path : None
         fortran_compiler : None
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
                  timeout : 60
                  nb_core : 32
                 run_mode : 2
bash-4.1$

Olivier Mattelaer (olivier-mattelaer)
Changed in mg5amcnlo:
importance: Undecided → Critical
assignee: nobody → marco zaro (marco-zaro)
importance: Critical → Undecided
Revision history for this message
marco zaro (marco-zaro) wrote : Re: [Bug 1307493] [NEW] Crash for VBF Higgs production in aMC@NLO 2.1.1
Download full text (6.6 KiB)

Hi,
in VBF this is expected, as loop diagrams involving pentagons (so with a
gluon exchanged among the quark lines) are not included in the computation.
These diagrams are numerically very small (less than one per mil) and it is
safe not to include them.
The non-inclusion triggers some poles miscancelation, which can be simply
ignored.
In order to have the code running, one should set
#IRPoleCheckThreshold
and
#PrecisionVirtualAtRunTime
to -1d0

inside the Cards/FKSparams.dat file.

Cheers,

Marco

Marco Zaro

On Mon, Apr 14, 2014 at 3:06 PM, Launchpad Bug Tracker <
<email address hidden>> wrote:

> Olivier Mattelaer (olivier-mattelaer) has assigned this bug to you for
> MadGraph5_aMC@NLO:
>
> I get a crash when trying to generate vbf higgs production at nlo.
>
> (same problem whether I use "generate p p > h j j [QCD]" or "generate p
> p > h j j $$ w+ w- z [QCD]" for the pure vbf contribution
>
> INFO: Checking test output:
> INFO: P0_uux_hssx
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20
> (tolerance=1.0e-05)
> INFO: P0_dxu_hcsx
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20
> (tolerance=1.0e-05)
> INFO: P0_uxc_hsdx
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20
> (tolerance=1.0e-05)
> INFO: P0_uux_huux
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> Error detected in "launch "
> write debug file
> /data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/run_01_tag_1_debug.log
> If you need help with this issue please contact us on
> https://answers.launchpad.net/madgraph5
> aMCatNLOError : Poles do not cancel, run cannot continue
> quit
> INFO:
>
> bash-4.1$ cat
> /data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/run_01_tag_1_debug.log
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 5.2.1.1 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* and *
> #* http://amcatnlo.cern.ch ...

Read more...

Revision history for this message
Josh Bendavid (joshbendavid) wrote :
Download full text (11.1 KiB)

Get a bit further now, but now I get

INFO: Output of the failing test:
 Enter xi_i, y_ij to be used in coll/soft tests
  Enter -2 to generate them randomly
 Enter number of tests for soft and collinear limits
 Sum over helicity (0), or random helicity (1)
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11799999999999999
 ==== STAND-ALONE LHAGLUE MODE TO USE LHAPDF ====
 *************************************
 * LHAPDF Version 5.9.1 *
 * Configured for the following: *
 * All PDFs *
 * FULL MEMORY option *
 * Maximum 3 concurrent set(s) *
 *************************************

 >>>>>> PDF description: <<<<<<
 CT10
 H.-L.Lai, M.Guzzi, J. Huston, Z.Li,
 P.M.Nadolsky, J.Pumplin and C.-P.Yuan.
 Reference: arXiv:1007.2241
 PDF fits using the standard CTEQ PDF evolution
 but using the HOPPET alphas_s running solution.
 Excluding the D0 Run-II W asymmetry data.
 This set has 52 member PDFs.
     mem=0 --> central value
     mem=1-52 --> eigenvector sets
 >>>>>> <<<<<<

 ==============================================
 PDFset name ../../lib/PDFsets/CT10.LHgrid
 with 52 members
 ==== initialized. ===========================
 Strong coupling at Mz for PDF is: 0.11798
 New value of alpha_s from PDF lhapdf : 0.11798150967350456
 using LHAPDF
 Give FKS configuration number ("0" loops over all)

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 1
 FKS partons are: i= 6 j= 1
 with PDGs: i= 21 j= -1

 Enter graph number (iconfig), '0' loops over all graphs
Using random seed offsets: 0 , 133 , 0
  with seed 32
 Ranmar initialization seeds 10729 9539
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.92456E-04 0.12500E+03 0.13000E+05
nFKSprocess: 1. Lower bound for tau is 0.12441E-03 0.14500E+03 0.13000E+05
nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.24969E-03 0.20542E+03 0.13000E+05
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.92456E-04 0.12500E+03 0.13000E+05
nFKSprocess: 2. Lower bound for tau is 0.12441E-03 0.14500E+03 0.13000E+05
nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.24969E-03 0.20542E+03 0.13000E+05
nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.12441E-03 0.14500E+03 0.13000E+05
nFKSprocess: 3. Lower bound for tau is (taking resonances into account) 0.24969E-03 0.20542E+03 0.13000E+05
nFKSprocess: 4. Absolute lower bound for tau at the Born is 0.12441E-03 0.14500E+03 0.13000E+05
nFKSprocess: 4. Lower bound for tau is (taking resonances into account) ...

Revision history for this message
marco zaro (marco-zaro) wrote : Re: [Bug 1307493] Re: Crash for VBF Higgs production in aMC@NLO 2.1.1
Download full text (17.6 KiB)

Hi,
which cuts have you put in the run_card for the jets pt?
Cheers,

Marco

On 14 Apr 2014, at 15:55, Josh Bendavid <email address hidden> wrote:

> Get a bit further now, but now I get
>
> INFO: Output of the failing test:
> Enter xi_i, y_ij to be used in coll/soft tests
> Enter -2 to generate them randomly
> Enter number of tests for soft and collinear limits
> Sum over helicity (0), or random helicity (1)
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11799999999999999
> ==== STAND-ALONE LHAGLUE MODE TO USE LHAPDF ====
> *************************************
> * LHAPDF Version 5.9.1 *
> * Configured for the following: *
> * All PDFs *
> * FULL MEMORY option *
> * Maximum 3 concurrent set(s) *
> *************************************
>
>>>>>>> PDF description: <<<<<<
> CT10
> H.-L.Lai, M.Guzzi, J. Huston, Z.Li,
> P.M.Nadolsky, J.Pumplin and C.-P.Yuan.
> Reference: arXiv:1007.2241
> PDF fits using the standard CTEQ PDF evolution
> but using the HOPPET alphas_s running solution.
> Excluding the D0 Run-II W asymmetry data.
> This set has 52 member PDFs.
> mem=0 --> central value
> mem=1-52 --> eigenvector sets
>>>>>>> <<<<<<
>
> ==============================================
> PDFset name ../../lib/PDFsets/CT10.LHgrid
> with 52 members
> ==== initialized. ===========================
> Strong coupling at Mz for PDF is: 0.11798
> New value of alpha_s from PDF lhapdf : 0.11798150967350456
> using LHAPDF
> Give FKS configuration number ("0" loops over all)
>
> =================================================
>
> NEW FKS CONFIGURATION:
> FKS configuration number is 1
> FKS partons are: i= 6 j= 1
> with PDGs: i= 21 j= -1
>
>
> Enter graph number (iconfig), '0' loops over all graphs
> Using random seed offsets: 0 , 133 , 0
> with seed 32
> Ranmar initialization seeds 10729 9539
> nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.92456E-04 0.12500E+03 0.13000E+05
> nFKSprocess: 1. Lower bound for tau is 0.12441E-03 0.14500E+03 0.13000E+05
> nFKSprocess: 1. Lower bound for tau is (taking resonances into account) 0.24969E-03 0.20542E+03 0.13000E+05
> nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.92456E-04 0.12500E+03 0.13000E+05
> nFKSprocess: 2. Lower bound for tau is 0.12441E-03 0.14500E+03 0.13000E+05
> nFKSprocess: 2. Lower bound for tau is (taking resonances into account) 0.24969E-03 0.20542E+03 0.13000E+05
> nFKSprocess: 3. Absolute lower bound for tau at the Born is 0.12441E-03 0.14500E+03 0.13000E+05
> nFKSprocess: 3. Lower bound for tau is (taking resonances into ...

Revision history for this message
Josh Bendavid (joshbendavid) wrote :

I'm using the default run cards.

To add more information, with the default loop_sm model and the changes you suggest for the FKSparams, it works. For the loop_sm-ckm model with flavour-changing contributions it fails.

(Whether the off-diagonal ckm contributions are relevant for this process is another question of course)

Revision history for this message
marco zaro (marco-zaro) wrote :
Download full text (6.8 KiB)

I see
can you try increasing the max_fail parameter in SubProcesses/symmetry_fks_symmetry_fks_test_ME.f
to something higher (0.6 should do)?
I am referring to these lines

c define here the maximum fraction of failures to consider the test
c passed
      double precision max_fail, fail_frac
      parameter (max_fail=0.3d0)

Cheers,

Marco

On 14 Apr 2014, at 17:30, Josh Bendavid <email address hidden> wrote:

> I'm using the default run cards.
>
> To add more information, with the default loop_sm model and the changes
> you suggest for the FKSparams, it works. For the loop_sm-ckm model with
> flavour-changing contributions it fails.
>
> (Whether the off-diagonal ckm contributions are relevant for this
> process is another question of course)
>
> --
> You received this bug notification because you are a bug assignee.
> https://bugs.launchpad.net/bugs/1307493
>
> Title:
> Crash for VBF Higgs production in aMC@NLO 2.1.1
>
> Status in MadGraph5_aMC@NLO Generator:
> New
>
> Bug description:
> I get a crash when trying to generate vbf higgs production at nlo.
>
> (same problem whether I use "generate p p > h j j [QCD]" or "generate
> p p > h j j $$ w+ w- z [QCD]" for the pure vbf contribution
>
> INFO: Checking test output:
> INFO: P0_uux_hssx
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> INFO: P0_dxu_hcsx
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> INFO: P0_uxc_hsdx
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> INFO: P0_uux_huux
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> Error detected in "launch "
> write debug file /data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/run_01_tag_1_debug.log
> If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
> aMCatNLOError : Poles do not cancel, run cannot continue
> quit
> INFO:
>
> bash-4.1$ cat /data/bendavid/mgstandalone/MG5_aMC_v2_1_1/vbfpuredefault/run_01_tag_1_debug.log
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 5...

Read more...

Revision history for this message
Josh Bendavid (joshbendavid) wrote :

Thanks, it works with 0.6.

Olivier Mattelaer (olivier-mattelaer)
Changed in mg5amcnlo:
status: New → Invalid
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.
Subscribing...

Other bug subscribers

Remote bug watches

Bug watches keep track of this bug in other bug trackers.