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', 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

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

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_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

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

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

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

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

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

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

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