ctools

#! /usr/bin/env python

#import os

#import sys

#import math

#import ctools

import gammalib

# ================================ #

# Compare counts map and model map #

# ================================ #

def compare(cname, mname):

    """

    Compare counts map and model map

    """

    # Print filenames

    print cname, mname

    # Load maps

    counts = gammalib.GCTAEventCube()

    model  = gammalib.GCTAEventCube()

    counts.load(cname)

    model.load(mname)

    # Retrieve total number of counts

    total_counts = 0.0

    total_model  = 0.0

    for imap in range(counts.ebins()):

        map_counts = 0.0

        map_model  = 0.0

        for index in range(counts.map().npix()):

            map_counts += counts.map()[index, imap]

            map_model  += model.map()[index, imap]

        total_counts += map_counts

        total_model  += map_model

        print imap, map_counts, map_model, map_counts-map_model

    print "Total: ", total_counts, total_model, total_counts-total_model

    print

    # Return

    return

# ================ #

# Main entry point #

# ================ #

if __name__ == '__main__':

    """

    Compare counts map and model map

    """

    # Compare maps

    #compare("cmap.fits", "model.fits")

    compare("cntmap_1.fits",  "modmap_1.fits")

    compare("cntmap_20.fits", "modmap_20.fits")

    compare("cntmap_40.fits", "modmap_40.fits")