From 537a0488d64da33231e88e39b31ff395f8f33016 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Fri, 13 May 2022 14:44:58 +0200 Subject: [PATCH 01/11] Add CDF transverse mass runcards --- .../CDF_WP_1960GEV_88FB_MT/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT/launch.txt | 34 ++++++++++ .../CDF_WP_1960GEV_88FB_MT/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT/output.txt | 7 +++ 4 files changed, 113 insertions(+) create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f new file mode 100644 index 00000000..3e85a152 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 80, 60d0, 100d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt new file mode 100644 index 00000000..2a605674 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt @@ -0,0 +1,34 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw @MW@ +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 60.0 +#user_defined_cut set mtwmax = 100.0 +set rphreco = 0.2 +set req_acc_FO 0.0005 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt new file mode 100644 index 00000000..c48507fe --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD] +output @OUTPUT@ +quit From 1895c0ff44952dcf92f629196802d4778a803fa4 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Fri, 20 May 2022 14:36:46 +0200 Subject: [PATCH 02/11] Add missing routines for custom cuts --- .../external/mg5/cuts_code/mtwmax.f | 22 +++++++++++++++++++ .../external/mg5/cuts_code/mtwmin.f | 22 +++++++++++++++++++ .../external/mg5/cuts_code/ptlmax.f | 11 ++++++++++ .../external/mg5/cuts_code/ptmissmax.f | 19 ++++++++++++++++ .../external/mg5/cuts_code/ptmissmin.f | 19 ++++++++++++++++ 5 files changed, 93 insertions(+) create mode 100644 runcardsrunner/external/mg5/cuts_code/mtwmax.f create mode 100644 runcardsrunner/external/mg5/cuts_code/mtwmin.f create mode 100644 runcardsrunner/external/mg5/cuts_code/ptlmax.f create mode 100644 runcardsrunner/external/mg5/cuts_code/ptmissmax.f create mode 100644 runcardsrunner/external/mg5/cuts_code/ptmissmin.f diff --git a/runcardsrunner/external/mg5/cuts_code/mtwmax.f b/runcardsrunner/external/mg5/cuts_code/mtwmax.f new file mode 100644 index 00000000..a43b5cdf --- /dev/null +++ b/runcardsrunner/external/mg5/cuts_code/mtwmax.f @@ -0,0 +1,22 @@ +c maximum transverse mass cut of W bosons + block + real*8 xmtw2 + + do i=3,nexternal + do j=i+1,nexternal + if (is_a_lm(i) .or. is_a_lp(i) .or. + & is_a_lm(j) .or. is_a_lp(j)) then + if (abs(ipdg_reco(i)+ipdg_reco(j)).eq.1) then + xmtw2=2d0*sqrt((p_reco(1,i)**2+p_reco(2,i)**2)* + & (p_reco(1,j)**2+p_reco(2,j)**2))- + & 2d0*(p_reco(1,i)*p_reco(1,j)+ + & p_reco(2,i)*p_reco(2,j)) + if (xmtw2 > (({})**2)) then + passcuts_user=.false. + return + endif + endif + endif + enddo + enddo + end block diff --git a/runcardsrunner/external/mg5/cuts_code/mtwmin.f b/runcardsrunner/external/mg5/cuts_code/mtwmin.f new file mode 100644 index 00000000..0c3c5269 --- /dev/null +++ b/runcardsrunner/external/mg5/cuts_code/mtwmin.f @@ -0,0 +1,22 @@ +c minimum transverse mass cut of W bosons + block + real*8 xmtw2 + + do i=3,nexternal + do j=i+1,nexternal + if (is_a_lm(i) .or. is_a_lp(i) .or. + & is_a_lm(j) .or. is_a_lp(j)) then + if (abs(ipdg_reco(i)+ipdg_reco(j)).eq.1) then + xmtw2=2d0*sqrt((p_reco(1,i)**2+p_reco(2,i)**2)* + & (p_reco(1,j)**2+p_reco(2,j)**2))- + & 2d0*(p_reco(1,i)*p_reco(1,j)+ + & p_reco(2,i)*p_reco(2,j)) + if (xmtw2 < (({})**2)) then + passcuts_user=.false. + return + endif + endif + endif + enddo + enddo + end block diff --git a/runcardsrunner/external/mg5/cuts_code/ptlmax.f b/runcardsrunner/external/mg5/cuts_code/ptlmax.f new file mode 100644 index 00000000..b337831f --- /dev/null +++ b/runcardsrunner/external/mg5/cuts_code/ptlmax.f @@ -0,0 +1,11 @@ +c maximum transverse momentum for all charged leptons + do i=3,nexternal + if (abs(ipdg_reco(i)) == 12 .or. + & abs(ipdg_reco(i)) == 14 .or. + & abs(ipdg_reco(i)) == 16) then + if (p_reco(1,i)**2+p_reco(2,i)**2 > (({})**2)) then + passcuts_user=.false. + return + endif + endif + enddo diff --git a/runcardsrunner/external/mg5/cuts_code/ptmissmax.f b/runcardsrunner/external/mg5/cuts_code/ptmissmax.f new file mode 100644 index 00000000..b16c7bee --- /dev/null +++ b/runcardsrunner/external/mg5/cuts_code/ptmissmax.f @@ -0,0 +1,19 @@ +c minimum missing transverse momentum + block + real*8 xptmiss(1:2) + xptmiss=0d0 + + do i=3,nexternal + if (abs(ipdg_reco(i))==12 .or. + & abs(ipdg_reco(i))==14 .or. + & abs(ipdg_reco(i))==16) then + xptmiss(1)=xptmiss(1)+p_reco(1,i) + xptmiss(2)=xptmiss(2)+p_reco(2,i) + endif + enddo + + if (xptmiss(1)**2+xptmiss(2)**2 > (({})**2)) then + passcuts_user=.false. + return + endif + end block diff --git a/runcardsrunner/external/mg5/cuts_code/ptmissmin.f b/runcardsrunner/external/mg5/cuts_code/ptmissmin.f new file mode 100644 index 00000000..97f59022 --- /dev/null +++ b/runcardsrunner/external/mg5/cuts_code/ptmissmin.f @@ -0,0 +1,19 @@ +c minimum missing transverse momentum + block + real*8 xptmiss(1:2) + xptmiss=0d0 + + do i=3,nexternal + if (abs(ipdg_reco(i))==12 .or. + & abs(ipdg_reco(i))==14 .or. + & abs(ipdg_reco(i))==16) then + xptmiss(1)=xptmiss(1)+p_reco(1,i) + xptmiss(2)=xptmiss(2)+p_reco(2,i) + endif + enddo + + if (xptmiss(1)**2+xptmiss(2)**2 < (({})**2)) then + passcuts_user=.false. + return + endif + end block From 92bd1e615d318b42bd5446482650f096031cfe70 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Fri, 20 May 2022 15:42:34 +0200 Subject: [PATCH 03/11] Add missing cut --- nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt | 1 + runcardsrunner/external/mg5/cuts_code/ptwmax.f | 12 ++++++++++++ 2 files changed, 13 insertions(+) create mode 100644 runcardsrunner/external/mg5/cuts_code/ptwmax.f diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt index 2a605674..efce994c 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt @@ -27,6 +27,7 @@ set etal = 1.0 #user_defined_cut set ptmissmax = 55.0 #user_defined_cut set mtwmin = 60.0 #user_defined_cut set mtwmax = 100.0 +#user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 set req_acc_FO 0.0005 set pineappl True diff --git a/runcardsrunner/external/mg5/cuts_code/ptwmax.f b/runcardsrunner/external/mg5/cuts_code/ptwmax.f new file mode 100644 index 00000000..01f17535 --- /dev/null +++ b/runcardsrunner/external/mg5/cuts_code/ptwmax.f @@ -0,0 +1,12 @@ +c maximum transverse mass cut of W bosons + do i=3,nexternal + do j=i+1,nexternal + if (abs(ipdg_reco(i)+ipdg_reco(j)).eq.1) then + if (((p_reco(1,i)+p_reco(1,j))**2+ + & (p_reco(2,i)+p_reco(2,j))**2) .gt. {}**2) then + passcuts_user=.false. + return + endif + endif + enddo + enddo From d856de4df0bf38c5e34f4197ec28e701ebd41e41 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Fri, 20 May 2022 15:42:45 +0200 Subject: [PATCH 04/11] Calculate EW corrections --- nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt index c48507fe..e33260d5 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/output.txt @@ -2,6 +2,6 @@ set complex_mass_scheme True import model loop_qcd_qed_sm_Gmu define p = p b b~ define j = p -generate p p > mu+ vm [QCD] +generate p p > mu+ vm [QCD QED] output @OUTPUT@ quit From b4da9bb8ddee30ed56d5ae38d81fae9450b81286 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Fri, 20 May 2022 15:44:09 +0200 Subject: [PATCH 05/11] Add PTL and PTNU observables --- .../CDF_WP_1960GEV_88FB_PTL/analysis.f | 56 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_PTL/launch.txt | 35 ++++++++++++ .../CDF_WP_1960GEV_88FB_PTL/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_PTL/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_PTNU/analysis.f | 56 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_PTNU/launch.txt | 35 ++++++++++++ .../CDF_WP_1960GEV_88FB_PTNU/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_PTNU/output.txt | 7 +++ 8 files changed, 214 insertions(+) create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/output.txt diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f new file mode 100644 index 00000000..154323f6 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f @@ -0,0 +1,56 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'ptl', 100, 30d0, 55d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision pplb(1:2) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(1:2)=p_reco(1:2,j) + enddo + + call HwU_fill(1,sqrt(pplb(1)**2+pplb(2)**2),wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt new file mode 100644 index 00000000..efce994c --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw @MW@ +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 60.0 +#user_defined_cut set mtwmax = 100.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.0005 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/metadata.txt new file mode 100644 index 00000000..978a50f8 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=ptl +x1_label_tex=$p_{\mathrm{T},\ell}$ +x1_unit=GeV +y_label=dsig/dptl +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}p_{\mathrm{T},\ell}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f new file mode 100644 index 00000000..3edbed19 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f @@ -0,0 +1,56 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'ptnu', 100, 30d0, 55d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(1:2) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.14) ppl(1:2)=p_reco(1:2,j) + enddo + + call HwU_fill(1,sqrt(ppl(1)**2+ppl(2)**2),wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt new file mode 100644 index 00000000..efce994c --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw @MW@ +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 60.0 +#user_defined_cut set mtwmax = 100.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.0005 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/metadata.txt new file mode 100644 index 00000000..5e971f7e --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=ptnu +x1_label_tex=$p_{\mathrm{T},\nu}$ +x1_unit=GeV +y_label=dsig/dptnu +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}p_{\mathrm{T},\nu}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit From bb3fb45fd640482a61f9f64b7773c601bb5dcea8 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Mon, 23 May 2022 11:11:56 +0200 Subject: [PATCH 06/11] Update pineappl dependency to fix anti-proton convolutions --- poetry.lock | 30 +++++++++++++++--------------- pyproject.toml | 2 +- 2 files changed, 16 insertions(+), 16 deletions(-) diff --git a/poetry.lock b/poetry.lock index 1a6718b9..bef9634c 100644 --- a/poetry.lock +++ b/poetry.lock @@ -589,7 +589,7 @@ python-versions = "*" [[package]] name = "pineappl" -version = "0.5.0" +version = "0.5.2" description = "Python bindings to PineAPPL" category = "main" optional = false @@ -1614,20 +1614,20 @@ pickleshare = [ {file = "pickleshare-0.7.5.tar.gz", hash = "sha256:87683d47965c1da65cdacaf31c8441d12b8044cdec9aca500cd78fc2c683afca"}, ] pineappl = [ - {file = "pineappl-0.5.0-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:01731da882e4624159c51c90c661ecf766eed4bb6fa78bb50ddba15829d66480"}, - {file = "pineappl-0.5.0-cp310-none-win_amd64.whl", hash = "sha256:6e072e18b1bcb105245ae2fa5d11a24fd00b540f4dc9b9bb6527c3b1c4d8792a"}, - {file = "pineappl-0.5.0-cp36-cp36m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:7691b5ccc5953c97062b1693733fbf5976d498b9bf8268b2dc0d35530ef82a62"}, - {file = "pineappl-0.5.0-cp36-none-win_amd64.whl", hash = "sha256:cedb1cf5af33ec32adf85a52b2815d6174ad03702a0de033f3273e250d8e03cf"}, - {file = "pineappl-0.5.0-cp37-cp37m-macosx_10_9_x86_64.macosx_11_0_arm64.macosx_10_9_universal2.whl", hash = "sha256:3ea10e89f52992f0cf6e4ccb467a316f5e6ef5496aff085608667fab69f28f93"}, - {file = "pineappl-0.5.0-cp37-cp37m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:218881a15b8c6f4cbc84847677004fa620ef7af210c5dddcf06627f4abffbbab"}, - {file = "pineappl-0.5.0-cp37-none-win_amd64.whl", hash = "sha256:47d968c9d5052cf62e864576081513d1cdc772bb756ac45a3aecb3be78e9c9bf"}, - {file = "pineappl-0.5.0-cp38-cp38-macosx_10_9_x86_64.macosx_11_0_arm64.macosx_10_9_universal2.whl", hash = "sha256:e6c931bc47619069a2029391c44a416d7d95ee340ee33ef8c03097ff422c218b"}, - {file = "pineappl-0.5.0-cp38-cp38-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:fc0c2696f2645f1876fa22b631fabeb10067067e3fc3a3fa1edc9ed9d53bdf50"}, - {file = "pineappl-0.5.0-cp38-none-win_amd64.whl", hash = "sha256:997822f02c3e2249bbcc215f80004e2c9c841b7879691ba4b20843e357680e87"}, - {file = "pineappl-0.5.0-cp39-cp39-macosx_10_9_x86_64.macosx_11_0_arm64.macosx_10_9_universal2.whl", hash = "sha256:883917faffddd6c077f64baf383682c779f4d3176a26b778a6879c680b16f6ad"}, - {file = "pineappl-0.5.0-cp39-cp39-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:0a447fb7d84de12af0416c07154fbb4234b5b6a7146f2ee8f54e70ca6c617309"}, - {file = "pineappl-0.5.0-cp39-none-win_amd64.whl", hash = "sha256:be9cc7c755b710d9113ad7bef13beea59ec32815d7e42b08622418e7cf4e3f5c"}, - {file = "pineappl-0.5.0.tar.gz", hash = "sha256:742a7f085db2a0c12dbf49750a1f1047e1aff4481627739cf0406b8c630d6fdd"}, + {file = "pineappl-0.5.2-cp310-cp310-macosx_10_9_x86_64.macosx_11_0_arm64.macosx_10_9_universal2.whl", hash = "sha256:d299e59c3ed65880db7b1b1cf693c737fb92544f9454d7bdb88198388d5ca380"}, + {file = "pineappl-0.5.2-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:8c81ad04a6785e584f22488b869327448f8ca9b73cab54714c62d0caac549bd5"}, + {file = "pineappl-0.5.2-cp310-none-win_amd64.whl", hash = "sha256:a22c7fa0bbd6fdcd174066c63f166dd107e4ad4e1301c57918d5895c87ab8b01"}, + {file = "pineappl-0.5.2-cp36-cp36m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:293a37bb0ec300ad2038270a723bb80c3f54803f38e41c026a8f5a56220759b5"}, + {file = "pineappl-0.5.2-cp37-cp37m-macosx_10_9_x86_64.macosx_11_0_arm64.macosx_10_9_universal2.whl", hash = "sha256:f16b2fea1b96cbd439f27f919a448d0e46d046fda94732e204ae726f2a0d22d0"}, + {file = "pineappl-0.5.2-cp37-cp37m-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:82e762eed440fc536ab82bb3b62ba81d77aaefb9c3dd04197c75e1d6fcb610e5"}, + {file = "pineappl-0.5.2-cp37-none-win_amd64.whl", hash = "sha256:106c1d5d91fae531cd6a7e5c1b9c926c61f1ef91a44285b3d053b9da830de452"}, + {file = "pineappl-0.5.2-cp38-cp38-macosx_10_9_x86_64.macosx_11_0_arm64.macosx_10_9_universal2.whl", hash = "sha256:e64f2cf7e9f26d5fcb12944a4bcd68a042c92ff7b5e813758d1490e117b8e966"}, + {file = "pineappl-0.5.2-cp38-cp38-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:1551154386e5b40dfb552aafe1e1c67c2aba0dbaf68c2a61cab945ab09546352"}, + {file = "pineappl-0.5.2-cp38-none-win_amd64.whl", hash = "sha256:0e7c96584c1564cc2c4b1ab3a714f75444458d9ff1f77d27d24191e1baef1710"}, + {file = "pineappl-0.5.2-cp39-cp39-macosx_10_9_x86_64.macosx_11_0_arm64.macosx_10_9_universal2.whl", hash = "sha256:7b285dfa3a0bc4b9a9665234667b3fe0f745fc16ab228cf4550b097eecf54b32"}, + {file = "pineappl-0.5.2-cp39-cp39-manylinux_2_17_x86_64.manylinux2014_x86_64.whl", hash = "sha256:77c45d043a49aa8f94fe9ed731b807ad95e9622afc1d93622ad7577d84f67a69"}, + {file = "pineappl-0.5.2-cp39-none-win_amd64.whl", hash = "sha256:cbc3e26f9aba7be7c04efc0b27db6ff9d97b3d3d1e9c309000d268be43d9ac9e"}, + {file = "pineappl-0.5.2.tar.gz", hash = "sha256:d8904c44edccad79a05d584743c4dfd33ab630a22a9d07c178c286c653134a1b"}, ] pkgconfig = [ {file = "pkgconfig-1.5.5-py3-none-any.whl", hash = "sha256:d20023bbeb42ee6d428a0fac6e0904631f545985a10cdd71a20aa58bc47a4209"}, diff --git a/pyproject.toml b/pyproject.toml index 133df5c9..1090da21 100644 --- a/pyproject.toml +++ b/pyproject.toml @@ -18,7 +18,7 @@ pygit2 = "^1.6.1" requests = "^2.26.0" a3b2bbc3ced97675ac3a71df45f55ba = "^6.4.0" lhapdf-management = "^0.2" -pineappl = "0.5.0" +pineappl = "0.5.2" # docs dependencies (for readthedocs) Sphinx = { version = "^4.2.0", optional = true } sphinx-rtd-theme = { version = "^1.0.0", optional = true } From 55282d03b1d4183bbb892956dbdc186aa8c04030 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Tue, 24 May 2022 11:07:45 +0200 Subject: [PATCH 07/11] Increase statistics --- nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt | 2 +- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt | 2 +- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt index efce994c..bd040c89 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt @@ -29,7 +29,7 @@ set etal = 1.0 #user_defined_cut set mtwmax = 100.0 #user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 -set req_acc_FO 0.0005 +set req_acc_FO 0.0001 set pineappl True done quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt index efce994c..bd040c89 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt @@ -29,7 +29,7 @@ set etal = 1.0 #user_defined_cut set mtwmax = 100.0 #user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 -set req_acc_FO 0.0005 +set req_acc_FO 0.0001 set pineappl True done quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt index efce994c..bd040c89 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt @@ -29,7 +29,7 @@ set etal = 1.0 #user_defined_cut set mtwmax = 100.0 #user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 -set req_acc_FO 0.0005 +set req_acc_FO 0.0001 set pineappl True done quit From 4107322444aba7a602d1994b00ebb456fe4e30c7 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Tue, 24 May 2022 13:19:52 +0200 Subject: [PATCH 08/11] Adjust ranges of distribution to reflect fit ranges --- nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f | 2 +- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f | 2 +- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f index 3e85a152..e6ed4029 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f @@ -7,7 +7,7 @@ subroutine analysis_begin(nwgt,weights_info) character*(*) weights_info(*) call HwU_inithist(nwgt,weights_info) - call HwU_book(1,'mt', 80, 60d0, 100d0) + call HwU_book(1,'mt', 50, 65d0, 90d0) return end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f index 154323f6..e7842fa7 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f @@ -7,7 +7,7 @@ subroutine analysis_begin(nwgt,weights_info) character*(*) weights_info(*) call HwU_inithist(nwgt,weights_info) - call HwU_book(1,'ptl', 100, 30d0, 55d0) + call HwU_book(1,'ptl', 64, 32d0, 48d0) return end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f index 3edbed19..70f8730f 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f @@ -7,7 +7,7 @@ subroutine analysis_begin(nwgt,weights_info) character*(*) weights_info(*) call HwU_inithist(nwgt,weights_info) - call HwU_book(1,'ptnu', 100, 30d0, 55d0) + call HwU_book(1,'ptnu', 64, 32d0, 48d0) return end From 723714a200911342e542940c5032067a9e9bb495 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Tue, 24 May 2022 21:03:56 +0200 Subject: [PATCH 09/11] Tighten cuts --- nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt | 4 ++-- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt | 4 ++-- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt | 4 ++-- 3 files changed, 6 insertions(+), 6 deletions(-) diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt index bd040c89..c8c1fa5e 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt @@ -25,8 +25,8 @@ set etal = 1.0 #user_defined_cut set ptlmax = 55.0 #user_defined_cut set ptmissmin = 30.0 #user_defined_cut set ptmissmax = 55.0 -#user_defined_cut set mtwmin = 60.0 -#user_defined_cut set mtwmax = 100.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 #user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 set req_acc_FO 0.0001 diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt index bd040c89..3978525a 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt @@ -20,9 +20,9 @@ set fixed_fac_scale True set mur_ref_fixed @MW@ set muf_ref_fixed @MW@ set reweight_scale True -set ptl = 30.0 +set ptl = 32.0 set etal = 1.0 -#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptlmax = 48.0 #user_defined_cut set ptmissmin = 30.0 #user_defined_cut set ptmissmax = 55.0 #user_defined_cut set mtwmin = 60.0 diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt index bd040c89..ecff3159 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt @@ -23,8 +23,8 @@ set reweight_scale True set ptl = 30.0 set etal = 1.0 #user_defined_cut set ptlmax = 55.0 -#user_defined_cut set ptmissmin = 30.0 -#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set ptmissmin = 32.0 +#user_defined_cut set ptmissmax = 48.0 #user_defined_cut set mtwmin = 60.0 #user_defined_cut set mtwmax = 100.0 #user_defined_cut set ptwmax = 15.0 From a7ccff44e8c95b754204e25f3579641f6e06cc36 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Tue, 24 May 2022 21:06:46 +0200 Subject: [PATCH 10/11] Increase statistics --- nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt | 2 +- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt | 2 +- nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt index c8c1fa5e..04c58de1 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt @@ -29,7 +29,7 @@ set etal = 1.0 #user_defined_cut set mtwmax = 90.0 #user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 -set req_acc_FO 0.0001 +set req_acc_FO 0.000025 set pineappl True done quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt index 3978525a..6886c219 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt @@ -29,7 +29,7 @@ set etal = 1.0 #user_defined_cut set mtwmax = 100.0 #user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 -set req_acc_FO 0.0001 +set req_acc_FO 0.000025 set pineappl True done quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt index ecff3159..682ba084 100644 --- a/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt @@ -29,7 +29,7 @@ set etal = 1.0 #user_defined_cut set mtwmax = 100.0 #user_defined_cut set ptwmax = 15.0 set rphreco = 0.2 -set req_acc_FO 0.0001 +set req_acc_FO 0.000025 set pineappl True done quit From 747e0f579a5d9fe203d21ff5acc4255d2640fcf2 Mon Sep 17 00:00:00 2001 From: Christopher Schwan Date: Thu, 9 Jun 2022 13:39:56 +0200 Subject: [PATCH 11/11] Add temples for MT distribution --- .../CDF_WP_1960GEV_88FB_MT_80354/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80354/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80354/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80354/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80358/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80358/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80358/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80358/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80362/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80362/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80362/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80362/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80366/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80366/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80366/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80366/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80370/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80370/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80370/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80370/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80374/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80374/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80374/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80374/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80378/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80378/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80378/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80378/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80382/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80382/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80382/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80382/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80386/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80386/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80386/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80386/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80390/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80390/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80390/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80390/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80394/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80394/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80394/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80394/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80398/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80398/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80398/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80398/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80402/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80402/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80402/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80402/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80406/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80406/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80406/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80406/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80410/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80410/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80410/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80410/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80414/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80414/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80414/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80414/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80418/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80418/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80418/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80418/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80422/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80422/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80422/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80422/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80426/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80426/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80426/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80426/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80430/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80430/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80430/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80430/output.txt | 7 +++ .../CDF_WP_1960GEV_88FB_MT_80434/analysis.f | 63 +++++++++++++++++++ .../CDF_WP_1960GEV_88FB_MT_80434/launch.txt | 35 +++++++++++ .../CDF_WP_1960GEV_88FB_MT_80434/metadata.txt | 9 +++ .../CDF_WP_1960GEV_88FB_MT_80434/output.txt | 7 +++ 84 files changed, 2394 insertions(+) create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/output.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/analysis.f create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/launch.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/metadata.txt create mode 100644 nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/output.txt diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/launch.txt new file mode 100644 index 00000000..877d7642 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.354 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80354/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/launch.txt new file mode 100644 index 00000000..2af12fa7 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.358 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80358/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/launch.txt new file mode 100644 index 00000000..a8c11146 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.362 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80362/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/launch.txt new file mode 100644 index 00000000..233b9ece --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.366 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80366/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/launch.txt new file mode 100644 index 00000000..7e21fa38 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.370 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80370/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/launch.txt new file mode 100644 index 00000000..2256600a --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.374 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80374/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/launch.txt new file mode 100644 index 00000000..245db628 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.378 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80378/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/launch.txt new file mode 100644 index 00000000..6de1da2e --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.382 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80382/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/launch.txt new file mode 100644 index 00000000..c2ad4211 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.386 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80386/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/launch.txt new file mode 100644 index 00000000..459adcb3 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.390 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80390/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/launch.txt new file mode 100644 index 00000000..7500921d --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.394 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80394/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/launch.txt new file mode 100644 index 00000000..79249ab8 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.398 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80398/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/launch.txt new file mode 100644 index 00000000..f2cc4323 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.402 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80402/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/launch.txt new file mode 100644 index 00000000..b6eb8e0a --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.406 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80406/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/launch.txt new file mode 100644 index 00000000..a3d97ae8 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.410 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80410/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/launch.txt new file mode 100644 index 00000000..7b69170e --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.414 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80414/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/launch.txt new file mode 100644 index 00000000..f573a409 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.418 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80418/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/launch.txt new file mode 100644 index 00000000..ae268ec1 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.422 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80422/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/launch.txt new file mode 100644 index 00000000..f4e1a9a0 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.426 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80426/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/launch.txt new file mode 100644 index 00000000..aad9e445 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.430 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80430/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/analysis.f b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/analysis.f new file mode 100644 index 00000000..e6ed4029 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/analysis.f @@ -0,0 +1,63 @@ +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_begin(nwgt,weights_info) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + integer nwgt + character*(*) weights_info(*) + + call HwU_inithist(nwgt,weights_info) + call HwU_book(1,'mt', 50, 65d0, 90d0) + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_end(dummy) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + double precision dummy + call HwU_write_file + return + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + subroutine analysis_fill(p,istatus,ipdg,wgts,ibody) +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + implicit none + include 'nexternal.inc' + include 'cuts.inc' + integer istatus(nexternal) + integer iPDG(nexternal) + integer ibody + integer i + integer j + double precision p(0:4,nexternal) + double precision wgts(*) + double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, + $ p, iPDG, is_nextph_iso, p_reco, iPDG_reco, + $ is_nextph_iso_reco) + + do j = nincoming+1, nexternal + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))- + & 2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2))) + + call HwU_fill(1,xmt,wgts) + end diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/launch.txt new file mode 100644 index 00000000..538eb5a9 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/launch.txt @@ -0,0 +1,35 @@ +launch @OUTPUT@ +fixed_order = ON +set maxjetflavor 5 +set gf @GF@ +set mh @MH@ +set mt @MT@ +set mw 80.434 +set mz @MZ@ +set wh @WH@ +set wt @WT@ +set ww @WW@ +set wz @WZ@ +set ebeam1 980 +set ebeam2 980 +set lpp2 -1 +set pdlabel lhapdf +set lhaid 324900 +set fixed_ren_scale True +set fixed_fac_scale True +set mur_ref_fixed @MW@ +set muf_ref_fixed @MW@ +set reweight_scale True +set ptl = 30.0 +set etal = 1.0 +#user_defined_cut set ptlmax = 55.0 +#user_defined_cut set ptmissmin = 30.0 +#user_defined_cut set ptmissmax = 55.0 +#user_defined_cut set mtwmin = 65.0 +#user_defined_cut set mtwmax = 90.0 +#user_defined_cut set ptwmax = 15.0 +set rphreco = 0.2 +set req_acc_FO 0.000025 +set pineappl True +done +quit diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/metadata.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/metadata.txt new file mode 100644 index 00000000..2979e5df --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=mt +x1_label_tex=$M_{\mathrm{T}}$ +x1_unit=GeV +y_label=dsig/dmt +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}M_{\mathrm{T}}}$ +y_unit=pb/GeV diff --git a/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/output.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/output.txt new file mode 100644 index 00000000..e33260d5 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT_80434/output.txt @@ -0,0 +1,7 @@ +set complex_mass_scheme True +import model loop_qcd_qed_sm_Gmu +define p = p b b~ +define j = p +generate p p > mu+ vm [QCD QED] +output @OUTPUT@ +quit