calculate differential cross section

Asked by abulikim2011

Hello
 I want to calculate the differential cross section of this process that p p > e- ve~ c [QCD] , in this calculation , i used the CT10nlo PDF , and make figure , but this figure is not same this ATLAS paper's figure.(ArXiv ePrint: 1402.6263 figure 13). on this calculation i applied same kinematic cuts like experimenters. this bis size same too, but this calculation gives me different results.
this kinematic cuts is :
p t > 20 GeV,
eta_l< 2.5,
pT_nu> 25 GeV,
mw_t > 40 GeV
for lepton
p T > 25 GeV,
 |η| < 2.5
R=0.4
for final jet
there some cuts that not able to modify on run_card.dat file , i modified cuts.f file.
and than i calculated the fix-order. the experimental results is:
19.86
,18.67
,18.38
,18.1
,16.0
,14.91
,14.4
,11.59
,9.43
,7.74
,6.75
my result is
1.8716418e+01
1.8448775e+01
1.7985600e+01
1.7212763e+01
1.6320482e+01
1.4880756e+01
1.3413127e+01
1.1852857e+01
1.0066175e+01
9.0576151e+00
6.6652629e+00

why this results is different. how i can same results like experimentalist result
thanks .

Question information

Language:
English Edit question
Status:
Expired
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

Hi,

What is the scale/pdf uncertainty on each number of that computation?
Have you check that you use the same scale/pdf as them?

Cheers,

Olivier

> On Nov 15, 2017, at 16:03, abulikim2011 <email address hidden> wrote:
>
> New question #660790 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/660790
>
> Hello
> I want to calculate the differential cross section of this process that p p > e- ve~ c [QCD] , in this calculation , i used the CT10nlo PDF , and make figure , but this figure is not same this ATLAS paper's figure.(ArXiv ePrint: 1402.6263 figure 13). on this calculation i applied same kinematic cuts like experimenters. this bis size same too, but this calculation gives me different results.
> this kinematic cuts is :
> p t > 20 GeV,
> eta_l< 2.5,
> pT_nu> 25 GeV,
> mw_t > 40 GeV
> for lepton
> p T > 25 GeV,
> |η| < 2.5
> R=0.4
> for final jet
> there some cuts that not able to modify on run_card.dat file , i modified cuts.f file.
> and than i calculated the fix-order. the experimental results is:
> 19.86
> ,18.67
> ,18.38
> ,18.1
> ,16.0
> ,14.91
> ,14.4
> ,11.59
> ,9.43
> ,7.74
> ,6.75
> my result is
> 1.8716418e+01
> 1.8448775e+01
> 1.7985600e+01
> 1.7212763e+01
> 1.6320482e+01
> 1.4880756e+01
> 1.3413127e+01
> 1.1852857e+01
> 1.0066175e+01
> 9.0576151e+00
> 6.6652629e+00
>
> why this results is different. how i can same results like experimentalist result
> thanks .
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

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

in this paper ,The dependence of the NLO prediction on the choice of renormalisation and factorisa-
tion scales is evaluated by independently halving and doubling their nominal value which
is chosen as the sum of the transverse mass of all final-state particles.so in run_card.dat file , i modify the dynamical scale , which equal to 2 . is it same?

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

Hi,

Their plot is actually NLO+PS accuracy, while your plot is NLO accuracy.
So you can indeed have difference in shape between those two precision (difference which can also depend of the parton-shower program that you use)
Additionally to that, the cut applied after parton-shower are not exactly the same as the cut before the parton-shower making even more plausible that you find such difference.

Cheers,

Olivier

> On Nov 16, 2017, at 01:38, abulikim2011 <email address hidden> wrote:
>
> Question #660790 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/660790
>
> Status: Answered => Open
>
> abulikim2011 is still having a problem:
> in this paper ,The dependence of the NLO prediction on the choice of renormalisation and factorisa-
> tion scales is evaluated by independently halving and doubling their nominal value which
> is chosen as the sum of the transverse mass of all final-state particles.so in run_card.dat file , i modify the dynamical scale , which equal to 2 . is it same?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

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

thanks your reply Olivier,
the Madgraph version which i used this process is MG5_aMC_v2_5_0_beta2, in this version i installed the pythia8 using install pythia8, and i applied those cuts and calculated, but final step there appear some error massage like not found -ld ....,
and i can't fixed this error , can you tell me how to fixed this ?

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

Please follow the suggestion of the following question:
https://answers.launchpad.net/mg5amcnlo/+question/659847

Cheers,

Olivier
> On Nov 16, 2017, at 11:44, abulikim2011 <email address hidden> wrote:
>
> Question #660790 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/660790
>
> Status: Answered => Open
>
> abulikim2011 is still having a problem:
> thanks your reply Olivier,
> the Madgraph version which i used this process is MG5_aMC_v2_5_0_beta2, in this version i installed the pythia8 using install pythia8, and i applied those cuts and calculated, but final step there appear some error massage like not found -ld ....,
> and i can't fixed this error , can you tell me how to fixed this ?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

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

thanks your help Olivier,
this problem is sloven , but i modify the run_card.dat shower=PYTHIA8 , and run it again, but there appear another error massage like that:
NFO: Idle: 0, Running: 0, Completed: 0 [ current time: 20h25 ]
Command "launch " interrupted with error:
IOError : [Errno 2] No such file or directory: '/home/ablikim/tools/MG5_aMC_v2_5_0_beta2_2/bin/PROCNLO_loop_sm_0/MCatNLO/RUN_PYTHIA8_1/events.hepmc.gz'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/ablikim/tools/MG5_aMC_v2_5_0_beta2_2/bin/PROCNLO_loop_sm_0/run_02_tag_1_debug.log'.
Please attach this file to your report.

what is happen here ?
thanks

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

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