本项目是一个使用 CMS Open Data 计算双μ子不变质量谱的简单分析样例。
本项目的分析代码基于 CERN Open Data 门户网站上的 Example code to produce the di-muon spectrum, 参考了另一个 GitHub 项目 DimuonSpectrum2011. 并加入了对蒙特卡洛样本的比较分析。分析流程参照 Searching in CMS Open Data for Dimuon Resonances with Substantial Transverse Momentum 一文,以下称为例文。
对原始代码的修改如下:
- 为了避免与工作区中任何现有的
DemoAnalyzer
插件冲突,将类名从DemoAnalyzer
更改为DimuonSpectrum2011MC
. - 文件路径已被修改为相对于配置文件中的相对路径,即它们指向
datasets
目录,该目录位于将要运行程序的目录下。 - 增加了例文中的选择条件,具体请参考
src/DimuonSpectrum2010.cc
文件
您可以在 CMS Open Data VM 上运行此代码。如果您尚未安装 CMSSW,请执行以下操作:
cmsrel CMSSW_5_3_32
您也可以在 CMS Open Data Docker 的容器上运行此代码。如执行以下操作:
docker run --name dimu2011mc -it cmsopendata/cmssw_5_3_32 bash
如果您已经安装了 CMSSW 或运行了 Docker 容器,请直接运行如下命令:
cd CMSSW_5_3_32/src
cmsenv
接下来,您需要创建一个工作目录,可以将其命名为 WorkDir
或任意其他名称。并将本项目下载到此目录下。
mkdir WorkDir
cd WorkDir
git clone git://github.com/PhyM73/DimuonSpectrum2011MC.git
(此处默认下载 master
分支,如需指定其他分支,添加参数 -b branchname
即可)
再转到本项目的目录,使用 scram b
进行编译。
cd DimuonSpectrum2011MC
scram b
其中输入文件已经在配置文件 demoanalyzer_cfg.py
中定义,并且放置在 datasets
目录下。接下来直接运行配置文件即可。
cmsRun demoanalyzer_cfg.py
输出文件是一个包含多个直方图的 ROOT 文件,默认情况下命名为 DoubleMuMC.root, 包含有10000个输入 event
。您可以使用 ROOT 查看这些内容。您还可以修改 demoanalyzer_cfg.py
中的相关部分以选择 datasets
中其他的输入文件,并重新运行和比较。或者使用 ROOT 软件将各个输出文件合并以获得更全面的分析样本(如 hadd
命令),同时也请注意调节 demoanalyzer_cfg.py
中的 process.maxEvent
参数。
有关更多详细信息,请参阅 src/DimuonSpectrum2011MC.cc
中的注释。
为了方便对全部文件进行分析,您还可以在高通量计算平台上运行本项目,请根据需要选择相应的分支(详见下方),并注意任务提交的次数。
为方便不同的分析样本和平台,本项目有数个分支如下:
master
用于个人笔记本上的 VM 或 Docker 平台上的测试和简单分析,默认分析单个蒙特卡洛样本文件中的10000个事例。datacondor
用于高通量计算平台,默认分析单个2011年 DoubleMu 数据文件全部事例。可配合多次提交以分析全部2011年数据。mccondor
用于高通量计算平台,默认分析单个2011年蒙特卡洛样本文件全部事例。可配合多次提交以分析全部蒙特卡洛样本文件。
datasets
目录下分别包含如下文件
-
datasets/double
DoubleMu primary dataset in AOD format from RunA of 2011 -
datasets/single
SingleMu primary dataset in AOD format from RunA of 2011 -
datasets/mc
-
CMS list of validated runs for primary datasets of 2011 data taking
如有需要,可以从 CMS luminosity information for 2011 CMS open data 中获取对撞数据的亮度信息;可以从 McM 上获取蒙特卡洛数据的反应截面。