-- JianqiaoDeng - 2022-11-10

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

Edit | Attach | Watch | Print version | History: r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r2 - 2022-11-17 - JianqiaoDeng
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2022 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback