print 'PL File Path:', filepath
print 'Select the system response file to use (if any)...'
sysres_filepath = open_dialog('*.*')
print 'System Response File Path:', sysres_filepath
measuredSpectrum = openMeasuredSpectrum(filepath, sysres_filepath)
energy = measuredSpectrum.getEnergy()
intensity = measuredSpectrum.getIntensity()
# Create some simulated spectrum
bg = ConstantSpectrum(energy=energy)
peak1 = GaussianSpectrum(energy=energy)
peak2 = GaussianSpectrum(energy=energy)
peak3 = GaussianSpectrum(energy=energy)
sum = SummedSpectrum(bg, peak1, peak2, peak3)
sum.energy = energy
bg.constant.max = peak1.amplitude.max = \
peak2.amplitude.max = peak3.amplitude.max = intensity.max()
peak1.center.min = peak2.center.min = peak3.center.min = energy.min()
peak1.center.max = peak2.center.max = peak3.center.max = energy.max()
peak1.center.value = peak2.center.value = \
peak3.center.value = energy[energy.size/2]
peak1.fwhm.min = peak2.fwhm.min = \
peak3.fwhm.min = (energy.max() - energy.min())/100.
peak1.fwhm.max = peak2.fwhm.max = \
peak3.fwhm.max = (energy.max() - energy.min())
peak1.fwhm.value = peak2.fwhm.value = \
peak3.fwhm.value = (energy.max() - energy.min())/10
parameters = []
parameters.extend(bg.parameters)
