MadGraph5_aMC@NLO

EPA with LHAPDF

Asked by Mamutjan

Dear experts,

I am trying to run two photon fusion process, say a a > w+ w-, with one beam type being EPA the other LHAPDF.
MadGraph is giving negative cross sections, something like -1.635e+12 +- 2.793e+09 pb.

I set lhapdf as pdlabel and chose the wanted lhaid then set lpp1=1,lpp2=2.
Then got the wrong answer. There is no error message, and it is running smoothly only with the wrong answer.

LHAPDF is working well when both beam types are PDF.

Could you help me with this?
Thanks.

Regards,
Mamut.

Question information

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

Did you use NLO PDF or LO ones?

Cheers,

Olivier
> On Apr 12, 2016, at 12:02, Mamutjan <email address hidden> wrote:
>
> New question #290726 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/290726
>
> Dear experts,
>
> I am trying to run two photon fusion process, say a a > w+ w-, with one beam type being EPA the other LHAPDF.
> MadGraph is giving negative cross sections, something like -1.635e+12 +- 2.793e+09 pb.
>
> I set lhapdf as pdlabel and chose the wanted lhaid then set lpp1=1,lpp2=2.
> Then got the wrong answer. There is no error message, and it is running smoothly only with the wrong answer.
>
> LHAPDF is working well when both beam types are PDF.
>
> Could you help me with this?
> Thanks.
>
> Regards,
> Mamut.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Mamutjan (mamutjan) said : 		#2

Dear Olivier

I have tried both MRST2004qed_proton(20463) and NNPDF23_nlo_as_0118_qed(244600), and didn't get the correct answer.

Regards
Mamut

At 2016-04-12 19:16:56, "Olivier Mattelaer" <email address hidden> wrote:
>Your question #290726 on MadGraph5_aMC@NLO changed:
>https://answers.launchpad.net/mg5amcnlo/+question/290726
>
> Status: Open => Answered
>
>Olivier Mattelaer proposed the following answer:
>Did you use NLO PDF or LO ones?
>
>Cheers,
>
>Olivier
>> On Apr 12, 2016, at 12:02, Mamutjan <email address hidden> wrote:
>>
>> New question #290726 on MadGraph5_aMC@NLO:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>
>> Dear experts,
>>
>> I am trying to run two photon fusion process, say a a > w+ w-, with one beam type being EPA the other LHAPDF.
>> MadGraph is giving negative cross sections, something like -1.635e+12 +- 2.793e+09 pb.
>>
>> I set lhapdf as pdlabel and chose the wanted lhaid then set lpp1=1,lpp2=2.
>> Then got the wrong answer. There is no error message, and it is running smoothly only with the wrong answer.
>>
>> LHAPDF is working well when both beam types are PDF.
>>
>> Could you help me with this?
>> Thanks.
>>
>> Regards,
>> Mamut.
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>
>--
>If this answers your question, please go to the following page to let us
>know that it is solved:
>https://answers.launchpad.net/mg5amcnlo/+question/290726/+confirm?answer_id=0
>
>If you still need help, you can reply to this email or go to the
>following page to enter your feedback:
>https://answers.launchpad.net/mg5amcnlo/+question/290726
>
>You received this question notification because you asked the question.

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

Hi,

I actually can not reproduce your result

I have run the following:
generate a a > w+ w-
output
launch
set lpp1 1
set lpp2 2
set pdlabel lhapdf
set lhaid 20463
set ptheavy 0
launch
set lhaid 244600

and get for
> MRST2004qed_proton -> 1.049 +- 0.002352 pb
 and for
> NNPDF23_nlo_as_0118_qed -> 0.8099 +- 0.001537 pb

I’m using the 2.3.4.beta version of the code.

Cheers,

Olivier

> On Apr 12, 2016, at 12:27, Mamutjan <email address hidden> wrote:
>
> Question #290726 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/290726
>
> Status: Answered => Open
>
> Mamutjan is still having a problem:
> Dear Olivier
>
> I have tried both MRST2004qed_proton(20463) and
> NNPDF23_nlo_as_0118_qed(244600), and didn't get the correct answer.
>
>
> Regards
> Mamut
>
>
>
>
> At 2016-04-12 19:16:56, "Olivier Mattelaer" <email address hidden> wrote:
>> Your question #290726 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Did you use NLO PDF or LO ones?
>>
>> Cheers,
>>
>> Olivier
>>> On Apr 12, 2016, at 12:02, Mamutjan <email address hidden> wrote:
>>>
>>> New question #290726 on MadGraph5_aMC@NLO:
>>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>>
>>> Dear experts,
>>>
>>> I am trying to run two photon fusion process, say a a > w+ w-, with one beam type being EPA the other LHAPDF.
>>> MadGraph is giving negative cross sections, something like -1.635e+12 +- 2.793e+09 pb.
>>>
>>> I set lhapdf as pdlabel and chose the wanted lhaid then set lpp1=1,lpp2=2.
>>> Then got the wrong answer. There is no error message, and it is running smoothly only with the wrong answer.
>>>
>>> LHAPDF is working well when both beam types are PDF.
>>>
>>> Could you help me with this?
>>> Thanks.
>>>
>>> Regards,
>>> Mamut.
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadGraph5_aMC@NLO.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726/+confirm?answer_id=0
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>
>> You received this question notification because you asked the question.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
abulikim2011 (abulikim2011) said : 		#4

Dear Olivier :

I have tried this process like you said , and I always have the negative Xcection value . But as long as both lpp set to 1 , it replays normal result .
I am using both the 2.3.3 beta version and 2.4.0 beta version .

best regards Nijat

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

Hi,

Did you run exactly the command that I have send?
Unfortunately, I will not be able to help you if I’m not able to reproduce this problem.

I have try on another computer (linux this time) and again, I have positive value.
I would suggest that you re-install the various tools (MG and lhapdf)

Cheers,

Olivier

> On Apr 13, 2016, at 07:27, abulikim2011 <email address hidden> wrote:
>
> Question #290726 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/290726
>
> abulikim2011 requested more information:
> Dear Olivier :
>
> I have tried this process like you said , and I always have the negative Xcection value . But as long as both lpp set to 1 , it replays normal result .
> I am using both the 2.3.3 beta version and 2.4.0 beta version .
>
> best regards Nijat
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Mamutjan (mamutjan) said : 		#6

Hi,

Thank you for your time.
I have reinstalled the LHAPDF6.1.6 and MadGraphv2_3_3 (I had tried MG5_aMC_v2.4.0.beta but had too many compile errors at launch command) and still my cross section is negative for a a > w+ w-. My OS is 14.04.1-Ubuntu x64.

However, I noticed that in output_folder/SubProcesses/P0_aa_wpwm/G1,1a0,1b0,.... there are run_01_app.log files read:

==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.1.6 loading /usr/local/share/LHAPDF/NNPDF23_nlo_as_0118_qed/NNPDF23_nlo_as_0118_qed_0000.dat
NNPDF23_nlo_as_0118_qed PDF set, member #0, version 4; LHAPDF ID = 244600
 Process in group number 0
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11799999999999999
 New value of alpha_s from PDF lhapdf : 0.11800232290472297
 using LHAPDF
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

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

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

 mdl_conjg__CKM3x3 = 1.0000000000000000
 mdl_CKM3x3 = 1.0000000000000000
 mdl_conjg__CKM1x1 = 1.0000000000000000
 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_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_vev = 246.21845810181637
 mdl_vev__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_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_ee__exp__2 = 9.4835522759998875E-002
 mdl_sw__exp__2 = 0.22224648578577769
 mdl_cw__exp__2 = 0.77775351421422245
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34351128074635334
 mdl_G__exp__2 = 1.4828317324943823
  Couplings of sm
  ---------------------------------

         GC_4 = 0.00000E+00 0.30795E+00
         GC_5 = 0.00000E+00 0.94836E-01

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

 Running at P a2 machine @ 13000.000000000000 GeV
 PDF set = lhapdf
 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

 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations 5000 9 3
Enter desired fractional accuracy:
 Desired fractional accuracy: -6.6925445054199998E-010
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 4
 Completed mapping 4
 about to integrate 4 5000 9 3 4 1
 Using non-zero grid deformation.
  4 dimensions 5000 events 4 invarients 9 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 fail to read data
 Error opening grid
 Using Uniform Grid! 16
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Set CM energy to 13000.00
 Mapping Graph 1 to config 1
Setting grid 1 0.38268E-04 1
 Transforming s_hat 1/s 3 1.5307020010344466E-004
 Masses: 0.000E+00 0.000E+00 0.804E+02 0.804E+02
   1 1 2 3 4
 Got moffset 2
Using random seed offsets 1 : 1
  with seed 27
 Ranmar initialization seeds 4635 15716

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

  Particle 3 4
      Et > 0.0 0.0
       E > 0.0 0.0
     Eta < -1.0 -1.0
   xqcut: 0.0 0.0
d R # 3 > -0.0 0.0
s min # 3> 0.0 0.0
xqcutij # 3> 0.0 0.0
 alpha_s for scale 94.377721023447251 is 0.11739861780502350
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 Added good helicity 1 5.1053409672572804E-006 in event 1
 Added good helicity 3 5.4781823738226461 in event 1
 Added good helicity 5 0.52500601144927206 in event 1
 Added good helicity 7 5.0341884841130294E-002 in event 1
 Added good helicity 9 5.1053409672572906E-006 in event 1
 Added good helicity 10 1.7115852508710567 in event 1
 Added good helicity 11 1.8575839065178188 in event 1
 Added good helicity 12 0.50400662717481126 in event 1
 Added good helicity 13 1.8575839065178188 in event 1
 Added good helicity 14 2.8669677293183371 in event 1
 Added good helicity 15 1.0492954964994250 in event 1
 Added good helicity 16 0.50400662717481171 in event 1
 Added good helicity 17 1.0492954964994252 in event 1
 Added good helicity 18 0.54613155896273602 in event 1
 Added good helicity 19 0.54613155896273613 in event 1
 Added good helicity 20 1.0492954964994252 in event 1
 Added good helicity 21 0.50400662717481193 in event 1
 Added good helicity 22 1.0492954964994250 in event 1
 Added good helicity 23 2.8669677293183375 in event 1
 Added good helicity 24 1.8575839065178190 in event 1
 Added good helicity 25 0.50400662717481115 in event 1
 Added good helicity 26 1.8575839065178190 in event 1
 Added good helicity 27 1.7115852508710574 in event 1
 Added good helicity 28 5.1053409672572788E-006 in event 1
 Added good helicity 30 5.0341884841130356E-002 in event 1
 Added good helicity 32 0.52500601144927295 in event 1
 Added good helicity 34 5.4781823738226461 in event 1
 Added good helicity 36 5.1053409672572754E-006 in event 1
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 Added good helicity 2 1.5433464669709732E-006 in event 4
 Added good helicity 4 1.0053686447502854E-006 in event 4
 Added good helicity 6 1.5433464669709063E-006 in event 4
 Added good helicity 8 1.0053686447503102E-006 in event 4
 Added good helicity 29 1.0053686447502651E-006 in event 4
 Added good helicity 31 1.5433464669709747E-006 in event 4
 Added good helicity 33 1.0053686447503004E-006 in event 4
 Added good helicity 35 1.5433464669708341E-006 in event 4
 beam type not supported in lhadpf
 beam type not supported in lhadpf
 beam type not supported in lhadpf
(millions of same line)
beam type not supported in lhadpf
 beam type not supported in lhadpf
 Iteration 9 Mean: -.1654E+13 Abs mean: 0.1654E+13 Fluctuation: 0.140E+11 0.145E+17 87.0%
  9 -0.1654E+13 0.1654E+13 +- 0.1398E+11 9.56
 Exceeded boost 1 94 213333.328
 Exceeded boost 1 95 213333.328
 Exceeded boost 1 96 1280000.00
 Relative summed weights:
         NaN NaN
         NaN NaN
 Relative number of events:
         NaN NaN
         NaN NaN
 Events:
           0 0
           0 0
 Found 128898 events.
 Wrote 1588 events.
 Actual xsec -1643598755412.9604
 Correct abs xsec 1643598755412.9604
 Event xsec 1643598755412.8787
 Events wgts > 1: 7
 % Cross section > 1: 15988371446.826084 0.97276609599340691
 Checking number of events 1099.9857819654158 1588 chi2: 1.4054530193604151

 -------------------------------------------------------------------------------
 Accumulated results: Integral = 0.1640E+13
                        Std dev = 0.9028E+10
                       Cross sec = -0.2344E+04
             Chi**2 per DoF. =
-------------------------------------------------
---------------------------
 Results Last 3 iters: Integral = -0.1644E+13
                     Abs integral = 0.1644E+13
                          Std dev = 0.9987E+10
                  Chi**2 per DoF. = 1.4055
-------------------------------------------------
---------------------------
 Status -6.6925445054199998E-010 10 9
 We found enough events 1588 1099.9857819654158

ls status:
events.lhe
ftn25
ftn26
input_sg.txt
moffset.dat
results.dat
run_01_app.log

==========================LOG FILE ENDS HERE==============================

I think the program is not working correctly when pdflabel is set to LHAPDF and lpp is set to EPA.

Regards
Mamut

At 2016-04-14 19:02:58, "Olivier Mattelaer" <email address hidden> wrote:
>Your question #290726 on MadGraph5_aMC@NLO changed:
>https://answers.launchpad.net/mg5amcnlo/+question/290726
>
>Olivier Mattelaer proposed the following answer:
>Hi,
>
>Did you run exactly the command that I have send?
>Unfortunately, I will not be able to help you if I’m not able to reproduce this problem.
>
>I have try on another computer (linux this time) and again, I have positive value.
>I would suggest that you re-install the various tools (MG and lhapdf)
>
>
>Cheers,
>
>Olivier
>
>
>> On Apr 13, 2016, at 07:27, abulikim2011 <email address hidden> wrote:
>>
>> Question #290726 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>
>> abulikim2011 requested more information:
>> Dear Olivier :
>>
>> I have tried this process like you said , and I always have the negative Xcection value . But as long as both lpp set to 1 , it replays normal result .
>> I am using both the 2.3.3 beta version and 2.4.0 beta version .
>>
>> best regards Nijat
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>
>--
>If this answers your question, please go to the following page to let us
>know that it is solved:
>https://answers.launchpad.net/mg5amcnlo/+question/290726/+confirm?answer_id=4
>
>If you still need help, you can reply to this email or go to the
>following page to enter your feedback:
>https://answers.launchpad.net/mg5amcnlo/+question/290726
>
>You received this question notification because you asked the question.

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

Ok,

This clarify why we do not have the same answer and the problem.
LHAPDF is indeed incompatible with EPA.
and the behavior differs if you use LHAPDF6 or LHAPDF5.

I have put in the development version a fix which I have tested with lhapdf6 (I do not have lhapdf5 anymore).
you can get it via the following command:
bzr branch lp:~maddevelopers/mg5amcnlo/2.4.1/

Could you do us a favor and create a bug report with your compilation error?
It is important for us to know what happens before 2.4.0 will be officially release.

Thanks,

Olivier

> On Apr 16, 2016, at 16:51, Mamutjan <email address hidden> wrote:
>
> Question #290726 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/290726
>
> Status: Answered => Open
>
> Mamutjan is still having a problem:
> Hi,
>
> Thank you for your time.
> I have reinstalled the LHAPDF6.1.6 and MadGraphv2_3_3 (I had tried MG5_aMC_v2.4.0.beta but had too many compile errors at launch command) and still my cross section is negative for a a > w+ w-. My OS is 14.04.1-Ubuntu x64.
>
> However, I noticed that in
> output_folder/SubProcesses/P0_aa_wpwm/G1,1a0,1b0,.... there are
> run_01_app.log files read:
>
> ==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
> LHAPDF 6.1.6 loading /usr/local/share/LHAPDF/NNPDF23_nlo_as_0118_qed/NNPDF23_nlo_as_0118_qed_0000.dat
> NNPDF23_nlo_as_0118_qed PDF set, member #0, version 4; LHAPDF ID = 244600
> Process in group number 0
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11799999999999999
> New value of alpha_s from PDF lhapdf : 0.11800232290472297
> using LHAPDF
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MZ = 91.188000000000002
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MH = 125.00000000000000
> mdl_MTA = 1.7769999999999999
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WT = 1.4915000000000000
> mdl_WH = 6.3823389999999999E-003
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM1x1 = 1.0000000000000000
> 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_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_vev = 246.21845810181637
> mdl_vev__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_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_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_cw__exp__2 = 0.77775351421422245
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> Couplings of sm
> ---------------------------------
>
> GC_4 = 0.00000E+00 0.30795E+00
> GC_5 = 0.00000E+00 0.94836E-01
>
> Collider parameters:
> --------------------
>
> Running at P a2 machine @ 13000.000000000000 GeV
> PDF set = lhapdf
> 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
>
>
> getting user params
> Enter number of events and max and min iterations:
> Number of events and iterations 5000 9 3
> Enter desired fractional accuracy:
> Desired fractional accuracy: -6.6925445054199998E-010
> 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 4
> Completed mapping 4
> about to integrate 4 5000 9 3 4 1
> Using non-zero grid deformation.
> 4 dimensions 5000 events 4 invarients 9 iterations 1 config(s), (0.99)
> Using h-tuple random number sequence.
> fail to read data
> Error opening grid
> Using Uniform Grid! 16
> Using uniform alpha 1.0000000000000000
> Grid defined OK
> Set CM energy to 13000.00
> Mapping Graph 1 to config 1
> Setting grid 1 0.38268E-04 1
> Transforming s_hat 1/s 3 1.5307020010344466E-004
> Masses: 0.000E+00 0.000E+00 0.804E+02 0.804E+02
> 1 1 2 3 4
> Got moffset 2
> Using random seed offsets 1 : 1
> with seed 27
> Ranmar initialization seeds 4635 15716
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Particle 3 4
> Et > 0.0 0.0
> E > 0.0 0.0
> Eta < -1.0 -1.0
> xqcut: 0.0 0.0
> d R # 3 > -0.0 0.0
> s min # 3> 0.0 0.0
> xqcutij # 3> 0.0 0.0
> alpha_s for scale 94.377721023447251 is 0.11739861780502350
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> Added good helicity 1 5.1053409672572804E-006 in event 1
> Added good helicity 3 5.4781823738226461 in event 1
> Added good helicity 5 0.52500601144927206 in event 1
> Added good helicity 7 5.0341884841130294E-002 in event 1
> Added good helicity 9 5.1053409672572906E-006 in event 1
> Added good helicity 10 1.7115852508710567 in event 1
> Added good helicity 11 1.8575839065178188 in event 1
> Added good helicity 12 0.50400662717481126 in event 1
> Added good helicity 13 1.8575839065178188 in event 1
> Added good helicity 14 2.8669677293183371 in event 1
> Added good helicity 15 1.0492954964994250 in event 1
> Added good helicity 16 0.50400662717481171 in event 1
> Added good helicity 17 1.0492954964994252 in event 1
> Added good helicity 18 0.54613155896273602 in event 1
> Added good helicity 19 0.54613155896273613 in event 1
> Added good helicity 20 1.0492954964994252 in event 1
> Added good helicity 21 0.50400662717481193 in event 1
> Added good helicity 22 1.0492954964994250 in event 1
> Added good helicity 23 2.8669677293183375 in event 1
> Added good helicity 24 1.8575839065178190 in event 1
> Added good helicity 25 0.50400662717481115 in event 1
> Added good helicity 26 1.8575839065178190 in event 1
> Added good helicity 27 1.7115852508710574 in event 1
> Added good helicity 28 5.1053409672572788E-006 in event 1
> Added good helicity 30 5.0341884841130356E-002 in event 1
> Added good helicity 32 0.52500601144927295 in event 1
> Added good helicity 34 5.4781823738226461 in event 1
> Added good helicity 36 5.1053409672572754E-006 in event 1
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> Added good helicity 2 1.5433464669709732E-006 in event 4
> Added good helicity 4 1.0053686447502854E-006 in event 4
> Added good helicity 6 1.5433464669709063E-006 in event 4
> Added good helicity 8 1.0053686447503102E-006 in event 4
> Added good helicity 29 1.0053686447502651E-006 in event 4
> Added good helicity 31 1.5433464669709747E-006 in event 4
> Added good helicity 33 1.0053686447503004E-006 in event 4
> Added good helicity 35 1.5433464669708341E-006 in event 4
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> (millions of same line)
> beam type not supported in lhadpf
> beam type not supported in lhadpf
> Iteration 9 Mean: -.1654E+13 Abs mean: 0.1654E+13 Fluctuation: 0.140E+11 0.145E+17 87.0%
> 9 -0.1654E+13 0.1654E+13 +- 0.1398E+11 9.56
> Exceeded boost 1 94 213333.328
> Exceeded boost 1 95 213333.328
> Exceeded boost 1 96 1280000.00
> Relative summed weights:
> NaN NaN
> NaN NaN
> Relative number of events:
> NaN NaN
> NaN NaN
> Events:
> 0 0
> 0 0
> Found 128898 events.
> Wrote 1588 events.
> Actual xsec -1643598755412.9604
> Correct abs xsec 1643598755412.9604
> Event xsec 1643598755412.8787
> Events wgts > 1: 7
> % Cross section > 1: 15988371446.826084 0.97276609599340691
> Checking number of events 1099.9857819654158 1588 chi2: 1.4054530193604151
>
> -------------------------------------------------------------------------------
> Accumulated results: Integral = 0.1640E+13
> Std dev = 0.9028E+10
> Cross sec = -0.2344E+04
> Chi**2 per DoF. =
> -------------------------------------------------
> ---------------------------
> Results Last 3 iters: Integral = -0.1644E+13
> Abs integral = 0.1644E+13
> Std dev = 0.9987E+10
> Chi**2 per DoF. = 1.4055
> -------------------------------------------------
> ---------------------------
> Status -6.6925445054199998E-010 10 9
> We found enough events 1588 1099.9857819654158
>
> ls status:
> events.lhe
> ftn25
> ftn26
> input_sg.txt
> moffset.dat
> results.dat
> run_01_app.log
>
> ==========================LOG FILE ENDS
> HERE==============================
>
> I think the program is not working correctly when pdflabel is set to
> LHAPDF and lpp is set to EPA.
>
>
> Regards
> Mamut
>
>
>
>
>
>
> At 2016-04-14 19:02:58, "Olivier Mattelaer" <email address hidden> wrote:
>> Your question #290726 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>
>> Olivier Mattelaer proposed the following answer:
>> Hi,
>>
>> Did you run exactly the command that I have send?
>> Unfortunately, I will not be able to help you if I’m not able to reproduce this problem.
>>
>> I have try on another computer (linux this time) and again, I have positive value.
>> I would suggest that you re-install the various tools (MG and lhapdf)
>>
>>
>> Cheers,
>>
>> Olivier
>>
>>
>>> On Apr 13, 2016, at 07:27, abulikim2011 <email address hidden> wrote:
>>>
>>> Question #290726 on MadGraph5_aMC@NLO changed:
>>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>>
>>> abulikim2011 requested more information:
>>> Dear Olivier :
>>>
>>> I have tried this process like you said , and I always have the negative Xcection value . But as long as both lpp set to 1 , it replays normal result .
>>> I am using both the 2.3.3 beta version and 2.4.0 beta version .
>>>
>>> best regards Nijat
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadGraph5_aMC@NLO.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726/+confirm?answer_id=4
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/mg5amcnlo/+question/290726
>>
>> You received this question notification because you asked the question.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Mamutjan (mamutjan) said : 		#8

Hi,

Thank you very much.

I have reported the error messages I had using version 2.4.1.

Mamut

Can you help with this problem?

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

To post a message you must log in.