From abba998579f7da419eaaf73f42abf734f982debf Mon Sep 17 00:00:00 2001 From: Jürgen Knödlseder Date: Fri, 16 Oct 2015 09:33:40 +0200 Subject: [PATCH] Add spatial models to science verification --- test/data/sciver/cena_lobes_parkes.fits | Bin 0 -> 267840 bytes test/data/sciver/crab_edisk.xml | 24 ++++++++++++++++++++++++ test/data/sciver/crab_egauss.xml | 24 ++++++++++++++++++++++++ test/data/sciver/crab_map.xml | 20 ++++++++++++++++++++ test/data/sciver/crab_ptsrc.xml | 21 +++++++++++++++++++++ test/data/sciver/crab_rdisk.xml | 22 ++++++++++++++++++++++ test/data/sciver/crab_rgauss.xml | 22 ++++++++++++++++++++++ test/data/sciver/crab_rshell.xml | 23 +++++++++++++++++++++++ test/science_verification.py | 123 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++---- 9 files changed, 275 insertions(+), 4 deletions(-) create mode 100644 test/data/sciver/cena_lobes_parkes.fits create mode 100644 test/data/sciver/crab_edisk.xml create mode 100644 test/data/sciver/crab_egauss.xml create mode 100644 test/data/sciver/crab_map.xml create mode 100644 test/data/sciver/crab_ptsrc.xml create mode 100644 test/data/sciver/crab_rdisk.xml create mode 100644 test/data/sciver/crab_rgauss.xml create mode 100644 test/data/sciver/crab_rshell.xml diff --git a/test/data/sciver/cena_lobes_parkes.fits b/test/data/sciver/cena_lobes_parkes.fits new file mode 100644 index 0000000..47d4721 Binary files /dev/null and b/test/data/sciver/cena_lobes_parkes.fits differ diff --git a/test/data/sciver/crab_edisk.xml b/test/data/sciver/crab_edisk.xml new file mode 100644 index 0000000..7bad7ca --- /dev/null +++ b/test/data/sciver/crab_edisk.xml @@ -0,0 +1,24 @@ + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/test/data/sciver/crab_egauss.xml b/test/data/sciver/crab_egauss.xml new file mode 100644 index 0000000..ee291d4 --- /dev/null +++ b/test/data/sciver/crab_egauss.xml @@ -0,0 +1,24 @@ + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/test/data/sciver/crab_map.xml b/test/data/sciver/crab_map.xml new file mode 100644 index 0000000..aff473c --- /dev/null +++ b/test/data/sciver/crab_map.xml @@ -0,0 +1,20 @@ + + + + + + + + + + + + + + + + + + + + diff --git a/test/data/sciver/crab_ptsrc.xml b/test/data/sciver/crab_ptsrc.xml new file mode 100644 index 0000000..7c4b378 --- /dev/null +++ b/test/data/sciver/crab_ptsrc.xml @@ -0,0 +1,21 @@ + + + + + + + + + + + + + + + + + + + + + diff --git a/test/data/sciver/crab_rdisk.xml b/test/data/sciver/crab_rdisk.xml new file mode 100644 index 0000000..33e66bd --- /dev/null +++ b/test/data/sciver/crab_rdisk.xml @@ -0,0 +1,22 @@ + + + + + + + + + + + + + + + + + + + + + + diff --git a/test/data/sciver/crab_rgauss.xml b/test/data/sciver/crab_rgauss.xml new file mode 100644 index 0000000..a78de6d --- /dev/null +++ b/test/data/sciver/crab_rgauss.xml @@ -0,0 +1,22 @@ + + + + + + + + + + + + + + + + + + + + + + diff --git a/test/data/sciver/crab_rshell.xml b/test/data/sciver/crab_rshell.xml new file mode 100644 index 0000000..3231819 --- /dev/null +++ b/test/data/sciver/crab_rshell.xml @@ -0,0 +1,23 @@ + + + + + + + + + + + + + + + + + + + + + + + diff --git a/test/science_verification.py b/test/science_verification.py index ac77928..e550b89 100755 --- a/test/science_verification.py +++ b/test/science_verification.py @@ -181,13 +181,22 @@ class sciver(gammalib.GPythonTestSuite): # Append spectral tests #self.append(self.spec_plaw, "Test power law model") - self.append(self.spec_plaw2, "Test power law 2 model") - self.append(self.spec_eplaw, "Test exponentially cut off power law model") + #self.append(self.spec_plaw2, "Test power law 2 model") + #self.append(self.spec_eplaw, "Test exponentially cut off power law model") #self.append(self.spec_supeplaw, "Test super exponentially cut off power law model") - self.append(self.spec_logparabola, "Test log parabola model") + #self.append(self.spec_logparabola, "Test log parabola model") #self.append(self.spec_gauss, "Test Gaussian model") #self.append(self.spec_filefct, "Test file function model") - self.append(self.spec_nodes, "Test nodes model") + #self.append(self.spec_nodes, "Test nodes model") + + # Append spatial tests + self.append(self.spec_ptsrc, "Test point source model") + self.append(self.spec_rdisk, "Test radial disk model") + self.append(self.spec_rgauss, "Test radial Gaussian model") + self.append(self.spec_rshell, "Test radial shell model") + self.append(self.spec_edisk, "Test elliptical disk model") + self.append(self.spec_egauss, "Test elliptical Gaussian model") + self.append(self.spec_map, "Test elliptical Gaussian model") # Return return @@ -329,6 +338,112 @@ class sciver(gammalib.GPythonTestSuite): self.test("Pull_Background_Index") return + # Test point source model + def spec_ptsrc(self): + """ + Test point source model. + """ + self.pull("data/sciver/crab_ptsrc") + self.test("Pull_Crab_Prefactor") + self.test("Pull_Crab_Index") + self.test("Pull_Crab_RA") + self.test("Pull_Crab_DEC") + self.test("Pull_Background_Prefactor") + self.test("Pull_Background_Index") + return + + # Test radial disk model + def spec_rdisk(self): + """ + Test radial disk model. + """ + self.pull("data/sciver/crab_rdisk") + self.test("Pull_Crab_Prefactor") + self.test("Pull_Crab_Index") + self.test("Pull_Crab_RA") + self.test("Pull_Crab_DEC") + self.test("Pull_Crab_Radius") + self.test("Pull_Background_Prefactor") + self.test("Pull_Background_Index") + return + + # Test radial Gaussian model + def spec_rgauss(self): + """ + Test radial Gaussian model. + """ + self.pull("data/sciver/crab_rgauss") + self.test("Pull_Crab_Prefactor") + self.test("Pull_Crab_Index") + self.test("Pull_Crab_RA") + self.test("Pull_Crab_DEC") + self.test("Pull_Crab_Sigma") + self.test("Pull_Background_Prefactor") + self.test("Pull_Background_Index") + return + + # Test radial shell model + def spec_rshell(self): + """ + Test radial shell model. + """ + self.pull("data/sciver/crab_rshell") + self.test("Pull_Crab_Prefactor") + self.test("Pull_Crab_Index") + self.test("Pull_Crab_RA") + self.test("Pull_Crab_DEC") + self.test("Pull_Crab_Radius") + self.test("Pull_Crab_Width") + self.test("Pull_Background_Prefactor") + self.test("Pull_Background_Index") + return + + # Test elliptical disk model + def spec_edisk(self): + """ + Test elliptical disk model. + """ + self.pull("data/sciver/crab_edisk") + self.test("Pull_Crab_Prefactor") + self.test("Pull_Crab_Index") + self.test("Pull_Crab_RA") + self.test("Pull_Crab_DEC") + self.test("Pull_Crab_PA") + self.test("Pull_Crab_MinorRadius") + self.test("Pull_Crab_MajorRadius") + self.test("Pull_Background_Prefactor") + self.test("Pull_Background_Index") + return + + # Test elliptical Gaussian model + def spec_egauss(self): + """ + Test elliptical Gaussian model. + """ + self.pull("data/sciver/crab_egauss") + self.test("Pull_Crab_Prefactor") + self.test("Pull_Crab_Index") + self.test("Pull_Crab_RA") + self.test("Pull_Crab_DEC") + self.test("Pull_Crab_PA") + self.test("Pull_Crab_MinorRadius") + self.test("Pull_Crab_MajorRadius") + self.test("Pull_Background_Prefactor") + self.test("Pull_Background_Index") + return + + # Test diffuse map model + def spec_map(self): + """ + Test diffuse map model. + """ + self.pull("data/sciver/crab_map") + self.test("Pull_Crab_Prefactor") + self.test("Pull_Crab_Index") + self.test("Pull_Background_Prefactor") + self.test("Pull_Background_Index") + return + # ======================== # # Main routine entry point # -- libgit2 0.21.4