#! /usr/bin/env python
# ==========================================================================
# This scripts performs tests processing pipelines.
#
# Copyright (C) 2014 Juergen Knoedlseder
#
# 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/licenses/>.
#
# ==========================================================================
import gammalib
import ctools


# ============================ #
# Test class for ctobssim tool #
# ============================ #
class Test(gammalib.GPythonTestSuite):
    """
    Test class for pipelines.
    """
    # Constructor
    def __init__(self):
        """
        Constructor.
        """
        # Call base class constructor
        gammalib.GPythonTestSuite.__init__(self)

        # Return
        return

    # Set test functions
    def set(self):
        """
        Set all test functions.
        """
        # Set test name
        self.name("pipelines")

        # Append tests
        self.append(self.test_unbinned_fits, "Test unbinned pipeline with FITS file saving")
        self.append(self.test_unbinned_mem, "Test unbinned in-memory pipeline")

        # Return
        return

    # Test unbinned pipeline with FITS file saving
    def test_unbinned_fits(self):
        """
        Test unbinned pipeline with FITS file saving.
        """
        # Set script parameters
        model_name           = "data/crab.xml"
        events_name          = "events.fits"
        selected_events_name = "selected_events.fits"
        result_name          = "results.xml"
        caldb                = "irf"
        irf                  = "cta_dummy_irf"
        ra                   =   83.63
        dec                  =   22.01
        rad_sim              =   10.0
        tstart               =    0.0
        tstop                = 1800.0
        emin                 =    0.1
        emax                 =  100.0
        rad_select           =    3.0
    
        # Simulate events
        sim = ctools.ctobssim()
        sim["infile"].filename(model_name)
        sim["outfile"].filename(events_name)
        sim["caldb"].string(caldb)
        sim["irf"].string(irf)
        sim["ra"].real(ra)
        sim["dec"].real(dec)
        sim["rad"].real(rad_sim)
        sim["tmin"].real(tstart)
        sim["tmax"].real(tstop)
        sim["emin"].real(emin)
        sim["emax"].real(emax)
        self.test_try("Execute ctobssim")
        try:
            sim.execute()
            self.test_try_success()
        except:
            self.test_try_failure("Exception occured in ctobssim.")
    
        # Select events
        select = ctools.ctselect()
        select["infile"].filename(events_name)
        select["outfile"].filename(selected_events_name)
        select["ra"].real(ra)
        select["dec"].real(dec)
        select["rad"].real(rad_select)
        select["tmin"].real(tstart)
        select["tmax"].real(tstop)
        select["emin"].real(emin)
        select["emax"].real(emax)
        self.test_try("Execute ctselect")
        try:
            select.execute()
            self.test_try_success()
        except:
            self.test_try_failure("Exception occured in ctselect.")
    
        # Perform maximum likelihood fitting
        like = ctools.ctlike()
        like["infile"].filename(selected_events_name)
        like["srcmdl"].filename(model_name)
        like["outmdl"].filename(result_name)
        like["caldb"].string(caldb)
        like["irf"].string(irf)
        self.test_try("Execute ctlike")
        try:
            like.execute()
            self.test_try_success()
        except:
            self.test_try_failure("Exception occured in ctlike.")


    # Test unbinned in-memory pipeline
    def test_unbinned_mem(self):
        """
        Test unbinned in-memory pipeline.
        """
        # Set script parameters
        model_name           = "data/crab.xml"
        caldb                = "irf"
        irf                  = "cta_dummy_irf"
        ra                   =   83.63
        dec                  =   22.01
        rad_sim              =   10.0
        tstart               =    0.0
        tstop                = 1800.0
        emin                 =    0.1
        emax                 =  100.0
        rad_select           =    3.0
    
        # Simulate events
        sim = ctools.ctobssim()
        sim["infile"].filename(model_name)
        sim["caldb"].string(caldb)
        sim["irf"].string(irf)
        sim["ra"].real(ra)
        sim["dec"].real(dec)
        sim["rad"].real(rad_sim)
        sim["tmin"].real(tstart)
        sim["tmax"].real(tstop)
        sim["emin"].real(emin)
        sim["emax"].real(emax)
        self.test_try("Run ctobssim")
        try:
            sim.run()
            self.test_try_success()
        except:
            self.test_try_failure("Exception occured in ctobssim.")
    
        # Select events
        select = ctools.ctselect(sim.obs())
        select["ra"].real(ra)
        select["dec"].real(dec)
        select["rad"].real(rad_select)
        select["tmin"].real(tstart)
        select["tmax"].real(tstop)
        select["emin"].real(emin)
        select["emax"].real(emax)
        self.test_try("Run ctselect")
        try:
            select.run()
            self.test_try_success()
        except:
            self.test_try_failure("Exception occured in ctselect.")
    
        # Perform maximum likelihood fitting
        like = ctools.ctlike(select.obs())
        self.test_try("Run ctlike")
        try:
            like.run()
            self.test_try_success()
        except:
            self.test_try_failure("Exception occured in ctlike.")