mirguest / MirguestIssueReport

5 stars 3 forks source link

一种检验信号是否含有本底的想法 #132

Open mirguest opened 11 years ago

mirguest commented 11 years ago

对IBD事例的观察

我们分别处理IBD和混合1.0MHz, 300ns本底的IBD事例,画出 edep_total, totalPE 及 edep_r3 的关系图:

我们看到,对于混合了本底的事例,有一些点落在了外面。 那么,能否用 edep_total, totalPE 及 edep_r3 三者的关系,作为判断某个信号是否包含本底的判断依据呢?

mirguest commented 11 years ago

图中包含了什么信息

对于一个特定信号的edep r,如果能量完全沉积,那么就可以给出: E_sig = k_sig(r_sig^3) * totalPE_sig 图中特定颜色就表示特定的edep r。 如果混合了本底,就会导致点的偏移。 即使某个本底没完全沉积,其edep_total与totalPE也会成一定关系: E_bkg = k_bkg(r_bkg^3) * totalPE_bkg 这样,我们观察的混合本底的信号,应该是 (E_sig+E_bkg) = k_new * (totalPE_sig+totalPE_bkg) 这样我们可以看到,k_new就会不一样,就是偏离了我们假设的点?

当然,这里只是草稿,需要另外计算。

yipeifeng commented 11 years ago

不同个数的Th混合之后的图:

yipeifeng commented 11 years ago
DingXuefeng commented 11 years ago

这些变量的定义有没有文档啊?比如edep_total的含义以及单位,还有得到edep_x,y,z的计算方法

mirguest commented 11 years ago

初步实现想法

我们可以寻找:

set k = edep_total / totalPE

寻找 k 和 edep_r^3 的关系

mirguest commented 11 years ago

下面是IBD eplus的结果: totalPE_div_edep_vs_edepr3

mirguest commented 11 years ago

下面是gamma 6130在玻璃罐厚度10cm下的情况 totalPE_div_edep_vs_edepr3_gamma6130

mirguest commented 11 years ago

由于edep_total可能比较小,结果导致totalPE/edep_total可能会很大。那么,我们换成使用k。 下面是 gamma 6130,30cm下的结果: edep_total_div_totalPE_vs_edepr3