[CSC 335] Check on Peaks

Andrew J. Pounds pounds_aj at mercer.edu
Tue Oct 29 02:36:28 EDT 2013


Gentlemen -- in class today Bryan asked me to check on my peak tops 
because they were an order or magnitude higher than the ones you were 
getting.   He is correct -- when I added the code to display the peak 
tops I copied a line in my code that included a Fortran edit descriptor 
"multiplier" that has to be there for the items displayed in scientific 
notation to show up properly - but it was not needed for the things 
shown in in floating format.

So for comparison all of my numbers in the "Top" column shown on the 
handout should be divided by 10.

Also, regarding the "extra peak" that some of you are finding, I think 
it is related to when and how you are filtering.  I replotted my 
filtered peaks and the baseline.  They are shown below.

nmr

In my results the peak at -3.6 is getting picked up as one big peak.  
Based on the plot, this makes sense because the minimum between the two 
peaks on the left is above the baseline.  If for any reason your 
filtering pulls that minimum down below the baseline then you should get 
another peak.

If you are curious, my order of operations is...

 1. Fix TMS location
 2. Adjust the baseline
 3. Filter
 4. Fit the spline
 5. Find the peaks
 6. Locate peak max values
 7. Integrate peaks
 8. Calculate hydrogens


If you are curious about your own code then I really encourage you to 
plot things in gnuplot.  You can create a datafile for checking your 
baseline.  If the baseling were 1400 then I would create a datafile with 
two points.

-5 1400
0 1400

If you plot this with your spectrum then you should see where your peaks 
are getting selected.  In the file above, which was after fixing the 
baseline and filtering, I set my file up as....

-5 0
0 0


-- 
Andrew J. Pounds, Ph.D.  (pounds_aj at mercer.edu)
Professor of Chemistry and Computer Science
Mercer University,  Macon, GA 31207   (478) 301-5627
http://faculty.mercer.edu/pounds_aj

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://theochem.mercer.edu/pipermail/csc335/attachments/20131029/fb9e4ed9/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot.png
Type: image/png
Size: 5504 bytes
Desc: not available
URL: <http://theochem.mercer.edu/pipermail/csc335/attachments/20131029/fb9e4ed9/attachment.png>


More information about the csc335 mailing list