在lxplus上产生相关过程的MC generator 事例

1)加载LHCb标准包,编译

注意替换为自己的用户名和目录,编译应该在5分钟左右完成。

首先需要登陆到 lxplus上,


cd /afs/cern.ch/work/h/hyin/workspace/

lb-set-platform x86_64-centos7-gcc9-opt

lb-dev Gauss/v54r2

cd ./GaussDev_v54r2

git lb-use Gauss

git lb-clone-pkg Gen/DecFiles

git lb-checkout Gauss/v54r2 Gen/GenCuts

cp /afs/cern.ch/user/h/hyin/public/For_YulongLiu/Xicc++_Xicpi_NoCUT.dec Gen/DecFiles/dkfiles/

cp /afs/cern.ch/user/h/hyin/public/For_YulongLiu/XiccDaughtersInLHCb.cpp Gen/GenCuts/src

make

Be careful here, GenCuts only needed for full phase-space Xicc++ study. Otherwise, just skip this step (GenCuts).

2)准备运行脚本


cd /afs/cern.ch/work/h/hyin/workspace/

mkdir Gauss_Run

cd Gauss_Run/

cp /afs/cern.ch/user/h/hyin/public/For_YulongLiu/Gauss-Job-1.py .

cp /afs/cern.ch/user/h/hyin/public/For_YulongLiu/run_gauss.sh .

cp /afs/cern.ch/user/h/hyin/public/For_YulongLiu/prepare.sh .

sh prepare.sh

这一步结束后,将会产生10个子目录,里面包含下一步运行所需要的文件。在转成ROOT文件之后,这些子目录就可以删掉,后面继续重复相关步骤。

注意,run_gauss.sh line 2, 根据自己的目录修改一下

MYDir=/afs/cern.ch/work/h/hyin/workspace/GaussDev_v54r2

3)运行脚本:

可以在华师服务器上运行screen命令,登陆到lxplus上

screen -S gauss_1

ssh XXX@lxplus.cern.ch

cd /afs/cern.ch/work/h/hyin/workspace/Gauss_Run/0

sh run_gauss.sh

...

到每个一个目录下,执行相关的脚本。大约3-4个子工作一次screen界面,然后每个子工作会运行1-2个小时左右。

(因为每个子工作,大约会产生2 G的文件,所以,不能一次提交太多子工作。每次10个左右正好)

产生的目标文件是每个子目录下面的 Gauss-26264056-10000ev-20201108.xgen 文件,以.xgen 结尾。

4)转换为root tuple:

下面将.xgen 转换为root可以分析的tree文件。


cd /afs/cern.ch/work/h/hyin/workspace/

mkdir DaVinci

cd DaVinci

cp /afs/cern.ch/user/h/hyin/public/For_YulongLiu/run_DaVinci.sh .

cp /afs/cern.ch/user/h/hyin/public/For_YulongLiu/Xicc_XicPi_DaVinci.py .

sh run_Davinci.sh

注意修改一下Xicc_XicPi_DaVinci.py, 把.xgen的路径修改正确,然后再执行。在line 19, as

datafile = "../Gauss_Run/0/Gauss-26264056-10000ev-20201108.xgen"

运行结束后,把产生的的root文件换一个名字,防止后面被新root文件冲掉。

此外,这里可以加多个.xgen 文件,比如

datafile = "../Gauss_Run/0/Gauss-26264056-10000ev-20201108.xgen, ../Gauss_Run/1/Gauss-26264056-10000ev-20201108.xgen"

重复2-3步骤,产生结束后,把所有产生的root tuple拿来做相关物理分析。

返回上一页

-- HangYin - 2020-11-08

Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r3 - 2020-11-09 - HangYin
 
    • 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-2023 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