Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add PDF4LHC21 cards #121

Closed
wants to merge 7 commits into from
Closed
Changes from 1 commit
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
16 changes: 8 additions & 8 deletions nnpdf31_proc/PDF4LHC_WP_13_TEV_21_PHENO/analysis.f
Original file line number Diff line number Diff line change
@@ -7,7 +7,7 @@ subroutine analysis_begin(nwgt,weights_info)
character*(*) weights_info(*)

call HwU_inithist(nwgt,weights_info)
call HwU_book(1,'mll', 6, 0d0, 6d0)
call HwU_book(1,'mll', 7, 0d0, 7d0)
return
end

@@ -55,19 +55,19 @@ subroutine analysis_fill(p,istatus,ipdg,wgts,ibody)
obs=-1d0

if (xmll.lt.1000d0) then
obs = -2d0
else if (xmll.lt.1500d0) then
obs=0.5d0
else if (xmll.lt.2000d0) then
else if (xmll.lt.1500d0) then
obs=1.5d0
else if (xmll.lt.3000d0) then
else if (xmll.lt.2000d0) then
obs=2.5d0
else if (xmll.lt.4000d0) then
else if (xmll.lt.3000d0) then
obs=3.5d0
else if (xmll.lt.5000d0) then
else if (xmll.lt.4000d0) then
obs=4.5d0
else if (xmll.lt.7000d0) then
else if (xmll.lt.5000d0) then
obs=5.5d0
else if (xmll.lt.7000d0) then
obs=6.5d0
else
obs=-3d0
endif
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
--- PDF4LHC_WP_13_TEV_21_PHENO/SubProcesses/cuts.f.old 2022-01-18 14:44:46.095088533 +0100
+++ PDF4LHC_WP_13_TEV_21_PHENO/SubProcesses/cuts.f 2022-01-18 14:49:52.421021878 +0100
@@ -975,25 +975,29 @@
@@ -975,25 +975,33 @@
c 'bias' = event_norm
c
implicit none
@@ -10,7 +10,7 @@
integer ipdg(nexternal),i
+ double precision ppl(0:3), pplb(0:3), ppv(0:3)
+ double precision p_reco(0:4,nexternal)
+ integer iPDG_reco(nexternal),grid
+ integer iPDG_reco(nexternal),grid,xmlnu
+

bias_wgt=1d0
@@ -39,7 +39,11 @@
+ ppv(i)=ppl(i)+pplb(i)
+ enddo
+
+ bias_wgt=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2)**3
+ xmlnu=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2)
+ if (xmlnu < 1000.0d0) then
+ bias_wgt=1d-10
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hi,
this is strange. How about using a smooth function, rather than a step function?
What if you set bias_wgt = xmlnu**2 for example?

+ endif
+ bias_wgt=xmlnu**3
return
end

Original file line number Diff line number Diff line change
@@ -45,8 +45,8 @@
+
+ xmll=sqrt(ppv(0)**2-ppv(1)**2-ppv(2)**2-ppv(3)**2)
+
+ if (xmll.lt.1000d0) then
+ write (*,*) "error: event outside bins", xmll
+ if (xmll.lt.0000d0) then
+ tmp=0.5*(0d0+1000d0)
+ else if (xmll.lt.1500d0) then
+ tmp=0.5*(1000d0+1500d0)
+ else if (xmll.lt.2000d0) then
10 changes: 9 additions & 1 deletion nnpdf31_proc/PDF4LHC_WP_13_TEV_21_PHENO/launch.txt
Original file line number Diff line number Diff line change
@@ -21,7 +21,15 @@ set ptl = 22.0
set etal = 2.4
#user_defined_cut set mmlnumax = 7000.0
set rphreco = 0.1
set req_acc_FO 0.001
# with a bias function we degrade the uncertainty of the integrated cross
# section, so don't optimize this quantity; instead, set the integration
# parameters explicitly
set req_acc_FO -1
set npoints_fo_grid 500000 # step 0 points
set niters_fo_grid 100 # step 0 iterations
set npoints_fo 100000000 # step 1 points
set niters_fo 50 # step 1 iterations
#
set pineappl True
done
quit
2 changes: 1 addition & 1 deletion nnpdf31_proc/PDF4LHC_WP_13_TEV_21_PHENO/postrun.sh
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
#!/bin/bash

pineappl remap "${GRID}" "${GRID}".tmp '1000,1500,2000,3000,4000,5000,7000'
pineappl remap "${GRID}" "${GRID}".tmp '0,1000,1500,2000,3000,4000,5000,7000'
mv "${GRID}".tmp "${GRID}"