在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