#!/usr/bin/python2 -u
# The -u option above turns off block buffering of python output. This 
# assures that each error message gets individually printed to the log file.
#
# Module: nodepowerAgent.py
#
# Description: This module acts as an agent between the mesh network and
# node power and enviromental sensors.  The agent periodically polls the
# sensors and processes the data returned from the sensors, including
#     - conversion of data items
#     - update a round robin (rrdtool) database with the sensor data
#     - periodically generate graphic charts for display in html documents
#     - write the processed node status to a JSON file for use by html
#       documents
#
# Copyright 2021 Jeff Owrey
#    This program is free software: you can redistribute it and/or modify
#    it under the terms of the GNU General Public License as published by
#    the Free Software Foundation, either version 3 of the License, or
#    (at your option) any later version.
#
#    This program is distributed in the hope that it will be useful,
#    but WITHOUT ANY WARRANTY; without even the implied warranty of
#    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
#    GNU General Public License for more details.
#
#    You should have received a copy of the GNU General Public License
#    along with this program.  If not, see http://www.gnu.org/license.
#
# Revision History
#   * v10 released 01 June 2021 by J L Owrey; first release
#
#2345678901234567890123456789012345678901234567890123456789012345678901234567890

# Import required python libraries.
import os
import sys
import signal
import subprocess
import multiprocessing
import time
import json

# Import sensor libraries.
import ina260 # power sensor
import tmp102 # temperature sensor

    ### ENVIRONMENT ###
_USER = os.environ['USER']

    ### SENSOR BUS ADDRESSES ###

# Set bus addresses of sensors.
_PWR_SENSOR_ADDR = 0X40
_BAT_TMP_SENSOR_ADDR = 0x48
_AMB_TMP_SENSOR_ADDR = 0x4B
# Set bus selector.
_BUS_SEL = 1

    ### FILE AND FOLDER LOCATIONS ###

# folder to contain html
_DOCROOT_PATH = "/home/%s/public_html/power/" % _USER
# folder to contain charts and output data file
_CHARTS_DIRECTORY = _DOCROOT_PATH + "dynamic/"
# location of JSON output data file
_OUTPUT_DATA_FILE = _DOCROOT_PATH + "dynamic/powerData.js"
# database that stores node data
_RRD_FILE = "/home/%s/database/powerData.rrd" % _USER

    ### GLOBAL CONSTANTS ###

# sensor data request interval in seconds
_DEFAULT_DATA_REQUEST_INTERVAL = 2
# rrdtool database update interval in seconds
_DATABASE_UPDATE_INTERVAL = 30
# chart update interval in seconds
_CHART_UPDATE_INTERVAL = 600
# standard chart width in pixels
_CHART_WIDTH = 600
# standard chart height in pixels
_CHART_HEIGHT = 150
# chart average line color
_AVERAGE_LINE_COLOR = '#006600'

   ### GLOBAL VARIABLES ###

# debug output options
debugOption = False
verboseDebug = False

# frequency of data requests to sensors
dataRequestInterval = _DEFAULT_DATA_REQUEST_INTERVAL
# how often charts get updated
chartUpdateInterval = _CHART_UPDATE_INTERVAL
# last node request time
lastDataPointTime = -1

# Define each sensor.  This also initialzes each sensor.
pwr = ina260.ina260(_PWR_SENSOR_ADDR, _BUS_SEL)
btmp = tmp102.tmp102(_BAT_TMP_SENSOR_ADDR, _BUS_SEL)
atmp = tmp102.tmp102(_AMB_TMP_SENSOR_ADDR, _BUS_SEL)

  ###  PRIVATE METHODS  ###

def getTimeStamp():
    """
    Get the local time and format as a text string.
    Parameters: none
    Returns: string containing the time stamp
    """
    return time.strftime( "%m/%d/%Y %T", time.localtime() )
## end def

def getEpochSeconds(sTime):
    """
    Convert the time stamp to seconds since 1/1/1970 00:00:00.
    Parameters: 
        sTime - the time stamp to be converted must be formatted
                   as %m/%d/%Y %H:%M:%S
    Returns: epoch seconds
    """
    try:
        t_sTime = time.strptime(sTime, '%m/%d/%Y %H:%M:%S')
    except Exception as exError:
        print('%s getEpochSeconds: %s' % (getTimeStamp(), exError))
        return None
    tSeconds = int(time.mktime(t_sTime))
    return tSeconds
## end def

def terminateAgentProcess(signal, frame):
    """
    Send a message to the log when the agent process gets killed
    by the operating system.  Inform downstream clients
    by removing output data files.
    Parameters:
        signal, frame - dummy parameters
    Returns: nothing
    """
    # Inform downstream clients by removing output data file.
    if os.path.exists(_OUTPUT_DATA_FILE):
       os.remove(_OUTPUT_DATA_FILE)
    print('%s terminating npw agent process' % getTimeStamp())
    sys.exit(0)
## end def

  ###  PUBLIC METHODS  ###

def getSensorData(dData):
    """
    Poll sensors for data. Store the data in a dictionary object for
    use by other subroutines.  The dictionary object passed in should
    an empty dictionary, i.e., dData = { }.
    Parameters: dData - a dictionary object to contain the sensor data
    Returns: True if successful, False otherwise
    """
    dData["time"] = getTimeStamp()
 
    try:
        dData["current"] = pwr.getCurrent()
        dData["voltage"] = pwr.getVoltage()
        dData["power"] = pwr.getPower()
        dData["battemp"] = btmp.getTempF()
        dData["ambtemp"] = atmp.getTempF()
    except Exception as exError:
        print("%s sensor error: %s" % (getTimeStamp(), exError))
        return False

    return True
## end def

def updateDatabase(dData):
    """
    Update the rrdtool database by executing an rrdtool system command.
    Format the command using the data extracted from the sensors.
    Parameters: dData - dictionary object containing data items to be
                        written to the rr database file
    Returns: True if successful, False otherwise
    """
 
    epochTime = getEpochSeconds(dData['time'])

    # Format the rrdtool update command.
    strFmt = "rrdtool update %s %s:%s:%s:%s:%s:%s"
    strCmd = strFmt % (_RRD_FILE, epochTime, dData['current'], \
             dData['voltage'], dData['power'], dData['battemp'], \
             dData['ambtemp'])

    if verboseDebug:
        print("%s" % strCmd) # DEBUG

    # Run the command as a subprocess.
    try:
        subprocess.check_output(strCmd, shell=True, \
            stderr=subprocess.STDOUT)
    except subprocess.CalledProcessError as exError:
        print("%s: rrdtool update failed: %s" % \
            (getTimeStamp(), exError.output))
        return False

    if debugOption and not verboseDebug:
        print("database updated")

    return True
## end def

def writeOutputFile(dData):
    """
    Write sensor data items to the output data file, formatted as 
    a Javascript file.  This file may then be requested and used by
    by downstream clients, for instance, an HTML document.
    Parameters:
        dData - a dictionary containing the data to be written
                   to the output data file
        Returns: True if successful, False otherwise
    """
    # Write a JSON formatted file for use by html clients.  The following
    # data items are sent to the client file.
    #    * The last database update date and time
    #    * The data request interval
    #    * The sensor values

    # Create a JSON formatted string from the sensor data.
    try:
        jsData = json.loads("{}")
        for key in dData:
            jsData.update({key:dData[key]})
        jsData.update({"chartUpdateInterval": chartUpdateInterval})
        sData = "[%s]" % json.dumps(jsData)
    except Exception as exError:
        print("%s writeOutputFile: %s" % (getTimeStamp(), exError))
        return False

    if verboseDebug:
        print(sData)

    # Write the JSON formatted data to the output data file.

    try:
        fc = open(_OUTPUT_DATA_FILE, "w")
        fc.write(sData)
        fc.close()
    except Exception as exError:
        print("%s write output file failed: %s" % \
              (getTimeStamp(), exError))
        return False

    return True
## end def

def createGraph(fileName, dataItem, gLabel, gTitle, gStart,
                lower=0, upper=0, trendLine=0, scaleFactor=1,
                autoScale=True, alertLine=""):
    """
    Uses rrdtool to create a graph of specified sensor data item.
    Parameters:
        fileName - name of file containing the graph
        dataItem - data item to be graphed
        gLabel - string containing a graph label for the data item
        gTitle - string containing a title for the graph
        gStart - beginning time of the graphed data
        lower - lower bound for graph ordinate #NOT USED
        upper - upper bound for graph ordinate #NOT USED
        trendLine 
            0, show only graph data
            1, show only a trend line
            2, show a trend line and the graph data
        scaleFactor - amount to pre-scale the data before charting
            the data [default=1]
        autoScale - if True, then use vertical axis auto scaling
            (lower and upper parameters must be zero)
        alertLine - value for which to print a critical
            low voltage alert line on the chart. If not provided
            alert line will not be printed.
    Returns: True if successful, False otherwise
    """
    gPath = _CHARTS_DIRECTORY + fileName + ".png"
    trendWindow = { 'end-1day': 7200,
                    'end-4weeks': 172800,
                    'end-12months': 604800 }
 
    # Format the rrdtool graph command.

    # Set chart start time, height, and width.
    strCmd = "rrdtool graph %s -a PNG -s %s -e now -w %s -h %s " \
             % (gPath, gStart, _CHART_WIDTH, _CHART_HEIGHT)
   
    # Set the range and scaling of the chart y-axis.
    if lower < upper:
        strCmd  +=  "-l %s -u %s -r " % (lower, upper)
    elif autoScale:
        strCmd += "-A "
    strCmd += "-Y "

    # Set the chart ordinate label and chart title. 
    strCmd += "-v %s -t %s " % (gLabel, gTitle)
 
    # Show the data, or a moving average trend line over
    # the data, or both.
    strCmd += "DEF:rSeries=%s:%s:LAST " % (_RRD_FILE, dataItem)
    strCmd += "CDEF:dSeries=rSeries,%s,/ " % (scaleFactor)

    if trendLine == 0:
        strCmd += "LINE1:dSeries#0400ff "
    elif trendLine == 1:
        strCmd += "CDEF:smoothed=dSeries,%s,TREND LINE2:smoothed%s " \
                  % (trendWindow[gStart], _AVERAGE_LINE_COLOR)
    elif trendLine == 2:
        strCmd += "LINE1:dSeries#0400ff "
        strCmd += "CDEF:smoothed=dSeries,%s,TREND LINE2:smoothed%s " \
                  % (trendWindow[gStart], _AVERAGE_LINE_COLOR)

    if alertLine != "":
        strCmd += "HRULE:%s#FF0000:Critical\ Low\ Voltage " % (alertLine)
     
    if verboseDebug:
        print("%s" % strCmd) # DEBUG
    
    # Run the formatted rrdtool command as a subprocess.
    try:
        result = subprocess.check_output(strCmd, \
                     stderr=subprocess.STDOUT,   \
                     shell=True)
    except subprocess.CalledProcessError as exError:
        print("rrdtool graph failed: %s" % (exError.output))
        return False

    if debugOption:
        print("rrdtool graph: %s" % result.decode('utf-8'))
    return True

## end def

def generateGraphs():
    """
    Generate graphs for display in html documents.
    Parameters: none
    Returns: nothing
    """

    # 24 hour stock charts

    createGraph('24hr_current', 'CUR', 'Amps',
                'Current\ -\ Last\ 24\ Hours', 'end-1day', \
                0, 0, 2, 1000)
    createGraph('24hr_voltage', 'VOLT', 'Volts',
                'Voltage\ -\ Last\ 24\ Hours', 'end-1day', \
                9, 15, 0, 1, True, 11)
    createGraph('24hr_power', 'PWR', 'Watts', 
                'Power\ -\ Last\ 24\ Hours', 'end-1day', \
                0, 0, 2, 1000)
    createGraph('24hr_battemp', 'BTMP', 'deg\ F', 
                'Battery\ Temperature\ -\ Last\ 24\ Hours', 'end-1day', \
                0, 0, 0)
    createGraph('24hr_ambtemp', 'ATMP', 'deg\ F', 
                'Ambient\ Temperature\ -\ Last\ 24\ Hours', 'end-1day', \
                0, 0, 0)

    # 4 week stock charts

    createGraph('4wk_current', 'CUR', 'Amps',
                'Current\ -\ Last\ 4\ Weeks', 'end-4weeks', \
                0, 0, 2, 1000)
    createGraph('4wk_voltage', 'VOLT', 'Volts',
                'Voltage\ -\ Last\ 4\ Weeks', 'end-4weeks', \
                9, 15, 0, 1, True, 11)
    createGraph('4wk_power', 'PWR', 'Watts', 
                'Power\ -\ Last\ 4\ Weeks', 'end-4weeks', \
                0, 0, 2, 1000)
    createGraph('4wk_battemp', 'BTMP', 'deg\ F', 
                'Battery\ Temperature\ -\ Last\ 4\ Weeks', 'end-4weeks', \
                0, 0, 2)
    createGraph('4wk_ambtemp', 'ATMP', 'deg\ F', 
                'Ambient\ Temperature\ -\ Last\ 4\ Weeks', 'end-4weeks', \
                0, 0, 2)

    # 12 month stock charts

    createGraph('12m_current', 'CUR', 'Amps',
                'Current\ -\ Past\ Year', 'end-12months', \
                0, 0, 2, 1000)
    createGraph('12m_voltage', 'VOLT', 'Volts',
                'Voltage\ -\ Past\ Year', 'end-12months', \
                9, 15, 0, 1, True, 11)
    createGraph('12m_power', 'PWR', 'Watts', 
                'Power\ -\ Past\ Year', 'end-12months', \
                0, 0, 2, 1000)
    createGraph('12m_battemp', 'BTMP', 'deg\ F', 
                'Battery\ Temperature\ -\ Past\ Year', 'end-12months', \
                0, 0, 2)
    createGraph('12m_ambtemp', 'ATMP', 'deg\ F', 
                'Ambient\ Temperature\ -\ Past\ Year', 'end-12months', \
                0, 0, 2)
## end def

def getCLarguments():
    """
    Get command line arguments.  There are three possible arguments
        -d turns on debug mode
        -v turns on verbose debug mode
        -p sets the sensor query period
        -c sets the chart update period
    Returns: nothing
    """
    global debugOption, verboseDebug, dataRequestInterval, chartUpdateInterval

    index = 1
    while index < len(sys.argv):
        if sys.argv[index] == '-d':
            debugOption = True
        elif sys.argv[index] == '-v':
            debugOption = True
            verboseDebug = True
        elif sys.argv[index] == '-p':
            try:
                dataRequestInterval = abs(int(sys.argv[index + 1]))
            except:
                print("invalid sensor query period")
                exit(-1)
            index += 1
        elif sys.argv[index] == '-c':
            try:
                chartUpdateInterval = abs(int(sys.argv[index + 1]))
            except:
                print("invalid chart update period")
                exit(-1)
            index += 1
        else:
            cmd_name = sys.argv[0].split('/')
            print("Usage: %s [-d | v] [-p seconds] [-c seconds]" \
                  % cmd_name[-1])
            exit(-1)
        index += 1
##end def

def main():
    """
    Handles timing of events and acts as executive routine managing
    all other functions.
    Parameters: none
    Returns: nothing
    """
    global dataRequestInterval

    signal.signal(signal.SIGTERM, terminateAgentProcess)
    signal.signal(signal.SIGINT, terminateAgentProcess)

    print('%s starting up node power agent process' % getTimeStamp())

    # last time output JSON file updated
    lastDataRequestTime = -1
    # last time charts generated
    lastChartUpdateTime = - 1
    # last time the rrdtool database updated
    lastDatabaseUpdateTime = -1

    ## Get command line arguments.
    getCLarguments()

    ## Exit with error if rrdtool database does not exist.
    if not os.path.exists(_RRD_FILE):
        print('rrdtool database does not exist\n' \
              'use createPowerRrd script to ' \
              'create rrdtool database\n')
        exit(1)
 
    ## main loop
    while True:

        currentTime = time.time() # get current time in seconds

        # Every data request interval read the sensors and process the
        # data from the sensors.
        if currentTime - lastDataRequestTime > dataRequestInterval:
            lastDataRequestTime = currentTime
            dData = {}
            result = True

            # Get the data from the sensors.
            result =getSensorData(dData)
 
            # If get data successful, write data to data files.
            if result:
                result = writeOutputFile(dData)

            # At the rrdtool database update interval, update the database.
            if result and (currentTime - lastDatabaseUpdateTime > \
                           _DATABASE_UPDATE_INTERVAL):   
                lastDatabaseUpdateTime = currentTime
                ## Update the round robin database with the parsed data.
                result = updateDatabase(dData)

        # At the chart generation interval, generate charts.
        if currentTime - lastChartUpdateTime > chartUpdateInterval:
            lastChartUpdateTime = currentTime
            p = multiprocessing.Process(target=generateGraphs, args=())
            p.start()
            
        # Relinquish processing back to the operating system until
        # the next update interval.

        elapsedTime = time.time() - currentTime
        if debugOption:
            if result:
                print("update successful: %6f sec\n"
                      % elapsedTime)
            else:
                print("update failed: %6f sec\n"
                      % elapsedTime)
        remainingTime = dataRequestInterval - elapsedTime
        if remainingTime > 0.0:
            time.sleep(remainingTime)
    ## end while
    return
## end def

if __name__ == '__main__':
    main()