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..e6ed4029 --- /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', 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/launch.txt b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/launch.txt new file mode 100644 index 00000000..04c58de1 --- /dev/null +++ b/nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/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 = 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/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..e33260d5 --- /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 QED] +output @OUTPUT@ +quit 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 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..e7842fa7 --- /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', 64, 32d0, 48d0) + 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..6886c219 --- /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 = 32.0 +set etal = 1.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 +#user_defined_cut set mtwmax = 100.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_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..70f8730f --- /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', 64, 32d0, 48d0) + 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..682ba084 --- /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 = 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 +set rphreco = 0.2 +set req_acc_FO 0.000025 +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 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 } 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 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