radiation 0.0.4 → 0.1.0

data/samples/efficiencies.rb

@@ -0,0 +1,10 @@
+# encoding: utf-8
+require 'radiation'
+
+file     = "./B0-Ra226.xml"
+source   = Radiation::Source.new(nuclide: "Ra-226", resource: "nucleide.org")
+
+puts ["E_ɣ", "I_ɣ", "ΔI_ɣ", "e", "Δe"].join("\t")
+Radiation::Spectrum.new(source: source).parse_hdtv(file).calibrate.efficiencies.select{|p| p[:intensity] > 0.3}.sort_by{|k| k[:energy]}.each do |p|
+	puts "#{p[:energy].to_f.round(1)}\t#{p[:intensity].value}\t#{p[:intensity].delta}\t#{p[:efficiency].value.round(1)}\t#{p[:efficiency].delta.round(1)}"
+end

data/spec/radiation_spec.rb

@@ -9,7 +9,7 @@ describe Radiation do
 		it "provides information about itself" do
 			source.nuclide.should eq(name)
 			#source.reference.should be_a_kind_of(String)
-			#source.halflive.should be_a_kind_of(Plusminus::PlusminusFloat)
+			#source.halflife.should be_a_kind_of(Plusminus::PlusminusFloat)
 		end
 
 		it "returns an array with gamma-ray energies for a given source" do
@@ -23,12 +23,13 @@ describe Radiation do
 			source.intensities.first.should be_a_kind_of(Hash)
 			source.intensities.length.should be > 0
 		end
+
 		it "can access different data resources" do
 			expect { Radiation::Source.new(nuclide: "Ra-226") }.to_not raise_error
 			expect { Radiation::Source.new(nuclide: "Ra-226", resource: "nucleide.org") }.to_not raise_error
 		end
 
-		it "cheks for valid nuclei" do
+		it "checks for valid nuclei" do
 			Radiation::Source.new.is_nuclide?("Nukular9000").should be false
 			Radiation::Source.new.is_nuclide?("226Ra-226").should be false
 			Radiation::Source.new.is_nuclide?("Ra-226").should be true
@@ -37,10 +38,15 @@ describe Radiation do
 			expect { Radiation::Source.new(nuclide: "Nukular9000") }.to raise_error
 			expect { Radiation::Source.new(nuclide: "Ra-226") }.to_not raise_error
 		end
+
 		it "gives an error if no record is found" do
 			expect { Radiation::Source.new(nuclide: "Ra-15") }.to raise_error
 			expect { Radiation::Source.new(nuclide: "Ra-15", resource: "nucleide.org") }.to raise_error
 		end
+
+		it "can access nuclide halflife" do
+			Radiation::Source.new(nuclide: "Ra-226", resource: "nucleide.org").halflife.to_f.should be_within(0.001).of(50.49E9)
+		end
 	end
 
 	describe Radiation::Spectrum do
@@ -63,8 +69,12 @@ describe Radiation do
 
 		it "can read peaks from hdtv xml data" do
 			expect { Radiation::Spectrum.new(source: source ).parse_hdtv("./samples/B0-Ra226.xml") }.to_not raise_error
-			Radiation::Spectrum.new(source: source ).parse_hdtv("./samples/B0-Ra226.xml").calibrate.calibration[1].round(3).should be == 0.100
+			Radiation::Spectrum.new(source: source ).parse_hdtv("./samples/B0-Ra226.xml").calibrate.calibration[1].to_f.should be_within(0.001).of(0.1)
+		end
+
+		it "can calculate detection efficiencies" do
+			Radiation::Spectrum.new(source: source).parse_hdtv("./samples/B0-Ra226.xml").calibrate.efficiencies.select{|p| p[:energy].to_f.to_i == 295}.first[:efficiency].should be_a_kind_of(Plusminus::PlusminusFloat)
 		end
 
 	end
-end
+end

metadata

@@ -1,97 +1,97 @@
 --- !ruby/object:Gem::Specification
 name: radiation
 version: !ruby/object:Gem::Version
-  version: 0.0.4
+  version: 0.1.0
 platform: ruby
 authors:
 - Jan Mayer
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2013-06-25 00:00:00.000000000 Z
+date: 2013-07-09 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: plusminus
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: open-uri-cached
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: combinatorics
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: linefit
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: xml-simple
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: thor
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
@@ -112,28 +112,28 @@ dependencies:
   name: rake
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 - !ruby/object:Gem::Dependency
   name: rspec
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
   type: :development
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - ! '>='
+    - - '>='
       - !ruby/object:Gem::Version
         version: '0'
 description: Decay Radiation
@@ -159,10 +159,10 @@ files:
 - data/Ra-226.csv
 - lib/radiation.rb
 - lib/radiation/cli.rb
+- lib/radiation/resource.rb
+- lib/radiation/resource/internal.rb
+- lib/radiation/resource/nucleideorg.rb
 - lib/radiation/source.rb
-- lib/radiation/source/resource.rb
-- lib/radiation/source/resource/internal.rb
-- lib/radiation/source/resource/nucleideorg.rb
 - lib/radiation/spectrum.rb
 - lib/radiation/version.rb
 - radiation.gemspec
@@ -170,6 +170,7 @@ files:
 - samples/B0-Ra226.xml
 - samples/calibrate_a_spectrum.rb
 - samples/calibration_from_hdtv_xml.rb
+- samples/efficiencies.rb
 - samples/get_some_data.rb
 - spec/radiation_spec.rb
 - spec/spec_helper.rb
@@ -183,12 +184,12 @@ require_paths:
 - lib
 required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
-  - - ! '>='
+  - - '>='
     - !ruby/object:Gem::Version
       version: '0'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
-  - - ! '>='
+  - - '>='
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []

data/lib/radiation/source/resource/internal.rb

@@ -1,17 +0,0 @@
-# encoding: utf-8
-require 'csv'
-
-class Radiation::Source
-	class Resource
-		private
-		def datasource_internal(nuclide)
-			begin
-				path = File.join(File.dirname(__FILE__), "../../../../data")
-				@data[:nuclide] = nuclide.to_s
-				@data[:lines] = CSV.read("#{path}/#{nuclide}.csv",  {col_sep: ';', converters: :numeric}).collect{|row| {:energy => row[0].pm(0), :intensity => row[1].pm(row[2])} }
-			rescue
-				raise "No Data for #{nuclide}"
-			end
-		end
-	end
-end

data/lib/radiation/source/resource/nucleideorg.rb

@@ -1,36 +0,0 @@
-# encoding: utf-8
-require 'plusminus'
-
-require 'open-uri'
-require 'open-uri/cached'
-require 'tmpdir'
-
-
-class Radiation::Source
-	class Resource
-
-		private
-		def datasource_nukleideorg(nuclide)
-			nuclide = "Ra-226D" if nuclide == "Ra-226" #Ra-226 in equilibrium with daughters
-			begin
-				@data[:nuclide] = nuclide.to_s
-				OpenURI::Cache.cache_path = "#{Dir.tmpdir}/radiation/"
-				uri = open("http://www.nucleide.org/DDEP_WG/Nuclides/#{nuclide}.lara.txt").readlines
-				start = 0
-				uri.each_with_index do |line, lineno|
-					if line.start_with?("------")
-						start = lineno + 2
-						break
-					end
-				end
-				return @data[:lines] = [] if start == 0 or uri.count < start
-				@data[:lines] = uri[start...-1].collect{|line| line.split(' ; ')}.select!{|row| row[4] == "g"}.collect do |row|
-				 {:energy => row[0].to_f.pm(row[1].to_f), :intensity => row[2].to_f.pm(row[3].to_f)} 
-				end
-			rescue
-				raise "No Data for #{nuclide}"
-			end
-		end
-
-	end
-end

data/lib/radiation/source/resource.rb

@@ -1,25 +0,0 @@
-# encoding: utf-8
-require 'radiation/source/resource/internal'
-require 'radiation/source/resource/nucleideorg'
-
-module Radiation
-class Source
-	class Resource
-		
-		def initialize(resource)
-			@resource = resource
-		end
-
-		def data(nuclide)
-			@data = {}
-			case @resource
-				when "internal" then datasource_internal(nuclide)
-				when "nucleide.org" then datasource_nukleideorg(nuclide)
-				else raise "Unknown Datasource"
-			end
-			return @data
-		end
-
-	end
-end
-end