Added pinecards for the normalising cross section of the CMS 2D mtt yt diff distr.

pinecards/CMS_TTB_8TEV_2D_TTM_TRAP_TOT/analysis.f

@@ -0,0 +1,74 @@
+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,'total rate      ', 1,0.5d0,1.5d0)
+
+      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'
+      integer istatus(nexternal)
+      integer iPDG(nexternal)
+      integer ibody
+      integer i
+      double precision p(0:4,nexternal)
+      double precision wgts(*)
+      double precision pttx(0:3), mtt, yt
+      double precision dot
+      double precision getabsy
+      external getabsy
+      double precision var
+
+      do i=0,3
+         pttx(i)=p(i,3)+p(i,4)
+      enddo
+
+      mtt = dsqrt(dot(pttx, pttx))
+      yt = abs(getabsy(p(0,3), p(3,3)))
+
+      var=1d0
+      if(mtt.gt.340d0.and.mtt.lt.1500d0)then
+         if(yt.lt.3.5d0)then
+            call HwU_fill(1,var,wgts)
+         endif
+      endif
+
+ 999  return
+      end
+
+ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+
+      function getabsy(en,pl)
+      implicit none
+      real*8 getabsy,en,pl,tmp
+      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

pinecards/CMS_TTB_8TEV_2D_TTM_TRAP_TOT/change_scale_to_et_avg.patch

@@ -0,0 +1,19 @@
+--- a/SubProcesses/setscales.f	2020-01-27 11:05:35.004477044 +0000
++++ b/SubProcesses/setscales.f	2020-01-27 11:06:31.970840053 +0000
+@@ -568,10 +568,12 @@
+ cc                 dynamical_scale_choice = 10                                   cc
+ cc      in the run_card (run_card.dat)                                           cc
+ ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+-         write(*,*) "User-defined scale not set"
+-         stop 1
+-         temp_scale_id='User-defined dynamical scale' ! use a meaningful string
+-         tmp = 0d0
++         tmp=0d0
++         do i=3,4
++           tmp=tmp+max(0d0,(pp(0,i)+pp(3,i))*(pp(0,i)-pp(3,i)))
++         enddo
++         tmp=dsqrt(tmp/2d0)
++         temp_scale_id='sqrt(sum_i mT(i)/2), i=final state'
+ ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
+ cc      USER-DEFINED SCALE: END OF USER CODE                                     cc
+ ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

pinecards/CMS_TTB_8TEV_2D_TTM_TRAP_TOT/launch.txt

@@ -0,0 +1,22 @@
+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 0.0
+set ww @WW@
+set wz @WZ@
+set ebeam1 4000
+set ebeam2 4000
+set pdlabel lhapdf
+set lhaid 324900
+set dynamical_scale_choice 10
+set reweight_scale True
+set req_acc_FO 0.0001
+set pineappl True
+done
+quit

pinecards/CMS_TTB_8TEV_2D_TTM_TRAP_TOT/metadata.txt

@@ -0,0 +1,7 @@
+arxiv=1703.01630
+description=CMS two-dimensional differential ttbar cross section at 8 TeV integrated over the fiducial phase space
+hepdata=10.17182/hepdata.77008.v1/t4
+nnpdf_id=CMS_TTBAR_2D_DIFF_MTT_TRAP_NORM_TOT
+y_label=sig
+y_label_tex=$\sigma$
+y_unit=pb

pinecards/CMS_TTB_8TEV_2D_TTM_TRAP_TOT/output.txt

@@ -0,0 +1,7 @@
+import model loop_qcd_qed_sm_Gmu
+set complex_mass_scheme True
+define p = p b b~
+define j = p
+generate p p > t t~ QCD^2<=4 QED^2<=2 [QCD]
+output @OUTPUT@
+quit