Install and test run:

https://gitlab.cern.ch/hyin/singlewz_generators/-/tree/master/MCFM#build

如果需要产生微分截面，可以在这里 MCFM-10.1/src/User/nplotter_Z_new.f90 (line 41-71, line 123)的分bin方式

以及W的 MCFM/MCFM-10.1/src/User/nplotter_W_new.f90 (line 41-62, line 113) （修改完注意需要编译）

产生的文件就在Bin目录下， MCFM-10.1/Bin/Z13TeV/Z_only_nnlo_CT18NNLO_1.00_1.00_6.0E-3_Z13TeV_total_cross.txt

How to produce W cross section at different eta bins?

add/change some lines in MCFM/MCFM-10.1/src/User/nplotter_W_new.f90

from"allocate(histos(4))" to "allocate(histos(5))"

histos(5) = plot_setup_uniform(2.0d0,4.5d0,0.25d0,'etalept_cms')

include 'src/Inc/nwz.f'

real(dp) :: etalept

if(nwz == +1) then
etalept = etarap(4,p)
else
etalept = etarap(3,p)
endif

from"

vals = [pt34,pt34,pt34,phistar]
wts = [wt*trans,wt*trans,wt*trans,wt*trans]"

to"

vals = [pt34,pt34,pt34,phistar,etalept]
wts = [wt*trans,wt*trans,wt*trans,wt*trans,wt*trans]"

if you made some changes and you need to make again:

cd MCFM-10.1

source ../setup_cvmfs.sh

cd Bin

cmake ..

make -j16

./mcfm input_LHCb_WW_LO.ini