Do not be overzelous with bin times, it produces unwanted side effects

StingraySoftware · Feb 2, 2024 · 32d2edb · 32d2edb
1 parent 06fc297
commit 32d2edb
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 6 deletions.
diff --git a/hendrics/fspec.py b/hendrics/fspec.py
@@ -317,7 +317,6 @@ def calc_pds(
 
     if pdsrebin is not None and pdsrebin != 1:
         pds = pds.rebin(pdsrebin)
-
     pds.instr = instr
     pds.fftlen = fftlen
     pds.back_phots = back_ctrate * fftlen
@@ -860,7 +859,7 @@ def main(args=None):
     log.setLevel(args.loglevel)
 
     with log.log_to_file("HENfspec.log"):
-        bintime = np.longdouble(interpret_bintime(args.bintime))
+        bintime = interpret_bintime(args.bintime)
 
         fftlen = args.fftlen
         pdsrebin = args.rebin

diff --git a/hendrics/lcurve.py b/hendrics/lcurve.py
@@ -354,7 +354,7 @@ def lcurve_from_events(
         evdata = deorbit_events(evdata, deorbit_par)
         deorbit_tag = "_deorb"
 
-    bintime = np.longdouble(interpret_bintime(bintime))
+    bintime = interpret_bintime(bintime)
 
     tag = ""
 

diff --git a/hendrics/tests/test_fspec.py b/hendrics/tests/test_fspec.py
@@ -145,7 +145,7 @@ def setup_class(cls):
         )
         hen.read_events.main(command.split())
 
-        command = "{0} {1}  --nproc 2 -b 0.5".format(cls.ev_fileA, cls.ev_fileB)
+        command = "{0} {1}  --nproc 2 -b -1".format(cls.ev_fileA, cls.ev_fileB)
         hen.lcurve.main(command.split())
         cls.lcA = cls.ev_fileA.replace("_ev", "_lc")
         cls.lcB = cls.ev_fileB.replace("_ev", "_lc")
@@ -219,14 +219,14 @@ def test_pds_leahy(self):
         evdata = self.ev_fileA
         lcdata = self.lcA
 
-        command = "{0} -f 128 -k PDS --norm leahy -b 0.5".format(evdata)
+        command = "{0} -f 128. -k PDS --norm leahy -b -1".format(evdata)
         hen.fspec.main(command.split())
         evout = evdata.replace("_ev", "_pds")
         assert os.path.exists(evout)
         evpds = hen.io.load_pds(evout)
         io.remove_pds(evout)
 
-        command = "{0} -f 128 -k PDS --save-all --norm leahy".format(lcdata)
+        command = "{0} -f 128. -k PDS --save-all --norm leahy".format(lcdata)
         hen.fspec.main(command.split())
         lcout = lcdata.replace("_lc", "_pds")
         assert os.path.exists(lcout)