From ac047732855fbaf26ea3029400184480fd7d4ae4 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Mon, 15 Dec 2014 14:46:12 +0100 Subject: Add doc/hitrate.html --- AUTHORS | 1 + Makefile.am | 3 ++- doc/hitrate.html | 80 +++++++++++++++++++++++++++++++++++++++++++++++++++++++ doc/hitrate.png | Bin 0 -> 48217 bytes 4 files changed, 83 insertions(+), 1 deletion(-) create mode 100644 doc/hitrate.html create mode 100644 doc/hitrate.png diff --git a/AUTHORS b/AUTHORS index b4f2af19..86dfdc40 100644 --- a/AUTHORS +++ b/AUTHORS @@ -58,6 +58,7 @@ * Takanori Nakane hdfsee ADSC export and scrollbars Some other bits and pieces + Hit rate calculator (doc/hitrate.{html,png}) * Keitaro Yamashita process_hkl bug fix diff --git a/Makefile.am b/Makefile.am index fc4e6df5..d18b4bd3 100644 --- a/Makefile.am +++ b/Makefile.am @@ -151,7 +151,8 @@ crystfeldoc_DATA = doc/twin-calculator.pdf doc/examples/lcls-dec.geom \ doc/examples/simple.geom \ doc/examples/cspad-feb2011.geom \ doc/examples/pilatus.geom \ - doc/examples/cell-example.cell + doc/examples/cell-example.cell \ + doc/hitrate.html doc/hitrate.png EXTRA_DIST += $(crystfeldoc_DATA) doc/twin-calculator.odt \ doc/reference/libcrystfel/xml/coding-standards.xml diff --git a/doc/hitrate.html b/doc/hitrate.html new file mode 100644 index 00000000..e8c95952 --- /dev/null +++ b/doc/hitrate.html @@ -0,0 +1,80 @@ + + + +SFX hitrate calculator + + + + +

Serial crystallography Hitrate calculator

+ +
+Mean number of crystals per shot: + +
+
+ +
+ +

Hitrate was modeled by Poisson distribution as in [1]. +The vertical line is at rate parameter = 1.

+ +
    +
  1. Hunter, Mark S., et +al. "Fixed-target +protein serial microcrystallography with an x-ray free electron +laser" Scientific reports 4 (2014).
    2. +
+ +

The plot was generated by the following R code.

+ +
+lambda <- seq(0, 3, by=0.05)
+plot(lambda, 100 - 100* dpois(0, lambda=lambda), type='l', 
+     ylab="Rate (%)", xlab="Mean(#Crystal/Shot)")
+lines(lambda, 100 * dpois(0, lambda=lambda), lty=2)
+lines(lambda, 100 * (1 - dpois(0, lambda=lambda) - dpois(1, lambda=lambda)), lty=3)
+lines(lambda, 100 * dpois(1, lambda=lambda), lty=1, col=2)
+lines(lambda, 100 * dpois(2, lambda=lambda), lty=1, col=3)
+lines(lambda, 100 * dpois(3, lambda=lambda), lty=1, col=4)
+abline(v=1)
+legend("topleft", lty=c(1, 2, 3, 1, 1, 1), col=c(1, 1, 1, 2, 3, 4),
+        c("Hit rate", "Blank frame", "Multiple hit", "Single lattice", 
+          "Two lattices", "Three lattices"))
+
+ + + + diff --git a/doc/hitrate.png b/doc/hitrate.png new file mode 100644 index 00000000..478a1b76 Binary files /dev/null and b/doc/hitrate.png differ -- cgit v1.2.3