Add CDF W mass determination runcards

nnpdf31_proc/CDF_WP_1960GEV_88FB_MT/analysis.f

@@ -0,0 +1,63 @@
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_begin(nwgt,weights_info)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      integer nwgt
+      character*(*) weights_info(*)
+
+      call HwU_inithist(nwgt,weights_info)
+      call HwU_book(1,'mt', 80, 60d0, 100d0)
+      return
+      end
+
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_end(dummy)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      double precision dummy
+      call HwU_write_file
+      return
+      end
+
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_fill(p,istatus,ipdg,wgts,ibody)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      include 'nexternal.inc'
+      include 'cuts.inc'
+      integer istatus(nexternal)
+      integer iPDG(nexternal)
+      integer ibody
+      integer i
+      integer j
+      double precision p(0:4,nexternal)
+      double precision wgts(*)
+      double precision ppl(0:3), pplb(0:3), ppv(0:3), xmt, getabsy
+      logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal)
+      external getabsy
+
+      double precision p_reco(0:4,nexternal)
+      integer iPDG_reco(nexternal)
+
+
+
+      call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco,
+     $                       p, iPDG, is_nextph_iso, p_reco, iPDG_reco,
+     $                       is_nextph_iso_reco)
+
+      do j = nincoming+1, nexternal
+        if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j)
+        if (iPDG_reco(j).eq.14) ppl(0:3)=p_reco(0:3,j)
+      enddo
+      do i=0,3
+        ppv(i)=ppl(i)+pplb(i)
+      enddo
+
+      xmt=sqrt(2d0*sqrt((ppl(1)**2+ppl(2)**2)*(pplb(1)**2+pplb(2)**2))-
+     &         2d0*(ppl(1)*pplb(1)+ppl(2)*pplb(2)))
+
+      call HwU_fill(1,xmt,wgts)
+      end

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 = 60.0
+#user_defined_cut set mtwmax = 100.0
+#user_defined_cut set ptwmax = 15.0
+set rphreco = 0.2
+set req_acc_FO 0.0005
+set pineappl True
+done
+quit

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

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

nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/analysis.f

@@ -0,0 +1,56 @@
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_begin(nwgt,weights_info)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      integer nwgt
+      character*(*) weights_info(*)
+
+      call HwU_inithist(nwgt,weights_info)
+      call HwU_book(1,'ptl', 100, 30d0, 55d0)
+      return
+      end
+
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_end(dummy)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      double precision dummy
+      call HwU_write_file
+      return
+      end
+
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_fill(p,istatus,ipdg,wgts,ibody)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      include 'nexternal.inc'
+      include 'cuts.inc'
+      integer istatus(nexternal)
+      integer iPDG(nexternal)
+      integer ibody
+      integer i
+      integer j
+      double precision p(0:4,nexternal)
+      double precision wgts(*)
+      double precision pplb(1:2)
+      logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal)
+      external getabsy
+
+      double precision p_reco(0:4,nexternal)
+      integer iPDG_reco(nexternal)
+
+
+
+      call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco,
+     $                       p, iPDG, is_nextph_iso, p_reco, iPDG_reco,
+     $                       is_nextph_iso_reco)
+
+      do j = nincoming+1, nexternal
+        if (iPDG_reco(j).eq.-13) pplb(1:2)=p_reco(1:2,j)
+      enddo
+
+      call HwU_fill(1,sqrt(pplb(1)**2+pplb(2)**2),wgts)
+      end

nnpdf31_proc/CDF_WP_1960GEV_88FB_PTL/launch.txt

@@ -0,0 +1,35 @@
+launch @OUTPUT@
+fixed_order = ON
+set maxjetflavor 5
+set gf @GF@
+set mh @MH@
+set mt @MT@
+set mw @MW@
+set mz @MZ@
+set wh @WH@
+set wt @WT@
+set ww @WW@
+set wz @WZ@
+set ebeam1 980
+set ebeam2 980
+set lpp2 -1
+set pdlabel lhapdf
+set lhaid 324900
+set fixed_ren_scale True
+set fixed_fac_scale True
+set mur_ref_fixed @MW@
+set muf_ref_fixed @MW@
+set reweight_scale True
+set ptl = 30.0
+set etal = 1.0
+#user_defined_cut set ptlmax = 55.0
+#user_defined_cut set ptmissmin = 30.0
+#user_defined_cut set ptmissmax = 55.0
+#user_defined_cut set mtwmin = 60.0
+#user_defined_cut set mtwmax = 100.0
+#user_defined_cut set ptwmax = 15.0
+set rphreco = 0.2
+set req_acc_FO 0.0005
+set pineappl True
+done
+quit

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

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

nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/analysis.f

@@ -0,0 +1,56 @@
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_begin(nwgt,weights_info)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      integer nwgt
+      character*(*) weights_info(*)
+
+      call HwU_inithist(nwgt,weights_info)
+      call HwU_book(1,'ptnu', 100, 30d0, 55d0)
+      return
+      end
+
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_end(dummy)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      double precision dummy
+      call HwU_write_file
+      return
+      end
+
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_fill(p,istatus,ipdg,wgts,ibody)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      include 'nexternal.inc'
+      include 'cuts.inc'
+      integer istatus(nexternal)
+      integer iPDG(nexternal)
+      integer ibody
+      integer i
+      integer j
+      double precision p(0:4,nexternal)
+      double precision wgts(*)
+      double precision ppl(1:2)
+      logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal)
+      external getabsy
+
+      double precision p_reco(0:4,nexternal)
+      integer iPDG_reco(nexternal)
+
+
+
+      call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco,
+     $                       p, iPDG, is_nextph_iso, p_reco, iPDG_reco,
+     $                       is_nextph_iso_reco)
+
+      do j = nincoming+1, nexternal
+        if (iPDG_reco(j).eq.14) ppl(1:2)=p_reco(1:2,j)
+      enddo
+
+      call HwU_fill(1,sqrt(ppl(1)**2+ppl(2)**2),wgts)
+      end

nnpdf31_proc/CDF_WP_1960GEV_88FB_PTNU/launch.txt

@@ -0,0 +1,35 @@
+launch @OUTPUT@
+fixed_order = ON
+set maxjetflavor 5
+set gf @GF@
+set mh @MH@
+set mt @MT@
+set mw @MW@
+set mz @MZ@
+set wh @WH@
+set wt @WT@
+set ww @WW@
+set wz @WZ@
+set ebeam1 980
+set ebeam2 980
+set lpp2 -1
+set pdlabel lhapdf
+set lhaid 324900
+set fixed_ren_scale True
+set fixed_fac_scale True
+set mur_ref_fixed @MW@
+set muf_ref_fixed @MW@
+set reweight_scale True
+set ptl = 30.0
+set etal = 1.0
+#user_defined_cut set ptlmax = 55.0
+#user_defined_cut set ptmissmin = 30.0
+#user_defined_cut set ptmissmax = 55.0
+#user_defined_cut set mtwmin = 60.0
+#user_defined_cut set mtwmax = 100.0
+#user_defined_cut set ptwmax = 15.0
+set rphreco = 0.2
+set req_acc_FO 0.0005
+set pineappl True
+done
+quit

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

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

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 }