Index: /Backup NB/Talks/MEMPHYS-Trigger/Binomial.nb
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/Binomial.nb	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/Binomial.nb	(revision 388)
@@ -0,0 +1,220 @@
+(************** Content-type: application/mathematica **************
+
+                    Mathematica-Compatible Notebook
+
+This notebook can be used with any Mathematica-compatible
+application, such as Mathematica, MathReader or Publicon. The data
+for the notebook starts with the line containing stars above.
+
+To get the notebook into a Mathematica-compatible application, do
+one of the following:
+
+* Save the data starting with the line of stars above into a file
+  with a name ending in .nb, then open the file inside the
+  application;
+
+* Copy the data starting with the line of stars above to the
+  clipboard, then use the Paste menu command inside the application.
+
+Data for notebooks contains only printable 7-bit ASCII and can be
+sent directly in email or through ftp in text mode.  Newlines can be
+CR, LF or CRLF (Unix, Macintosh or MS-DOS style).
+
+NOTE: If you modify the data for this notebook not in a Mathematica-
+compatible application, you must delete the line below containing
+the word CacheID, otherwise Mathematica-compatible applications may
+try to use invalid cache data.
+
+For more information on notebooks and Mathematica-compatible 
+applications, contact Wolfram Research:
+  web: http://www.wolfram.com
+  email: info@wolfram.com
+  phone: +1-217-398-0700 (U.S.)
+
+Notebook reader applications are available free of charge from 
+Wolfram Research.
+*******************************************************************)
+
+(*CacheID: 232*)
+
+
+(*NotebookFileLineBreakTest
+NotebookFileLineBreakTest*)
+(*NotebookOptionsPosition[      4421,        144]*)
+(*NotebookOutlinePosition[      5546,        180]*)
+(*  CellTagsIndexPosition[      5418,        173]*)
+(*WindowFrame->Normal*)
+
+
+
+Notebook[{
+Cell[BoxData[
+    \(BinomPDF[k_, N_, dt_, f_] := 
+      PDF[BinomialDistribution[N, 1 - Exp[\(-f\)*dt]], k]\)], "Input"],
+
+Cell[BoxData[
+    \(Trigger1[nTh_, N_, dt_, f_] := 
+      Evaluate[Sum[
+            BinomPDF[k, N, 5/10^9, f]*
+              BinomPDF[nTh - k, N, dt, 
+                f]*\((1 - Sum[BinomPDF[j, N, 5/10^9, f], {j, 0, k}])\), {k, 
+              0, nTh}]]*\((10^9/5)\)\)], "Input"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(\(?BinomPDF\)\)], "Input"],
+
+Cell["Global`BinomPDF", "Print",
+  CellTags->"Info3362315879-8409761"],
+
+Cell[BoxData[
+    InterpretationBox[GridBox[{
+          {GridBox[{
+                {\(BinomPDF[k_, N_, dt_, f_] := 
+                    PDF[BinomialDistribution[N, 
+                        1 - \[ExponentialE]\^\(\(-f\)\ dt\)], k]\)}
+                },
+              GridBaseline->{Baseline, {1, 1}},
+              ColumnWidths->0.999,
+              ColumnAlignments->{Left}]}
+          },
+        GridBaseline->{Baseline, {1, 1}},
+        ColumnAlignments->{Left}],
+      Definition[ "BinomPDF"],
+      Editable->False]], "Print",
+  CellTags->"Info3362315879-8409761"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(\(?Trigger1\)\)], "Input"],
+
+Cell[BoxData[
+    \(Information::"notfound" \(\(:\)\(\ \)\) 
+      "Symbol \!\(\"Trigger1\"\) not found."\)], "Message"]
+}, Open  ]],
+
+Cell[BoxData[
+    \(Trigger1[nTh_, N_, dt_, f_] := 
+      Evaluate[Sum[
+            BinomPDF[k, N, 0.0001/10^9, f]*
+              BinomPDF[nTh - k, N, dt, 
+                f]*\((1 - 
+                  Sum[BinomPDF[j, N, 0.0001/10^9, f], {j, 0, k}])\), {k, 0, 
+              nTh}]]*\((10^9/0.0001)\)\)], "Input"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(\(?Trigger1\)\)], "Input"],
+
+Cell["Global`Trigger1", "Print",
+  CellTags->"Info3362315904-1598254"],
+
+Cell[BoxData[
+    InterpretationBox[GridBox[{
+          {GridBox[{
+                {\(Trigger1[nTh_, N_, dt_, f_] := 
+                    1\/5\ Evaluate[\[Sum]\+\(k = 0\)\%nTh BinomPDF[k, N, 
+                              5\/10\^9, f]\ BinomPDF[nTh - k, N, dt, 
+                              f]\ \((1 - \[Sum]\+\(j = 0\)\%k BinomPDF[j, N, 
+                                    5\/10\^9, f])\)]\ 10\^9\)}
+                },
+              GridBaseline->{Baseline, {1, 1}},
+              ColumnWidths->0.999,
+              ColumnAlignments->{Left}]}
+          },
+        GridBaseline->{Baseline, {1, 1}},
+        ColumnAlignments->{Left}],
+      Definition[ "Trigger1"],
+      Editable->False]], "Print",
+  CellTags->"Info3362315904-1598254"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(N[Trigger1[249, 81000, 425/10^9, 6\ 10^3]]\)], "Input"],
+
+Cell[BoxData[
+    \(192707.40988527454`\)], "Output"]
+}, Open  ]]
+},
+FrontEndVersion->"4.1 for Microsoft Windows",
+ScreenRectangle->{{0, 1024}, {0, 685}},
+WindowSize->{495, 527},
+WindowMargins->{{0, Automatic}, {Automatic, 0}}
+]
+
+(*******************************************************************
+Cached data follows.  If you edit this Notebook file directly, not
+using Mathematica, you must remove the line containing CacheID at
+the top of  the file.  The cache data will then be recreated when
+you save this file from within Mathematica.
+*******************************************************************)
+
+(*CellTagsOutline
+CellTagsIndex->{
+  "Info3362315879-8409761"->{
+    Cell[2179, 67, 70, 1, 25, "Print",
+      CellTags->"Info3362315879-8409761"],
+    Cell[2252, 70, 568, 15, 44, "Print",
+      CellTags->"Info3362315879-8409761"]},
+  "Info3362315904-1598254"->{
+    Cell[3426, 112, 70, 1, 25, "Print",
+      CellTags->"Info3362315904-1598254"],
+    Cell[3499, 115, 738, 17, 77, "Print",
+      CellTags->"Info3362315904-1598254"]}
+  }
+*)
+
+(*CellTagsIndex
+CellTagsIndex->{
+  {"Info3362315879-8409761", 5027, 159},
+  {"Info3362315904-1598254", 5224, 164}
+  }
+*)
+
+(*NotebookFileOutline
+Notebook[{
+Cell[1705, 50, 120, 2, 50, "Input"],
+Cell[1828, 54, 277, 6, 130, "Input"],
+
+Cell[CellGroupData[{
+Cell[2130, 64, 46, 1, 30, "Input"],
+Cell[2179, 67, 70, 1, 25, "Print",
+  CellTags->"Info3362315879-8409761"],
+Cell[2252, 70, 568, 15, 44, "Print",
+  CellTags->"Info3362315879-8409761"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+Cell[2857, 90, 46, 1, 30, "Input"],
+Cell[2906, 93, 120, 2, 24, "Message"]
+}, Open  ]],
+Cell[3041, 98, 311, 7, 130, "Input"],
+
+Cell[CellGroupData[{
+Cell[3377, 109, 46, 1, 30, "Input"],
+Cell[3426, 112, 70, 1, 25, "Print",
+  CellTags->"Info3362315904-1598254"],
+Cell[3499, 115, 738, 17, 77, "Print",
+  CellTags->"Info3362315904-1598254"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+Cell[4274, 137, 75, 1, 30, "Input"],
+Cell[4352, 140, 53, 1, 29, "Output"]
+}, Open  ]]
+}
+]
+*)
+
+
+
+(*******************************************************************
+End of Mathematica Notebook file.
+*******************************************************************)
+
Index: /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.aux
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.aux	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.aux	(revision 388)
@@ -0,0 +1,14 @@
+\relax 
+\@writefile{toc}{\contentsline {section}{\numberline {1}Motivations}{1}}
+\@writefile{toc}{\contentsline {section}{\numberline {2}A toy MC}{1}}
+\@writefile{toc}{\contentsline {subsection}{\numberline {2.1}Dark current generation}{2}}
+\newlabel{eq:singleProba}{{1}{2}}
+\newlabel{eq:poissonLimit}{{2}{2}}
+\@writefile{lof}{\contentsline {figure}{\numberline {1}{\ignorespaces  {\it  left panel}: exponential generation of $\delta t$ with $1/\tau = 1$\nobreakspace  {}MHz dark current noise; {\it  right panel}: number of PMTs fired among 50 PMTs, per $\Delta t = 100$\nobreakspace  {}ns windows. One expects a Poisson distribution with a mean of $N\Delta t/\tau = 5$ as it is shown. }}{2}}
+\newlabel{fig:PMTPoisson}{{1}{2}}
+\newlabel{eq:BinomialProb}{{3}{2}}
+\@writefile{toc}{\contentsline {subsection}{\numberline {2.2}Trigger simulation}{2}}
+\@writefile{lof}{\contentsline {figure}{\numberline {2}{\ignorespaces  {\it  left panel}: Digital sum of the time evolution of $N=81,000$ PMTs with a dark current of $6$\nobreakspace  {}kHz, computed for a time sliding window of $425$\nobreakspace  {}ns every $5$\nobreakspace  {}ns; {\it  right panel}: zoom around the location of a trigger threshold at $250$ PMTs. }}{3}}
+\newlabel{fig:DigitalSum}{{2}{3}}
+\@writefile{lot}{\contentsline {table}{\numberline {1}{\ignorespaces  Mean number of PMTs ($m$)fired among $N=81,000$ PMTs on a time sliding window of $\Delta t=425$\nobreakspace  {}ns corresponding to the dark current (DC) rate ($f$) per PMT (this is $N\Delta t * f$). Then, setting a threshold to $m+3\sqrt  m$, one reads the trigger rate (see text for the definition of a trigger).}}{4}}
+\newlabel{tab:TriggerRate}{{1}{4}}
Index: /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.blg
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.blg	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.blg	(revision 388)
@@ -0,0 +1,5 @@
+This is BibTeX, Version 0.99cThe top-level auxiliary file: C:\Presentation\NuFactSB\MEMPHYS-Trigger\MEMPHYS-Trigger.aux
+I found no \citation commands---while reading file C:\Presentation\NuFactSB\MEMPHYS-Trigger\MEMPHYS-Trigger.aux
+I found no \bibdata command---while reading file C:\Presentation\NuFactSB\MEMPHYS-Trigger\MEMPHYS-Trigger.aux
+I found no \bibstyle command---while reading file C:\Presentation\NuFactSB\MEMPHYS-Trigger\MEMPHYS-Trigger.aux
+(There were 3 error messages)
Index: /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.log
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.log	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.log	(revision 388)
@@ -0,0 +1,194 @@
+This is e-TeX, Version 3.141592-2.1 (MiKTeX 2.4) (preloaded format=latex 2004.7.19)  18 JUL 2006 20:59
+entering extended mode
+**MEMPHYS-Trigger.tex
+(MEMPHYS-Trigger.tex
+LaTeX2e <2003/12/01>
+Babel <v3.8a> and hyphenation patterns for english, french, german, ngerman, du
+mylang, nohyphenation, loaded.
+(C:\Program Files\MikTex\tex\latex\base\article.cls
+Document Class: article 2004/02/16 v1.4f Standard LaTeX document class
+(C:\Program Files\MikTex\tex\latex\base\size12.clo
+File: size12.clo 2004/02/16 v1.4f Standard LaTeX file (size option)
+)
+\c@part=\count79
+\c@section=\count80
+\c@subsection=\count81
+\c@subsubsection=\count82
+\c@paragraph=\count83
+\c@subparagraph=\count84
+\c@figure=\count85
+\c@table=\count86
+\abovecaptionskip=\skip41
+\belowcaptionskip=\skip42
+\bibindent=\dimen102
+)
+(C:\Program Files\MikTex\tex\latex\graphics\graphicx.sty
+Package: graphicx 1999/02/16 v1.0f Enhanced LaTeX Graphics (DPC,SPQR)
+
+(C:\Program Files\MikTex\tex\latex\graphics\keyval.sty
+Package: keyval 1999/03/16 v1.13 key=value parser (DPC)
+\KV@toks@=\toks14
+)
+(C:\Program Files\MikTex\tex\latex\graphics\graphics.sty
+Package: graphics 2001/07/07 v1.0n Standard LaTeX Graphics (DPC,SPQR)
+
+(C:\Program Files\MikTex\tex\latex\graphics\trig.sty
+Package: trig 1999/03/16 v1.09 sin cos tan (DPC)
+)
+(C:\Program Files\MikTex\tex\latex\00miktex\graphics.cfg
+File: graphics.cfg 2003/03/12 v1.1 MiKTeX 'graphics' configuration
+)
+Package graphics Info: Driver file: dvips.def on input line 80.
+
+(C:\Program Files\MikTex\tex\latex\graphics\dvips.def
+File: dvips.def 1999/02/16 v3.0i Driver-dependant file (DPC,SPQR)
+))
+\Gin@req@height=\dimen103
+\Gin@req@width=\dimen104
+)
+(C:\Program Files\MikTex\tex\latex\graphics\epsfig.sty
+Package: epsfig 1999/02/16 v1.7a (e)psfig emulation (SPQR)
+\epsfxsize=\dimen105
+\epsfysize=\dimen106
+)
+(C:\Program Files\MikTex\tex\latex\amsfonts\amssymb.sty
+Package: amssymb 2002/01/22 v2.2d
+
+(C:\Program Files\MikTex\tex\latex\amsfonts\amsfonts.sty
+Package: amsfonts 2001/10/25 v2.2f
+\@emptytoks=\toks15
+\symAMSa=\mathgroup4
+\symAMSb=\mathgroup5
+LaTeX Font Info:    Overwriting math alphabet `\mathfrak' in version `bold'
+(Font)                  U/euf/m/n --> U/euf/b/n on input line 132.
+))
+(C:\Program Files\MikTex\tex\latex\amsmath\amsmath.sty
+Package: amsmath 2000/07/18 v2.13 AMS math features
+\@mathmargin=\skip43
+
+For additional information on amsmath, use the `?' option.
+(C:\Program Files\MikTex\tex\latex\amsmath\amstext.sty
+Package: amstext 2000/06/29 v2.01
+
+(C:\Program Files\MikTex\tex\latex\amsmath\amsgen.sty
+File: amsgen.sty 1999/11/30 v2.0
+\@emptytoks=\toks16
+\ex@=\dimen107
+))
+(C:\Program Files\MikTex\tex\latex\amsmath\amsbsy.sty
+Package: amsbsy 1999/11/29 v1.2d
+\pmbraise@=\dimen108
+)
+(C:\Program Files\MikTex\tex\latex\amsmath\amsopn.sty
+Package: amsopn 1999/12/14 v2.01 operator names
+)
+\inf@bad=\count87
+LaTeX Info: Redefining \frac on input line 211.
+\uproot@=\count88
+\leftroot@=\count89
+LaTeX Info: Redefining \overline on input line 307.
+\classnum@=\count90
+\DOTSCASE@=\count91
+LaTeX Info: Redefining \ldots on input line 379.
+LaTeX Info: Redefining \dots on input line 382.
+LaTeX Info: Redefining \cdots on input line 467.
+\Mathstrutbox@=\box26
+\strutbox@=\box27
+\big@size=\dimen109
+LaTeX Font Info:    Redeclaring font encoding OML on input line 567.
+LaTeX Font Info:    Redeclaring font encoding OMS on input line 568.
+\macc@depth=\count92
+\c@MaxMatrixCols=\count93
+\dotsspace@=\muskip10
+\c@parentequation=\count94
+\dspbrk@lvl=\count95
+\tag@help=\toks17
+\row@=\count96
+\column@=\count97
+\maxfields@=\count98
+\andhelp@=\toks18
+\eqnshift@=\dimen110
+\alignsep@=\dimen111
+\tagshift@=\dimen112
+\tagwidth@=\dimen113
+\totwidth@=\dimen114
+\lineht@=\dimen115
+\@envbody=\toks19
+\multlinegap=\skip44
+\multlinetaggap=\skip45
+\mathdisplay@stack=\toks20
+LaTeX Info: Redefining \[ on input line 2666.
+LaTeX Info: Redefining \] on input line 2667.
+)
+(C:\Program Files\MikTex\tex\latex\base\latexsym.sty
+Package: latexsym 1998/08/17 v2.2e Standard LaTeX package (lasy symbols)
+\symlasy=\mathgroup6
+LaTeX Font Info:    Overwriting symbol font `lasy' in version `bold'
+(Font)                  U/lasy/m/n --> U/lasy/b/n on input line 47.
+) (MEMPHYS-Trigger.aux)
+LaTeX Font Info:    Checking defaults for OML/cmm/m/it on input line 29.
+LaTeX Font Info:    ... okay on input line 29.
+LaTeX Font Info:    Checking defaults for T1/cmr/m/n on input line 29.
+LaTeX Font Info:    ... okay on input line 29.
+LaTeX Font Info:    Checking defaults for OT1/cmr/m/n on input line 29.
+LaTeX Font Info:    ... okay on input line 29.
+LaTeX Font Info:    Checking defaults for OMS/cmsy/m/n on input line 29.
+LaTeX Font Info:    ... okay on input line 29.
+LaTeX Font Info:    Checking defaults for OMX/cmex/m/n on input line 29.
+LaTeX Font Info:    ... okay on input line 29.
+LaTeX Font Info:    Checking defaults for U/cmr/m/n on input line 29.
+LaTeX Font Info:    ... okay on input line 29.
+LaTeX Font Info:    Try loading font information for U+msa on input line 34.
+
+(C:\Program Files\MikTex\tex\latex\amsfonts\umsa.fd
+File: umsa.fd 2002/01/19 v2.2g AMS font definitions
+)
+LaTeX Font Info:    Try loading font information for U+msb on input line 34.
+
+(C:\Program Files\MikTex\tex\latex\amsfonts\umsb.fd
+File: umsb.fd 2002/01/19 v2.2g AMS font definitions
+)
+LaTeX Font Info:    Try loading font information for U+lasy on input line 34.
+
+(C:\Program Files\MikTex\tex\latex\base\ulasy.fd
+File: ulasy.fd 1998/08/17 v2.2e LaTeX symbol font definitions
+)
+Overfull \hbox (1.70615pt too wide) in paragraph at lines 58--58
+[]   \OT1/cmtt/m/n/12 PulseWidth = 5.0*ns this is the time quantum (a pulse wid
+th equivalent)[] 
+ []
+
+
+Overfull \hbox (51.10617pt too wide) in paragraph at lines 58--58
+[]   \OT1/cmtt/m/n/12 CountThreshold = number of PMTs at least to be fired in a
+ TriggerSlidingGate to[] 
+ []
+
+
+Overfull \hbox (2.13571pt too wide) in paragraph at lines 62--63
+\OT1/cmr/m/n/12 there is no cor-re-la-tion be-tween two noise hit (at least on 
+the scale of []\OT1/cmtt/m/n/12 TriggerSlidingGate\OT1/cmr/m/n/12 ),
+ []
+
+[1
+
+]
+File: deltaT.eps Graphic file (type eps)
+ <deltaT.eps>
+File: NumberPMTs.eps Graphic file (type eps)
+ <NumberPMTs.eps>
+File: SumBaseline.eps Graphic file (type eps)
+ <SumBaseline.eps>
+File: SumBaselineZoom.eps Graphic file (type eps)
+ <SumBaselineZoom.eps>
+[2] [3] [4] (MEMPHYS-Trigger.aux) ) 
+Here is how much of TeX's memory you used:
+ 1604 strings out of 95897
+ 17114 string characters out of 1194805
+ 66793 words of memory out of 1066474
+ 4561 multiletter control sequences out of 35000
+ 15969 words of font info for 67 fonts, out of 500000 for 1000
+ 14 hyphenation exceptions out of 607
+ 29i,9n,24p,1116b,299s stack positions out of 1500i,500n,5000p,200000b,32768s
+
+Output written on MEMPHYS-Trigger.dvi (4 pages, 10248 bytes).
Index: /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.tex
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.tex	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/MEMPHYS-Trigger.tex	(revision 388)
@@ -0,0 +1,138 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\documentclass[12pt,a4paper]{article}
+
+%-- used packages ------------------------------------------------------
+
+%%%JEC 3/2/06 does not compile at LAL \usepackage{cite}
+%\usepackage[T1]{fontenc}
+%\usepackage[latin1]{inputenc}
+\usepackage{graphicx}
+\usepackage{epsfig}
+\usepackage{amssymb}
+\usepackage{amsmath}
+\usepackage{latexsym}
+
+\setlength{\textheight}{23cm}
+\setlength{\textwidth}{16cm}
+\setlength{\topmargin}{0cm}
+\setlength{\oddsidemargin}{0cm}
+\setlength{\evensidemargin}{0cm}
+\setlength{\marginparwidth}{0pt}
+
+\def\beq{\begin{equation}}
+\def\eeq{\end{equation}}
+\def\bea{\begin{eqnarray}}
+\def\eea{\end{eqnarray}}
+
+
+
+\begin{document}
+
+\title{Dark Current induced trigger in MEMPHYS\thanks{Preliminary}}
+\author{J.E Campagne - LAL - Orsay}
+\date{\today}
+\maketitle
+
+\begin{abstract}
+Blabla
+\end{abstract}
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\section{Motivations}
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+On 5th of July, J.~Bouchez had demonstrated the need to investigate deeper than expected the dark current induced trigger for a MEMPHYS tank. As a matter of fact, a 3~MeV reactor neutrino produce around 15~p.e, so basically 15 PMTs on a Cerenkov cone are hit, while in the same time at 2~kHz dark current noise per PMT in a 400~ns time window one can expect around 160 PMT hit randomly distributed. Although, it may be not a problem (not yet proved) to recover the PMTs of the signal, the main concern is the trigger.  J.~Bouchez has made a preliminary numerical analysis to show that the question is relevant. 
+A long this line, I have investigated a toy MC and analytical analysis.
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\section{A toy MC}
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+The toy MC consists simply of the generation of dark current random noise on each PMT baseline, then it simulates a digital sum of all the PMT baselines, and applies a threshold on a total number of PMT hit during a sliding window. This toy MC can be adapted according to: 
+\begin{verbatim}
+   NumberOfPMT  = 81,000 for MEMPHYS
+   DarkCurrentRate = 6.0*kHz  more likely for 12"
+   TriggerSlidingGate = 425.0*ns for MEMPHYS
+   PulseWidth = 5.0*ns this is the time quantum (a pulse width equivalent)
+   TimeDurationOfExp = 6000.0*micros (+TriggerSlidingGate) 
+                       as a total time for investigation  
+   CountThreshold = number of PMTs at least to be fired in a TriggerSlidingGate to 
+                    produce a Trigger
+\end{verbatim}
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\subsection{Dark current generation}
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+According to the fact that dark current is a random noise, if ones make the assumption that there is no correlation between two noise hit (at least on the scale of \verb|TriggerSlidingGate|), then the time interval between two hit $\delta t$ follow an exponential probability with parameter $\tau = 1/\verb|DarkCurrentRate|$. So, the probability to get 1 hit in the time window $[t_0,t_0+\Delta t]$ reads:
+\bea
+P(1;[t_0,t_0+\Delta t]) & = & \int_{t_0}^{t_0+\Delta t} e^{(t-t_0)/\tau} \frac{dt}{\tau} \nonumber \\
+	& = & 1 - e^{-\Delta t/\tau} = \frac{\Delta t}{\tau} - \frac{1}{2} \left( \frac{\Delta t}{\tau}\right)^2 \dots
+\label{eq:singleProba}
+\eea
+and is independent of $t_0$. 
+
+If one piles up $N$ PMTs and if $\Delta t/\tau \rightarrow 0$, then it is well known that the probability to get $n$ PMTs hit in a time window of $\Delta t$ tends to a Poisson distribution:
+\beq
+P(n|N;[t,t+\Delta t]) \stackrel{\Delta t/\tau \rightarrow 0}{\longrightarrow}\mathrm{Poisson}(n;N\Delta t/\tau)
+\label{eq:poissonLimit}
+\eeq
+One can see these properties reproduced by the toy MC on Fig.~\ref{fig:PMTPoisson} for some parameters values used on purpose.
+\begin{figure}[htb]
+ \begin{minipage}[c]{0.44\textwidth}
+\includegraphics[width=\textwidth]{deltaT.eps}
+\end{minipage}
+ \begin{minipage}[c]{0.05\textwidth}
+\end{minipage}
+ \begin{minipage}[c]{0.44\textwidth}
+\includegraphics[width=\textwidth]{NumberPMTs.eps}
+\end{minipage}
+\caption{\label{fig:PMTPoisson} {\it left panel}: exponential generation of $\delta t$ with $1/\tau = 1$~MHz dark current noise; {\it right panel}: number of PMTs fired among 50 PMTs, per $\Delta t = 100$~ns windows. One expects a Poisson distribution with a mean of $N\Delta t/\tau = 5$ as it is shown.
+}
+\end{figure}
+
+In fact, as it is mentioned, the Poisson distribution is an asymptotic distribution. The real probability to get exactly $n$ PMTs fired among a collection of $N$ PMTs in a $\Delta t$ time window, with the hypothesis that each PMT  dark current rate has a $1/\tau$ frequency, reads (using Eq.~\ref{eq:singleProba}):
+\bea
+P(n|N;[t,t+\Delta t]) & = & C_N^n P(1;[t,t+\Delta t])^n \left(1-P(1;[t,t+\Delta t])\right)^{N-n} \nonumber \\
+											& = &\mathrm{Binomial}(n,N;P(1,[t,t+\Delta t]))
+\label{eq:BinomialProb}
+\eea
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\subsection{Trigger simulation}
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+One may wish to know the occurrence of "how many times per second, $n_{Th}$ PMTs are fired among $N$ PMTs?". In fact, J.~Bouchez noticed that this question is not complete as one may trigger only if the number of PMTs fired in a $\Delta t$ sliding windows just exceeds $n_{Th}$. That is to say, one should consider the stepping function of the number of PMTs  fired in the sliding window $\Delta t$ examined at $t$, noted $n(t)$, and consider as "trigger" the passage $n(t)=n_{Th} \rightarrow n(t+dt)=n_{Th}+1$ in an infinitesimal time duration $dt$. In fact, looking at Fig.~\ref{fig:DigitalSum}, where is displayed the digital sum of $N=81,000$ PMts with 6~kHz dark current noise rate in sliding window of $425$~ns (computed every $dt=5$~ns), one notices that it is true that the threshold\footnote{mean number of PMTS is 206 and the threshold is set 3 sigma above} at $250$  is reached by steps, but one may also specify that if there is a trigger at $t_1$, before to re-trigger the $n(t>t_1)$ function may pass through a sequence $n(t)=n_{Th}+1 \rightarrow n(t+dt)=n_{Th}$. 
+\begin{figure}[htb]
+ \begin{minipage}[c]{0.44\textwidth}
+\includegraphics[width=\textwidth]{SumBaseline.eps}
+\end{minipage}
+ \begin{minipage}[c]{0.05\textwidth}
+\end{minipage}
+ \begin{minipage}[c]{0.44\textwidth}
+\includegraphics[width=\textwidth]{SumBaselineZoom.eps}
+\end{minipage}
+\caption{\label{fig:DigitalSum} {\it left panel}: Digital sum of the time evolution of $N=81,000$ PMTs with a dark current of $6$~kHz, computed for a time sliding window of $425$~ns every $5$~ns; {\it right panel}: zoom around the location of a trigger threshold at $250$ PMTs.
+}
+\end{figure}
+
+So in practice, one considers as  a trigger the passage $n(t)=n_{Th} \rightarrow n(t+dt)=n_{Th}+1$ and groups together the following $n(t)$ values contiguous in time for which $n(t)>=n_{Th}$. Using this definition, one gets the results for one MEMPHYS tank (Tab.\ref{tab:TriggerRate}): $N=81,000$ PMTs, a sliding window of $425$~ns which corresponds to the maximum time of flight of \v{C}erenkov light in the $\phi: 65~\mathrm{m} \times H:65~\mathrm{m}$ cylinder water volume ($v = 22,6~\mathrm{cm/ns}$). 
+\begin{table}
+\begin{center}
+\begin{tabular}{cccc}\hline\hline
+DC per PMT (kHz) & Mean of PMTs fired & Threshold (at 3$\sigma$)  & Trigger (kHz) & Th (kHz)\\
+1                 &      34           &  52                     & 412             & 57\\
+2                 &      69           &  93                     & 378              & 64.9\\
+3                 &      103           &  133                     & 198            & 2.8 10^3\\
+4                 &      138           &  172                     & 378            & 7.1\\
+5                 &      172           &  211                     & 268            & 3.1\\
+6                 &      206          &  249                     & 416             & 1.5\\ \hline
+DC per PMT (kHz) & Mean of PMTs fired & Threshold  & Trigger \& rel. err (kHz) \\
+1                 &      34           & 40 (1$\sigma$)                 &   30 811 (0.7\%) & 1475.5 \\
+1                 &      34           & 46 (2$\sigma$)                 &    4580 (2\%)    & 363.4\\
+1                 &      34           & 52 (3$\sigma$)                 &     412 (7\%)     & 41.2   \\
+1                 &      34           & 57 (4$\sigma$)                 &      28 (12\%)    & 3.9\\
+\hline\hline
+\end{tabular}
+\caption{\label{tab:TriggerRate} Mean number of PMTs ($m$)fired among $N=81,000$ PMTs on a time sliding window of $\Delta t=425$~ns corresponding to the dark current (DC) rate ($f$) per PMT (this is $N\Delta t * f$). Then, setting a threshold to $m+3\sqrt m$, one reads the trigger rate (see text for the definition of a trigger). The relative errors is also reported as it is not negligeable.}
+\end{center}
+\end{table} 
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+\subsection{Analytical formula}
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+
+\end{document}
Index: /Backup NB/Talks/MEMPHYS-Trigger/PoissApprox.nb
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/PoissApprox.nb	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/PoissApprox.nb	(revision 388)
@@ -0,0 +1,227 @@
+(************** Content-type: application/mathematica **************
+
+                    Mathematica-Compatible Notebook
+
+This notebook can be used with any Mathematica-compatible
+application, such as Mathematica, MathReader or Publicon. The data
+for the notebook starts with the line containing stars above.
+
+To get the notebook into a Mathematica-compatible application, do
+one of the following:
+
+* Save the data starting with the line of stars above into a file
+  with a name ending in .nb, then open the file inside the
+  application;
+
+* Copy the data starting with the line of stars above to the
+  clipboard, then use the Paste menu command inside the application.
+
+Data for notebooks contains only printable 7-bit ASCII and can be
+sent directly in email or through ftp in text mode.  Newlines can be
+CR, LF or CRLF (Unix, Macintosh or MS-DOS style).
+
+NOTE: If you modify the data for this notebook not in a Mathematica-
+compatible application, you must delete the line below containing
+the word CacheID, otherwise Mathematica-compatible applications may
+try to use invalid cache data.
+
+For more information on notebooks and Mathematica-compatible 
+applications, contact Wolfram Research:
+  web: http://www.wolfram.com
+  email: info@wolfram.com
+  phone: +1-217-398-0700 (U.S.)
+
+Notebook reader applications are available free of charge from 
+Wolfram Research.
+*******************************************************************)
+
+(*CacheID: 232*)
+
+
+(*NotebookFileLineBreakTest
+NotebookFileLineBreakTest*)
+(*NotebookOptionsPosition[      4243,        146]*)
+(*NotebookOutlinePosition[      5366,        182]*)
+(*  CellTagsIndexPosition[      5238,        175]*)
+(*WindowFrame->Normal*)
+
+
+
+Notebook[{
+Cell[BoxData[
+    \(\(\(<< Statistics`DiscreteDistributions`\)\(\[IndentingNewLine]\)
+    \)\)], "Input"],
+
+Cell[BoxData[
+    \(PoissPDF[k_, N_, dt_, f_] := 
+      PDF[PoissonDistribution[N*f*dt], k]\)], "Input"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(\(?PoissPDF\)\)], "Input"],
+
+Cell["Global`PoissPDF", "Print",
+  CellTags->"Info3362307184-6683970"],
+
+Cell[BoxData[
+    InterpretationBox[GridBox[{
+          {GridBox[{
+                {\(PoissPDF[k_, N_, dt_, f_] := 
+                    PDF[PoissonDistribution[N\ f\ dt], k]\)}
+                },
+              GridBaseline->{Baseline, {1, 1}},
+              ColumnWidths->0.999,
+              ColumnAlignments->{Left}]}
+          },
+        GridBaseline->{Baseline, {1, 1}},
+        ColumnAlignments->{Left}],
+      Definition[ "PoissPDF"],
+      Editable->False]], "Print",
+  CellTags->"Info3362307184-6683970"]
+}, Open  ]],
+
+Cell[BoxData[
+    \(Trigger[nTh_, N_, dt_, f_] := 
+      Evaluate[Sum[
+            PoissPDF[k, N, 5/10^9, f]*
+              PoissPDF[nTh - k, N, dt, 
+                f]*\((1 - Sum[PoissPDF[j, N, 5/10^9, f], {j, 0, k}])\), {k, 
+              0, nTh}]]*\((10^9/5)\)\)], "Input"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(\(?Trigger\)\)], "Input"],
+
+Cell["Global`Trigger", "Print",
+  CellTags->"Info3362307731-1140847"],
+
+Cell[BoxData[
+    InterpretationBox[GridBox[{
+          {GridBox[{
+                {\(Trigger[nTh_, N_, dt_, f_] := 
+                    1\/5\ Evaluate[\[Sum]\+\(k = 0\)\%nTh PoissPDF[k, N, 
+                              5\/10\^9, f]\ PoissPDF[nTh - k, N, dt, 
+                              f]\ \((1 - \[Sum]\+\(j = 0\)\%k PoissPDF[j, N, 
+                                    5\/10\^9, f])\)]\ 10\^9\)}
+                },
+              GridBaseline->{Baseline, {1, 1}},
+              ColumnWidths->0.999,
+              ColumnAlignments->{Left}]}
+          },
+        GridBaseline->{Baseline, {1, 1}},
+        ColumnAlignments->{Left}],
+      Definition[ "Trigger"],
+      Editable->False]], "Print",
+  CellTags->"Info3362307731-1140847"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(N[Trigger[133, 81000, 425/10^9, 3\ 10^3]]\)], "Input"],
+
+Cell[BoxData[
+    \(57679.52805210537`\)], "Output"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(N[%]\)], "Input"],
+
+Cell[BoxData[
+    \(0.000285805480422623`\)], "Output"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(N[PoissPDF[249, 81000, 425/10^9, 6\ 10^3]*81000*6*10^3]\)], "Input"],
+
+Cell[BoxData[
+    \(205465.67261408278`\)], "Output"]
+}, Open  ]]
+},
+FrontEndVersion->"4.1 for Microsoft Windows",
+ScreenRectangle->{{0, 1024}, {0, 685}},
+WindowSize->{495, 527},
+WindowMargins->{{0, Automatic}, {Automatic, 0}}
+]
+
+(*******************************************************************
+Cached data follows.  If you edit this Notebook file directly, not
+using Mathematica, you must remove the line containing CacheID at
+the top of  the file.  The cache data will then be recreated when
+you save this file from within Mathematica.
+*******************************************************************)
+
+(*CellTagsOutline
+CellTagsIndex->{
+  "Info3362307184-6683970"->{
+    Cell[1991, 63, 70, 1, 25, "Print",
+      CellTags->"Info3362307184-6683970"],
+    Cell[2064, 66, 512, 14, 44, "Print",
+      CellTags->"Info3362307184-6683970"]},
+  "Info3362307731-1140847"->{
+    Cell[2940, 96, 69, 1, 25, "Print",
+      CellTags->"Info3362307731-1140847"],
+    Cell[3012, 99, 736, 17, 77, "Print",
+      CellTags->"Info3362307731-1140847"]}
+  }
+*)
+
+(*CellTagsIndex
+CellTagsIndex->{
+  {"Info3362307184-6683970", 4849, 161},
+  {"Info3362307731-1140847", 5046, 166}
+  }
+*)
+
+(*NotebookFileOutline
+Notebook[{
+Cell[1705, 50, 105, 2, 50, "Input"],
+Cell[1813, 54, 104, 2, 50, "Input"],
+
+Cell[CellGroupData[{
+Cell[1942, 60, 46, 1, 30, "Input"],
+Cell[1991, 63, 70, 1, 25, "Print",
+  CellTags->"Info3362307184-6683970"],
+Cell[2064, 66, 512, 14, 44, "Print",
+  CellTags->"Info3362307184-6683970"]
+}, Open  ]],
+Cell[2591, 83, 276, 6, 130, "Input"],
+
+Cell[CellGroupData[{
+Cell[2892, 93, 45, 1, 30, "Input"],
+Cell[2940, 96, 69, 1, 25, "Print",
+  CellTags->"Info3362307731-1140847"],
+Cell[3012, 99, 736, 17, 77, "Print",
+  CellTags->"Info3362307731-1140847"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+Cell[3785, 121, 74, 1, 30, "Input"],
+Cell[3862, 124, 52, 1, 29, "Output"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+Cell[3951, 130, 37, 1, 30, "Input"],
+Cell[3991, 133, 55, 1, 29, "Output"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+Cell[4083, 139, 88, 1, 30, "Input"],
+Cell[4174, 142, 53, 1, 29, "Output"]
+}, Open  ]]
+}
+]
+*)
+
+
+
+(*******************************************************************
+End of Mathematica Notebook file.
+*******************************************************************)
+
Index: /Backup NB/Talks/MEMPHYS-Trigger/PoissonApp2.nb
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/PoissonApp2.nb	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/PoissonApp2.nb	(revision 388)
@@ -0,0 +1,207 @@
+(************** Content-type: application/mathematica **************
+
+                    Mathematica-Compatible Notebook
+
+This notebook can be used with any Mathematica-compatible
+application, such as Mathematica, MathReader or Publicon. The data
+for the notebook starts with the line containing stars above.
+
+To get the notebook into a Mathematica-compatible application, do
+one of the following:
+
+* Save the data starting with the line of stars above into a file
+  with a name ending in .nb, then open the file inside the
+  application;
+
+* Copy the data starting with the line of stars above to the
+  clipboard, then use the Paste menu command inside the application.
+
+Data for notebooks contains only printable 7-bit ASCII and can be
+sent directly in email or through ftp in text mode.  Newlines can be
+CR, LF or CRLF (Unix, Macintosh or MS-DOS style).
+
+NOTE: If you modify the data for this notebook not in a Mathematica-
+compatible application, you must delete the line below containing
+the word CacheID, otherwise Mathematica-compatible applications may
+try to use invalid cache data.
+
+For more information on notebooks and Mathematica-compatible 
+applications, contact Wolfram Research:
+  web: http://www.wolfram.com
+  email: info@wolfram.com
+  phone: +1-217-398-0700 (U.S.)
+
+Notebook reader applications are available free of charge from 
+Wolfram Research.
+*******************************************************************)
+
+(*CacheID: 232*)
+
+
+(*NotebookFileLineBreakTest
+NotebookFileLineBreakTest*)
+(*NotebookOptionsPosition[      4203,        143]*)
+(*NotebookOutlinePosition[      5092,        173]*)
+(*  CellTagsIndexPosition[      5005,        167]*)
+(*WindowFrame->Normal*)
+
+
+
+Notebook[{
+Cell[BoxData[
+    \(<< Statistics`DiscreteDistributions`\)], "Input",
+  FontFamily->"Arial TUR",
+  FontSize->16],
+
+Cell[BoxData[
+    \(PoissPDF[k_, N_, dt_, f_] := 
+      PDF[PoissonDistribution[N*f*dt], k]\)], "Input",
+  FontFamily->"Arial TUR",
+  FontSize->16],
+
+Cell[BoxData[
+    \(Trigger0[nTh_, N_, dt_, f_] := 
+      Evaluate[Sum[
+            PoissPDF[k, N, 5/10^9, f]*
+              PoissPDF[nTh - k, N, dt, 
+                f]*\((1 - Sum[PoissPDF[j, N, 5/10^9, f], {j, 0, k}])\), {k, 
+              0, nTh}]]*\((10^9/5)\)\)], "Input",
+  FontFamily->"Arial TUR",
+  FontSize->16],
+
+Cell[BoxData[
+    \(Trigger[nTh_, N_, dt_, f_] := 
+      Evaluate[Sum[
+            PoissPDF[k, N, 0.0001/10^9, f]*
+              PoissPDF[nm - k, N, dt - 0.0001/10^9, 
+                f]*\((1 - 
+                  Sum[PoissPDF[j, N, 0.0001/10^9, f], {j, 0, 
+                      k + nTh - nm}])\), {nm, 0, nTh}, {k, 0, nm}]]*\((10^9/
+            0.0001)\)\)], "Input",
+  FontFamily->"Arial TUR",
+  FontSize->16],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(\(?Trigger\)\)], "Input",
+  FontFamily->"Arial TUR",
+  FontSize->16],
+
+Cell["Global`Trigger", "Print",
+  FontFamily->"Arial TUR",
+  FontSize->16,
+  CellTags->"Info3362411253-9063905"],
+
+Cell[BoxData[
+    InterpretationBox[GridBox[{
+          {GridBox[{
+                {\(Trigger[nTh_, N_, dt_, f_] := 
+                    1\/5\ Evaluate[\[Sum]\+\(nm = 0\)\%nTh\(\[Sum]\+\(k = 0\)\
+\%nm PoissPDF[k, N, 5\/10\^9, f]\ PoissPDF[nm - k, N, dt - 5\/10\^9, 
+                                f]\ \((1 - \[Sum]\+\(j = 0\)\%\(k + nTh - \
+nm\)PoissPDF[j, N, 5\/10\^9, f])\)\)]\ 10\^9\)}
+                },
+              GridBaseline->{Baseline, {1, 1}},
+              ColumnWidths->0.999,
+              ColumnAlignments->{Left}]}
+          },
+        GridBaseline->{Baseline, {1, 1}},
+        ColumnAlignments->{Left}],
+      Definition[ "Trigger"],
+      Editable->False]], "Print",
+  FontFamily->"Arial TUR",
+  FontSize->16,
+  CellTags->"Info3362411253-9063905"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(N[Trigger0[52, 81000, 425/10^9, 10^3]]\)], "Input",
+  FontFamily->"Arial TUR",
+  FontSize->16],
+
+Cell[BoxData[
+    \(57161.0960845246`\)], "Output",
+  FontFamily->"Arial TUR",
+  FontSize->16]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+    \(N[Trigger[52, 81000, 425/10^9, 1\ 10^3]]\)], "Input",
+  FontFamily->"Arial TUR",
+  FontSize->16],
+
+Cell[BoxData[
+    \(93044.35202434182`\)], "Output"]
+}, Open  ]],
+
+Cell[BoxData[""], "Input"]
+},
+FrontEndVersion->"4.1 for Microsoft Windows",
+ScreenRectangle->{{0, 1024}, {0, 685}},
+WindowSize->{805, 527},
+WindowMargins->{{50, Automatic}, {Automatic, -8}}
+]
+
+(*******************************************************************
+Cached data follows.  If you edit this Notebook file directly, not
+using Mathematica, you must remove the line containing CacheID at
+the top of  the file.  The cache data will then be recreated when
+you save this file from within Mathematica.
+*******************************************************************)
+
+(*CellTagsOutline
+CellTagsIndex->{
+  "Info3362411253-9063905"->{
+    Cell[2820, 90, 112, 3, 28, "Print",
+      CellTags->"Info3362411253-9063905"],
+    Cell[2935, 95, 767, 19, 128, "Print",
+      CellTags->"Info3362411253-9063905"]}
+  }
+*)
+
+(*CellTagsIndex
+CellTagsIndex->{
+  {"Info3362411253-9063905", 4811, 158}
+  }
+*)
+
+(*NotebookFileOutline
+Notebook[{
+Cell[1705, 50, 112, 3, 33, "Input"],
+Cell[1820, 55, 147, 4, 33, "Input"],
+Cell[1970, 61, 320, 8, 81, "Input"],
+Cell[2293, 71, 411, 10, 105, "Input"],
+
+Cell[CellGroupData[{
+Cell[2729, 85, 88, 3, 33, "Input"],
+Cell[2820, 90, 112, 3, 28, "Print",
+  CellTags->"Info3362411253-9063905"],
+Cell[2935, 95, 767, 19, 128, "Print",
+  CellTags->"Info3362411253-9063905"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+Cell[3739, 119, 114, 3, 33, "Input"],
+Cell[3856, 124, 94, 3, 32, "Output"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+Cell[3987, 132, 116, 3, 33, "Input"],
+Cell[4106, 137, 52, 1, 29, "Output"]
+}, Open  ]],
+Cell[4173, 141, 26, 0, 30, "Input"]
+}
+]
+*)
+
+
+
+(*******************************************************************
+End of Mathematica Notebook file.
+*******************************************************************)
+
Index: /Backup NB/Talks/MEMPHYS-Trigger/Trigger.nb
===================================================================
--- /Backup NB/Talks/MEMPHYS-Trigger/Trigger.nb	(revision 388)
+++ /Backup NB/Talks/MEMPHYS-Trigger/Trigger.nb	(revision 388)
@@ -0,0 +1,9 @@
+Trigger[n_, N_, dt_, f_] := P[n, N, dt, f]*P[0, N, 5/10^9, f]*
+     P[1, N, 5/10^9, f]*(10^9/5)
+ 
+P[n_, N_, dt_, f_] := PDF[BinomialDistribution[N, p[dt, f]], n]
+ 
+PDF /: PDF::usage = "PDF[distribution, x] gives the probability density \
+function of the specified statistical distribution evaluated at x."
+ 
+p[dt_, f_] := 1 - Exp[(-dt)*f]
