Changeset 18291 for branches/trigger_burst_research

Timestamp:

08/25/15 12:08:59 (9 years ago)

Author:

dneise

Message:

some mods, so the stuff can maybe also be executed by somebody else :-|

File:

• branches/trigger_burst_research/sketch.py (modified) (3 diffs)

branches/trigger_burst_research/sketch.py

@@ -1 +1 @@
 # coding: utf-8
-from collections import namedtuple
-from sqlalchemy import create_engine
-import pandas as pd
+import calendar
 import numpy as np
 from astropy.io import fits
-import fact
+import pandas as pd
+from sqlalchemy import create_engine
 
 night_int = 20150721
+
+
+def fjd(datetime_inst):
+    """ convert datetime instance to FJD
+    """
+    return calendar.timegm(datetime_inst.utctimetuple()) / (24.*3600.)
 
 def create_DB_connection():
@@ -52 +57 @@
 
 
-def get_trigger_rates(night_int, base_path='other_aux_data/'):
+def get_trigger_rates(night_int, base_path='/fact/aux/'):
+    night_string = str(night_int)
     fits_file = fits.open(
-        base_path+'{0}.FTM_CONTROL_TRIGGER_RATES.fits'.format(night_int)
+        base_path+'{y}/{m}/{d}/{n}.FTM_CONTROL_TRIGGER_RATES.fits'.format(
+            ni=night_string,
+            y=night_string[0:4],
+            m=night_string[4:6],
+            d=night_string[6:8],
+            )
         )
     trigger_rates = fits_file[1].data
@@ -71 +82 @@
         this_run = data_runs.ix[dataframe_index]
         mask = (
-            (trigger_rates['Time'] > fact.fjd(this_run['fRunStart'])) *
-            (trigger_rates['Time'] < fact.fjd(this_run['fRunStop']))
+            (trigger_rates['Time'] > fjd(this_run['fRunStart'])) *
+            (trigger_rates['Time'] < fjd(this_run['fRunStop']))
             )
         this_board_rates = trigger_rates['BoardRate'][mask]