From 9dfda7afd4372d679ca5534ab6f04020d7c64fbc Mon Sep 17 00:00:00 2001
From: Michael Mayer <michael.mayer@desy.de>
Date: Wed, 6 Apr 2016 15:58:31 +0200
Subject: [PATCH] Added seed parameter to cspull

---
 cscripts/cspull.par | 2 ++
 cscripts/cspull.py  | 6 ++++--
 2 files changed, 6 insertions(+), 2 deletions(-)

diff --git a/cscripts/cspull.par b/cscripts/cspull.par
index 638886f..0968d4b 100644
--- a/cscripts/cspull.par
+++ b/cscripts/cspull.par
@@ -44,6 +44,8 @@ rad,      r, h, 5.0,0,180, "Radius of ROI (deg)"
 pattern,  s, h, single,single|four,, "Observation pattern"
 offset,   r, h, 1.5,,, "Observation pattern offset (deg)"
 profile,  b, h, no,,, "Use likelihood profile method for errors?"
+seed,     i, h, 1,,, "Initial random number generator seed"
+
 
 #
 # Standard parameters
diff --git a/cscripts/cspull.py b/cscripts/cspull.py
index ff8bbfe..ab02524 100755
--- a/cscripts/cspull.py
+++ b/cscripts/cspull.py
@@ -64,6 +64,7 @@ class cspull(ctools.cscript):
         self._binsz       = 0.0
         self._pattern     = "single"
         self._enumbins    = 0
+        self._seed        = 1
         self._chatter     = 2
 
         # Initialise observation container from constructor arguments.
@@ -116,6 +117,7 @@ class cspull(ctools.cscript):
         self._ntrials = self["ntrials"].integer()
         self._edisp   = self["edisp"].boolean()
         self._offset  = self["offset"].real()
+        self._seed    = self["seed"].integer()
         self._chatter = self["chatter"].integer()
 
         # Query some parameters
@@ -410,8 +412,8 @@ class cspull(ctools.cscript):
         # Loop over trials
         for seed in range(self._ntrials):
 
-            # Make a trial
-            result = self._trial(seed)
+            # Make a trial and add initial seed
+            result = self._trial(seed + self._seed)
 
             # Write out result immediately
             if seed == 0:
--
libgit2 0.21.4