xiaoguangbin/AnalysisSystemForRadionuclide
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

feat：beta energy calibration通过b_energy数据重新拟合一条初始beta曲线

Browse Source
This commit is contained in:
xiaoguangbin 2024-07-29 19:34:25 +08:00
parent cf369218da
commit f905f9951d
1 changed files with 31 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

52
jeecg-module-spectrum-analysis/src/main/java/org/jeecg/modules/service/impl/SelfStationServiceImpl.java
View File

@ -1778,7 +1778,7 @@ public class SelfStationServiceImpl implements ISelfStationService {
} }
map.put("CToE", fittingParaStr); map.put("CToE", fittingParaStr);
List<Double> channel = new LinkedList<>(); List<Double> channel = new LinkedList<>();
for (int i=0; i<255; ++i) { for (int i=0; i<511; ++i) {
channel.add(Double.valueOf(i)); channel.add(Double.valueOf(i));
} }
List<Double> energy = EnergySpectrumHandler.GetFileFittingData(channel, fittingPara); List<Double> energy = EnergySpectrumHandler.GetFileFittingData(channel, fittingPara);
@ -1804,28 +1804,38 @@ public class SelfStationServiceImpl implements ISelfStationService {
// 高斯函数拟合得到拟合的中心道址，显示Energy vs. Channel // 高斯函数拟合得到拟合的中心道址，显示Energy vs. Channel
//页面展示的表单数据数组 //页面展示的表单数据数组
List<TableWidget> table = new LinkedList<>(); List<TableWidget> table = new LinkedList<>();
for (int i = 0; i < oldScatterSeries.size(); i++) { List<List<Double>> gammaEnergy = (List<List<Double>>) map.get("gammaEnergy");
SeriseData data = oldScatterSeries.get(i); for (int i = 0; i < bElectronEnergy.size(); i++) {
// 获取beta energy
Double be = bElectronEnergy.get(i);
// energy 转化为 channel
Double bc = fittingParaToUi.get(0) + fittingParaToUi.get(1) * be + fittingParaToUi.get(2) * Math.sqrt(be);
// 计算得到gamma energy
Double ge = 661.656 - be;
// todo 需要换成gamma公式
Double gc = 1.51421142654972 + 3.018099280336795 * ge + 2.6690095368611742E-6 * Math.sqrt(ge);
// 通过gamma channel获取gammaGate数据
List<SeriseData> seriseDatas = getGateGamma(5, 4096, gc.intValue(), betaDataFile);
// 开始计算中心道址， 范围：左[C_β-3√(C_β ), C_β+3√(C_β )]右
WeightedObservedPoints points = new WeightedObservedPoints(); WeightedObservedPoints points = new WeightedObservedPoints();
// 通过散点图gamma的channel拿到当前beta横向数据 double left = bc - Math.cbrt(bc);
List<SeriseData> seriseDatas = this.getGateGamma(5, 256, (int) data.getX(), betaDataFile); double right = bc + Math.cbrt(bc);
seriseDatas.forEach(f->{ int sum = 0;
// x = beta channel, y = beta energy for (SeriseData f : seriseDatas) {
double x = f.getX(); if(left < f.getX() && right > f.getX()){
double left = x - Math.cbrt(x); points.add(f.getX(), f.getY());
double right = x + Math.cbrt(x); sum += f.getX();
points.add(x, left, right); }
}); }
// 使用高斯曲线拟合观测点 // GaussianCurveFitter fitter = GaussianCurveFitter.create();
GaussianCurveFitter fitter = GaussianCurveFitter.create(); // fitter.fit(points.toList());
// todo 高斯拟合有问题，暂时使用平均值
double[] parameters = fitter.fit(points.toList()); double v = sum / (right - left);
//表单数据信息 //表单数据信息
TableWidget tableWidget = new TableWidget(); TableWidget tableWidget = new TableWidget();
tableWidget.setRowCount(i+1); tableWidget.setRowCount(i+1);
tableWidget.setChannel(parameters[0]); tableWidget.setChannel(v);
tableWidget.setEnergy(parameters[0]); tableWidget.setEnergy(be);
System.out.println("table:" + parameters[0]);
table.add(tableWidget); table.add(tableWidget);
} }
map.put("tableWidgets", table); map.put("tableWidgets", table);
@ -3336,7 +3346,7 @@ public class SelfStationServiceImpl implements ISelfStationService {
} }
map.put("CToE", fittingParaStr); map.put("CToE", fittingParaStr);
//填充道值数组的数据 //填充道值数组的数据
for (int i=0; i<255; i++) { for (int i=0; i<512; i++) {
channels.add(Double.valueOf(i)); channels.add(Double.valueOf(i));
} }
//计算道值的能量 //计算道值的能量
@ -3556,7 +3566,7 @@ public class SelfStationServiceImpl implements ISelfStationService {
private List<SeriseData> getGateGamma(int channelWidth, int chartHeight, int gammaChannel, SelfStationData betaDataFile){ private List<SeriseData> getGateGamma(int channelWidth, int chartHeight, int gammaChannel, SelfStationData betaDataFile){
List<SeriseData> serise_data = new LinkedList<>(); List<SeriseData> serise_data = new LinkedList<>();
//选择矩形框高度 //选择矩形框高度
Integer flagHeight = channelWidth * (chartHeight/256); Integer flagHeight = channelWidth * (chartHeight/512);
int value = Double.valueOf(flagHeight / 2).intValue(); int value = Double.valueOf(flagHeight / 2).intValue();
//计算得到最高值 //计算得到最高值
int up = gammaChannel - value; int up = gammaChannel - value;