Epsilom的博客

原创 处理TTree的几种方法

TTree是ROOT的一个存储数据的类，大概可以理解成一个二维数据表格。每一个Branch代表一种属性（或者一列），Entry代表一个事例（或者一行）。用户可以通过读取每一Entry中的特定Branch，来统计一个样本的某个属性的分布。读取TTree数据的方法有很多，这里总结几个比较常用的，他们有的比较复杂但比较灵活，有的简单快捷，但是面对更复杂的数据结构或者操作就无能为力。不同的场景适用不同的方法，但个人觉得RDataFrame在各种情况都相对好用一些，写这篇博客也是起因于此。利用TTree本身的成

原创 图像重建过程中由于分立的传感器导致的周期性偏差

不想学习，随便记录一个偶尔想到的小问题。宏观世界是连续的，但测量事物的工具是不连续的。比如数码照相机以及人的视网膜等等，都是通过有限个灵敏单元组成一个阵列来探测事物的。这样的话在重建物体的形状或者位置的时候就会产生一定的偏差。最简单地，可以在一维空间举个例子：假设一维空间有一列传感器，每个传感器有一定的尺寸a，当一个传感器所在的[x, x+a)范围内有物体（或者说信号）存在，该传感器就返回1，否则返回0。如果一个尺寸和传感器尺寸相等的物体恰好与一个传感器重合的时候，重建出的位置恰好在这个传感器的中

原创 Valgrind: 这个代码诊断工具有点东西 : )

今天在网上逛街，发现了这样一个小玩意儿，感觉很有意思。起因大概是，同学短短几行代码，不知道为什么在运行的时候会突破服务器大约4G的memory限制而被系统挂起。我一直觉得这件事情很难解决，bug只要涉及到内存管理，出了毛病就很难找到。我想这样的事情可能有望借助今天找到的Valgrind工具解决。Valgrind: 这个诊断工具有点东西有很多博客讲这个工具，比如说知乎，CSDN，对于工具包massif，还有对应的可视化工具可以用massif-visualizer，等等。这里主要记一些我喜欢的认为比较常

原创 关于Docker中出现Failed to get D-Bus connection且--privileged也不管用的问题

简单讲就是不要用centos:7而要用centos:centos7.9.2009一些小众的科研软件非常古老，兼容性极低，在较常用的ubuntu、manjaro、centos8下很难装，只能生活在docker中。但他们经常在服务器中开发更新的版本，有的时候服务器系统升级了，新版本的软件就不兼容原来旧的docker。就很烦。最近想把官方的docker从sl6改到centos7，从而适应更新一点的软件。遇到了这样的问题：Failed to get D-Bus connection: Operation n

原创 ROOT叠加TH1并保留所有statistics box

详见ROOT官方文档https://root.cern/doc/master/classTPaveStats.html以及论坛帖子：https://root-forum.cern.ch/t/how-to-move-statistics-box-when-using-draw-sames/5139在调用TH1::Draw()方法时，可以使用“same”选项将两个TH1叠加在同一个TCanvas上，但这样做会导致图中只保留一个直方图的统计信息。例如：{ TH1D* h1 = new TH1D("h1"

原创 C++ Note: 字符串和容器

string转char：C风格字符串和C++风格字符串C++98和C++11的区别：定义 & 遍历

原创 Pandoc Note: markdown to PDF & templates

不想学习……最近绝对typora越来越不称手了，导出pdf的表格和图片总是会给分页造成麻烦，然后动不动就处半页纸都是白的。想把markdown转换成顺眼一点的pdf，或者用markdown写ppt什么的。然后就引出了什么beamer、marple、org-mode等等乱七八糟的工具。每个工具都好用的样子，但做不出选择是最烦的。所以干脆就先看看pandoc好了。pandoc可以完成各种文本格式的转换，当然包括markdown到pdf（latex）。重要的一点是pandoc有一套约定，用这套约定可以写好模

原创 Geant4:调整可视化界面的视角

刚开始赶工G4的作业，好不容易建好一个小人，塞到B1里面去发现怎么也找不到一个合适的视角，小人以各种崎岖的姿势站在坐标系里。很头疼，后来终于找到几条有用的命令。赶紧记一下。tldr版本几个指令：这些命令可以写在vis.mac脚本里，也可以编译之后通过Idle输入/vis/viewer/set/upVector [x] [y] [z]:在坐标系中，设置一个up-vector，调整viewe...

原创 第一次尝试RooFit

我想开学……课程作业要求做个拟合的小练习，其中TH1::Fit()已经用过了，试一试RooFit。我的理解上RooFit大概就是另外一票人基于ROOT做的另一个工具，在使用方法上，我感觉和ROOT本身有点格格不入。ROOT的手册里都用指针的，但是RooFit手册里都不用指针，就很烦。而且两边的数据类型都没有很圆滑地接口，各自用自己的一套。也没有时间详细写什么，先把码儿们放下，以后看心情吧。...

原创 使用zotero管理文献和书籍的一种同步方式

使用Zotero管理文献和书籍Zotero是类似于Endnote的一款文献管理软件。由于它是开源的、跨平台的软件，所以我在没有获得正版Endnote、并且偶尔使用linux的情况下，选用了这一款文献管理工具。到现在我断断续续地使用了将近两年的时间。这里并不介绍它的基本功能，而是记录一下，在短期轻度的使用过程中，我摸索到的一种自认为比较好用的同步方案。对尝试过的几种方案发些牢骚网络上有丰富...

原创 一个菜鸟的Geant4入门之路：alpha粒子轰击金箔的例子

一个菜鸟的Geant4入门之路：α\alphaα粒子轰击金箔的例子文章目录一个菜鸟的Geant4入门之路：α\alphaα粒子轰击金箔的例子前言去哪里找资料：几个重要的类一个活的G4程序需要哪几个类runManagerDetectorConstruction、PhysicsList和PrimaryGenerator一个“不平庸”的G4程序额外需要的类visManagerOptional User...

转载 转载：最舒服的色彩搭配RGB值配色 （赞、实用）

原文出处：https://www.icoa.cn/a/512.html最舒服的色彩搭配RGB值配色来源：逍遥峡谷 发布时间：2011-06-20 更新时间：2018...

转载 常用RGB

<div...

原创 Eorror: dyld: Symbol not found: __cg_jpeg_resync_to_restart

在MBP上编译安装完root，运行时发现报错：dyld: Symbol not found: __cg_jpeg_resync_to_restart重新将系统库中的库文件映射到/usr/local/lib下：cd /System/Library/Frameworks/ApplicationServices.framework/Versions/A/Frameworks/ImageIO.fr...

转载 转载：Ceph论文译文

译者注：本文是出于作者对于ceph的兴趣，在开源中国上关注ceph翻译，没有看到ceph论文的相关翻译，索性在阅读过程中把它翻译了出来，花费了几个周末时间，翻译过程中收获颇多，现把译文分享出来，如对您有益则倍感荣幸，肯定有很多不足之处。如有纰漏之处，欢迎指...

原创 Xrootd

Xrootd-NoteXrootd是高能物理领域通用的数据访问协议。起初是root中的一个模块，后来有了独立的版本。基本使用元数据处理：xrdfs对远端数据实现基本的ls,cd,mv,cat等操作：交互式：xrdfs root://ceph.xx.xx.cn/这样就进入了交互式页面，可以使用类似于bash的语法。命令行使用：直接把命令跟在上述指令后面：xrdfs root:/...

原创 MacOS Catalina词典参考源勾选失败的解决方法

macOS Catalina下的词典突然不好用了，总是说查不到此条。进入偏好设置之后取消相应词典的勾选，重新勾选的时候发现勾选的蓝色对勾总是闪一下就没了，没办法勾选。打了客服，各种重置还覆盖安装了系统还是不管用…技术支持也是不容易，估计也是本科毕业生，接受了一些基本的培训就开始工作了，这种冷门的问题他们可能也只会重启重装重买233。但网上看到一个解决方法十分简单：Finder中查找：.d...

原创 ROOT (a Data analysis Framework) TMatrix类的简单操作

简单使用TMatrix类ROOT中的TMatrix类可以实现矩阵运算。详细介绍见ROOTUsersGuide-14 Linear Algebra in ROOT一个内含Double_t元素的矩阵可以这样定义：TMatrixD x = TMatrix(nRow, nCol,)...

转载 Bug: version `GLIBCXX_3.4.26' not found (required by /usr/lib/root/libPyROOT.so)

ROOT --notebook的小问题今天打开jupyter发现root内核打不开了（说不定是manjaro滚挂了），主要报错是这样：Traceback (most recent call last): File "/opt/miniconda3/lib/python3.7/runpy.py", line 183, in _run_module_as_main mod_name, ...

原创 Excel2013: VLOOKUP函数

如何批量查找内容并填写相应的数据问题：在处理Excel表的过程中，经常会遇到要从一个现有的表格中挑选数据填入另一个表格中的问题。例如，有一个名单，填写成员的性别、身高等信息。另一个名单列举了部分需要被调用的成员名单，用户要从原名单中找出相应的成员，并将其性别身高等信息填入表格中。如果被调用的信息很多，反复手动查找就会很复杂。这时就需要一个函数来批量处理：VLOOKUP函数可以完成这样的工作。...

原创 C++参数传递：传值和传引用

C++参数传递：传值和传引用C++不允许函数return数组或多个对象，但很多时候，确实需要一个函数传回多个值。这时就希望函数直接对外部的变量或指针做操作，绕过return来达到目的。引用的概念曾经一直是盲区，觉得&仅可能是取地址符，&a就是取a的地址的意思。由此导致很多错误。引用在《C++ Primer》（中文第五版，45～46页）中讲到，引用其实是给变量取了一个别名：i...

原创 ROOT(a Data analysis Framework)-Note6: iSTEP day3-Random

Note5: Random3其实每天进展都还是听慢的，今天也只考虑了一个小问题，就是生成随机数。Monte-Carlo方法产生特定分布的随机数给定一个分布，或许不一定归一化：P(x=x′)=f(x)P(x = x&#x27;) = f(x)P(x=x′)=f(x)方法一：生成一个平均分布的随机数x∈[a1,a2]x \in [a_1,a_2]x∈[a1​,a2​]生成一...

原创 ROOT(a Data analysis Framework)-Note5: iSTEP day2-TH1::Fit()

Note5 ：数据拟合今天还是关于TH1和TF1的内容，使用Fit方法，使用不同的函数对数据进行拟合。下面以问题的形式做个小的总结记录。Section 1目标： Read histogram “\h1” from file “\hists.root” that I give you. And Draw it.读取一个.root文件：基于TFile类,有这样几种方法：TFile f(...

原创 ROOT(a Data analysis Framework)-Note4: iSTEP day1-TH1&TF1

ROOT(a Data analysis Framework)-Note4: iSTEP day1最近参加了一个叫做iSTEP的暑期学校，这里教授一些基本的ROOT使用方法。曾经用root的过程忘的差不多了，也没有记录。现在正好记录一下。TF1类TF1类是一个数学函数的类，基本的用法就是画一些函数的图像。root [0] TF1* f = new TF1("f","sin(x)",-10,...

原创 ROOT (a Data analysis Framework) - Note2：关于编译和运行

Note2：关于编译运行写在前面ROOT是基于C++的，即使有了cling的帮忙，也还是不能像python那样方便地交互运行。怎么说呢，就比如你写了一个脚本，在交互式界面中直接.x运行，ROOT还是把它当做一个C++脚本去看待的，诸如.L、.help之类的命令是不能出现在脚本中的。这样如果你自己写了一个类，在主函数中调用它之前，一般是include这个类的头文件，但是这在ROOT上并不好用，R...

原创 ROOT (a Data analysis Framework) - Note3：TCanvas类

TCanvas类TCanvas是一个画布的类，任何一张图，无论是直方图、函数曲线还是其他图，最终显示，都需要在这样一个画布上完成。实际上，如上一节那样，当你没有建立TCanvas实例就使用h.Draw()这个方法时，系统总会默认给这个Draw()方法创建一个画布，并且默认给一个名字叫c1。Info in <TCanvas::MakeDefCanvas>: created defa...

原创 ROOT (a Data analysis Framework) - Note1：先定一个小目标

ROOT-NoteROOT是一个基于C++的软件包，它被广泛地应用于高能物理实验的数据处理中。为了应付毕设， 近期开始学习ROOT的基本操作方法。打开ROOTRoot软件是基于C++的，但它嵌入了cling，因此可以通过命令行交互，还可以运行macro文件（实际上是一个C++脚本），可以解释运行也可以编译运行，这就很方便。在设置好环境变量的前提下，打开ROOT交互界面可以输入命令：roo...

原创 Linux连接投影仪

命令：xrandrxrandrSet the size, orientation and/or reflection of the outputs for a screen. - Display the current state of the system (known screens, resolutions, ...): xrandr --query - Disable d...

原创 ROOT类初始化后成员指针变成野指针的问题

ROOT类初始化后成员指针变成野指针的问题现在对面向对象和面向过程理解的不够透彻，只是希望用面向对象的形式写出一个类，把一套计算工作写成一个类似与工具包的样子。写出一个简易的雏形，大概构造是这样的：class A: public TObject{private: TH1D* a; TH2D* b; TH3D* data;public: A(); ~A(); void APl...

原创 远程使用Jupyter Notebook

远程使用Jupyter Notebook生成密码用来登录python终端下：In [1]: from IPython.lib import passwdIn [2]: passwd()Enter password: ...............Verify password: ................Out[2]: '....................' 创建配...

原创 关于SSH配置

关于SSH配置今天又想搞ssh玩，折腾了好久。主要问题是发现只有本机可以ssh自己，其他的网络就不行。本机和另一个终端用热点连接也是可以的，但是热点外的设备依然ssh不通。到目前只找到了解决办法，但是没有找到原因。主要步骤下载一个openssh：sudo apt-get install openssh-*然后配置sshd_config文件（ssh_config管客户端，就是如何连其他主...

原创 Geant4程序 T1program

Geant4程序 T1program基本内容本项目是初学Geant4的第一次尝试，希望完成一个满足基本运行要求的Geant4程序，支持可视化交互模式。物理过程是一束α\alphaα射线入射到金箔上。包括DetectorConstruction的构造（一片金箔），以及一支粒子枪，距离金箔数厘米，垂直射向金箔。文件结构：Program/T1program.ccCMakeLists.t...

原创 面向对象(Object-Oriented Programming) Note

面向对象(Object-Oriented Programming)What is different?// Ctypedef struct pont3d { float x; float y; float z;} Point3d;void Point3d_print(const Point3d* pd);Point3d a;a.x = 1; a.y = ...

原创 重新挂载目录的一个简便方法

重新挂载目录的一个简便方法deepin分区分的不好感觉固态40g很快就用完了，软件都装到/usr、/opt里去感觉不行。想办法把/opt迁移到/home下边。找到了个简便方法，就是在/home下弄一个/opt，然后把/opt挂到/home/opt上，设置下自动挂载：代码： 全选cd / ; sudo tar cf - usr opt | (cd /home ; sudo tar xfv ...

原创 日常消遣——改了一下WSL配色

日常消遣——改了一下WSL配色首先，我发现这个WSL实际上并不必须通过Ubuntu图标进入，而可以直接Win+R输入bash进入命令行界面，也可以开启多终端。cmd &amp; bash配色管理bash的配色是跟随cmd终端界面的设置而变的。有一个小软件ColorTool可以简便地切换配色方案：Github: https://github.com/microsoft/console下载后...

原创 Geant4 Note

Geant4 使用（1）一般的Program目录结构+- Program/ | CMakeLists.txt //CMake编译时需要的文件 | Program.cc // main()函数 | +- include | | ...include.hh... //放置头文件，声明需要的类 | +- src | | .....

原创 geant4安装以及Root配置

geant4安装以及Root配置先贴一个豆瓣链接:https://www.douban.com/note/431788626/ROOT按照cern root 安装：https://leeyeel.github.io/2017/07/01/root-install/简要来说就是如下几步：安装依赖关系： sudo apt-get install git dpkg-dev cmake ...

原创 摸索摸索Python3_一个巨大的空气劈_新浪博客

原创 第一次练习使用markdown_一个巨大的空气劈_新浪博客

---style: plain---# 关于蒙特卡洛模拟求解二维Ising模型的笔记—— markdown笔记学习## 基本原理假设有 $L \times L = N$大小的Ising模型1. 周期性边界条件：$$A(i,j) =A(i+m_{i}L,j+m_{j}L)$$（其中，$m_i,m_j,...

原创 MATLAB-钠原子光谱数据处理_一个巨大的空气劈_新浪博客

最近强迫症很严重，花好多工夫做了一些没什么用的事情。做完了感觉想留下来，然后又觉得占地方。突然就想出来还有博客这么一个东西，想想也是很好玩的。如果以后还有强迫症的时候，就把这些东西堆在这。如果遇上高手指点，还可以顺带学点东西……关于钠原子发射光谱实验，实在感觉数据分析比较繁杂。其它近代物理实验大多有编好的程序可以高效分析数据（虽然有些是老师或学长随手编出来的MathCAD小脚本，还...