在二维拟合脚本里面最后输出一个datamc.dat文件 第二行第二个是包含了三部分sig+peaking+nonpeaking,计算纯度我们需要把第二行第一个 除以 第二个值 但是用sideband的时候去估计peaking本底时需要我们去用第一个值,在up down left right 四个部分分别进行二维拟合得到第二行第一个值
这样我们用二维拟合得到的第二行第二个值 - sideband区域的四个第二行第一个之和 / 4
第三行是误差
double Err = (nsig.getError())*fr;
out1<<Err<<endl;
要在二维拟合程序里面加上这两段