diff --git a/nnpdf31_proc/CDF_DY_1960GEV_21FB/analysis.f b/nnpdf31_proc/CDF_DY_1960GEV_21FB/analysis.f
new file mode 100644
index 00000000..4c999020
--- /dev/null
+++ b/nnpdf31_proc/CDF_DY_1960GEV_21FB/analysis.f
@@ -0,0 +1,84 @@
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+      subroutine analysis_begin(nwgt,weights_info)
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      implicit none
+      integer nwgt
+      character*(*) weights_info(*)
+
+      call set_error_estimation(1)
+      call HwU_inithist(nwgt,weights_info)
+      call HwU_book(1,'rap', 27, 0.0d0, 2.7d0)
+      call HwU_book(2,'rap',  1, 2.7d0, 2.9d0)
+      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), xyll, 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.11) ppl(0:3)=p_reco(0:3,j)
+        if (iPDG_reco(j).eq.-11) pplb(0:3)=p_reco(0:3,j)
+      enddo
+      do i=0,3
+        ppv(i)=ppl(i)+pplb(i)
+      enddo
+
+      xyll=getabsy(ppv(0),ppv(3))
+
+      call HwU_fill(1,xyll,wgts)
+      call HwU_fill(2,xyll,wgts)
+
+ 999  return
+      end
+
+cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      function getabsy(en,pl)
+      implicit none
+      real*8 getabsy,en,pl,tmp
+c
+      tmp=pl/en
+      if(abs(tmp).lt.1d0)then
+        tmp=abs(atanh(tmp))
+      else
+        write(*,*)'Attempt to compute atanh(x) with x > 1'
+        stop
+      endif
+      getabsy=tmp
+      return
+      end
diff --git a/nnpdf31_proc/CDF_DY_1960GEV_21FB/launch.txt b/nnpdf31_proc/CDF_DY_1960GEV_21FB/launch.txt
new file mode 100644
index 00000000..947a5af4
--- /dev/null
+++ b/nnpdf31_proc/CDF_DY_1960GEV_21FB/launch.txt
@@ -0,0 +1,32 @@
+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 @MZ@
+set muf_ref_fixed @MZ@
+set reweight_scale True
+set ptl = 0.0
+set etal = -1.0
+#user_defined_cut set yz = 2.9
+set mll_sf = 66
+#user_defined_cut set mmllmax = 116.0
+set rphreco = 0.2
+set req_acc_FO 0.0001
+set pineappl True
+done
+quit
diff --git a/nnpdf31_proc/CDF_DY_1960GEV_21FB/metadata.txt b/nnpdf31_proc/CDF_DY_1960GEV_21FB/metadata.txt
new file mode 100644
index 00000000..a8a6a682
--- /dev/null
+++ b/nnpdf31_proc/CDF_DY_1960GEV_21FB/metadata.txt
@@ -0,0 +1,9 @@
+arxiv=0908.3914
+description=CDF differential Drell–Yan cross section at 1.96 TeV
+hepdata=10.17182/hepdata.52674.v1/t2
+x1_label=yll
+x1_label_tex=$y_{\ell\bar{\ell}}$
+x1_unit=GeV
+y_label=dsig/dyll
+y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}y_{\ell\bar{\ell}}}$
+y_unit=pb
diff --git a/nnpdf31_proc/CDF_DY_1960GEV_21FB/output.txt b/nnpdf31_proc/CDF_DY_1960GEV_21FB/output.txt
new file mode 100644
index 00000000..2f087638
--- /dev/null
+++ b/nnpdf31_proc/CDF_DY_1960GEV_21FB/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 > e+ e- [QCD QED]
+output @OUTPUT@
+quit
diff --git a/run.sh b/run.sh
index 2752d32b..1c375f44 100755
--- a/run.sh
+++ b/run.sh
@@ -328,6 +328,11 @@ merge() {
     fi
     mv "${grid}".tmp "${grid}"
 
+    if [[ -x ../nnpdf31_proc/"${dataset}"/precheck.sh ]]; then
+        cp ../nnpdf31_proc/"${dataset}"/precheck.sh .
+        GRID=$grid ./precheck.sh
+    fi
+
     # find out which PDF set was used to generate the predictions
     pdfstring=$(grep "set lhaid" "${launch_file}" | sed 's/set lhaid \([0-9]\+\)/\1/')