---primary

Using the offline software at BESIII
---primary
LI Gang
Group meeting, 2006-06-20
1

Content
• BOSS– short introduction
• Starts the jobs of simulation,
reconstruction and analysis
• Using of ROOT
• Popular fit package --- roofit
• Using your generator
All based on BOSS 6.0.0
2

BOSS--- BESIII Offline Software System
Machine: lxplus.ihep.ac.cn, lxslc.ihep.ac.cn
OS : Scientific Linux Cern (SLC)
Compiler : gcc-3.2.3
Framework :Based on GAUDI_v16r4
External lib: lcg/external; geant4 and etc.
Developing language : C++,some FORTRAN
Database : PostgreSql, MySql
Version management tool : CVS
Browsing files in CVS repository : ViewCVS
Configuration management tool : CMT
3

BOSS
Structure(link)
BOOST integrated with BOSS
4

Starts the jobs
• 可 以 按 照 网 页 上 的 说 明 来 建 立 自 己 的 工
作 环 境 。
• 也 可 以 采 用 网 上 给 出 的 脚 本 自 动 建 立 。
• 但 是 , 前 提 是 必 须 有 lxplus 和 lxslc 上 的 帐
号 , 以 及 koala 上 的 cvs 帐 号 ( 用 来 从 cvs
checkout 相 关 的 文 件 )
5

Setting up your CMT environment
• 首 先 , 连 接 CMT
– $mkdir cmthome
– $cp /afs/ihep.ac.cn/bes3/offline/template/cmthome/setupCMT.csh ~/cmthome
注 : 如 果 你 使 用 的 是 bash, 则 copy setupCMT.sh 到 你 的
cmthome 目 录 下 。
• 然 后 , 修 改 login requirements, 设 置 符 合 自 己 要 求 的 配 置
– $cp /afs/ihep.ac.cn/bes3/offline/template/cmthome/requirements ~/cmthome
– 编 辑 requirements 文 件 , 使 其 符 合 你 的 要 求
$vi requirements
注 : 需 要 修 改 的 可 能 有 两 个 地 方 :
set BES_RELEASE 6.0.0
macro WorkArea "/ihepbatch/bes/lig/workarea"
6

author ligang@ihep.ac.cn
# set the SITEROOT
set SITEROOT "/afs/ihep.ac.cn/bes3/offline"
…… … …
# Set release number
macro BES_RELEASE "6.0.0"
#macro BES_RELEASE "nightly_build"
# Setup the CMT search path
set BesArea "$(BES_DIST_AREA)/$(BES_RELEASE)"
path_remove CMTPATH "/sw/dist/"
path_prepend CMTPATH "${BesArea}"
# Setup the CMT search path
set BesArea "$(BES_DIST_AREA)/$(BES_RELEASE)"
path_remove CMTPATH "dist"
path_prepend CMTPATH "${BesArea}"
#Add your worarea to CMTPATH
macro WorkArea "/ihepbatch/bes/lig/workarea"
…… … …
7

• 最 后 , 设 置 环 境 完 成
– source setupCMT.csh
– cmt config
– source setup.csh
– mkdir workarea; cd workarea; cmt co TestRealse
注 : 最 好 把 以 下 语 句 放 到 .tcshrc 文 件 中 :
– source /ihepbatch/bes/lig/cmthome/setupCMT.csh
– source /ihepbatch/bes/lig/cmthome/setupCVS.csh
– source /ihepbatch/bes/lig/cmthome/setup.csh
– source /ihepbatch/bes/lig/workarea/TestRelease/TestRelease-00-00-
14/cmt/setup.csh
8

How to run simulation job
• Assuming you have done:
– Setup your environment in ~/cmthome
– Checkout package TestRelease to your $workarea
• (1) cd TestRelease/*/cmt
• (2) vi requirements
– #======== for Generator ============
– use PartPropSvc PartPropSvc-* Generator
– use Tester Tester-* Generator/BesGenInterface
– use Howl Howl-* Generator/BesGenInterface
– use Radee Radee-* Generator/BesGenInterface
– ……………………………………………………………..
– #========= for Simulation ==============
– use BesSim BesSim-* Simulation/BOOST
9

How to run simulation job
• (3) cmt config
• (4) cmt broadcast gmake
• (5) source setup.csh
• (6) cd ../run
• (7) edit jobOption file: jobOptions_sim.txt
10

jobOptions_sim.txt
– ApplicationMgr.DLLs += { "Tester" , "BesSim" };
– ApplicationMgr.TopAlg += { "BesGenModule" , "BesSim"};
– #include "Bes_Gen.txt"
– #include "PartPropSvc.txt“
– BesSim.Mdc = 1; (0:no construction 1; construction)
– BesSim.Tof = 1;
– BesSim.Emc = 1;
– BesSim.Muc = 1;
– BesSim.Field = true;
– BesSim.PhysicsList = 4;
• 1:BesPhysicsList 2:LHEP_GN
• 3:QGSP_GN 4:QGSP
for BOSS 6.0.0
11

jobOptions_sim.txt
• ApplicationMgr.DLLs += { "G4Svc" };
• ApplicationMgr.ExtSvc += { "G4Svc" };
• // interactive mode
• G4Svc.InteractiveG4 = false;
• //Geant4 run macro, to set directory and file name of output data
• G4Svc.FADSMacro = "run.mac";
• // geant4 verbosity
• G4Svc.RunVerbosity = 1;
• G4Svc.EventVerbosity = 0;
• G4Svc.TrackingVerbosity = 0;
for BOSS 6.0.0
• //number of events to simulate
• ApplicationMgr.EvtMax = 10;
12

How to run simulation job
• (8) edit generator card (run.cards)
– same as genbes.cards in BESII
– the generator name in run.cards must match with algorithm
name in jobOptions_sim.txt
• ApplicationMgr.DLLs += { “Rhopi" , "BesSim" };
• ApplicationMgr.TopAlg += { "BesGenModule" , "BesSim"};
• generator_name 'rhopi‘
• rhopi_ipar 1
– Ineffective control in run.cards
• Event range (controlled in jobOption)
• Random seed (will be added in the future by yugw)
for BOSS 6.0.0
13

Particle ID in run.cards
same as BESII
14

• (9) edit run.mac
How to run simulation job
– /run/runID -9
– /runAction/ascii 1 1 1 1 1 1 1 1 /home/lig/boost.dat
for BOSS 6.0.0
MDCTRUTH TOFTRUTH EMCTRUTH MUCTRUT
H
MDCDIGI
TOFDIG
I
EMCDIGI
MUCDIGI
– you must specify a full directory name, such as
• /ihepbatch/besdata/public/lig/boost.dat
• (10) boss.exe jobOptions_sim.txt
15

• In brief
How to run simulation job
– specify generator algorithm in jobOptions_sim.txt
• ApplicationMgr.DLLs += { “Rhopi" , "BesSim" };
• G4Svc.FADSMacro = "run.mac";
• ApplicationMgr.EvtMax = 10;
– open generator in run.cards
• generator_name ‘rhopi’
– specify output data file in run.mac
• /runAction/ascii 1 1 1 1 1 1 1 1 /home/lig/boost.dat
– using all the other default settings
for BOSS 6.0.0
16

How to run reconstruction job
boss -q jobOptions_rec.txt
#include "Ascii.txt“
#include "PartPropSvc.txt"
//zhangy
#include "$MDCTRKRECONROOT/share/jobOptions_MdcPatRec.txt“
//zangsl
//#include "$TRKRECOROOT/share/jobOptions_TrkReco.txt"
#include "$MDCDEDXALGROOT/share/job_dedx_all.txt"
#include "$EMCRECROOT/share/EmcRecOptions.txt"
#include "$TRKEXTALGROOT/share/TrkExtAlgOption.txt"
#include "$TOFRECROOT/share/jobOptions_TofRec.txt"
#include "$MUCRECALGROOT/share/jobOptions_MucRecTrkExt.txt"
#include "$ROOTIOROOT/share/jobOptions_Save2Root.txt"
17

How to run reconstruction job
// Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING,
5=ERROR, 6=FATAL )
MessageSvc.OutputLevel = 5;
// Number of events to be processed (default is 10)
ApplicationMgr.EvtMax = 200;
EventSelector.InputFiles="AsciiDmp_out.dat";
ApplicationMgr.HistogramPersistency = "ROOT";
RootCnvSvc.digiRootOutputFile =“ggchic.root";
NTupleSvc.Output = { "FILE101 DATAFILE='Mdc.root' OPT='NEW'
TYP='ROOT'", "FILE104 DATAFILE='kalfit.root' OPT='NEW'
TYP='ROOT'", "FILE103 DATAFILE='DeDx.root' OPT='NEW'
TYP='ROOT'", "FILE201 DATAFILE='Tof.root' OPT='NEW'
TYP='ROOT'", "FILE301 DATAFILE='Emc.root' OPT='NEW'
TYP='ROOT'", "FILE401 DATAFILE='Muc.root' OPT='NEW'
TYP='ROOT'" };
18

How to run analysis job
#include "$ROOTIOROOT/share/jobOptions_ReadDst.txt"
#include "$DSTEVENTASSEMBLYALGROOT/share/jobOptions_DstEventAssembly.txt"
#include "jobOptions_BParticles_true.txt"
#include "/ihepbatch/bes/lig/workarea/Analysis/Analysis-00-00-3/Physics/JpsiTagAlg/share/jobOptions_JPsiTag.txt"
// Input DST file name
EventCnvSvc.digiRootInputFile={
"/besdata/public/lig/data/zhangy/psip_kk2f_8105.root",
"/besdata/public/lig/data/zhangy/psip_kk2f_8106.root"};
boss -q jobOptions_ana.txt
G2MuMu.CmsEnergy = 3.686;
RootEvtSelector.Input = "ROOT";
// Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
MessageSvc.OutputLevel = 4;
// Number of events to be processed (default is 10)
ApplicationMgr.EvtMax = 500000;
// Histogram and Ntuple
ApplicationMgr.DLLs += { "RootHistCnv" };
ApplicationMgr.DLLs += { "NTuples" };
ApplicationMgr.HistogramPersistency = "ROOT";
NTupleSvc.Output = { "FILE1 DATAFILE='tag_zhangy_true.root' OPT='NEW' TYP='ROOT'"};
19

Something important
• /Event/AnalEvent/ : Bparticle
• /Event/EventModel
• /Event/McTruth
• External Package: CLHEP
20

Using of ROOT
•Batch mode
•Interactive mode
•Graphical mode
21

Popular fit package --- roofit
• Roofit 的 功 能 很 强 , 比 如 单 纯 的 拟 合 , 卷
积 , 考 虑 效 率 , 可 扩 展 性 等 。
• 先 安 装 roofit, 如 果 root 的 版 本 比 较 低 。
• 高 版 本 的 root 已 经 把 roofit 集 成 进 来 了 。
•epcfarm 上 , 我 已 经 安 装 了 roofit, 不 过 ,
我 建 议 大 家 自 己 下 载 一 个 新 版 本 的 root 试
着 安 装 一 下 , 这 样 最 简 单 。
• 下 面 看 个 例 子 :
22

#define ETA 0
void eta(){
TFile *f;TTree *t;
// define ppm
f = new TFile("etajuu.root");
RooRealVar ppm("ppm", "two gamma invariant mass", 0.3500, 0.7500,"GeV/c2");
RooRealVar mass("mass0","mass",0.54700,0.3000,0.6000);
//
RooRealVar mass1("mass1","mass",0.13580,0.1000,0.1500);
RooRealVar massresn("mres0","mass Resolution",4.00e-03,0.,0.040);
RooRealVar massresn1("mres1","mass Resolution1",1.00e-02,0.,0.040);
RooRealVar alpha("malpha","mass CB alpha ",1.65e-03,0.,3.000);
RooRealVar dmass("dmass","dmass CB n",1.73645,0.0,30.0);
// get data
t = (TTree *)f->Get("3");
RooDataSet data = RooDataSet("data","data",t,RooArgSet(ppm));
RooPlot* xframe1 = ppm.frame() ;
data.plotOn(xframe1) ; xframe1->Draw() ;
23

}
RooRealVar poly0("poly0","poly0",-1.0e10,1.0e10);
RooRealVar poly1("poly1","poly1",-1.0e10,1.0e10);
RooRealVar poly2("poly2","poly2",-1.0e10,1.0e10);
RooRealVar g1frac("g1frac","fraction of gauss",9.0e-01,0.0,1.0) ;
RooGaussian gauss("gauss", "Signal mass Distribution", ppm, mass1, massresn1);
RooCBShape cbs("cbs","cbs",ppm,mass, massresn,alpha,dmass);
RooPolynomial poly("poly", "Polynomial background ",ppm, RooArgList(poly0, poly1));
RooAddPdf sigDM("CBS","just CBS",RooArgList(cbs),RooArgList(g1frac));
sigDM.fitTo(data,"mhr") ;
RooPlot* xframe2 = ppm.frame() ;
sigDM.paramOn(xframe2,&data,"",1,"NELUC");
data.plotOn(xframe2) ;
sigDM.plotOn(xframe2,Components(cbs),LineStyle(kDashed));
xframe2->getAttLine()->SetLineColor(4);
sigDM.plotOn(xframe2);
TPostScript psout("./eps/etajuu_ggm.eps",113);psout.Off();
if(!psout){cout

Usung your own generator
• Check out the genbes package
• Put your own source code in it and rewrite
the requirement file, add code of
yours;compiling
• Write the interface subroutine (just
modifying a little words)
• Compile them
• Use it as usual
25

Bear in mind these, you can find
help everywhere!
• View the code:
http://koala.ihep.ac.cn/cgi-bin/viewcvs.cgi/BossCvs/
• find help and tutorial :
http://boss.ihep.ac.cn/GuideAndTraining/SetEnv-6.0.2.htm
• You’d better learn some shell or python program languages for
batch mode of simulation, reconstruction and analysis!
Thanks!
26