diff --git a/test/data/sciver/crab_smoothbplaw.xml b/test/data/sciver/crab_smoothbplaw.xml new file mode 100644 index 0000000..b502d0b --- /dev/null +++ b/test/data/sciver/crab_smoothbplaw.xml @@ -0,0 +1,24 @@ + + + + + + + + + + + + + + + + + + + + + + + + diff --git a/test/science_verification.py b/test/science_verification.py index 77c8faf..586896d 100755 --- a/test/science_verification.py +++ b/test/science_verification.py @@ -245,6 +245,7 @@ class sciver(gammalib.GPythonTestSuite): self.append(self.spec_plaw_edisp, 'Test power law model with energy dispersion') self.append(self.spec_plaw_stacked, 'Test power law model with stacked analysis') self.append(self.spec_plaw2, 'Test power law 2 model') + self.append(self.spec_smoothbplaw, 'Test smoothly broken power law 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') @@ -420,6 +421,20 @@ class sciver(gammalib.GPythonTestSuite): self.test('Pull_Background_Index') return + # Test smoothly broken power law model + def spec_smoothbplaw(self): + """ + Test smoothly broken power law model + """ + self.pull('data/sciver/crab_smoothbplaw') + self.test('Pull_Crab_Prefactor') + self.test('Pull_Crab_Index1') + self.test('Pull_Crab_Index2') + self.test('Pull_Crab_BreakEnergy') + self.test('Pull_Background_Prefactor') + self.test('Pull_Background_Index') + return + # Test exponentially cut off power law model def spec_eplaw(self): """