diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/analysis.f b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/analysis.f new file mode 100644 index 00000000..2c5c2ecd --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/analysis.f @@ -0,0 +1,88 @@ +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,'lmlp yrap',50, -1d0, 1d0) + + 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) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + double precision xcos, getcostheta + external getcostheta + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + is_nextph_iso(:) = .false. + 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) ppl(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xcos=getcostheta(ppl(0),ppl(1),ppl(2),ppl(3),pplb(0),pplb(1), + & pplb(2),pplb(3)) + + call HwU_fill(1,xcos,wgts) + end + + function getcostheta(enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal) + implicit none + real*8 getcostheta,enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal + real*8 pzll,mll,p1p,p1m,p2p,p2m,pt2ll +c implementation of first formula on page 6 of https://arxiv.org/abs/1710.05167 + p1p = enl + pzl + p1m = enl - pzl + p2p = enal + pzal + p2m = enal - pzal + pzll = pzl + pzal + pt2ll = (ptxl + ptxal) * (ptxl + ptxal) + + & (ptyl + ptyal) * (ptyl + ptyal) + mll = sqrt((enl + enal) * (enl + enal) - (pt2ll + pzll * pzll)) + getcostheta = (p1p*p2m-p1m*p2p)/sqrt(mll*mll+pt2ll)/mll* + & sign(1d0,pzll) + + if (abs(getcostheta) > 1d0) then + print *, 'Calculating cos with cos > 1' + stop + endif + end diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/launch.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/launch.txt new file mode 100644 index 00000000..b10b2cec --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/launch.txt @@ -0,0 +1,26 @@ +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 7000 +set ebeam2 7000 +set pdlabel lhapdf +set lhaid @LHAPDF_ID@ +set dynamical_scale_choice 10 +set reweight_scale True +set ptl = 10.0 +set etal = 2.4 +set mll_sf = 60.0 +#user_defined_cut set mmllmax = 120.0 +set req_acc_FO 0.00001 +set pineappl True +done +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/metadata.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/metadata.txt new file mode 100644 index 00000000..8a1ddd7d --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=costh +x1_label_tex=$\cos \theta^*$ +x1_unit= +y_label=dsig/dcosth +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}\cos \theta^*}$ +y_unit=pb diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/output.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/output.txt new file mode 100644 index 00000000..a64cd92a --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/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+ mu- [QCD] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/set_scale_to_mll.patch b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/set_scale_to_mll.patch new file mode 100644 index 00000000..10c218d7 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_060_120/set_scale_to_mll.patch @@ -0,0 +1,52 @@ +--- NLO/SubProcesses/setscales.f 2020-05-21 17:23:55.126143088 +0200 ++++ NLO/SubProcesses/setscales.f.new 2020-05-21 17:27:26.262700419 +0200 +@@ -527,6 +527,18 @@ + integer i,j + character*80 temp_scale_id + common/ctemp_scale_id/temp_scale_id ++ integer iPDG_reco(nexternal) ++ double precision ppl(0:3), pplb(0:3), ppv(0:3) ++ logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) ++ double precision p_reco(0:4,nexternal), p_in(0:4,nexternal) ++c les houches accord stuff to identify particles ++c ++ integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc), ++ & icolup(2,nexternal,maxflow),niprocs ++ common /c_leshouche_inc/idup,mothup,icolup,niprocs ++c Masses of external particles ++ double precision pmass(nexternal) ++ common/to_mass/pmass + c + tmp=0 + if(ickkw.eq.-1)then +@@ -568,10 +579,26 @@ + 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 ++ temp_scale_id='FOOBAR' ! use a meaningful string ++ tmp = -1d0 ++ do i=1,nexternal ++ p_in(0:3,i) = pp(0:3,i) ++ p_in(4,i) = pmass(i) ++ enddo ++ is_nextph_iso(:) = .false. ++ call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, ++ $ p_in, idup(1,1), is_nextph_iso, p_reco, ++ $ iPDG_reco, is_nextph_iso_reco) ++ ++ do j = nincoming+1, nexternal ++ if (iPDG_reco(j).eq.13) ppl(0:3)=p_reco(0:3,j) ++ if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) ++ enddo ++ do i=0,3 ++ ppv(i)=ppl(i)+pplb(i) ++ enddo ++ ++ tmp=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2) + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + cc USER-DEFINED SCALE: END OF USER CODE cc + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/analysis.f b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/analysis.f new file mode 100644 index 00000000..2c5c2ecd --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/analysis.f @@ -0,0 +1,88 @@ +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,'lmlp yrap',50, -1d0, 1d0) + + 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) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + double precision xcos, getcostheta + external getcostheta + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + is_nextph_iso(:) = .false. + 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) ppl(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xcos=getcostheta(ppl(0),ppl(1),ppl(2),ppl(3),pplb(0),pplb(1), + & pplb(2),pplb(3)) + + call HwU_fill(1,xcos,wgts) + end + + function getcostheta(enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal) + implicit none + real*8 getcostheta,enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal + real*8 pzll,mll,p1p,p1m,p2p,p2m,pt2ll +c implementation of first formula on page 6 of https://arxiv.org/abs/1710.05167 + p1p = enl + pzl + p1m = enl - pzl + p2p = enal + pzal + p2m = enal - pzal + pzll = pzl + pzal + pt2ll = (ptxl + ptxal) * (ptxl + ptxal) + + & (ptyl + ptyal) * (ptyl + ptyal) + mll = sqrt((enl + enal) * (enl + enal) - (pt2ll + pzll * pzll)) + getcostheta = (p1p*p2m-p1m*p2p)/sqrt(mll*mll+pt2ll)/mll* + & sign(1d0,pzll) + + if (abs(getcostheta) > 1d0) then + print *, 'Calculating cos with cos > 1' + stop + endif + end diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/launch.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/launch.txt new file mode 100644 index 00000000..5fca57b5 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/launch.txt @@ -0,0 +1,25 @@ +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 7000 +set ebeam2 7000 +set pdlabel lhapdf +set lhaid @LHAPDF_ID@ +set dynamical_scale_choice 10 +set reweight_scale True +set ptl = 10.0 +set etal = 2.4 +set mll_sf = 3000.0 +set req_acc_FO 0.0001 +set pineappl True +done +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/metadata.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/metadata.txt new file mode 100644 index 00000000..8a1ddd7d --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=costh +x1_label_tex=$\cos \theta^*$ +x1_unit= +y_label=dsig/dcosth +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}\cos \theta^*}$ +y_unit=pb diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/output.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/output.txt new file mode 100644 index 00000000..a64cd92a --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/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+ mu- [QCD] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/set_scale_to_mll.patch b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/set_scale_to_mll.patch new file mode 100644 index 00000000..10c218d7 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_3000/set_scale_to_mll.patch @@ -0,0 +1,52 @@ +--- NLO/SubProcesses/setscales.f 2020-05-21 17:23:55.126143088 +0200 ++++ NLO/SubProcesses/setscales.f.new 2020-05-21 17:27:26.262700419 +0200 +@@ -527,6 +527,18 @@ + integer i,j + character*80 temp_scale_id + common/ctemp_scale_id/temp_scale_id ++ integer iPDG_reco(nexternal) ++ double precision ppl(0:3), pplb(0:3), ppv(0:3) ++ logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) ++ double precision p_reco(0:4,nexternal), p_in(0:4,nexternal) ++c les houches accord stuff to identify particles ++c ++ integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc), ++ & icolup(2,nexternal,maxflow),niprocs ++ common /c_leshouche_inc/idup,mothup,icolup,niprocs ++c Masses of external particles ++ double precision pmass(nexternal) ++ common/to_mass/pmass + c + tmp=0 + if(ickkw.eq.-1)then +@@ -568,10 +579,26 @@ + 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 ++ temp_scale_id='FOOBAR' ! use a meaningful string ++ tmp = -1d0 ++ do i=1,nexternal ++ p_in(0:3,i) = pp(0:3,i) ++ p_in(4,i) = pmass(i) ++ enddo ++ is_nextph_iso(:) = .false. ++ call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, ++ $ p_in, idup(1,1), is_nextph_iso, p_reco, ++ $ iPDG_reco, is_nextph_iso_reco) ++ ++ do j = nincoming+1, nexternal ++ if (iPDG_reco(j).eq.13) ppl(0:3)=p_reco(0:3,j) ++ if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) ++ enddo ++ do i=0,3 ++ ppv(i)=ppl(i)+pplb(i) ++ enddo ++ ++ tmp=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2) + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + cc USER-DEFINED SCALE: END OF USER CODE cc + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/analysis.f b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/analysis.f new file mode 100644 index 00000000..2c5c2ecd --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/analysis.f @@ -0,0 +1,88 @@ +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,'lmlp yrap',50, -1d0, 1d0) + + 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) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + double precision xcos, getcostheta + external getcostheta + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + is_nextph_iso(:) = .false. + 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) ppl(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xcos=getcostheta(ppl(0),ppl(1),ppl(2),ppl(3),pplb(0),pplb(1), + & pplb(2),pplb(3)) + + call HwU_fill(1,xcos,wgts) + end + + function getcostheta(enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal) + implicit none + real*8 getcostheta,enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal + real*8 pzll,mll,p1p,p1m,p2p,p2m,pt2ll +c implementation of first formula on page 6 of https://arxiv.org/abs/1710.05167 + p1p = enl + pzl + p1m = enl - pzl + p2p = enal + pzal + p2m = enal - pzal + pzll = pzl + pzal + pt2ll = (ptxl + ptxal) * (ptxl + ptxal) + + & (ptyl + ptyal) * (ptyl + ptyal) + mll = sqrt((enl + enal) * (enl + enal) - (pt2ll + pzll * pzll)) + getcostheta = (p1p*p2m-p1m*p2p)/sqrt(mll*mll+pt2ll)/mll* + & sign(1d0,pzll) + + if (abs(getcostheta) > 1d0) then + print *, 'Calculating cos with cos > 1' + stop + endif + end diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/launch.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/launch.txt new file mode 100644 index 00000000..eae149c9 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/launch.txt @@ -0,0 +1,25 @@ +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 7000 +set ebeam2 7000 +set pdlabel lhapdf +set lhaid @LHAPDF_ID@ +set dynamical_scale_choice 10 +set reweight_scale True +set ptl = 10.0 +set etal = 2.4 +set mll_sf = 4000.0 +set req_acc_FO 0.0001 +set pineappl True +done +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/metadata.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/metadata.txt new file mode 100644 index 00000000..8a1ddd7d --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=costh +x1_label_tex=$\cos \theta^*$ +x1_unit= +y_label=dsig/dcosth +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}\cos \theta^*}$ +y_unit=pb diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/output.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/output.txt new file mode 100644 index 00000000..a64cd92a --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/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+ mu- [QCD] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/set_scale_to_mll.patch b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/set_scale_to_mll.patch new file mode 100644 index 00000000..10c218d7 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_4000/set_scale_to_mll.patch @@ -0,0 +1,52 @@ +--- NLO/SubProcesses/setscales.f 2020-05-21 17:23:55.126143088 +0200 ++++ NLO/SubProcesses/setscales.f.new 2020-05-21 17:27:26.262700419 +0200 +@@ -527,6 +527,18 @@ + integer i,j + character*80 temp_scale_id + common/ctemp_scale_id/temp_scale_id ++ integer iPDG_reco(nexternal) ++ double precision ppl(0:3), pplb(0:3), ppv(0:3) ++ logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) ++ double precision p_reco(0:4,nexternal), p_in(0:4,nexternal) ++c les houches accord stuff to identify particles ++c ++ integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc), ++ & icolup(2,nexternal,maxflow),niprocs ++ common /c_leshouche_inc/idup,mothup,icolup,niprocs ++c Masses of external particles ++ double precision pmass(nexternal) ++ common/to_mass/pmass + c + tmp=0 + if(ickkw.eq.-1)then +@@ -568,10 +579,26 @@ + 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 ++ temp_scale_id='FOOBAR' ! use a meaningful string ++ tmp = -1d0 ++ do i=1,nexternal ++ p_in(0:3,i) = pp(0:3,i) ++ p_in(4,i) = pmass(i) ++ enddo ++ is_nextph_iso(:) = .false. ++ call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, ++ $ p_in, idup(1,1), is_nextph_iso, p_reco, ++ $ iPDG_reco, is_nextph_iso_reco) ++ ++ do j = nincoming+1, nexternal ++ if (iPDG_reco(j).eq.13) ppl(0:3)=p_reco(0:3,j) ++ if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) ++ enddo ++ do i=0,3 ++ ppv(i)=ppl(i)+pplb(i) ++ enddo ++ ++ tmp=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2) + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + cc USER-DEFINED SCALE: END OF USER CODE cc + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/analysis.f b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/analysis.f new file mode 100644 index 00000000..2c5c2ecd --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/analysis.f @@ -0,0 +1,88 @@ +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,'lmlp yrap',50, -1d0, 1d0) + + 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) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + double precision xcos, getcostheta + external getcostheta + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + is_nextph_iso(:) = .false. + 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) ppl(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xcos=getcostheta(ppl(0),ppl(1),ppl(2),ppl(3),pplb(0),pplb(1), + & pplb(2),pplb(3)) + + call HwU_fill(1,xcos,wgts) + end + + function getcostheta(enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal) + implicit none + real*8 getcostheta,enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal + real*8 pzll,mll,p1p,p1m,p2p,p2m,pt2ll +c implementation of first formula on page 6 of https://arxiv.org/abs/1710.05167 + p1p = enl + pzl + p1m = enl - pzl + p2p = enal + pzal + p2m = enal - pzal + pzll = pzl + pzal + pt2ll = (ptxl + ptxal) * (ptxl + ptxal) + + & (ptyl + ptyal) * (ptyl + ptyal) + mll = sqrt((enl + enal) * (enl + enal) - (pt2ll + pzll * pzll)) + getcostheta = (p1p*p2m-p1m*p2p)/sqrt(mll*mll+pt2ll)/mll* + & sign(1d0,pzll) + + if (abs(getcostheta) > 1d0) then + print *, 'Calculating cos with cos > 1' + stop + endif + end diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/launch.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/launch.txt new file mode 100644 index 00000000..5e341583 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/launch.txt @@ -0,0 +1,25 @@ +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 7000 +set ebeam2 7000 +set pdlabel lhapdf +set lhaid @LHAPDF_ID@ +set dynamical_scale_choice 10 +set reweight_scale True +set ptl = 10.0 +set etal = 2.4 +set mll_sf = 5000.0 +set req_acc_FO 0.0001 +set pineappl True +done +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/metadata.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/metadata.txt new file mode 100644 index 00000000..8a1ddd7d --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=costh +x1_label_tex=$\cos \theta^*$ +x1_unit= +y_label=dsig/dcosth +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}\cos \theta^*}$ +y_unit=pb diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/output.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/output.txt new file mode 100644 index 00000000..a64cd92a --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/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+ mu- [QCD] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/set_scale_to_mll.patch b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/set_scale_to_mll.patch new file mode 100644 index 00000000..10c218d7 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_5000/set_scale_to_mll.patch @@ -0,0 +1,52 @@ +--- NLO/SubProcesses/setscales.f 2020-05-21 17:23:55.126143088 +0200 ++++ NLO/SubProcesses/setscales.f.new 2020-05-21 17:27:26.262700419 +0200 +@@ -527,6 +527,18 @@ + integer i,j + character*80 temp_scale_id + common/ctemp_scale_id/temp_scale_id ++ integer iPDG_reco(nexternal) ++ double precision ppl(0:3), pplb(0:3), ppv(0:3) ++ logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) ++ double precision p_reco(0:4,nexternal), p_in(0:4,nexternal) ++c les houches accord stuff to identify particles ++c ++ integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc), ++ & icolup(2,nexternal,maxflow),niprocs ++ common /c_leshouche_inc/idup,mothup,icolup,niprocs ++c Masses of external particles ++ double precision pmass(nexternal) ++ common/to_mass/pmass + c + tmp=0 + if(ickkw.eq.-1)then +@@ -568,10 +579,26 @@ + 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 ++ temp_scale_id='FOOBAR' ! use a meaningful string ++ tmp = -1d0 ++ do i=1,nexternal ++ p_in(0:3,i) = pp(0:3,i) ++ p_in(4,i) = pmass(i) ++ enddo ++ is_nextph_iso(:) = .false. ++ call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, ++ $ p_in, idup(1,1), is_nextph_iso, p_reco, ++ $ iPDG_reco, is_nextph_iso_reco) ++ ++ do j = nincoming+1, nexternal ++ if (iPDG_reco(j).eq.13) ppl(0:3)=p_reco(0:3,j) ++ if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) ++ enddo ++ do i=0,3 ++ ppv(i)=ppl(i)+pplb(i) ++ enddo ++ ++ tmp=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2) + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + cc USER-DEFINED SCALE: END OF USER CODE cc + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/analysis.f b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/analysis.f new file mode 100644 index 00000000..2c5c2ecd --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/analysis.f @@ -0,0 +1,88 @@ +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,'lmlp yrap',50, -1d0, 1d0) + + 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) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + double precision xcos, getcostheta + external getcostheta + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + is_nextph_iso(:) = .false. + 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) ppl(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) + enddo + do i=0,3 + ppv(i)=ppl(i)+pplb(i) + enddo + + xcos=getcostheta(ppl(0),ppl(1),ppl(2),ppl(3),pplb(0),pplb(1), + & pplb(2),pplb(3)) + + call HwU_fill(1,xcos,wgts) + end + + function getcostheta(enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal) + implicit none + real*8 getcostheta,enl,ptxl,ptyl,pzl,enal,ptxal,ptyal,pzal + real*8 pzll,mll,p1p,p1m,p2p,p2m,pt2ll +c implementation of first formula on page 6 of https://arxiv.org/abs/1710.05167 + p1p = enl + pzl + p1m = enl - pzl + p2p = enal + pzal + p2m = enal - pzal + pzll = pzl + pzal + pt2ll = (ptxl + ptxal) * (ptxl + ptxal) + + & (ptyl + ptyal) * (ptyl + ptyal) + mll = sqrt((enl + enal) * (enl + enal) - (pt2ll + pzll * pzll)) + getcostheta = (p1p*p2m-p1m*p2p)/sqrt(mll*mll+pt2ll)/mll* + & sign(1d0,pzll) + + if (abs(getcostheta) > 1d0) then + print *, 'Calculating cos with cos > 1' + stop + endif + end diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/launch.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/launch.txt new file mode 100644 index 00000000..4f84828a --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/launch.txt @@ -0,0 +1,25 @@ +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 7000 +set ebeam2 7000 +set pdlabel lhapdf +set lhaid @LHAPDF_ID@ +set dynamical_scale_choice 10 +set reweight_scale True +set ptl = 10.0 +set etal = 2.4 +set mll_sf = 6000.0 +set req_acc_FO 0.0001 +set pineappl True +done +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/metadata.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/metadata.txt new file mode 100644 index 00000000..8a1ddd7d --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=costh +x1_label_tex=$\cos \theta^*$ +x1_unit= +y_label=dsig/dcosth +y_label_tex=$\frac{\mathrm{d}\sigma}{\mathrm{d}\cos \theta^*}$ +y_unit=pb diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/output.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/output.txt new file mode 100644 index 00000000..a64cd92a --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/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+ mu- [QCD] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/set_scale_to_mll.patch b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/set_scale_to_mll.patch new file mode 100644 index 00000000..10c218d7 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_COS_6000/set_scale_to_mll.patch @@ -0,0 +1,52 @@ +--- NLO/SubProcesses/setscales.f 2020-05-21 17:23:55.126143088 +0200 ++++ NLO/SubProcesses/setscales.f.new 2020-05-21 17:27:26.262700419 +0200 +@@ -527,6 +527,18 @@ + integer i,j + character*80 temp_scale_id + common/ctemp_scale_id/temp_scale_id ++ integer iPDG_reco(nexternal) ++ double precision ppl(0:3), pplb(0:3), ppv(0:3) ++ logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) ++ double precision p_reco(0:4,nexternal), p_in(0:4,nexternal) ++c les houches accord stuff to identify particles ++c ++ integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc), ++ & icolup(2,nexternal,maxflow),niprocs ++ common /c_leshouche_inc/idup,mothup,icolup,niprocs ++c Masses of external particles ++ double precision pmass(nexternal) ++ common/to_mass/pmass + c + tmp=0 + if(ickkw.eq.-1)then +@@ -568,10 +579,26 @@ + 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 ++ temp_scale_id='FOOBAR' ! use a meaningful string ++ tmp = -1d0 ++ do i=1,nexternal ++ p_in(0:3,i) = pp(0:3,i) ++ p_in(4,i) = pmass(i) ++ enddo ++ is_nextph_iso(:) = .false. ++ call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, ++ $ p_in, idup(1,1), is_nextph_iso, p_reco, ++ $ iPDG_reco, is_nextph_iso_reco) ++ ++ do j = nincoming+1, nexternal ++ if (iPDG_reco(j).eq.13) ppl(0:3)=p_reco(0:3,j) ++ if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) ++ enddo ++ do i=0,3 ++ ppv(i)=ppl(i)+pplb(i) ++ enddo ++ ++ tmp=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2) + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + cc USER-DEFINED SCALE: END OF USER CODE cc + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/analysis.f b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/analysis.f new file mode 100644 index 00000000..48675833 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/analysis.f @@ -0,0 +1,87 @@ +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,'lmlp yrap',25, 0d0, 1d0) + + 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) + logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) + double precision xyll, getabsy + external getabsy + + double precision p_reco(0:4,nexternal) + integer iPDG_reco(nexternal) + + + + is_nextph_iso(:) = .false. + 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) ppl(0:3)=p_reco(0:3,j) + if (iPDG_reco(j).eq.-13) 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)) + + !if (xyll > 1.0d0) then + ! write (*,*) "error: event outside bin", xyll + ! stop 1 + !endif + + call HwU_fill(1,xyll,wgts) + end + +cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + + function getabsy(en,pl) + implicit none + real*8 getabsy,en,pl,tmp +c + tmp=pl/en + if (abs(tmp) < 1d0) then + tmp=abs(atanh(tmp)) + else + print *, 'Attempt to compute atanh(x) with x > 1' + stop + endif + getabsy=tmp + return + end diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/launch.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/launch.txt new file mode 100644 index 00000000..5e341583 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/launch.txt @@ -0,0 +1,25 @@ +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 7000 +set ebeam2 7000 +set pdlabel lhapdf +set lhaid @LHAPDF_ID@ +set dynamical_scale_choice 10 +set reweight_scale True +set ptl = 10.0 +set etal = 2.4 +set mll_sf = 5000.0 +set req_acc_FO 0.0001 +set pineappl True +done +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/metadata.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/metadata.txt new file mode 100644 index 00000000..6bd2c102 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/metadata.txt @@ -0,0 +1,9 @@ +arxiv= +description= +hepdata= +x1_label=yll +x1_label_tex=$|y_{\ell\bar{\ell}}|$ +x1_unit= +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/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/output.txt b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/output.txt new file mode 100644 index 00000000..a64cd92a --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/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+ mu- [QCD] +output @OUTPUT@ +quit diff --git a/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/set_scale_to_mll.patch b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/set_scale_to_mll.patch new file mode 100644 index 00000000..10c218d7 --- /dev/null +++ b/nnpdf31_proc/NNPDF_DY_14TEV_BSM_AFB_YLL_5000/set_scale_to_mll.patch @@ -0,0 +1,52 @@ +--- NLO/SubProcesses/setscales.f 2020-05-21 17:23:55.126143088 +0200 ++++ NLO/SubProcesses/setscales.f.new 2020-05-21 17:27:26.262700419 +0200 +@@ -527,6 +527,18 @@ + integer i,j + character*80 temp_scale_id + common/ctemp_scale_id/temp_scale_id ++ integer iPDG_reco(nexternal) ++ double precision ppl(0:3), pplb(0:3), ppv(0:3) ++ logical is_nextph_iso(nexternal),is_nextph_iso_reco(nexternal) ++ double precision p_reco(0:4,nexternal), p_in(0:4,nexternal) ++c les houches accord stuff to identify particles ++c ++ integer idup(nexternal,maxproc),mothup(2,nexternal,maxproc), ++ & icolup(2,nexternal,maxflow),niprocs ++ common /c_leshouche_inc/idup,mothup,icolup,niprocs ++c Masses of external particles ++ double precision pmass(nexternal) ++ common/to_mass/pmass + c + tmp=0 + if(ickkw.eq.-1)then +@@ -568,10 +579,26 @@ + 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 ++ temp_scale_id='FOOBAR' ! use a meaningful string ++ tmp = -1d0 ++ do i=1,nexternal ++ p_in(0:3,i) = pp(0:3,i) ++ p_in(4,i) = pmass(i) ++ enddo ++ is_nextph_iso(:) = .false. ++ call recombine_momenta(rphreco, etaphreco, lepphreco, quarkphreco, ++ $ p_in, idup(1,1), is_nextph_iso, p_reco, ++ $ iPDG_reco, is_nextph_iso_reco) ++ ++ do j = nincoming+1, nexternal ++ if (iPDG_reco(j).eq.13) ppl(0:3)=p_reco(0:3,j) ++ if (iPDG_reco(j).eq.-13) pplb(0:3)=p_reco(0:3,j) ++ enddo ++ do i=0,3 ++ ppv(i)=ppl(i)+pplb(i) ++ enddo ++ ++ tmp=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2) + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc + cc USER-DEFINED SCALE: END OF USER CODE cc + ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc diff --git a/nnpdf31_proc/NNPDF_TTB_14TEV_40_PHENO/launch.txt b/nnpdf31_proc/NNPDF_TTB_14TEV_40_PHENO/launch.txt index 142c2b23..d9e3098c 100644 --- a/nnpdf31_proc/NNPDF_TTB_14TEV_40_PHENO/launch.txt +++ b/nnpdf31_proc/NNPDF_TTB_14TEV_40_PHENO/launch.txt @@ -16,8 +16,6 @@ set pdlabel lhapdf set lhaid @LHAPDF_ID@ set dynamical_scale_choice 10 set reweight_scale True -#user_defined_tau_min 340.0 -#enable_patch no_pole_cancellation_checks set req_acc_FO 0.0001 set pineappl True done