Karim Sabri / ctools

Commit faaaab17a11849dc65be893a1dfeb0bdb9f0cad5

Authored by Jürgen Knödlseder
1 parent 3d1b77e4

Update initial scaling factors for Bremsstrahlung and Inverse Compton maps

Inline Side-by-side
Showing 1 changed file with 8 additions and 8 deletions
modules/comscripts/comobsmodel.py
  View file @faaaab1
@@ -334,10 +334,10 @@ class comobsmodel(ctools.csobservation): @@ -334,10 +334,10 @@ class comobsmodel(ctools.csobservation):
334 # Set spectral model using initial scaling factors that correspond 334 # Set spectral model using initial scaling factors that correspond
335 # to the results of Bloemen et al. (1999) 335 # to the results of Bloemen et al. (1999)
336 spectral = gammalib.GModelSpectralNodes() 336 spectral = gammalib.GModelSpectralNodes()
337 - spectral.append(gammalib.GEnergy(0.87, 'MeV'), 6.475e-3)  
338 - spectral.append(gammalib.GEnergy(1.73, 'MeV'), 1.403e-3)  
339 - spectral.append(gammalib.GEnergy(5.48, 'MeV'), 1.279e-4)  
340 - spectral.append(gammalib.GEnergy(17.32, 'MeV'), 1.552e-5) 337 + spectral.append(gammalib.GEnergy(0.87, 'MeV'), 0.93)
  338 + spectral.append(gammalib.GEnergy(1.73, 'MeV'), 0.20)
  339 + spectral.append(gammalib.GEnergy(5.48, 'MeV'), 0.018)
  340 + spectral.append(gammalib.GEnergy(17.32, 'MeV'), 0.0022)
341 341
342 # Recast model 342 # Recast model
343 spectral = self._recast_model(spectral, ebounds, -2.0, 1.0e-10, 1.0) 343 spectral = self._recast_model(spectral, ebounds, -2.0, 1.0e-10, 1.0)
@@ -393,10 +393,10 @@ class comobsmodel(ctools.csobservation): @@ -393,10 +393,10 @@ class comobsmodel(ctools.csobservation):
393 # Set spectral model using initial scaling factors that correspond 393 # Set spectral model using initial scaling factors that correspond
394 # to the results of Bloemen et al. (1999) 394 # to the results of Bloemen et al. (1999)
395 spectral = gammalib.GModelSpectralNodes() 395 spectral = gammalib.GModelSpectralNodes()
396 - spectral.append(gammalib.GEnergy(0.87, 'MeV'), 1.069e-2)  
397 - spectral.append(gammalib.GEnergy(1.73, 'MeV'), 3.043e-3)  
398 - spectral.append(gammalib.GEnergy(5.48, 'MeV'), 3.776e-4)  
399 - spectral.append(gammalib.GEnergy(17.32, 'MeV'), 5.175e-5) 396 + spectral.append(gammalib.GEnergy(0.87, 'MeV'), 2.23)
  397 + spectral.append(gammalib.GEnergy(1.73, 'MeV'), 0.64)
  398 + spectral.append(gammalib.GEnergy(5.48, 'MeV'), 0.079)
  399 + spectral.append(gammalib.GEnergy(17.32, 'MeV'), 0.011)
400 400
401 # Recast model 401 # Recast model
402 spectral = self._recast_model(spectral, ebounds, -2.0, 1.0e-10, 1.0) 402 spectral = self._recast_model(spectral, ebounds, -2.0, 1.0e-10, 1.0)