Merge branch 'julia'

1 files changed, 26 insertions, 0 deletions

diff --git a/julia/make-indexed.jl b/julia/make-indexed.jl
new file mode 100644
index 00000000..7a24f606
--- /dev/null
+++ b/julia/make-indexed.jl
@@ -0,0 +1,26 @@
+using CrystFEL
+using Random
+
+# "Simulate" a diffraction pattern from the reflections
+function sketch_pattern(image, cr)
+    reflist = predictreflections(cr, image)
+    peaklist = PeakList()
+    for refl in reflist
+        if randn() > 0
+            let dpos = refl.detectorposition
+                push!(peaklist, dpos.fs, dpos.ss, dpos.panelnumber, 100.0)
+            end
+        end
+    end
+    return peaklist
+end
+
+
+dtempl = loaddatatemplate("julia/alignment-test.geom")
+cell = UnitCell(MonoclinicLattice, PrimitiveCell, 123, 45, 80, 90, 97, 90)
+indexer = Indexer("asdf", dtempl, cell, retry=false, multilattice=false, refine=true)
+
+image = Image(dtempl)
+cr = Crystal(rotatecell(cell))
+image.peaklist = sketch_pattern(image, cr)
+index(image, indexer)