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 --- doc/hitrate.html | 80 +++++++++++++++++++++++++++++++++++++++++++++++++++++++ doc/hitrate.png | Bin 0 -> 48217 bytes 2 files changed, 80 insertions(+) create mode 100644 doc/hitrate.html create mode 100644 doc/hitrate.png (limited to 'doc') 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