private void MakeED()
{
//Fill Rho array
try
{
RhoCalc.IsOneSigma = Holdsigma.Checked;
RhoCalc.ZOffset = Zoffset.ToDouble();
RhoCalc.GetSubRoughness = SubRough.ToDouble();
RhoCalc.BoxCount = BoxCount.ToInt();
RhoCalc.SubphaseSLD = SubphaseSLD.ToDouble();
RhoCalc.SuperphaseSLD = SupSLDTB.ToDouble();
//Blank our Rho data from the previous iteration
RhoCalc.RhoArray.Clear();
RhoCalc.LengthArray.Clear();
RhoCalc.SigmaArray.Clear();
BoxRhoArray.ForEach(p => RhoCalc.RhoArray.Add(p.ToDouble()));
BoxLengthArray.ForEach(p => RhoCalc.LengthArray.Add(p.ToDouble()));
BoxSigmaArray.ForEach(p => RhoCalc.SigmaArray.Add(p.ToDouble()));
if (Holdsigma.Checked)
{
ChangeRoughnessArray();
}
GreyFields();
RhoCalc.UpdateProfile();
}
catch (Exception ex)
{
MessageBox.Show(ex.Message.ToString());
}
}
//Update our variables after validation
private void Variable_Changed(object sender, EventArgs e)
{
try
{
if (m_isupdating == false)
{
RhoGraphing.SupSLD = ReflGraphing.SupSLD = SupSLDTB.ToDouble();
RhoGraphing.SubSLD = ReflGraphing.SubSLD = SubphaseSLD.ToDouble();
ReflGraphing.SetGraphType(Properties.Settings.Default.ForceRQ4, DBFCB.Checked);
ReflGraphing.SetGraphType(Properties.Settings.Default.ForceRQ4, DBFCB.Checked);
MakeReflectivity();
}
}
catch { }
}
private void MakeReflectivity()
{
//Fill Rho array
try
{
ReflCalc.IsOneSigma = Holdsigma.Checked;
ReflCalc.ZOffset = 25;
ReflCalc.GetSubRoughness = SubRough.ToDouble();
ReflCalc.BoxCount = BoxCount.ToInt();
ReflCalc.SubphaseSLD = SubphaseSLD.ToDouble();
ReflCalc.SuperphaseSLD = SupSLDTB.ToDouble();
ReflCalc.HighQOffset = Rightoffset.ToInt();
ReflCalc.LowQOffset = CritOffset.ToInt();
ReflCalc.NormalizationFactor = NormCorrectTB.ToDouble();
ReflCalc.ImpNormCB = ImpNormCB.Checked;
ReflCalc.QSpreadTB = QSpreadTB.ToDouble();
//Blank our Rho data from the previous iteration
ReflCalc.RhoArray.Clear();
ReflCalc.LengthArray.Clear();
ReflCalc.SigmaArray.Clear();
RhoArray.ForEach(p => ReflCalc.RhoArray.Add(p.ToDouble()));
LengthArray.ForEach(p => ReflCalc.LengthArray.Add(p.ToDouble()));
SigmaArray.ForEach(p => ReflCalc.SigmaArray.Add(p.ToDouble()));
if (Holdsigma.Checked)
{
ChangeRoughnessArray();
}
ReflCalc.UpdateProfile();
}
catch (Exception ex)
{
MessageBox.Show(ex.Message.ToString());
}
}
private void FillSettingsStruct(ref ModelSettings settings)
{
settings.Title = TitleTB.Text;
settings.CritEdgeOffset = int.Parse(critedgeoffTB.Text);
settings.HighQOffset = int.Parse(HQoffsetTB.Text);
settings.SurflayerSLD = Rholipid.ToDouble();
settings.SubSLD = SubSLDTB.ToDouble();
settings.SupSLD = SupSLDTB.ToDouble();
settings.Boxes = Boxlayers.ToInt();
settings.Surflayerlength = layerlength.ToDouble();
settings.SurflayerAbs = SurfAbs.ToDouble();
settings.Wavelength = wavelength.ToDouble();
settings.SubAbs = SubAbs.ToDouble();
settings.SupAbs = SupAbsTB.ToDouble();
settings.UseAbs = UseAbsCB.Checked;
settings.SupOffset = Double.Parse(SupoffsetTB.Text);
settings.Percerror = Double.Parse(QErrTB.Text);
settings.Forcenorm = ForceNormCB.Checked;
settings.Forcesig = Double.Parse(SigTSTB.Text);
if (double.Parse(SigTSTB.Text) == 0)
{
settings.SigmaSearchPerc = int.Parse(SigmaSearchTB.Text);
}
else
{
settings.SigmaSearchPerc = 0;
}
settings.Debug = debugToolStripMenuItem.Checked;
settings.ForceXR = forceXRToolStripMenuItem1.Checked;
settings.Resolution = int.Parse(ResolutionTB.Text);
settings.Totallength = Double.Parse(TotlengthTB.Text);
settings.ImpNorm = ImpNormCB.Checked;
settings.FitFunc = objectiveCB.SelectedIndex;
settings.ParamTemp = Double.Parse(ParamTempTB.Text);
if (UseAbsCB.Checked)
{
settings.AbsorptionSearchPerc = int.Parse(AbsorptionSearchTB.Text);
}
else
{
settings.AbsorptionSearchPerc = 0;
}
if (ImpNormCB.Checked)
{
settings.NormalizationSearchPerc = int.Parse(NormSearchTB.Text);
}
else
{
settings.NormalizationSearchPerc = 0;
}
settings.Algorithm = AlgorithmCB.SelectedIndex;
settings.AnnealInitTemp = m_dAnnealtemp;
settings.AnnealTempPlat = m_iAnnealplat;
settings.AnnealSlope = m_dAnnealslope;
settings.AnnealGamma = m_dSTUNGamma;
settings.STUNfunc = m_iSTUNfunc;
settings.STUNAdaptive = m_bSTUNAdaptive;
settings.STUNtempiter = m_iSTUNtempiter;
settings.STUNdeciter = m_iSTUNdeciter;
settings.STUNgammadec = m_dSTUNgammadec;
}
private void LoadDataFile(string origreflfilename)
{
try
{
if (!ReflData.Instance.SetReflData(origreflfilename, !errorsAreInVarianceToolStripMenuItem.Checked))
{
throw new Exception("Could not load file");
}
FileNameTB.Text = origreflfilename;
settingsfile = origreflfilename + "settings.xml";
reflgraphobject.SetAllFonts("Garamond", 20, 18);
reflgraphobject.SubSLD = SubSLDTB.ToDouble();
reflgraphobject.SupSLD = SupSLDTB.ToDouble();
reflgraphobject.GetHighQOffset = ReflData.Instance.GetNumberDataPoints;
reflgraphobject.GetLowQOffset = 0;
reflgraphobject.SetGraphType(forceRQ4GraphingToolStripMenuItem.Checked, fresnelcb.Checked);
//Load the modeled files if they are available
if (File.Exists(settingsfile) && !m_bloadfromcommandline)
{
if (MessageBox.Show("Do you want to load the previous run?", "Resume?", MessageBoxButtons.YesNo) == DialogResult.Yes)
{
rhographobject.Clear();
reflgraphobject.Clear();
if (!LoadSettings())
{
return;
}
string tempfile = ReflData.Instance.GetWorkingDirectory + "\\rho.dat";
if (File.Exists(tempfile))
{
LoadZ(tempfile);
rhographobject.SubSLD = double.Parse(SubSLDTB.Text);
rhographobject.UseSLD = UseSLDToolStripMenuItem.Checked;
if (UseSLDToolStripMenuItem.Checked)
{
rhographobject.SetAxisTitles("Z", "SLD");
rhographobject.Title = "Model Independent SLD Fit";
}
else
{
rhographobject.SetAxisTitles("Z", "Normalized Electron Density");
rhographobject.Title = "Model Independent Electron Density Fit";
}
rhographobject.LoadFiletoGraph(tempfile.ToString(), rhomodelname, "Model Independent Electron Density Fit", Color.Tomato, SymbolType.None, 0, true);
}
tempfile = ReflData.Instance.GetWorkingDirectory + "\\rf.dat";
if (File.Exists(tempfile))
{
reflgraphobject.SubSLD = double.Parse(SubSLDTB.Text);
reflgraphobject.SupSLD = double.Parse(SupSLDTB.Text);
reflgraphobject.SetGraphType(forceRQ4GraphingToolStripMenuItem.Checked, fresnelcb.Checked);
reflgraphobject.GetLowQOffset = int.Parse(critedgeoffTB.Text);
reflgraphobject.GetHighQOffset = ReflData.Instance.GetNumberDataPoints - int.Parse(HQoffsetTB.Text);
//Load the data file to the graph
reflgraphobject.LoadDatawithErrorstoGraph("Reflectivity Data", Color.Black, SymbolType.Circle, 5, ReflData.Instance.GetQData, ReflData.Instance.GetReflData);
reflgraphobject.LoadFiletoGraph(tempfile.ToString(), modelreflname, "Model Independent Reflectivity", Color.Tomato, SymbolType.Square, 2, true);
}
GraphCollection.Instance.MainReflGraph = reflgraphobject;
GraphCollection.Instance.MainRhoGraph = rhographobject;
UpdateReportParameters();
}
else
{
//Clear the report
ReportGenerator.Instance.ClearAll();
//Load the data file to the graph
reflgraphobject.LoadDatawithErrorstoGraph("Reflectivity Data", Color.Black, SymbolType.Circle, 5, ReflData.Instance.GetQData, ReflData.Instance.GetReflData);
FileInfo info = new FileInfo(settingsfile);
if (File.Exists(info.DirectoryName + "\\pop.dat"))
{
//Create a backup folder if we don't have one
if (!System.IO.Directory.Exists(info.DirectoryName + "\\FitBackUp"))
{
System.IO.Directory.CreateDirectory(info.DirectoryName + "\\FitBackUp");
}
int index = 0;
for (; ;)
{
if (!System.IO.Directory.Exists(info.DirectoryName + "\\FitBackUp\\Fit" + index.ToString()))
{
System.IO.Directory.CreateDirectory(info.DirectoryName + "\\FitBackUp\\Fit" + index.ToString());
break;
}
index++;
}
string fileloc = info.DirectoryName + "\\FitBackUp\\Fit" + index.ToString();
File.Move(info.DirectoryName + "\\pop.dat", fileloc + "\\pop.dat");
File.Move(settingsfile, fileloc + "\\" + info.Name);
if (File.Exists(info.DirectoryName + "\\rf.dat"))
{
File.Move(info.DirectoryName + "\\rf.dat", fileloc + "\\rf.dat");
}
if (File.Exists(info.DirectoryName + "\\rho.dat"))
{
File.Move(info.DirectoryName + "\\rho.dat", fileloc + "\\rho.dat");
}
if (File.Exists(info.DirectoryName + "\\reflfile.dat"))
{
File.Move(info.DirectoryName + "\\reflfile.dat", fileloc + "\\reflfile.dat");
}
if (File.Exists(info.DirectoryName + "\\reflrhofit.dat"))
{
File.Move(info.DirectoryName + "\\reflrhofit.dat", fileloc + "\\reflrhofit.dat");
}
}
}
}
else
{
//Load the data file to the graph
reflgraphobject.LoadDatawithErrorstoGraph("Reflectivity Data", Color.Black, SymbolType.Circle, 5, ReflData.Instance.GetQData, ReflData.Instance.GetReflData);
}
DisableInterface(false);
Cancelbutton.Enabled = false;
reflgraphobject.SetBounds();
SubAbs.Enabled = SupAbsTB.Enabled = SurfAbs.Enabled = UseAbsCB.Checked;
}
catch (Exception ex)
{
origreflfilename = string.Empty;
FileNameTB.Clear();
DisableInterface(true);
MessageBox.Show(ex.Message);
}
}
public ConstrainedReflmodeling(RhoFit RhoCalc)
{
InitializeComponent();
ReflCalc = new ReflFit(RhoCalc as BoxReflFitBase);
m_bUseSLD = Properties.Settings.Default.UseSLDSingleSession;
SigmaArray = new List <TextBox>(6);
LengthArray = new List <TextBox>(6);
RhoArray = new List <TextBox>(6);
HoldSigmaArray = new List <CheckBox>(6);
HoldLengthArray = new List <CheckBox>(6);
HoldRhoArray = new List <CheckBox>(6);
if (m_bUseSLD)
{
tabControl1.TabPages[1].Text = "SLD";
RhoLabel.Text = "SLD";
}
else
{
tabControl1.TabPages[1].Text = "Electron Density";
}
WavelengthTB.Text = ((double)(1.24)).ToString();
//Initialize the arrays
MakeArrays();
//Copy over the values
BoxCount.Text = RhoCalc.BoxCount.ToString();
SubRough.Text = RhoCalc.GetSubRoughness.ToString();
SubphaseSLD.Text = RhoCalc.SubphaseSLD.ToString();
SupSLDTB.Text = RhoCalc.SuperphaseSLD.ToString();
Holdsigma.Checked = RhoCalc.IsOneSigma;
//Setup the Graph
ReflGraphing = new Graphing(string.Empty);
ReflGraphing.SetGraphType(Properties.Settings.Default.ForceRQ4, DBFCB.Checked);
ReflGraphing.SubSLD = SubphaseSLD.ToDouble();
ReflGraphing.SupSLD = SupSLDTB.ToDouble();
ReflGraphing.CreateGraph(RhoGraph, "Reflectivity", "Q/Qc", "Intensity / Fresnel", AxisType.Log);
ReflGraphing.LoadDatawithErrorstoGraph("Reflectivity Data", System.Drawing.Color.Black, SymbolType.Circle, 5, ReflData.Instance.GetQData, ReflData.Instance.GetReflData);
//Set up ED Graph
RhoGraphing = new Graphing(string.Empty);
RhoGraphing.SubSLD = SubphaseSLD.ToDouble();
RhoGraphing.UseSLD = m_bUseSLD;
RhoGraphing.SetGraphType(false, false);
if (m_bUseSLD == false)
{
RhoGraphing.CreateGraph(EDzedGraphControl1, "Electron Density Profile", "Z", "Normalized Electron Density", AxisType.Linear);
}
else
{
RhoGraphing.CreateGraph(EDzedGraphControl1, "SLD Profile", "Z", "SLD", AxisType.Linear);
}
RhoGraphing.Pane.XAxis.Scale.Min = 0;
RhoGraphing.Pane.XAxis.Scale.Max = ReflCalc.Z[ReflCalc.Z.Length - 1];
loadingCircle2.Visible = loadingCircle1.Visible = false;
loadingCircle2.NumberSpoke = loadingCircle1.NumberSpoke = 25;
loadingCircle2.InnerCircleRadius = loadingCircle1.InnerCircleRadius = 60;
loadingCircle2.OuterCircleRadius = loadingCircle1.OuterCircleRadius = 120;
loadingCircle2.RotationSpeed = loadingCircle1.RotationSpeed = 150;
loadingCircle2.SpokeThickness = loadingCircle1.SpokeThickness = 3;
//Initialize constrain form
ConstrForm = new Constraints(6);
GreyFields();
//Setup the callback if the graph updates the bounds
ReflGraphing.ChangedBounds += new Graphing.ChangedEventHandler(PointChanged);
//Setup the callback to update the frontend with new information
ReflCalc.Update += new BoxReflFitBase.UpdateProfileHandler(ReflCalc_Update);
ReflCalc_Update(null, null);
MakeReflectivity();
}