MadGraph5_aMC@NLO

e+e- at NLO

Asked by Konrad Tywoniuk

Launching the process e+ e- > b b~ [QCD] crashes (while e.g. p p > b b~ [QCD] works fine). I have explicitly set the pdf set options to 0 to generate lepton lepton collisions. I'm running without the parton shower (just interested in the lhe events).

I'm on a Mac, the OSX is Mavericks, XCode is installed together with command line tools, and I'm using gfortran/g++ 4.8.2.

The launch terminates with the following message:

------------------------------------------------------------------------------------------------------------------------
INFO: Starting run
INFO: Compiling the code
INFO: Lepton-Lepton collision: Ignoring 'pdlabel' and 'lhaid' in the run_card.
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 4 cores
INFO: Compiling P0_epem_bbx...
INFO: P0_epem_bbx done.
INFO: Checking test output:
INFO: P0_epem_bbx
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: Starting run
INFO: Using 4 cores
INFO: Cleaning previous results
INFO: Generating events without running the shower.
INFO: Setting up grid
INFO: Idle: 0, Running: 0, Completed: 2 [ current time: 11h43 ]
grep: P*/GF*/res_0: No such file or directory
Traceback (most recent call last):
  File "./sumres.py", line 130, in <module>
    (totABS, math.sqrt(errABS), math.sqrt(errABS)/totABS *100.,tot, math.sqrt(err), math.sqrt(err)/tot *100.)
ZeroDivisionError: float division by zero
INFO: Return code of the event collection: None
INFO: Output of the event collection:
N of directories: 2
Determining the number of unweighted events per channel
2 jobs did not terminate correctly
P0_epem_bbx/GF1
P0_epem_bbx/GF2
Found 0 correctly terminated jobs
0 2
Integrated abs(cross-section)
Integrated cross-section

Error detected in "launch auto "
write debug file /Users/konradtywoniuk/phys/MG5_aMC_v2_1_1/ee_bbar_NLO_1/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
aMCatNLOError : An error occurred during the collection of results.
 Please check the .log files inside the directories which failed.
quit
INFO:
------------------------------------------------------------------------------------------------------------------------

and the debug file contains the following information:

------------------------------------------------------------------------------------------------------------------------
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.1.1 2014-03-31 *
#* *
#* 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 auto
Traceback (most recent call last):
  File "/Users/konradtywoniuk/phys/MG5_aMC_v2_1_1/madgraph/interface/extended_cmd.py", line 872, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/konradtywoniuk/phys/MG5_aMC_v2_1_1/madgraph/interface/extended_cmd.py", line 865, in onecmd_orig
    return func(arg, **opt)
  File "/Users/konradtywoniuk/phys/MG5_aMC_v2_1_1/madgraph/interface/amcatnlo_run_interface.py", line 1185, in do_launch
    evt_file = self.run(mode, options)
  File "/Users/konradtywoniuk/phys/MG5_aMC_v2_1_1/madgraph/interface/amcatnlo_run_interface.py", line 1506, in run
    self.cross_sect_dict = self.read_results(output, mode)
  File "/Users/konradtywoniuk/phys/MG5_aMC_v2_1_1/madgraph/interface/amcatnlo_run_interface.py", line 1600, in read_results
    'Please check the .log files inside the directories which failed.')
aMCatNLOError: An error occurred during the collection of results.
Please check the .log files inside the directories which failed.
Value of current Options:
              text_editor : None
              web_browser : Safari
        cluster_temp_path : None
       group_subprocesses : False
                  timeout : 60
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : lhapdf-config
ignore_six_quark_processes : False
             stdout_level : 20
    loop_optimized_output : True
    cluster_status_update : (600, 30)
               hepmc_path : None
             pythia8_path : None
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
      output_dependencies : external
          pythia-pgs_path : None
                 mg5_path : /Users/konradtywoniuk/phys/MG5_aMC_v2_1_1
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
      exrootanalysis_path : None
         fortran_compiler : None
                      OLP : MadLoop
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : None
                  fastjet : fastjet-config
                    gauge : unitary
      complex_mass_scheme : False
                  nb_core : 4
                 run_mode : 2
 ------------------------------------------------------------------------------------------------------------------------

Thanks for any help! ;)

Konrad

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said : 		#1

Dear Konrad,
I managed to have the same process running on my mac.
Which center-of-mass energy are you running at?
Can you please paste the content of the proc/run/param cards (attaching them won't work), and also of the log file inside P0_epem_bbx/GF1?

That wopuld help to understand what is going wrong...
Cheers,

Marco

Revision history for this message
Konrad Tywoniuk (konrad-tywoniuk) said : 		#2

Dear Marco,

I was generating events at very high com energy, (6500 + 6500 GeV). When I
reduced the energy to 45.8 + 45.8 GeV the program generated output.

Thanks a lot for the rapid answer!

Just for reference, here is the output of the file log_MINT0.txt (it seems
the other log files in that directory were ok):

 ===============================================================
 INFO: MadFKS read these parameters from FKS_params.dat

 ===============================================================
  > IRPoleCheckThreshold = 1.0000000000000001E-005
  > PrecisionVirtualAtRunTime = 1.0000000000000000E-003
  > NHelForMCoverHels = 4
  > VirtualFraction = 1.0000000000000000
  > MinVirtualFraction = 5.0000000000000001E-003
 ===============================================================
 Process in group number 0
 No PDF is used, alpha_s(MZ) from param_card is used
 Value of alpha_s from param_card: 0.11799999999999999
 The default order of alpha_s running is fixed to 2
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

 MU_R = 91.188000000000002
 aEWM1 = 132.50700000000001
 mdl_Gf = 1.1663900000000000E-005
 aS = 0.11799999999999999
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 173.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MT = 173.00000000000000
 mdl_MB = 4.7000000000000002
 mdl_MZ = 91.188000000000002
 mdl_MH = 125.00000000000000
 mdl_MTA = 1.7769999999999999
 mdl_WT = 1.4915000000000000
 mdl_WZ = 2.4414039999999999
 mdl_WW = 2.0476000000000001
 mdl_WH = 6.3823389999999999E-003
  Internal Params
  ---------------------------------

 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_CKM22 = 1.0000000000000000
 mdl_lhv = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM22 = 1.0000000000000000
 mdl_conjg__CKM33 = 1.0000000000000000
 mdl_Ncol = 3.0000000000000000
 mdl_CA = 3.0000000000000000
 mdl_TF = 0.50000000000000000
 mdl_CF = 1.3333333333333333
 mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
 mdl_MZ__exp__2 = 8315.2513440000002
 mdl_MZ__exp__4 = 69143404.913893804
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15625.000000000000
 mdl_Ncol__exp__2_m_1 = 8.0000000000000000
 mdl_MB__exp__2 = 22.090000000000003
 mdl_MT__exp__2 = 29929.000000000000
 mdl_Ncol__exp__2_m_1_0 = 8.0000000000000000
 mdl_aEW = 7.5467711139788835E-003
 mdl_MW = 80.419002445756163
 mdl_sqrt__aEW = 8.6872153846781555E-002
 mdl_ee = 0.30795376724436879
 mdl_MW__exp__2 = 6467.2159543705357
 mdl_sw2 = 0.22224648578577766
 mdl_cw = 0.88190334743339216
 mdl_sqrt__sw2 = 0.47143025548407230
 mdl_sw = 0.47143025548407230
 mdl_g1 = 0.34919219678733299
 mdl_gw = 0.65323293034757990
 mdl_v = 246.21845810181637
 mdl_v__exp__2 = 60623.529110035903
 mdl_lam = 0.12886910601690263
 mdl_yb = 2.6995554250465490E-002
 mdl_yt = 0.99366614581500623
 mdl_ytau = 1.0206617000654717E-002
 mdl_muH = 88.388347648318430
 mdl_AxialZUp = -0.18517701861793787
 mdl_AxialZDown = 0.18517701861793787
 mdl_VectorZUp = 7.5430507588273299E-002
 mdl_VectorZDown = -0.13030376310310560
 mdl_VectorAUp = 0.20530251149624587
 mdl_VectorADown = -0.10265125574812294
 mdl_VectorWmDxU = 0.23095271737156670
 mdl_AxialWmDxU = -0.23095271737156670
 mdl_VectorWpUxD = 0.23095271737156670
 mdl_AxialWpUxD = -0.23095271737156670
 mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
 mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
 mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
 mdl_Vector_tbGp = (-0.96667059156454072 , 0.0000000000000000 )
 mdl_Axial_tbGp = ( -1.0206617000654716 , -0.0000000000000000 )
 mdl_Vector_tbGm = ( 0.96667059156454072 , 0.0000000000000000 )
 mdl_Axial_tbGm = ( -1.0206617000654716 , -0.0000000000000000 )
 mdl_gw__exp__2 = 0.42671326129048615
 mdl_cw__exp__2 = 0.77775351421422245
 mdl_ee__exp__2 = 9.4835522759998875E-002
 mdl_sw__exp__2 = 0.22224648578577769
 mdl_yb__exp__2 = 7.2875994928982540E-004
 mdl_yt__exp__2 = 0.98737240933884918
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_MU_R__exp__2 = 8315.2513440000002
 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
 mdl_G__exp__3 = 1.8056676068262196
 mdl_G__exp__4 = 2.1987899468922913
  Couplings of loop_sm
  ---------------------------------

   UVWfct_b_0 = -0.13642E+00 -0.00000E+00
       R2_DDA = 0.00000E+00 0.25704E-02
    R2_DDZ_V2 = 0.00000E+00 0.72127E-02
    R2_UUZ_V5 = 0.00000E+00 0.68702E-03
         GC_5 = 0.00000E+00 0.12177E+01
 UVWfct_b_0_1eps -0.18780E-01 0.00000E+00
         GC_1 = -0.00000E+00 -0.10265E+00
        GC_21 = -0.00000E+00 -0.28804E+00
        GC_23 = -0.00000E+00 -0.27437E-01
        GC_24 = 0.00000E+00 0.82310E-01
      GC_1003 = -0.00000E+00 -0.30795E+00

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

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

 Diagram information for clustering has been set-up for nFKSprocess
    1
 Diagram information for clustering has been set-up for nFKSprocess
    2
 getting user params
Enter number of events and iterations:
 Number of events and iterations -1 12
Enter desired fractional accuracy:
 Desired fractional accuracy: 2.9999999999999999E-002
 Enter alpha, beta for G_soft
   Enter alpha<0 to set G_soft=1 (no ME soft)
 for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001

 Enter alpha, beta for G_azi
   Enter alpha>0 to set G_azi=0 (no azi corr)
 for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001

 Doing the S and H events together
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Summing over 1 helicities/event for virt
Enter Configuration Number:
Running Configuration Number: 1
Enter running mode for MINT:
0 to set-up grids, 1 to integrate, 2 to generate events
 MINT running mode: 0
Set the three folding parameters for MINT
xi_i, phi_i, y_ij
           1 1 1
 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
 Enter 'born0' or 'virt0' to perform
  a pure n-body integration (no S functions)
 doing the all of this channel
 Normal integration (Sfunction != 1)
 Not subdividing B.W.
 about to integrate 5 -1 12 1
 imode is 0
 ------- iteration 1
 Update # PS points (even): 400 --> 384
Using random seed offsets: 1 , 1 , 0
  with seed 34
 Ranmar initialization seeds 13168 9409
 Total number of FKS directories is 2
 FKS process map (sum= 3 ) :
           1 --> 1 : 1
           2 --> 1 : 2
 ================================
 process combination map (specified per FKS dir):
  1 map 1
  1 inv. map 1
  2 map 1
  2 inv. map 1
 ================================
nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.52284E-06
 0.94000E+01 0.13000E+05
nFKSprocess: 1. Lower bound for tau is (taking resonances into account)
0.52284E-06 0.94000E+01 0.13000E+05
nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.52284E-06
 0.94000E+01 0.13000E+05
nFKSprocess: 2. Lower bound for tau is (taking resonances into account)
0.52284E-06 0.94000E+01 0.13000E+05
 bpower is 0.0000000000000000
 Scale values (may change event by event):
 muR, muR_reference: 0.580612D+04 0.580612D+04 1.00
 muF1, muF1_reference: 0.580612D+04 0.580612D+04 1.00
 muF2, muF2_reference: 0.580612D+04 0.580612D+04 1.00
 QES, QES_reference: 0.580612D+04 0.580612D+04 1.00

 muR_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF1_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 muF2_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state
 QES_reference [functional form]:
    H_T/2 := sum_i mT(i)/2, i=final state

 alpha_s= 7.2835541857615441E-002
 alpha_s value used for the virtuals is (for the first PS point):
 7.2835541857615441E-002
 ===============================================================
 INFO: MadLoop read these parameters from MadLoopParams.dat

 ===============================================================
  > CTModeRun = -1
  > MLStabThres = 1.0000000000000000E-003
  > NRotations_DP = 1
  > NRotations_QP = 0
  > CTStabThres = 1.0000000000000000E-002
  > CTLoopLibrary = 2
  > CTModeInit = 1
  > CheckCycle = 3
  > MaxAttempts = 10
  > UseLoopFilter = F
  > ImprovePSPoint = 2
  > DoubleCheckHelicityFilter = T
  > LoopInitStartOver = F
  > HelInitStartOver = F
  > ZeroThres = 1.0000000000000001E-009
 ===============================================================

------------------------------------------------------------------------
| You are using CutTools - Version 1.9.2 |
| Authors: G. Ossola, C. Papadopoulos, R. Pittau |
| Published in JHEP 0803:042,2008 |
| http://www.ugr.es/~pittau/CutTools |
| |
| Compiler with 34 significant digits detetected |
  ----------------------------------------------------------------------

########################################################################
# #
# You are using OneLOop-3.4 #
# #
# for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
# #
# author: Andreas van Hameren <email address hidden> #
# date: 02-01-2014 #
# #
# Please cite #
# A. van Hameren, #
# Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
# A. van Hameren, C.G. Papadopoulos and R. Pittau, #
# JHEP 0909:106,2009, arXiv:0903.4665 #
# in publications with results obtained with the help of this program. #
# #
########################################################################
 alpha_s value used for the virtuals is (for the first PS point):
 7.3761629254844654E-002
 POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
  COEFFICIENT DOUBLE POLE:
        MadFKS: 0.0000000000000000 OLP:
 0.0000000000000000
  COEFFICIENT SINGLE POLE:
        MadFKS: 1.2937557745598925E-003 OLP:
 1.2937871394031789E-003
  FINITE:
           OLP: 9.5925705431879789E-003

  MOMENTA (Exyzm):
           1 6500.0000000000000 0.0000000000000000
 0.0000000000000000 6500.0000000000000 0.0000000000000000
           2 6500.0000000000000 -0.0000000000000000
-0.0000000000000000 -6500.0000000000000 0.0000000000000000
           3 6500.0000000000000 -5045.8607560508080
-398.61297129080060 4078.0356594401046 4.7000000000000002
           4 6500.0000000000000 5045.8607560508080
 398.61297129080060 -4078.0356594401046 4.7000000000000002

  SCALE**2: 25619625.172964834
  Trying another PS point
 alpha_s value used for the virtuals is (for the first PS point):
 7.2361406150092192E-002
 POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
  COEFFICIENT DOUBLE POLE:
        MadFKS: 0.0000000000000000 OLP:
 0.0000000000000000
  COEFFICIENT SINGLE POLE:
        MadFKS: 2.4012150042950218E-003 OLP:
 2.4012732175481047E-003
  FINITE:
           OLP: 1.8806964942823662E-002

  MOMENTA (Exyzm):
           1 6500.0000000000000 0.0000000000000000
 0.0000000000000000 6500.0000000000000 0.0000000000000000
           2 6500.0000000000000 -0.0000000000000000
-0.0000000000000000 -6500.0000000000000 0.0000000000000000
           3 6500.0000000000000 4311.8298899974525
-4506.9266181820312 1829.1291283140260 4.7000000000000002
           4 6500.0000000000000 -4311.8298899974525
 4506.9266181820312 -1829.1291283140260 4.7000000000000002

  SCALE**2: 38904286.633089729
  Trying another PS point
 alpha_s value used for the virtuals is (for the first PS point):
 7.2140979913912648E-002
 POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
  COEFFICIENT DOUBLE POLE:
        MadFKS: 0.0000000000000000 OLP:
 0.0000000000000000
  COEFFICIENT SINGLE POLE:
        MadFKS: 3.2320728212418097E-003 OLP:
 3.2321511771877111E-003
  FINITE:
           OLP: 2.5531865799127904E-002

  MOMENTA (Exyzm):
           1 6500.0000000000000 0.0000000000000000
 0.0000000000000000 6500.0000000000000 0.0000000000000000
           2 6500.0000000000000 -0.0000000000000000
-0.0000000000000000 -6500.0000000000000 0.0000000000000000
           3 6500.0000000000000 1479.8299673302513
-6278.6237650420508 799.35273493027876 4.7000000000000002
           4 6500.0000000000000 -1479.8299673302513
 6278.6237650420508 -799.35273493027876 4.7000000000000002

  SCALE**2: 41611035.201112099
  Trying another PS point
 alpha_s value used for the virtuals is (for the first PS point):
 7.2168874147859094E-002
 POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
  COEFFICIENT DOUBLE POLE:
        MadFKS: 0.0000000000000000 OLP:
 0.0000000000000000
  COEFFICIENT SINGLE POLE:
        MadFKS: 3.0573482959829969E-003 OLP:
 3.0574224160387804E-003
  FINITE:
           OLP: 2.4125529468230239E-002

  MOMENTA (Exyzm):
           1 6500.0000000000000 0.0000000000000000
 0.0000000000000000 6500.0000000000000 0.0000000000000000
           2 6500.0000000000000 -0.0000000000000000
-0.0000000000000000 -6500.0000000000000 0.0000000000000000
           3 6500.0000000000000 3727.7785418867043
-5230.7795121781883 996.28848614542653 4.7000000000000002
           4 6500.0000000000000 -3727.7785418867043
 5230.7795121781883 -996.28848614542653 4.7000000000000002

  SCALE**2: 41257409.249448217
  Trying another PS point
 alpha_s value used for the virtuals is (for the first PS point):
 7.2108014426281147E-002
 POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
  COEFFICIENT DOUBLE POLE:
        MadFKS: 0.0000000000000000 OLP:
 0.0000000000000000
  COEFFICIENT SINGLE POLE:
        MadFKS: 3.5454333441789459E-003 OLP:
 3.5455192970020259E-003
  FINITE:
           OLP: 2.8043059503156151E-002

  MOMENTA (Exyzm):
           1 6500.0000000000000 0.0000000000000000
 0.0000000000000000 6500.0000000000000 0.0000000000000000
           2 6500.0000000000000 -0.0000000000000000
-0.0000000000000000 -6500.0000000000000 0.0000000000000000
           3 6500.0000000000000 5210.0000992210571
-3858.6397608945545 465.59217321620292 4.7000000000000002
           4 6500.0000000000000 -5210.0000992210571
 3858.6397608945545 -465.59217321620292 4.7000000000000002

  SCALE**2: 42033223.937885426
  Trying another PS point
 alpha_s value used for the virtuals is (for the first PS point):
 7.2775156870115043E-002
 POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
  COEFFICIENT DOUBLE POLE:
        MadFKS: 0.0000000000000000 OLP:
 0.0000000000000000
  COEFFICIENT SINGLE POLE:
        MadFKS: 1.7965130048913504E-003 OLP:
 1.7965565582033017E-003
  FINITE:
           OLP: 1.3845960589353980E-002

  MOMENTA (Exyzm):
           1 6500.0000000000000 0.0000000000000000
 0.0000000000000000 6500.0000000000000 0.0000000000000000
           2 6500.0000000000000 -0.0000000000000000
-0.0000000000000000 -6500.0000000000000 0.0000000000000000
           3 6500.0000000000000 5726.5722391184545
-1238.7804314620125 2814.5641479249725 4.7000000000000002
           4 6500.0000000000000 -5726.5722391184545
 1238.7804314620125 -2814.5641479249725 4.7000000000000002

  SCALE**2: 34328228.643745333
  TOO MANY FAILURES, QUITTING
Time in seconds: 0

On Tue, May 13, 2014 at 2:17 PM, marco zaro <
<email address hidden>> wrote:

> Your question #248615 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/248615
>
> Status: Open => Needs information
>
> marco zaro requested more information:
> Dear Konrad,
> I managed to have the same process running on my mac.
> Which center-of-mass energy are you running at?
> Can you please paste the content of the proc/run/param cards (attaching
> them won't work), and also of the log file inside P0_epem_bbx/GF1?
>
> That wopuld help to understand what is going wrong...
> Cheers,
>
> Marco
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https://answers.launchpad.net/mg5amcnlo/+question/248615
>
> You received this question notification because you asked the question.
>

Revision history for this message
marco zaro (marco-zaro) said : 		#3

indeed, it may be that you got instabilities for such a high center of mass energy.
Cheers,

Marco

On 13 May 2014, at 14:46, Konrad Tywoniuk <email address hidden> wrote:

> Question #248615 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/248615
>
> Status: Needs information => Open
>
> Konrad Tywoniuk gave more information on the question:
> Dear Marco,
>
> I was generating events at very high com energy, (6500 + 6500 GeV). When I
> reduced the energy to 45.8 + 45.8 GeV the program generated output.
>
> Thanks a lot for the rapid answer!
>
> Just for reference, here is the output of the file log_MINT0.txt (it seems
> the other log files in that directory were ok):
>
> ===============================================================
> INFO: MadFKS read these parameters from FKS_params.dat
>
> ===============================================================
>> IRPoleCheckThreshold = 1.0000000000000001E-005
>> PrecisionVirtualAtRunTime = 1.0000000000000000E-003
>> NHelForMCoverHels = 4
>> VirtualFraction = 1.0000000000000000
>> MinVirtualFraction = 5.0000000000000001E-003
> ===============================================================
> Process in group number 0
> No PDF is used, alpha_s(MZ) from param_card is used
> Value of alpha_s from param_card: 0.11799999999999999
> The default order of alpha_s running is fixed to 2
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> MU_R = 91.188000000000002
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MZ = 91.188000000000002
> mdl_MH = 125.00000000000000
> mdl_MTA = 1.7769999999999999
> mdl_WT = 1.4915000000000000
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WH = 6.3823389999999999E-003
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_CKM22 = 1.0000000000000000
> mdl_lhv = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM22 = 1.0000000000000000
> mdl_conjg__CKM33 = 1.0000000000000000
> mdl_Ncol = 3.0000000000000000
> mdl_CA = 3.0000000000000000
> mdl_TF = 0.50000000000000000
> mdl_CF = 1.3333333333333333
> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__2 = 15625.000000000000
> mdl_Ncol__exp__2_m_1 = 8.0000000000000000
> mdl_MB__exp__2 = 22.090000000000003
> mdl_MT__exp__2 = 29929.000000000000
> mdl_Ncol__exp__2_m_1_0 = 8.0000000000000000
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_v = 246.21845810181637
> mdl_v__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_AxialZUp = -0.18517701861793787
> mdl_AxialZDown = 0.18517701861793787
> mdl_VectorZUp = 7.5430507588273299E-002
> mdl_VectorZDown = -0.13030376310310560
> mdl_VectorAUp = 0.20530251149624587
> mdl_VectorADown = -0.10265125574812294
> mdl_VectorWmDxU = 0.23095271737156670
> mdl_AxialWmDxU = -0.23095271737156670
> mdl_VectorWpUxD = 0.23095271737156670
> mdl_AxialWpUxD = -0.23095271737156670
> mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_Vector_tbGp = (-0.96667059156454072 , 0.0000000000000000 )
> mdl_Axial_tbGp = ( -1.0206617000654716 , -0.0000000000000000 )
> mdl_Vector_tbGm = ( 0.96667059156454072 , 0.0000000000000000 )
> mdl_Axial_tbGm = ( -1.0206617000654716 , -0.0000000000000000 )
> mdl_gw__exp__2 = 0.42671326129048615
> mdl_cw__exp__2 = 0.77775351421422245
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_yb__exp__2 = 7.2875994928982540E-004
> mdl_yt__exp__2 = 0.98737240933884918
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_MU_R__exp__2 = 8315.2513440000002
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> mdl_G__exp__3 = 1.8056676068262196
> mdl_G__exp__4 = 2.1987899468922913
> Couplings of loop_sm
> ---------------------------------
>
> UVWfct_b_0 = -0.13642E+00 -0.00000E+00
> R2_DDA = 0.00000E+00 0.25704E-02
> R2_DDZ_V2 = 0.00000E+00 0.72127E-02
> R2_UUZ_V5 = 0.00000E+00 0.68702E-03
> GC_5 = 0.00000E+00 0.12177E+01
> UVWfct_b_0_1eps -0.18780E-01 0.00000E+00
> GC_1 = -0.00000E+00 -0.10265E+00
> GC_21 = -0.00000E+00 -0.28804E+00
> GC_23 = -0.00000E+00 -0.27437E-01
> GC_24 = 0.00000E+00 0.82310E-01
> GC_1003 = -0.00000E+00 -0.30795E+00
>
> Collider parameters:
> --------------------
>
> Running at e e machine @ 13000.000000000000 GeV
> PDF set = none
> alpha_s(Mz)= 0.1180 running at 2 loops.
> alpha_s(Mz)= 0.1180 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
>
>
> Diagram information for clustering has been set-up for nFKSprocess
> 1
> Diagram information for clustering has been set-up for nFKSprocess
> 2
> getting user params
> Enter number of events and iterations:
> Number of events and iterations -1 12
> Enter desired fractional accuracy:
> Desired fractional accuracy: 2.9999999999999999E-002
> Enter alpha, beta for G_soft
> Enter alpha<0 to set G_soft=1 (no ME soft)
> for G_soft: alpha= 1.0000000000000000 , beta= -0.10000000000000001
>
> Enter alpha, beta for G_azi
> Enter alpha>0 to set G_azi=0 (no azi corr)
> for G_azi: alpha= -1.0000000000000000 , beta= -0.10000000000000001
>
> Doing the S and H events together
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Summing over 1 helicities/event for virt
> Enter Configuration Number:
> Running Configuration Number: 1
> Enter running mode for MINT:
> 0 to set-up grids, 1 to integrate, 2 to generate events
> MINT running mode: 0
> Set the three folding parameters for MINT
> xi_i, phi_i, y_ij
> 1 1 1
> 'all ', 'born', 'real', 'virt', 'novi' or 'grid'?
> Enter 'born0' or 'virt0' to perform
> a pure n-body integration (no S functions)
> doing the all of this channel
> Normal integration (Sfunction != 1)
> Not subdividing B.W.
> about to integrate 5 -1 12 1
> imode is 0
> ------- iteration 1
> Update # PS points (even): 400 --> 384
> Using random seed offsets: 1 , 1 , 0
> with seed 34
> Ranmar initialization seeds 13168 9409
> Total number of FKS directories is 2
> FKS process map (sum= 3 ) :
> 1 --> 1 : 1
> 2 --> 1 : 2
> ================================
> process combination map (specified per FKS dir):
> 1 map 1
> 1 inv. map 1
> 2 map 1
> 2 inv. map 1
> ================================
> nFKSprocess: 1. Absolute lower bound for tau at the Born is 0.52284E-06
> 0.94000E+01 0.13000E+05
> nFKSprocess: 1. Lower bound for tau is (taking resonances into account)
> 0.52284E-06 0.94000E+01 0.13000E+05
> nFKSprocess: 2. Absolute lower bound for tau at the Born is 0.52284E-06
> 0.94000E+01 0.13000E+05
> nFKSprocess: 2. Lower bound for tau is (taking resonances into account)
> 0.52284E-06 0.94000E+01 0.13000E+05
> bpower is 0.0000000000000000
> Scale values (may change event by event):
> muR, muR_reference: 0.580612D+04 0.580612D+04 1.00
> muF1, muF1_reference: 0.580612D+04 0.580612D+04 1.00
> muF2, muF2_reference: 0.580612D+04 0.580612D+04 1.00
> QES, QES_reference: 0.580612D+04 0.580612D+04 1.00
>
> muR_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF1_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> muF2_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
> QES_reference [functional form]:
> H_T/2 := sum_i mT(i)/2, i=final state
>
> alpha_s= 7.2835541857615441E-002
> alpha_s value used for the virtuals is (for the first PS point):
> 7.2835541857615441E-002
> ===============================================================
> INFO: MadLoop read these parameters from MadLoopParams.dat
>
> ===============================================================
>> CTModeRun = -1
>> MLStabThres = 1.0000000000000000E-003
>> NRotations_DP = 1
>> NRotations_QP = 0
>> CTStabThres = 1.0000000000000000E-002
>> CTLoopLibrary = 2
>> CTModeInit = 1
>> CheckCycle = 3
>> MaxAttempts = 10
>> UseLoopFilter = F
>> ImprovePSPoint = 2
>> DoubleCheckHelicityFilter = T
>> LoopInitStartOver = F
>> HelInitStartOver = F
>> ZeroThres = 1.0000000000000001E-009
> ===============================================================
>
> ------------------------------------------------------------------------
> | You are using CutTools - Version 1.9.2 |
> | Authors: G. Ossola, C. Papadopoulos, R. Pittau |
> | Published in JHEP 0803:042,2008 |
> | http://www.ugr.es/~pittau/CutTools |
> | |
> | Compiler with 34 significant digits detetected |
> ----------------------------------------------------------------------
>
> ########################################################################
> # #
> # You are using OneLOop-3.4 #
> # #
> # for the evaluation of 1-loop scalar 1-, 2-, 3- and 4-point functions #
> # #
> # author: Andreas van Hameren <email address hidden> #
> # date: 02-01-2014 #
> # #
> # Please cite #
> # A. van Hameren, #
> # Comput.Phys.Commun. 182 (2011) 2427-2438, arXiv:1007.4716 #
> # A. van Hameren, C.G. Papadopoulos and R. Pittau, #
> # JHEP 0909:106,2009, arXiv:0903.4665 #
> # in publications with results obtained with the help of this program. #
> # #
> ########################################################################
> alpha_s value used for the virtuals is (for the first PS point):
> 7.3761629254844654E-002
> POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
> COEFFICIENT DOUBLE POLE:
> MadFKS: 0.0000000000000000 OLP:
> 0.0000000000000000
> COEFFICIENT SINGLE POLE:
> MadFKS: 1.2937557745598925E-003 OLP:
> 1.2937871394031789E-003
> FINITE:
> OLP: 9.5925705431879789E-003
>
> MOMENTA (Exyzm):
> 1 6500.0000000000000 0.0000000000000000
> 0.0000000000000000 6500.0000000000000 0.0000000000000000
> 2 6500.0000000000000 -0.0000000000000000
> -0.0000000000000000 -6500.0000000000000 0.0000000000000000
> 3 6500.0000000000000 -5045.8607560508080
> -398.61297129080060 4078.0356594401046 4.7000000000000002
> 4 6500.0000000000000 5045.8607560508080
> 398.61297129080060 -4078.0356594401046 4.7000000000000002
>
> SCALE**2: 25619625.172964834
> Trying another PS point
> alpha_s value used for the virtuals is (for the first PS point):
> 7.2361406150092192E-002
> POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
> COEFFICIENT DOUBLE POLE:
> MadFKS: 0.0000000000000000 OLP:
> 0.0000000000000000
> COEFFICIENT SINGLE POLE:
> MadFKS: 2.4012150042950218E-003 OLP:
> 2.4012732175481047E-003
> FINITE:
> OLP: 1.8806964942823662E-002
>
> MOMENTA (Exyzm):
> 1 6500.0000000000000 0.0000000000000000
> 0.0000000000000000 6500.0000000000000 0.0000000000000000
> 2 6500.0000000000000 -0.0000000000000000
> -0.0000000000000000 -6500.0000000000000 0.0000000000000000
> 3 6500.0000000000000 4311.8298899974525
> -4506.9266181820312 1829.1291283140260 4.7000000000000002
> 4 6500.0000000000000 -4311.8298899974525
> 4506.9266181820312 -1829.1291283140260 4.7000000000000002
>
> SCALE**2: 38904286.633089729
> Trying another PS point
> alpha_s value used for the virtuals is (for the first PS point):
> 7.2140979913912648E-002
> POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
> COEFFICIENT DOUBLE POLE:
> MadFKS: 0.0000000000000000 OLP:
> 0.0000000000000000
> COEFFICIENT SINGLE POLE:
> MadFKS: 3.2320728212418097E-003 OLP:
> 3.2321511771877111E-003
> FINITE:
> OLP: 2.5531865799127904E-002
>
> MOMENTA (Exyzm):
> 1 6500.0000000000000 0.0000000000000000
> 0.0000000000000000 6500.0000000000000 0.0000000000000000
> 2 6500.0000000000000 -0.0000000000000000
> -0.0000000000000000 -6500.0000000000000 0.0000000000000000
> 3 6500.0000000000000 1479.8299673302513
> -6278.6237650420508 799.35273493027876 4.7000000000000002
> 4 6500.0000000000000 -1479.8299673302513
> 6278.6237650420508 -799.35273493027876 4.7000000000000002
>
> SCALE**2: 41611035.201112099
> Trying another PS point
> alpha_s value used for the virtuals is (for the first PS point):
> 7.2168874147859094E-002
> POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
> COEFFICIENT DOUBLE POLE:
> MadFKS: 0.0000000000000000 OLP:
> 0.0000000000000000
> COEFFICIENT SINGLE POLE:
> MadFKS: 3.0573482959829969E-003 OLP:
> 3.0574224160387804E-003
> FINITE:
> OLP: 2.4125529468230239E-002
>
> MOMENTA (Exyzm):
> 1 6500.0000000000000 0.0000000000000000
> 0.0000000000000000 6500.0000000000000 0.0000000000000000
> 2 6500.0000000000000 -0.0000000000000000
> -0.0000000000000000 -6500.0000000000000 0.0000000000000000
> 3 6500.0000000000000 3727.7785418867043
> -5230.7795121781883 996.28848614542653 4.7000000000000002
> 4 6500.0000000000000 -3727.7785418867043
> 5230.7795121781883 -996.28848614542653 4.7000000000000002
>
> SCALE**2: 41257409.249448217
> Trying another PS point
> alpha_s value used for the virtuals is (for the first PS point):
> 7.2108014426281147E-002
> POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
> COEFFICIENT DOUBLE POLE:
> MadFKS: 0.0000000000000000 OLP:
> 0.0000000000000000
> COEFFICIENT SINGLE POLE:
> MadFKS: 3.5454333441789459E-003 OLP:
> 3.5455192970020259E-003
> FINITE:
> OLP: 2.8043059503156151E-002
>
> MOMENTA (Exyzm):
> 1 6500.0000000000000 0.0000000000000000
> 0.0000000000000000 6500.0000000000000 0.0000000000000000
> 2 6500.0000000000000 -0.0000000000000000
> -0.0000000000000000 -6500.0000000000000 0.0000000000000000
> 3 6500.0000000000000 5210.0000992210571
> -3858.6397608945545 465.59217321620292 4.7000000000000002
> 4 6500.0000000000000 -5210.0000992210571
> 3858.6397608945545 -465.59217321620292 4.7000000000000002
>
> SCALE**2: 42033223.937885426
> Trying another PS point
> alpha_s value used for the virtuals is (for the first PS point):
> 7.2775156870115043E-002
> POLES MISCANCELLATION, DIFFERENCE > 1.0000000000000003E-005
> COEFFICIENT DOUBLE POLE:
> MadFKS: 0.0000000000000000 OLP:
> 0.0000000000000000
> COEFFICIENT SINGLE POLE:
> MadFKS: 1.7965130048913504E-003 OLP:
> 1.7965565582033017E-003
> FINITE:
> OLP: 1.3845960589353980E-002
>
> MOMENTA (Exyzm):
> 1 6500.0000000000000 0.0000000000000000
> 0.0000000000000000 6500.0000000000000 0.0000000000000000
> 2 6500.0000000000000 -0.0000000000000000
> -0.0000000000000000 -6500.0000000000000 0.0000000000000000
> 3 6500.0000000000000 5726.5722391184545
> -1238.7804314620125 2814.5641479249725 4.7000000000000002
> 4 6500.0000000000000 -5726.5722391184545
> 1238.7804314620125 -2814.5641479249725 4.7000000000000002
>
> SCALE**2: 34328228.643745333
> TOO MANY FAILURES, QUITTING
> Time in seconds: 0
>
>
> On Tue, May 13, 2014 at 2:17 PM, marco zaro <
> <email address hidden>> wrote:
>
>> Your question #248615 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/248615
>>
>> Status: Open => Needs information
>>
>> marco zaro requested more information:
>> Dear Konrad,
>> I managed to have the same process running on my mac.
>> Which center-of-mass energy are you running at?
>> Can you please paste the content of the proc/run/param cards (attaching
>> them won't work), and also of the log file inside P0_epem_bbx/GF1?
>>
>> That wopuld help to understand what is going wrong...
>> Cheers,
>>
>> Marco
>>
>> --
>> To answer this request for more information, you can either reply to
>> this email or enter your reply at the following page:
>> https://answers.launchpad.net/mg5amcnlo/+question/248615
>>
>> You received this question notification because you asked the question.
>>
>
> --
> You received this question notification because you are a direct
> subscriber of the question.

Can you help with this problem?

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

To post a message you must log in.