static internal OperationResult GetWarningFromString(this OperationResult result, string msg)
{
result.AddWarningMessage(msg);
result.Status = StatusResult.Warning;
Logger.Instance.WriteLog.Warn(msg);
return result;
}
static internal OperationResult СombineResult(this OperationResult result, OperationResult newResult)
{
if (newResult == null)
{
return result;
}
if (newResult.Status == StatusResult.Warning)
{
result.AddWarningMessage(newResult.WarningMessages);
result.Status = newResult.Status;
}
if (newResult.Status == StatusResult.Fail)
{
result.ErrorMessage = newResult.ErrorMessage;
result.Status = newResult.Status;
}
return result;
}
static void CalculateVerificationResults(this Measurement meas, Multiplicity mkey, MultiplicityCountingRes results)
{
Tuple normal_mass = new Tuple(-1, 0), backup_mass = new Tuple(-1, 0);
try
{
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.CalibrationCurve))
{
// get the current results_cal_curve_rec and cal_curve params
//dev note: the rates as computed by the first and second phases are not yet on ccres, because they exist soley on the counting results MultiplicityCountingRes instance
INCCMethodResults.results_cal_curve_rec ccres = (INCCMethodResults.results_cal_curve_rec)
meas.INCCAnalysisResults.LookupMethodResults(mkey, meas.INCCAnalysisState.Methods.selector, AnalysisMethod.CalibrationCurve, true);
INCCAnalysisParams.cal_curve_rec cal_curve = (INCCAnalysisParams.cal_curve_rec)meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.CalibrationCurve);
INCCAnalysisParams.CalCurveResult status = INCCAnalysisParams.CalCurveResult.Unknown;
if (cal_curve == null)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 10199, "No " + AnalysisMethod.KnownA.FullName() + " method parameters found");
return;
}
if (cal_curve.CalCurveType != INCCAnalysisParams.CalCurveType.HM)
{
Tuple pu240e = new Tuple();
Tuple doubles;
if (cal_curve.cev.useSingles) // the 2009 MTS hack
doubles = new Tuple(results.rates.GetDTCRates(RatesAdjustments.DeadtimeCorrected).Singles);
else
doubles = new Tuple(results.rates.GetDTCRates(RatesAdjustments.DeadtimeCorrected).Doubles);
status = INCCAnalysis.CalculateCalibrationCurveOnly(cal_curve.cev, out pu240e, results.rates, doubles, RatesAdjustments.DeadtimeCorrected); // rates (triples) not used
ccres.pu240e_mass = pu240e;
}
else
{
// get the item id from the acquire record or the ItemId on the measurement itself
// if there is no item id use the empty default item id
// dev note: at some point the acquire record item id becomes a full ItemId record on the measurmeent
//if not NC.App.DB.ItemIdSet.Contains AcquireState.item_id then
// get the default empty one
//end
// from HEAVY_M.cpp
INCCAnalysis.calc_heavy_metal(
cal_curve.heavy_metal_corr_factor,
cal_curve.heavy_metal_reference,
results.rates.DTCRates.Singles,
results.rates.DTCRates.Doubles,
ref ccres.heavy_metal_content,
ref ccres.heavy_metal_correction,
ref ccres.heavy_metal_corr_singles,
ref ccres.heavy_metal_corr_doubles, meas);
status = INCCAnalysis.CalculateCalibrationCurveOnly(cal_curve.cev,
out ccres.pu240e_mass,
results.rates, ccres.heavy_metal_corr_doubles,
RatesAdjustments.DeadtimeCorrected);
ccres.pu240e_mass.v *= meas.MeasurementId.Item.length;
ccres.pu240e_mass.err *= meas.MeasurementId.Item.length;
}
if (status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
string msg = String.Format("Passive calibration curve failed mass limits of {0} and {1}", cal_curve.cev.lower_mass_limit, cal_curve.cev.upper_mass_limit);
meas.AddErrorMessage(msg, 10196, mkey);
}
else if (status != INCCAnalysisParams.CalCurveResult.Success)
{
meas.AddErrorMessage("Passive calibration curve analysis error", 10197, mkey);
}
if (status == INCCAnalysisParams.CalCurveResult.Success || status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
ccres.dcl_pu_mass = meas.AcquireState.mass; // another requirement for the acquire state
meas.Logger.TraceEvent(NCCReporter.LogLevels.Verbose, 10133, "calc_mass/calc_u235_mass are called next");
if (cal_curve.CalCurveType != INCCAnalysisParams.CalCurveType.U)
{
INCCAnalysis.calc_mass(ccres.pu240e_mass,
ref ccres.pu_mass, ref ccres.dcl_pu_mass, ref ccres.dcl_pu240e_mass, ref ccres.dcl_minus_asy_pu_mass, ref ccres.dcl_minus_asy_pu_mass_pct, ref ccres.pass,
meas);
}
else
{
INCCAnalysis.calc_u235_mass(cal_curve.percent_u235, ccres.pu240e_mass,
ref ccres.pu_mass, ref ccres.dcl_pu_mass, ref ccres.dcl_minus_asy_pu_mass, ref ccres.dcl_minus_asy_pu_mass_pct, ref ccres.pass,
meas);
}
if (!ccres.pass)
{
meas.AddWarningMessage("Passive calibration curve: failed stratum rejection limits", 10198, mkey);
}
else
{
meas.AddWarningMessage("Passive calibration curve: passed stratum rejection limits", 10200, mkey);
}
if (ccres.pu240e_mass.v > ccres.methodParams.cev.upper_mass_limit)
{
meas.AddWarningMessage("Passive calibration curve: upper Pu240e mass limit exceeded.", 10210, mkey);
}
if (ccres.pu240e_mass.v < ccres.methodParams.cev.lower_mass_limit)
{
meas.AddWarningMessage("Passive calibration curve: lower Pu240e mass limit exceeded.", 10211, mkey);
}
}
// normal and backup retention
if (meas.INCCAnalysisState.Methods.Normal == AnalysisMethod.CalibrationCurve)
{
normal_mass.CopyFrom(ccres.pu240e_mass);
}
if (meas.INCCAnalysisState.Methods.Backup == AnalysisMethod.CalibrationCurve)
{
backup_mass.CopyFrom(ccres.pu240e_mass);
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.KnownA))
{
INCCMethodResults.results_known_alpha_rec kares = (INCCMethodResults.results_known_alpha_rec)
meas.INCCAnalysisResults.LookupMethodResults(mkey, meas.INCCAnalysisState.Methods.selector, AnalysisMethod.KnownA, true);
INCCAnalysisParams.known_alpha_rec ka_params = (INCCAnalysisParams.known_alpha_rec)meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.KnownA);
if (ka_params == null)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 10199, "No Known alpha method parameters found");
return;
}
bool success = false;
kares.dcl_pu_mass = meas.AcquireState.mass; // dev note: another use of acq, a requirement, here
// copy the input calibration params to the copy on the results rec, to be saved with the KA results
kares.methodParams = new INCCAnalysisParams.known_alpha_rec(ka_params);
if (ka_params.known_alpha_type == INCCAnalysisParams.KnownAlphaVariant.Conventional)
{
INCCMethodResults.results_known_alpha_rec karesdup = INCCAnalysis.CalculateKnownAlpha(mkey, results.rates, meas, RatesAdjustments.DeadtimeCorrected); // rates (triples) not used
if (karesdup != null) // we have the new mass results, and they are preserved in the results map
{
success = true;
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 240, "Known alpha results for pu240E {0} +- {1}", karesdup.pu240e_mass.v, karesdup.pu240e_mass.err);
}
}
else if (ka_params.known_alpha_type == INCCAnalysisParams.KnownAlphaVariant.HeavyMetalCorrection)
{
INCCAnalysis.calc_heavy_metal(
ka_params.heavy_metal_corr_factor,
ka_params.heavy_metal_reference,
results.rates.DTCRates.Singles,
results.rates.DTCRates.Doubles,
ref kares.heavy_metal_content,
ref kares.heavy_metal_correction,
ref kares.corr_singles,
ref kares.corr_doubles, meas);
Rates HMSDRates = new Rates();
HMSDRates.DeadtimeCorrectedRates.Singles.CopyFrom(kares.corr_singles);
HMSDRates.DeadtimeCorrectedRates.Doubles.CopyFrom(kares.corr_doubles);
INCCMethodResults.results_known_alpha_rec karesdup = INCCAnalysis.CalculateKnownAlpha(mkey, HMSDRates, meas, RatesAdjustments.DeadtimeCorrected);
if (karesdup != null) // we have the new mass results, and they are preserved in the results map
{
success = true;
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 240, "Known alpha HM results for pu240E {0} +- {1}", karesdup.pu240e_mass.v, karesdup.pu240e_mass.err);
}
kares.pu240e_mass.v *= meas.MeasurementId.Item.length;
kares.pu240e_mass.err *= meas.MeasurementId.Item.length;
}
else if (ka_params.known_alpha_type == INCCAnalysisParams.KnownAlphaVariant.MoistureCorrAppliedToDryAlpha)
{
success =
INCCAnalysis.calc_known_alpha_moisture_corr(
results.rates.DTCRates.Singles,
results.rates.DTCRates.Doubles,
results.Scaler1, results.Scaler2,
ref kares.corr_singles, /* ring ratio */
ref kares.corr_factor,
ref kares.dry_alpha_or_mult_dbls, /* dry alpha */
ref kares.mult_corr_doubles,
ref kares.mult,
ref kares.alphaK,
ref kares.pu240e_mass, ka_params, meas, mkey);
if (success)
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 240, "Known alpha MoistureCorrAppliedToDryAlpha results for pu240E {0} +- {1}", kares.pu240e_mass.v, kares.pu240e_mass.err);
}
else if (ka_params.known_alpha_type == INCCAnalysisParams.KnownAlphaVariant.MoistureCorrAppliedToMultCorrDoubles)
{
success =
INCCAnalysis.calc_known_alpha_moisture_corr_mult_doubles(
results.rates.DTCRates.Singles,
results.rates.DTCRates.Doubles,
results.Scaler1, results.Scaler2,
ref kares.corr_singles, /* ring ratio */
ref kares.corr_factor,
ref kares.dry_alpha_or_mult_dbls, /* moist mult_corr_doubles */
ref kares.mult_corr_doubles,
ref kares.mult,
ref kares.alphaK,
ref kares.pu240e_mass, ka_params, meas, mkey);
if (success)
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 240, "Known alpha MoistureCorrAppliedToMultCorrDoubles results for pu240E {0} +- {1}", kares.pu240e_mass.v, kares.pu240e_mass.err);
}
if (success)
{
INCCAnalysis.calc_mass(kares.pu240e_mass,
ref kares.pu_mass, ref kares.dcl_pu_mass, ref kares.dcl_pu240e_mass, ref kares.dcl_minus_asy_pu_mass, ref kares.dcl_minus_asy_pu_mass_pct, ref kares.pass,
meas);
}
else
{
meas.AddErrorMessage("Known alpha analysis error", 10199, mkey);
}
if (!kares.pass)
{
meas.AddWarningMessage("Known alpha: failed stratum rejection limits", 10198, mkey);
}
else
{
meas.AddWarningMessage("Known alpha: passed stratum rejection limits", 10200, mkey);
}
if (kares.pu240e_mass.v > kares.methodParams.cev.upper_mass_limit)
{
meas.AddWarningMessage("Known alpha: upper Pu240e mass limit exceeded.", 10210, mkey);
}
if (kares.pu240e_mass.v < kares.methodParams.cev.lower_mass_limit)
{
meas.AddWarningMessage("Known alpha: lower Pu240e mass limit exceeded.", 10211, mkey);
}
// normal and backup retention
if (meas.INCCAnalysisState.Methods.Normal == AnalysisMethod.KnownA)
{
normal_mass.CopyFrom(kares.pu240e_mass);
}
if (meas.INCCAnalysisState.Methods.Backup == AnalysisMethod.KnownA)
{
backup_mass.CopyFrom(kares.pu240e_mass);
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.Multiplicity))
{
double error = 0.0;
INCCAnalysisParams.multiplicity_rec mul_param = (INCCAnalysisParams.multiplicity_rec)meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.Multiplicity);
if (mul_param == null)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 10198, "No Multiplicity method parameters found");
return;
}
INCCMethodResults.results_multiplicity_rec mmres = (INCCMethodResults.results_multiplicity_rec)meas.INCCAnalysisResults.LookupMethodResults(
mkey, meas.INCCAnalysisState.Methods.selector, AnalysisMethod.Multiplicity, true);
// weird rates (triples) used, but they are wrong (see note line 568 avg_sums.cs)
INCCMethodResults.results_multiplicity_rec mmresdup = INCCAnalysis.CalculateMultiplicity(mkey, results.covariance_matrix, results.rates.GetDTCRates(RatesAdjustments.DeadtimeCorrected), meas, RatesAdjustments.DeadtimeCorrected);
if (mmresdup != null) // we have the new mass results, and they are preserved in the results map
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 240, "Multiplicity results for pu240E {0} +- {1}", mmres.pu240e_mass.v, mmres.pu240e_mass.err);
if (meas.AcquireState.acquire_type == AcquireConvergence.Pu240EffPrecision)
{
if (mmres.pu240e_mass.v != 0.0)
{
error = mmres.pu240e_mass.err /
mmres.pu240e_mass.v * 100.0;
if (error > meas.AcquireState.meas_precision)
{
meas.AddWarningMessage(String.Format("Multiplicity: Pu240e error = {0}%", error), 10198, mkey);
}
}
}
mmres.dcl_pu_mass = meas.AcquireState.mass; // another use of acq, a requirement, here
INCCAnalysis.calc_mass(mmres.pu240e_mass,
ref mmres.pu_mass, ref mmres.dcl_pu_mass, ref mmres.dcl_pu240e_mass, ref mmres.dcl_minus_asy_pu_mass, ref mmres.dcl_minus_asy_pu_mass_pct, ref mmres.pass,
meas);
if (!mmres.pass)
{
meas.AddWarningMessage("Multiplicity: failed stratum rejection limits", 10198, mkey);
}
else
{
meas.AddWarningMessage("Multiplicity: passed stratum rejection limits", 10200, mkey);
}
if (mul_param.solve_efficiency == INCCAnalysisParams.MultChoice.CONVENTIONAL_MULT_WEIGHTED) // todo: implement Weighted
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 36010, "CONVENTIONAL_MULT_WEIGHTED Multiplicity measurement results");
}
else if (mul_param.solve_efficiency == INCCAnalysisParams.MultChoice.MULT_DUAL_ENERGY_MODEL) // todo: implement DE
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 36010, "MULT_DUAL_ENERGY_MODEL Multiplicity measurement results");
}
// normal and backup retention
if (meas.INCCAnalysisState.Methods.Normal == AnalysisMethod.Multiplicity)
{
normal_mass.CopyFrom(mmres.pu240e_mass);
}
if (meas.INCCAnalysisState.Methods.Backup == AnalysisMethod.Multiplicity)
{
backup_mass.CopyFrom(mmres.pu240e_mass);
}
}
else
{
meas.AddErrorMessage("Multiplicity analysis error", 10198, mkey);
}
// copy the input calib to the results rec
mmres.methodParams = new INCCAnalysisParams.multiplicity_rec(mul_param);
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.KnownM))
{
INCCMethodResults.results_known_m_rec kmres = INCCAnalysis.CalculateKnownM(mkey, results, meas, RatesAdjustments.DeadtimeCorrected);
// calc mass
if (kmres != null) // you have calculated well my child
{
INCCAnalysis.calc_mass(kmres.pu240e_mass,
ref kmres.pu_mass, ref kmres.dcl_pu_mass, ref kmres.dcl_pu240e_mass,
ref kmres.dcl_minus_asy_pu_mass, ref kmres.dcl_minus_asy_pu_mass_pct, ref kmres.pass, meas);
if (!kmres.pass)
{
meas.AddWarningMessage("Known M: failed stratum rejection limits", 10198, mkey);
}
else
{
meas.AddWarningMessage("Known M: passed stratum rejection limits", 10200, mkey);
}
if (kmres.pu240e_mass.v > kmres.methodParams.upper_mass_limit)
{
meas.AddWarningMessage("Known M: upper Pu240e mass limit exceeded.", 10210, mkey);
}
if (kmres.pu240e_mass.v < kmres.methodParams.lower_mass_limit)
{
meas.AddWarningMessage("Known M: lower Pu240e mass limit exceeded.", 10211, mkey);
}
// normal and backup retention
if (meas.INCCAnalysisState.Methods.Normal == AnalysisMethod.KnownM)
{
normal_mass.CopyFrom(kmres.pu240e_mass);
}
if (meas.INCCAnalysisState.Methods.Backup == AnalysisMethod.KnownM)
{
backup_mass.CopyFrom(kmres.pu240e_mass);
}
}
else
{
meas.AddErrorMessage("Known M: analysis error", 10198, mkey);
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.Active))
{
INCCAnalysisParams.active_rec act_param = (INCCAnalysisParams.active_rec)meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.Active);
if (act_param == null)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 10198, "No Active method parameters found");
return;
}
INCCMethodResults.results_active_rec actres = (INCCMethodResults.results_active_rec)meas.INCCAnalysisResults.LookupMethodResults(
mkey, meas.INCCAnalysisState.Methods.selector, AnalysisMethod.Active, true);
/* calculate active doubles rate corrected for source yield factor */
//line 331 if calc_asy.cpp
Measurement.SourceYieldFactoredRates syfr = new Measurement.SourceYieldFactoredRates(results, meas);
// line 1267 of calc_asy.cpp
//Martyn says we need stuff here to deal with Cf active measurements HN 7.23.2015
actres.k0.v = syfr.source_yield_factor;
actres.k = new Tuple(syfr.total_corr_fact);
actres.k1 = new Tuple(meas.Norm.currNormalizationConstant);
// This stays the same for Cf. HN 7.23.2015
INCCAnalysisParams.CalCurveResult status = INCCAnalysis.CalculateCalibrationCurveOnly(act_param.cev,
out actres.u235_mass, results.rates, syfr.corrected_doubles,
RatesAdjustments.DeadtimeCorrected);
if (status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
string msg = String.Format("Active calibration curve failed mass limits of {0} and {1}", act_param.cev.lower_mass_limit, act_param.cev.upper_mass_limit);
meas.AddErrorMessage(msg, 10196, mkey);
}
else if (status != INCCAnalysisParams.CalCurveResult.Success)
{
meas.AddErrorMessage("Active calibration curve analysis error", 10197, mkey);
}
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 240, "Active results for U235 {0} +- {1}", actres.u235_mass.v, actres.u235_mass.err);
if (status == INCCAnalysisParams.CalCurveResult.Success || status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
actres.dcl_u235_mass = meas.AcquireState.mass;
INCCAnalysis.calc_decl_minus_assay_u235(actres.u235_mass, actres.dcl_u235_mass, ref actres.dcl_minus_asy_u235_mass, ref actres.dcl_minus_asy_u235_mass_pct, ref actres.pass, meas);
if (!actres.pass)
{
meas.AddWarningMessage("Active calibration curve: failed stratum rejection limits", 10198, mkey);
}
else if (!meas.Stratum.Unset)
{
meas.AddWarningMessage("Active calibration curve: passed stratum rejection limits", 10200, mkey);
}
if (actres.u235_mass.v > actres.methodParams.cev.upper_mass_limit)
{
meas.AddWarningMessage("Active calibration curve: upper U235 mass limit exceeded.", 10210, mkey);
}
if (actres.u235_mass.v < actres.methodParams.cev.lower_mass_limit)
{
meas.AddWarningMessage("Active calibration curve: lower U235 mass limit exceeded.", 10211, mkey);
}
}
// normal and backup retention not performed for active
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.ActiveMultiplicity))
{
// sets results class' mult v,err values at end of calculation
INCCMethodResults.results_active_mult_rec res = INCCAnalysis.CalculateActiveMultiplicity(mkey, results, meas, RatesAdjustments.DeadtimeCorrected);
if (res == null)
{
meas.AddErrorMessage("Active multiplicity analysis error", 10152, mkey);
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.TruncatedMultiplicity))
{
INCCMethodResults.results_truncated_mult_rec res = INCCAnalysis.CalculateTruncatedMult(mkey, results, meas, RatesAdjustments.DeadtimeCorrected);
// normal and backup retention
if (meas.INCCAnalysisState.Methods.Normal == AnalysisMethod.TruncatedMultiplicity)
{
if (res.methodParams.known_eff)
normal_mass.CopyFrom(res.k.pu240e_mass);
else
normal_mass.CopyFrom(res.s.pu240e_mass);
}
if (meas.INCCAnalysisState.Methods.Backup == AnalysisMethod.TruncatedMultiplicity)
{
if (res.methodParams.known_eff)
backup_mass.CopyFrom(res.k.pu240e_mass);
else
backup_mass.CopyFrom(res.s.pu240e_mass);
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.Collar))
{
meas.AddWarningMessage("Collar mass results", 10153, mkey); // NEXT: Collar is incomplete, new design from IAEA is pending, this is a big task
INCCAnalysis.CalculateCollar(mkey, results, meas, RatesAdjustments.DeadtimeCorrected);
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.ActivePassive))
{
INCCAnalysisParams.active_passive_rec act_param = (INCCAnalysisParams.active_passive_rec)meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.ActivePassive);
if (act_param == null)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 10198, "No active/passive method parameters found");
return;
}
INCCMethodResults.results_active_passive_rec actres = (INCCMethodResults.results_active_passive_rec)meas.INCCAnalysisResults.LookupMethodResults(
mkey, meas.INCCAnalysisState.Methods.selector, AnalysisMethod.ActivePassive, true);
INCCAnalysis.CalculateActivePassive(mkey, results, meas, RatesAdjustments.DeadtimeCorrected);
//line 1134 if calc_asy.cpp
// calculate delta doubles and error from passive and active doubles
Measurement.SourceYieldFactoredRates syfr = new Measurement.SourceYieldFactoredRates(results, meas);
actres.k0.v = syfr.source_yield_factor;
actres.k = new Tuple(syfr.total_corr_fact);
actres.k1 = new Tuple(meas.Norm.currNormalizationConstant);
actres.delta_doubles.v = syfr.corrected_doubles.v - results.DeadtimeCorrectedDoublesRate.v;
actres.delta_doubles.err = Math.Sqrt((syfr.corrected_doubles.v * results.DeadtimeCorrectedDoublesRate.err) + (syfr.corrected_doubles.v * results.DeadtimeCorrectedDoublesRate.err));
INCCAnalysisParams.CalCurveResult status = INCCAnalysis.CalculateCalibrationCurveOnly(act_param.cev,
out actres.u235_mass, results.rates, actres.delta_doubles,
RatesAdjustments.DeadtimeCorrected);
if (status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
string msg = String.Format("Active/passive calibration curve failed mass limits of {0} and {1}", act_param.cev.lower_mass_limit, act_param.cev.upper_mass_limit);
meas.AddErrorMessage(msg, 10196, mkey);
}
else if (status != INCCAnalysisParams.CalCurveResult.Success)
{
meas.AddErrorMessage("Active/passive calibration curve analysis error", 10197, mkey);
}
if (status == INCCAnalysisParams.CalCurveResult.Success || status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
actres.dcl_u235_mass = meas.AcquireState.mass;
INCCAnalysis.calc_decl_minus_assay_u235(actres.u235_mass, actres.dcl_u235_mass, ref actres.dcl_minus_asy_u235_mass, ref actres.dcl_minus_asy_u235_mass_pct, ref actres.pass, meas);
if (!actres.pass)
{
meas.AddWarningMessage("Active/passive: failed stratum rejection limits", 10198, mkey);
}
else if (!meas.Stratum.Unset)
{
meas.AddWarningMessage("Active/passive: passed stratum rejection limits", 10200, mkey);
}
if (actres.u235_mass.v > actres.methodParams.cev.upper_mass_limit)
{
meas.AddWarningMessage("Active/passive: upper U235 mass limit exceeded.", 10210, mkey);
}
if (actres.u235_mass.v < actres.methodParams.cev.lower_mass_limit)
{
meas.AddWarningMessage("Active/passive: lower U235 mass limit exceeded.", 10211, mkey);
}
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.AddASource))
{
INCCAnalysisParams.CalCurveResult status;
INCCMethodResults.results_add_a_source_rec res = null;
// gotta do a sanity check for the AAS cycles, they may not be there due to unfinished processing in the overall code
if (meas.CFCycles != null)
{
res = INCCAnalysis.CalculateAddASource(mkey, results, meas, RatesAdjustments.DeadtimeCorrected, out status);
}
else
{
meas.AddErrorMessage("Add-a-source calibration curve cycles not present error", 10197, mkey);
status = INCCAnalysisParams.CalCurveResult.EpicFailLOL;
}
if (status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
string msg = String.Format("Add-a-source calibration curve failed mass limits of {0} and {1}", res.methodParams.cev.lower_mass_limit, res.methodParams.cev.upper_mass_limit);
meas.AddErrorMessage(msg, 10196, mkey);
}
else if (status != INCCAnalysisParams.CalCurveResult.Success)
{
meas.AddErrorMessage("Add-a-source calibration curve analysis error", 10197, mkey);
}
if (status == INCCAnalysisParams.CalCurveResult.Success || status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
res.dcl_pu_mass = meas.AcquireState.mass;
INCCAnalysis.calc_mass(res.pu240e_mass,
ref res.pu_mass, ref res.dcl_pu_mass, ref res.dcl_pu240e_mass,
ref res.dcl_minus_asy_pu_mass, ref res.dcl_minus_asy_pu_mass_pct, ref res.pass, meas);
if (!res.pass)
{
meas.AddWarningMessage("Add-a-source calibration curve: failed stratum rejection limits", 10198, mkey);
}
else
{
meas.AddWarningMessage("Add-a-source calibration curve: passed stratum rejection limits", 10200, mkey);
}
if (res.pu240e_mass.v > res.methodParams.cev.upper_mass_limit)
{
meas.AddWarningMessage("Add-a-source calibration curve: upper Pu240e mass limit exceeded.", 10210, mkey);
}
if (res.pu240e_mass.v < res.methodParams.cev.lower_mass_limit)
{
meas.AddWarningMessage("Add-a-source calibration curve: lower Pu240e mass limit exceeded.", 10211, mkey);
}
// normal and backup retention
if (meas.INCCAnalysisState.Methods.Normal == AnalysisMethod.AddASource)
{
normal_mass.CopyFrom(res.pu240e_mass);
}
if (meas.INCCAnalysisState.Methods.Backup == AnalysisMethod.AddASource)
{
backup_mass.CopyFrom(res.pu240e_mass);
}
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.CuriumRatio))
{
INCCAnalysisParams.CalCurveResult status;
INCCMethodResults.results_curium_ratio_rec res = INCCAnalysis.CalculateCuriumRatio(mkey, results, meas, RatesAdjustments.DeadtimeCorrected, out status);
if (status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
string msg = String.Format("Curium ratio calibration curve failed mass limits of {0} and {1}", res.methodParams.cev.lower_mass_limit, res.methodParams.cev.upper_mass_limit);
meas.AddErrorMessage(msg, 10196, mkey);
}
else if (status != INCCAnalysisParams.CalCurveResult.Success)
{
meas.AddErrorMessage("Curium ratio calibration curve analysis error", 10197, mkey);
}
if (status == INCCAnalysisParams.CalCurveResult.Success || status == INCCAnalysisParams.CalCurveResult.FailedOnMassLimit)
{
INCCAnalysisParams.cm_pu_ratio_rec cm_pu_ratio = NC.App.DB.Cm_Pu_RatioParameters.Get(); // load from DB, just like test params,
// dev note: better not to ref DB here, because this is a one-off state retrieval and no other DB access occurs during mass calc processing, but that is how it works this morning
//calc curium mass
INCCAnalysis.calc_curium_mass(res, cm_pu_ratio, meas);
res.u.dcl_mass = cm_pu_ratio.cm_dcl_u_mass;
res.u235.dcl_mass = cm_pu_ratio.cm_dcl_u235_mass;
if (!res.pu.pass)
{
meas.AddWarningMessage("Curium ratio: Pu failed stratum rejection limits", 10198, mkey);
}
else
{
meas.AddWarningMessage("Curium ratio: Pu passed stratum rejection limits", 10200, mkey);
}
if (!res.u.pass)
{
meas.AddWarningMessage("Curium ratio: U failed stratum rejection limits", 10198, mkey);
}
else
{
meas.AddWarningMessage("Curium ratio: U passed stratum rejection limits", 10200, mkey);
}
// normal and backup retention
if (meas.INCCAnalysisState.Methods.Normal == AnalysisMethod.CuriumRatio)
{
normal_mass.CopyFrom(res.pu.pu240e_mass);
}
if (meas.INCCAnalysisState.Methods.Backup == AnalysisMethod.CuriumRatio)
{
backup_mass.CopyFrom(res.pu.pu240e_mass);
}
}
}
// annotate the final results method marker on the INCC results instance.
INCCMethodResults imr = null;
if (normal_mass.v != -1.0)
{
bool got = meas.INCCAnalysisResults.TryGetINCCResults(mkey, out imr);
if (got)
imr.primaryMethod = meas.INCCAnalysisState.Methods.Normal;
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10100, "Verification primary method {0} with mass {1} is from the normal method", imr.primaryMethod.ToString(), normal_mass.v);
}
else if (backup_mass.v != -1.0)
{
bool got = meas.INCCAnalysisResults.TryGetINCCResults(mkey, out imr);
if (got)
imr.primaryMethod = meas.INCCAnalysisState.Methods.Backup;
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10101, "Verification primary method {0} with mass {1} is from the backup method", imr.primaryMethod.ToString(), backup_mass.v);
}
if ((normal_mass.v != -1.0) && (backup_mass.v != -1.0))
{
if (meas.INCCAnalysisState.Methods.Backup.IsNone())
{
double delta = Math.Abs(normal_mass.v - backup_mass.v);
double delta_error = Math.Sqrt(normal_mass.err * normal_mass.err +
backup_mass.err * backup_mass.err);
bool got = meas.INCCAnalysisResults.TryGetINCCResults(mkey, out imr);
if (delta <= (delta_error * meas.Tests.normalBackupAssayTestLimit))
{
if (got) imr.primaryMethod = meas.INCCAnalysisState.Methods.Normal;
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10102, "Verification primary method {0} with masses {1} and {2} is from the normal method", imr.primaryMethod.ToString(), normal_mass.v, backup_mass.v);
}
else
{
if (got) imr.primaryMethod = meas.INCCAnalysisState.Methods.Backup;
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10103, "Verification primary method {0} with masses {1} and {2} is from the backup method", imr.primaryMethod.ToString(), normal_mass.v, backup_mass.v);
}
}
}
}
catch (Exception e)
{
meas.Logger.TraceException(e);
}
}
/// <summary>
/// after each cycle, run this cycle limit test and exit code
/// </summary>
/// <param name="seq"></param>
/// <param name="mkey"></param>
public static void CycleStatusTerminationCheck(this Measurement meas, Cycle cc)
{
bool stopAndComputeResults = meas.AcquireState.lm.SaveOnTerminate;
IEnumerator iter = cc.CountingAnalysisResults.GetMultiplicityEnumerator();
while (iter.MoveNext())
{
Multiplicity mkey = (Multiplicity)((KeyValuePair<SpecificCountingAnalyzerParams, object>)(iter.Current)).Key;
MultiplicityCountingRes mcr = (MultiplicityCountingRes)((KeyValuePair<SpecificCountingAnalyzerParams, object>)(iter.Current)).Value;
if (meas.MeasCycleStatus.num_checksum_failures > meas.Tests.maxNumFailures)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 7071, "Maximum checksum failure count met, cycle {0} {1}", cc.seq, mkey);
if (stopAndComputeResults)
{
// URGENT: set state to end of measurement so that state moves out of cycles and skips forward to CalculateMeasurementResults
}
else
{
// exit processing completely, do not compute or save results
NC.App.Opstate.Abort();
}
}
if (meas.MeasCycleStatus.num_acc_sngl_failures > meas.Tests.maxNumFailures)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 7072, "Maximum A/S failure count met, cycle {0} {1}", cc.seq, mkey);
// ditto above comment block
}
if (meas.MeasCycleStatus.num_high_voltage_failures > meas.Tests.maxNumFailures)
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Warning, 7073, "Maximum checksum failure count met, cycle {0} {1}", cc.seq, mkey);
// ditto above comment block
}
if ((meas.MeasCycleStatus.num_acc_sngl_failures >= NUM_ACC_SNGL_WARNING) &&
(!meas.MeasCycleStatus.acc_sngl_warning_sent))
{
meas.AddWarningMessage("At least " + meas.MeasCycleStatus.num_acc_sngl_failures + " A/S test failures", 7074, mkey);
meas.MeasCycleStatus.acc_sngl_warning_sent = true;
}
if ((meas.MeasCycleStatus.num_checksum_failures >= NUM_CHECKSUM_WARNING) &&
(!meas.MeasCycleStatus.checksum_warning_sent))
{
meas.AddWarningMessage("At least " + meas.MeasCycleStatus.num_checksum_failures + " checksum failures", 7075, mkey);
meas.MeasCycleStatus.checksum_warning_sent = true;
}
if ((meas.MeasCycleStatus.num_high_voltage_failures >= NUM_HIGH_VOLTAGE_WARNING) &&
(!meas.MeasCycleStatus.high_voltage_warning_sent))
{
meas.AddWarningMessage("At least " + meas.MeasCycleStatus.num_high_voltage_failures + " high voltage failures", 7076, mkey);
meas.MeasCycleStatus.high_voltage_warning_sent = true;
}
if (meas.AcquireState.acquire_type == AcquireConvergence.CycleCount)
{
meas.MeasCycleStatus.acquire_num_runs = meas.MeasCycleStatus.initial_num_runs +
meas.MeasCycleStatus.num_outlier_failures + meas.MeasCycleStatus.num_acc_sngl_failures +
meas.MeasCycleStatus.num_checksum_failures;
}
else if (meas.Cycles.GetValidCycleCountForThisKey(mkey) < 100) // an arbitrary limit, make it a config item
{
// URGENT: set state to end of measurement so that state moves out of cycles and skips forward to CalculateMeasurementResults
// todo: see take_data.cpp line 733 for mass check processing steps
if (meas.Cycles.GetValidCycleCountForThisKey(mkey) > 1)
{
INCCResult results;
MeasOptionSelector ar = new MeasOptionSelector(meas.MeasOption, mkey);
bool found = meas.INCCAnalysisResults.TryGetValue(ar, out results);
double error = 0;
if (meas.AcquireState.acquire_type == AcquireConvergence.DoublesPrecision)
{
if (results.rates.DTCRates.Doubles.v != 0.0)
{
error = results.rates.DTCRates.Doubles.err / results.rates.DTCRates.Doubles.v * 100.0;
}
}
else if (meas.AcquireState.acquire_type == AcquireConvergence.TriplesPrecision)
{
if (results.rates.DTCRates.Triples.v != 0.0)
{
error = results.rates.DTCRates.Triples.err / results.rates.DTCRates.Triples.v * 100.0;
}
}
else if (meas.AcquireState.acquire_type == AcquireConvergence.Pu240EffPrecision)
{
RatesAdjustments dtchoice = RatesAdjustments.Raw; // NEXT: design out this choice everywhere, it shoud be a flag on the current computational state somewhere
INCCMethodResults.results_multiplicity_rec res = INCCAnalysis.CalculateMultiplicity(mkey, results.covariance_matrix, results.rates.GetDTCRates(dtchoice), meas, dtchoice);
if (res != null)
{
// the results obtained within the method are already stored on the result and they do not need to be copied here (VERIFY THIS)
if (res.pu240e_mass.v != 0.0)
{
error = res.pu240e_mass.err /
res.pu240e_mass.v * 100.0;
}
else
error = 0.0;
}
}
meas.MeasCycleStatus.acquire_num_runs = (uint)((double)cc.seq * (error / meas.AcquireState.meas_precision) * (error / meas.AcquireState.meas_precision)) + 1;
meas.MeasCycleStatus.acquire_num_runs += meas.MeasCycleStatus.num_outlier_failures +
meas.MeasCycleStatus.num_acc_sngl_failures + meas.MeasCycleStatus.num_checksum_failures;
}
else
meas.MeasCycleStatus.acquire_num_runs = meas.AcquireState.max_num_runs;
if (meas.MeasCycleStatus.acquire_num_runs > meas.AcquireState.max_num_runs)
meas.MeasCycleStatus.acquire_num_runs = meas.AcquireState.max_num_runs;
if (meas.MeasCycleStatus.acquire_num_runs < meas.AcquireState.min_num_runs)
meas.MeasCycleStatus.acquire_num_runs = meas.AcquireState.min_num_runs;
}
}
}
/// <summary>
/// Do the third phase . . .
/// From calc_res.cpp,
/// </summary>
public static void CalculateResults(this Measurement meas)
{
IEnumerator iter = meas.CountingAnalysisResults.GetMultiplicityEnumerator();
while (iter.MoveNext())
{
Multiplicity mkey = (Multiplicity)((KeyValuePair<SpecificCountingAnalyzerParams, object>)(iter.Current)).Key;
if (NC.App.Opstate.IsAbortRequested)
return;
INCCResult results;
MeasOptionSelector ar = new MeasOptionSelector(meas.MeasOption, mkey);
bool found = meas.INCCAnalysisResults.TryGetValue(ar, out results);
if (!found)
{
meas.AddErrorMessage("No results available", 10151, mkey);
return;
}
/* if using measure to precision, and max # runs reached, then add warning message indicating actual precision reached. */
if (meas.AcquireState.acquire_type == AcquireConvergence.DoublesPrecision)
{
if (results.rates.DTCRates.DoublesRate != 0.0)
{
double error = results.rates.DTCRates.Doubles.err / results.rates.DTCRates.DoublesRate * 100.0;
if (error > meas.AcquireState.meas_precision)
{
meas.AddWarningMessage("Measurement doubles error = " + error.ToString("F2"), 10126, mkey);
}
}
}
else if (meas.AcquireState.acquire_type == AcquireConvergence.TriplesPrecision)
{
if (results.rates.DTCRates.TriplesRate != 0.0)
{
double error = results.rates.DTCRates.Triples.err / results.rates.DTCRates.TriplesRate * 100.0;
if (error > meas.AcquireState.meas_precision)
{
meas.AddWarningMessage("Measurement triples error = " + error.ToString("F2"), 10127, mkey);
}
}
}
if (meas.MeasOption == AssaySelector.MeasurementOption.rates) // Doug Requirement #2
{
// for a rates only measurement, all done!
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10180, "Rates Only measurement complete");
continue;
}
try
{
switch (meas.MeasOption)
{
case AssaySelector.MeasurementOption.background: // Doug Requirement #3
// a bkg is an average over a bunch of cycles with only the deadtime correction applied,
// but note the variant for Truncated Mult background, where an additional calculation is made
if (results.rates.DTCRates.DoublesRate > meas.Tests.bkgDoublesRateLimit)
meas.AddWarningMessage("Background doubles rate " + results.rates.DTCRates.DoublesRate + " greater than " + meas.Tests.bkgDoublesRateLimit, 10141, mkey);
if (results.rates.DTCRates.TriplesRate > meas.Tests.bkgTriplesRateLimit)
meas.AddWarningMessage("Background triples rate " + results.rates.DTCRates.TriplesRate + " greater than " + meas.Tests.bkgTriplesRateLimit, 10142, mkey);
// dev note: this is a background measurement per se, so we copy the results to the Background class and store the entire thing that way
// NEXT: Need to account for Active bkg beginning here.
if (Integ.GetCurrentAcquireParams().well_config == WellConfiguration.Active)// Is an active acquisition.
{
meas.Background.INCCActive.CopyFrom(results.rates.DeadtimeCorrectedRates);
}
else
meas.Background.CopyFrom(results.rates);
// maybe if (INCCAnalysisState.Methods.Has(AnalysisMethod.TruncatedMultiplicity))
if (meas.Background.TMBkgParams.ComputeTMBkg)
// Trunc Mult Bkg step, calc_tm_rates, sets TM bkg rates on Measurement.Background
INCCAnalysis.calc_tm_rates(mkey, results, meas, meas.Background.TMBkgParams, meas.Detector.Id.SRType);
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10181, "Background measurement complete");
break;
case AssaySelector.MeasurementOption.initial:
INCCResults.results_init_src_rec results_init_src = (INCCResults.results_init_src_rec)results;
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10182, "Calculating Initial source measurement results");
bool funhappy = INCCAnalysis.initial_source_meas(meas, mkey, RatesAdjustments.DeadtimeCorrected);
if (!funhappy || !results_init_src.pass)
{
meas.AddWarningMessage("Initial source measurement failed", 10123, mkey);
}
// on fail, only the relevant results_init_src_rec is saved
// on pass, the normalization parameters are modified with the results_init_src_rec results, and so both are udpated.
break;
case AssaySelector.MeasurementOption.normalization:
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10183, "Calculating Normalization measurement results");
bool happyfun = INCCAnalysis.bias_test(meas, mkey, RatesAdjustments.DeadtimeCorrected);
if (!happyfun)
{
meas.AddWarningMessage("Normalization test -- data quality is inadequate", 10124, mkey);
}
break;
case AssaySelector.MeasurementOption.precision:
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10184, "Calculating Precision measurement results");
bool charmyfun = INCCAnalysis.precision_test(meas, mkey, RatesAdjustments.DeadtimeCorrected);
if (!charmyfun)
{
meas.AddWarningMessage("Precision test failed", 10125, mkey);
}
break;
case AssaySelector.MeasurementOption.calibration: // from calc_res.cpp
if (meas.INCCAnalysisState.Methods.CalibrationAnalysisSelected())
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10185, "Calculating Calibration measurement results");
// dev note: since the analysis routines have similar signatures, design a class OK?
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.CalibrationCurve))
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10191, "Calculating " + AnalysisMethod.CalibrationCurve.FullName() + " measurement results");
// get the current results_cal_curve_rec
INCCMethodResults.results_cal_curve_rec ccres = (INCCMethodResults.results_cal_curve_rec)
meas.INCCAnalysisResults.LookupMethodResults(mkey, meas.INCCAnalysisState.Methods.selector, AnalysisMethod.CalibrationCurve, true);
INCCAnalysisParams.cal_curve_rec cal_curve = (INCCAnalysisParams.cal_curve_rec)meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.CalibrationCurve);
ccres.methodParams = new INCCAnalysisParams.cal_curve_rec(cal_curve);
double pu_mass = ccres.pu_mass.v;
meas.Isotopics.UpdateDeclaredPuMass(meas.MeasDate, ref pu_mass, INCCAnalysis.INCCParity);
ccres.pu_mass.v = pu_mass;
ccres.pu_mass.err = 0.0;
ccres.pu240e_mass.err = 0.0;
double pu240e_mass = ccres.pu240e_mass.v;
meas.Isotopics.calc_pu240e(pu_mass, out pu240e_mass, meas);
ccres.pu240e_mass.v = pu240e_mass;
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.KnownA))
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10192, "Calculating " + AnalysisMethod.KnownA.FullName() + " calibration results");
INCCMethodResults.results_known_alpha_rec kares = INCCAnalysis.CalculateKnownAlpha(mkey, results.rates, meas, RatesAdjustments.DeadtimeCorrected);
if (kares == null)
{
meas.AddErrorMessage("Known alpha analysis error", 10153, mkey);
}
else
{
kares.pu_mass.v = meas.AcquireState.mass;
double pu_mass = kares.pu_mass.v;
meas.Isotopics.UpdateDeclaredPuMass(meas.MeasDate, ref pu_mass, INCCAnalysis.INCCParity);
kares.pu_mass.v = pu_mass;
kares.pu_mass.err = 0.0;
kares.pu240e_mass.err = 0.0;
double pu240e_mass = kares.pu240e_mass.v;
meas.Isotopics.calc_pu240e(pu_mass, out pu240e_mass, meas);
kares.pu240e_mass.v = pu240e_mass;
}
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.Active))
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10192, "Arranging " + AnalysisMethod.Active.FullName() + " calibration results");
INCCSelector sel = new INCCSelector(meas.INCCAnalysisState.Methods.selector);
INCCAnalysisParams.active_rec act;
INCCAnalysisParams.INCCMethodDescriptor surr = meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.Active);
if (surr == null)
{
act = new INCCAnalysisParams.active_rec();
meas.INCCAnalysisState.Methods.AddMethod(AnalysisMethod.Active, act);
}
else
{
act = (INCCAnalysisParams.active_rec)surr;
}
INCCMethodResults.results_active_rec res;
res = (INCCMethodResults.results_active_rec)meas.INCCAnalysisResults.LookupMethodResults(mkey, sel, AnalysisMethod.Active, true);
// all this does is copy the declared mass over to the results, WTF
res.u235_mass.v = meas.AcquireState.mass;
}
if (meas.INCCAnalysisState.Methods.Has(AnalysisMethod.AddASource))
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10152, "Calculating " + AnalysisMethod.AddASource.FullName() + " calibration results");
INCCSelector sel = new INCCSelector(meas.INCCAnalysisState.Methods.selector);
INCCAnalysisParams.add_a_source_rec aas;
INCCMethodResults.results_add_a_source_rec acres;
INCCAnalysisParams.INCCMethodDescriptor surr = meas.INCCAnalysisState.Methods.GetMethodParameters(AnalysisMethod.AddASource);
if (surr == null)
{
aas = new INCCAnalysisParams.add_a_source_rec();
meas.INCCAnalysisState.Methods.AddMethod(AnalysisMethod.AddASource, aas);
}
else
{
aas = (INCCAnalysisParams.add_a_source_rec)surr;
}
acres = (INCCMethodResults.results_add_a_source_rec)meas.INCCAnalysisResults.LookupMethodResults(mkey, sel, AnalysisMethod.AddASource, true);
acres.pu_mass.v = meas.AcquireState.mass;
// update_declared_pu_mass
double pu_mass = acres.pu_mass.v;
meas.Isotopics.UpdateDeclaredPuMass(meas.MeasDate, ref pu_mass, INCCAnalysis.INCCParity);
acres.pu_mass.v = pu_mass;
acres.pu_mass.err = 0.0;
acres.pu240e_mass.err = 0.0;
double pu240e_mass = acres.pu240e_mass.v;
meas.Isotopics.calc_pu240e(pu_mass, out pu240e_mass, meas);
acres.pu240e_mass.v = pu240e_mass;
}
}
else
{
meas.AddErrorMessage("No Calibration analysis methods selected", 10153, mkey);
}
break;
case AssaySelector.MeasurementOption.verification:
// see INCC calc_asy.cpp
// dev note: check for item in the item table, make sure to place this item id on the MeasurementId.item property
if (!string.IsNullOrEmpty(meas.AcquireState.item_id))
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10194, "Using item id '{0}'", meas.AcquireState.item_id);
}
else
meas.Logger.TraceEvent(NCCReporter.LogLevels.Verbose, 10194, "No item id");
if (meas.INCCAnalysisState.Methods.VerificationAnalysisSelected())
{
meas.Logger.TraceEvent(NCCReporter.LogLevels.Info, 10186, "Calculating {0} measurement results", meas.MeasOption.PrintName());
meas.CalculateVerificationResults(mkey, results);
}
else
{
meas.AddErrorMessage("No analysis methods selected", 10156, mkey);
}
break;
case AssaySelector.MeasurementOption.holdup: // NEXT: Hold-up held up, implement it #35
meas.Logger.TraceEvent(NCCReporter.LogLevels.Error, 10187, "Holdup analysis unsupported");
break;
}
}
catch (Exception e)
{
meas.Logger.TraceException(e);
}
}
}
