diff --git a/doc/source/users/tutorials/howto/comptel/howto_comptel_butterfly.png b/doc/source/users/tutorials/howto/comptel/howto_comptel_butterfly.png index 521cdae..c432723 100755 --- a/doc/source/users/tutorials/howto/comptel/howto_comptel_butterfly.png +++ b/doc/source/users/tutorials/howto/comptel/howto_comptel_butterfly.png diff --git a/doc/source/users/tutorials/howto/comptel/howto_comptel_dri.rst b/doc/source/users/tutorials/howto/comptel/howto_comptel_dri.rst index c2edcdd..be0f4e9 100644 --- a/doc/source/users/tutorials/howto/comptel/howto_comptel_dri.rst +++ b/doc/source/users/tutorials/howto/comptel/howto_comptel_dri.rst @@ -161,8 +161,9 @@ that should look as follows: not ``file`` since the ``IAQ`` parameter is not necessarily a file but can be also a response name of the calibration database. -Before doing a model fit you need a model. Create the following -:ref:`model definition file `: +Before doing a model fit you need a model. We included one in the ctools package +that you can find at ``$CTOOLS/share/models/comptel_howto_gc.xml``. Here the +first few lines of this model: .. code-block:: xml @@ -182,10 +183,6 @@ Before doing a model fit you need a model. Create the following - - - - ... @@ -212,15 +209,15 @@ map of the region around the Galactic centre as follows: $ cttsmap Input event list, counts cube or observation definition XML file [events.fits] obs_dri.xml Test source name [Crab] GC - Input model definition XML file [$CTOOLS/share/models/crab.xml] models.xml + Input model definition XML file [$CTOOLS/share/models/crab.xml] $CTOOLS/share/models/comptel_howto_gc.xml + Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL] GAL + Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR] First coordinate of image center in degrees (RA or galactic l) (0-360) [83.63] 0.0 Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] 0.0 - Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR] - Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL] GAL Image scale (in degrees/pixel) [0.02] 1.0 Size of the X axis in pixels [200] 50 Size of the Y axis in pixels [200] 30 - Output Test Statistic map file [tsmap.fits] tsmap.fits + Output Test Statistic map file [tsmap.fits] Below is the resulting Test Statistic map that shows 1.8 MeV emission following the Galactic plane and peaking near the Galactic centre. diff --git a/doc/source/users/tutorials/howto/comptel/howto_comptel_fitting.rst b/doc/source/users/tutorials/howto/comptel/howto_comptel_fitting.rst index 07d56b3..6292102 100755 --- a/doc/source/users/tutorials/howto/comptel/howto_comptel_fitting.rst +++ b/doc/source/users/tutorials/howto/comptel/howto_comptel_fitting.rst @@ -187,78 +187,78 @@ log file ``ctlike.log`` that was also created by :ref:`ctlike`: .. code-block:: none - 2017-08-26T22:05:05: +=================================+ - 2017-08-26T22:05:05: | Maximum likelihood optimisation | - 2017-08-26T22:05:05: +=================================+ - 2017-08-26T22:05:06: >Iteration 0: -logL=988086.091, Lambda=1.0e-03 - 2017-08-26T22:05:07: >Iteration 1: -logL=670534.792, Lambda=1.0e-03, delta=317551.299, step=1.0e+00, max(|grad|)=86219.859480 [Index:3] - 2017-08-26T22:05:08: >Iteration 2: -logL=387046.862, Lambda=1.0e-04, delta=283487.930, step=1.0e+00, max(|grad|)=69886.869623 [Index:3] - 2017-08-26T22:05:08: >Iteration 3: -logL=154426.362, Lambda=1.0e-05, delta=232620.499, step=1.0e+00, max(|grad|)=47994.333526 [Index:3] - 2017-08-26T22:05:09: >Iteration 4: -logL=-12393.885, Lambda=1.0e-06, delta=166820.247, step=1.0e+00, max(|grad|)=24763.811238 [Index:3] - 2017-08-26T22:05:10: >Iteration 5: -logL=-108271.655, Lambda=1.0e-07, delta=95877.770, step=1.0e+00, max(|grad|)=7544.090872 [Index:3] - 2017-08-26T22:05:11: >Iteration 6: -logL=-145906.957, Lambda=1.0e-08, delta=37635.302, step=1.0e+00, max(|grad|)=1167.483336 [Index:3] - 2017-08-26T22:05:12: >Iteration 7: -logL=-152898.212, Lambda=1.0e-09, delta=6991.255, step=1.0e+00, max(|grad|)=110.879761 [Index:3] - 2017-08-26T22:05:12: >Iteration 8: -logL=-153217.134, Lambda=1.0e-10, delta=318.921, step=1.0e+00, max(|grad|)=-2.969681 [Scale factor 2:111] - 2017-08-26T22:05:13: >Iteration 9: -logL=-153218.344, Lambda=1.0e-11, delta=1.210, step=1.0e+00, max(|grad|)=-0.047390 [Index:3] - 2017-08-26T22:05:14: >Iteration 10: -logL=-153218.344, Lambda=1.0e-12, delta=0.000, step=1.0e+00, max(|grad|)=-0.001891 [Index:3] - 2017-08-26T22:05:15: - 2017-08-26T22:05:15: +====================================+ - 2017-08-26T22:05:15: | Maximum likelihood re-optimisation | - 2017-08-26T22:05:15: +====================================+ - 2017-08-26T22:05:15: >Iteration 0: -logL=-152378.830, Lambda=1.0e-03 - 2017-08-26T22:05:16: >Iteration 1: -logL=-152759.255, Lambda=1.0e-03, delta=380.425, step=1.0e+00, max(|grad|)=-2.480864 [Scale factor 10:21] - 2017-08-26T22:05:16: >Iteration 2: -logL=-152760.518, Lambda=1.0e-04, delta=1.264, step=1.0e+00, max(|grad|)=-0.019031 [Scale factor 10:21] - 2017-08-26T22:05:16: >Iteration 3: -logL=-152760.518, Lambda=1.0e-05, delta=0.000, step=1.0e+00, max(|grad|)=-0.000001 [Scale factor 10:21] - 2017-08-26T22:05:16: - 2017-08-26T22:05:16: +============================================+ - 2017-08-26T22:05:16: | Maximum likelihood re-optimisation results | - 2017-08-26T22:05:16: +============================================+ - 2017-08-26T22:05:16: === GOptimizerLM === - 2017-08-26T22:05:16: Optimized function value ..: -152760.518 - 2017-08-26T22:05:16: Absolute precision ........: 0.005 - 2017-08-26T22:05:16: Acceptable value decrease .: 2 - 2017-08-26T22:05:16: Optimization status .......: converged - 2017-08-26T22:05:16: Number of parameters ......: 200 - 2017-08-26T22:05:16: Number of free parameters .: 87 - 2017-08-26T22:05:16: Number of iterations ......: 3 - 2017-08-26T22:05:16: Lambda ....................: 1e-06 - 2017-08-26T22:05:16: - 2017-08-26T22:05:16: +=========================================+ - 2017-08-26T22:05:16: | Maximum likelihood optimisation results | - 2017-08-26T22:05:16: +=========================================+ - 2017-08-26T22:05:16: === GOptimizerLM === - 2017-08-26T22:05:16: Optimized function value ..: -153218.344 - 2017-08-26T22:05:16: Absolute precision ........: 0.005 - 2017-08-26T22:05:16: Acceptable value decrease .: 2 - 2017-08-26T22:05:16: Optimization status .......: converged - 2017-08-26T22:05:16: Number of parameters ......: 206 - 2017-08-26T22:05:16: Number of free parameters .: 89 - 2017-08-26T22:05:16: Number of iterations ......: 10 - 2017-08-26T22:05:16: Lambda ....................: 1e-13 - 2017-08-26T22:05:16: Maximum log likelihood ....: 153218.344 - 2017-08-26T22:05:16: Observed events (Nobs) ...: 527595.000 - 2017-08-26T22:05:16: Predicted events (Npred) ..: 527594.000 (Nobs - Npred = 1.00004221417475) - 2017-08-26T22:05:16: === GModels === - 2017-08-26T22:05:16: Number of models ..........: 5 - 2017-08-26T22:05:16: Number of parameters ......: 206 - 2017-08-26T22:05:16: === GModelSky === - 2017-08-26T22:05:16: Name ......................: Crab - 2017-08-26T22:05:16: Instruments ...............: all - 2017-08-26T22:05:16: Test Statistic ............: 915.650908945536 - 2017-08-26T22:05:16: Instrument scale factors ..: unity - 2017-08-26T22:05:16: Observation identifiers ...: all - 2017-08-26T22:05:16: Model type ................: PointSource - 2017-08-26T22:05:16: Model components ..........: "PointSource" * "PowerLaw" * "Constant" - 2017-08-26T22:05:16: Number of parameters ......: 6 - 2017-08-26T22:05:16: Number of spatial par's ...: 2 - 2017-08-26T22:05:16: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) - 2017-08-26T22:05:16: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) - 2017-08-26T22:05:16: Number of spectral par's ..: 3 - 2017-08-26T22:05:16: Prefactor ................: 0.00237517119227384 +/- 0.000107673211964802 [-1,1] ph/cm2/s/MeV (free,scale=0.001,gradient) - 2017-08-26T22:05:16: Index ....................: -2.47153532398462 +/- 0.0500838601327289 [5,-5] (free,scale=-1,gradient) - 2017-08-26T22:05:16: PivotEnergy ..............: 1 [0.01,1000] MeV (fixed,scale=1,gradient) - 2017-08-26T22:05:16: Number of temporal par's ..: 1 - 2017-08-26T22:05:16: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) + 2019-04-04T15:11:31: +=================================+ + 2019-04-04T15:11:31: | Maximum likelihood optimisation | + 2019-04-04T15:11:31: +=================================+ + 2019-04-04T15:11:32: >Iteration 0: -logL=988079.891, Lambda=1.0e-03 + 2019-04-04T15:11:32: >Iteration 1: -logL=670528.685, Lambda=1.0e-03, delta=317551.206, step=1.0e+00, max(|grad|)=86223.268050 [Index:3] + 2019-04-04T15:11:33: >Iteration 2: -logL=387041.036, Lambda=1.0e-04, delta=283487.648, step=1.0e+00, max(|grad|)=69889.358952 [Index:3] + 2019-04-04T15:11:34: >Iteration 3: -logL=154421.136, Lambda=1.0e-05, delta=232619.900, step=1.0e+00, max(|grad|)=47995.573198 [Index:3] + 2019-04-04T15:11:35: >Iteration 4: -logL=-12398.051, Lambda=1.0e-06, delta=166819.188, step=1.0e+00, max(|grad|)=24764.034734 [Index:3] + 2019-04-04T15:11:35: >Iteration 5: -logL=-108274.155, Lambda=1.0e-07, delta=95876.104, step=1.0e+00, max(|grad|)=7544.055951 [Index:3] + 2019-04-04T15:11:36: >Iteration 6: -logL=-145907.646, Lambda=1.0e-08, delta=37633.490, step=1.0e+00, max(|grad|)=1167.458562 [Index:3] + 2019-04-04T15:11:36: >Iteration 7: -logL=-152898.181, Lambda=1.0e-09, delta=6990.535, step=1.0e+00, max(|grad|)=110.883829 [Index:3] + 2019-04-04T15:11:37: >Iteration 8: -logL=-153217.048, Lambda=1.0e-10, delta=318.867, step=1.0e+00, max(|grad|)=-2.969678 [Scale factor 2:111] + 2019-04-04T15:11:37: >Iteration 9: -logL=-153218.258, Lambda=1.0e-11, delta=1.210, step=1.0e+00, max(|grad|)=-0.047105 [Index:3] + 2019-04-04T15:11:38: >Iteration 10: -logL=-153218.258, Lambda=1.0e-12, delta=0.000, step=1.0e+00, max(|grad|)=-0.001877 [Index:3] + 2019-04-04T15:11:38: + 2019-04-04T15:11:38: +====================================+ + 2019-04-04T15:11:38: | Maximum likelihood re-optimisation | + 2019-04-04T15:11:38: +====================================+ + 2019-04-04T15:11:39: >Iteration 0: -logL=-152378.895, Lambda=1.0e-03 + 2019-04-04T15:11:39: >Iteration 1: -logL=-152759.255, Lambda=1.0e-03, delta=380.360, step=1.0e+00, max(|grad|)=-2.480233 [Scale factor 10:21] + 2019-04-04T15:11:39: >Iteration 2: -logL=-152760.518, Lambda=1.0e-04, delta=1.263, step=1.0e+00, max(|grad|)=-0.019021 [Scale factor 10:21] + 2019-04-04T15:11:40: >Iteration 3: -logL=-152760.518, Lambda=1.0e-05, delta=0.000, step=1.0e+00, max(|grad|)=-0.000001 [Scale factor 10:21] + 2019-04-04T15:11:40: + 2019-04-04T15:11:40: +============================================+ + 2019-04-04T15:11:40: | Maximum likelihood re-optimisation results | + 2019-04-04T15:11:40: +============================================+ + 2019-04-04T15:11:40: === GOptimizerLM === + 2019-04-04T15:11:40: Optimized function value ..: -152760.518 + 2019-04-04T15:11:40: Absolute precision ........: 0.005 + 2019-04-04T15:11:40: Acceptable value decrease .: 2 + 2019-04-04T15:11:40: Optimization status .......: converged + 2019-04-04T15:11:40: Number of parameters ......: 200 + 2019-04-04T15:11:40: Number of free parameters .: 87 + 2019-04-04T15:11:40: Number of iterations ......: 3 + 2019-04-04T15:11:40: Lambda ....................: 1e-06 + 2019-04-04T15:11:40: + 2019-04-04T15:11:40: +=========================================+ + 2019-04-04T15:11:40: | Maximum likelihood optimisation results | + 2019-04-04T15:11:40: +=========================================+ + 2019-04-04T15:11:40: === GOptimizerLM === + 2019-04-04T15:11:40: Optimized function value ..: -153218.258 + 2019-04-04T15:11:40: Absolute precision ........: 0.005 + 2019-04-04T15:11:40: Acceptable value decrease .: 2 + 2019-04-04T15:11:40: Optimization status .......: converged + 2019-04-04T15:11:40: Number of parameters ......: 206 + 2019-04-04T15:11:40: Number of free parameters .: 89 + 2019-04-04T15:11:40: Number of iterations ......: 10 + 2019-04-04T15:11:40: Lambda ....................: 1e-13 + 2019-04-04T15:11:40: Maximum log likelihood ....: 153218.258 + 2019-04-04T15:11:40: Observed events (Nobs) ...: 527595.000 + 2019-04-04T15:11:40: Predicted events (Npred) ..: 527594.000 (Nobs - Npred = 1.00004248635378) + 2019-04-04T15:11:40: === GModels === + 2019-04-04T15:11:40: Number of models ..........: 5 + 2019-04-04T15:11:40: Number of parameters ......: 206 + 2019-04-04T15:11:40: === GModelSky === + 2019-04-04T15:11:40: Name ......................: Crab + 2019-04-04T15:11:40: Instruments ...............: all + 2019-04-04T15:11:40: Test Statistic ............: 915.478960024426 + 2019-04-04T15:11:40: Instrument scale factors ..: unity + 2019-04-04T15:11:40: Observation identifiers ...: all + 2019-04-04T15:11:40: Model type ................: PointSource + 2019-04-04T15:11:40: Model components ..........: "PointSource" * "PowerLaw" * "Constant" + 2019-04-04T15:11:40: Number of parameters ......: 6 + 2019-04-04T15:11:40: Number of spatial par's ...: 2 + 2019-04-04T15:11:40: RA .......................: 83.6331 [-360,360] deg (fixed,scale=1) + 2019-04-04T15:11:40: DEC ......................: 22.0145 [-90,90] deg (fixed,scale=1) + 2019-04-04T15:11:40: Number of spectral par's ..: 3 + 2019-04-04T15:11:40: Prefactor ................: 0.00237485707890478 +/- 0.000107670251181817 [-1,1] ph/cm2/s/MeV (free,scale=0.001,gradient) + 2019-04-04T15:11:40: Index ....................: -2.47146337740947 +/- 0.0500850436160353 [5,-5] (free,scale=-1,gradient) + 2019-04-04T15:11:40: PivotEnergy ..............: 1 [0.01,1000] MeV (fixed,scale=1,gradient) + 2019-04-04T15:11:40: Number of temporal par's ..: 1 + 2019-04-04T15:11:40: Normalization ............: 1 (relative value) (fixed,scale=1,gradient) The fit converged after 10 iterations and did a second fit without the source component to determine the Test Statistic value of the Crab. diff --git a/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.png b/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.png index 7fdbd3b..e5c4dc3 100755 --- a/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.png +++ b/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.png diff --git a/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.rst b/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.rst index 04cf1ad..c3ecb2f 100755 --- a/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.rst +++ b/doc/source/users/tutorials/howto/comptel/howto_comptel_spectrum.rst @@ -8,65 +8,30 @@ Generate a spectral energy distribution for a source You will learn how to use :ref:`ctlike` to derive the spectral energy density of a source from COMPTEL data. -So far the :ref:`csspec` script does not work for COMPTEL data, but you -will learn here a method that you can use to derive the spectrum of any -source using COMPTEL data. For this you have to modify the -:ref:`model definition file ` -and replace the spectral component of the Crab by a node function: - -.. code-block:: xml - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - ... - - -Instead of the power law the model now contains a node function with four -intensity nodes between 0.75 MeV and 30 MeV, and the :ref:`ctlike` tool will -now adjust the intensities of each of the nodes which will result in a -spectrum. To do this type now: +To determine a spectrum of the Crab nebula you can use the :ref:`csspec` +script. To do this, type .. code-block:: bash - $ ctlike - Input event list, counts cube or observation definition XML file [obs.xml] - Input model definition XML file [models.xml] model_nodes.xml - Output model definition XML file [crab_results.xml] crab_nodes_results.xml - -This will produce the -:ref:`model definition file ` -``crab_nodes_results.xml`` with the intensities of the fitted source at -each energy. The results are displayed graphically below. The red data points -show the combined Crab nebular and pulsar spectrum from + $ csspec + Input event list, counts cube, or observation definition XML file [events.fits] obs.xml + Input model definition XML file [$CTOOLS/share/models/crab.xml] models.xml + Source name [Crab] + Spectrum generation method (SLICE|NODES|AUTO) [AUTO] + Algorithm for defining energy bins (FILE|LIN|LOG) [LOG] + Start value for first energy bin in TeV [0.1] 0.00000075 + Stop value for last energy bin in TeV [100.0] 0.00003 + Number of energy bins (1-200) [20] 4 + Output spectrum file [spectrum.fits] + +This will generate a logarithmically spaced spectrum composed of 4 energy bins +comprised within 0.75 MeV and 30 GeV. The :ref:`csspec` tool is run in the ``AUTO`` +mode, which for COMPTEL data corresponds to the ``NODES`` method. :ref:`csspec` +will convert the Crab spectrum in the input model into a node function and fit +all intensity nodes using :ref:`ctlike`. The resulting node intensities are +then interpreted as the spectral point, and are displayed as red data points +in the figure below. The black data points show the combined Crab nebular and +pulsar spectrum from `Kuiper et al. 2001, A&A, 378, 918 `_ for comparison. @@ -75,7 +40,3 @@ for comparison. :align: center *Spectrum of the Crab derived using a node function* - -.. note:: - The ``NodeFunction`` spectral models implements a piece-wise power law - that is defined by the intensity values at each node. diff --git a/doc/source/users/tutorials/howto/comptel/howto_comptel_tsmap.rst b/doc/source/users/tutorials/howto/comptel/howto_comptel_tsmap.rst index 5fac357..1c2ddb6 100755 --- a/doc/source/users/tutorials/howto/comptel/howto_comptel_tsmap.rst +++ b/doc/source/users/tutorials/howto/comptel/howto_comptel_tsmap.rst @@ -20,14 +20,14 @@ source. You do this by using the :ref:`cttsmap` tool as follows: Input event list, counts cube or observation definition XML file [events.fits] obs.xml Test source name [Crab] Input model definition XML file [$CTOOLS/share/models/crab.xml] models.xml + Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL] + Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR] First coordinate of image center in degrees (RA or galactic l) (0-360) [83.63] Second coordinate of image center in degrees (DEC or galactic b) (-90-90) [22.01] - Projection method (AIT|AZP|CAR|GLS|MER|MOL|SFL|SIN|STG|TAN) [CAR] - Coordinate system (CEL - celestial, GAL - galactic) (CEL|GAL) [CEL] Image scale (in degrees/pixel) [0.02] 1.0 Size of the X axis in pixels [200] 10 Size of the Y axis in pixels [200] 10 - Output Test Statistic map file [tsmap.fits] + Output Test Statistic map file [tsmap.fits] :ref:`cttsmap` created a 10 x 10 pixels map with a pixel size of 1 degree around the location of the Crab and saved the result in the FITS file diff --git a/doc/source/users/tutorials/howto/comptel/howto_comptel_ulimit.rst b/doc/source/users/tutorials/howto/comptel/howto_comptel_ulimit.rst index 2e8e5ab..766efc3 100755 --- a/doc/source/users/tutorials/howto/comptel/howto_comptel_ulimit.rst +++ b/doc/source/users/tutorials/howto/comptel/howto_comptel_ulimit.rst @@ -29,28 +29,28 @@ create the following - + - + - + - + @@ -98,9 +98,9 @@ The results of the upper limit computation can then be extracted from the .. code-block:: bash - 2017-08-27T08:13:44: +=====================+ - 2017-08-27T08:13:44: | Upper limit results | - 2017-08-27T08:13:44: +=====================+ - 2017-08-27T08:13:44: Differential flux limit ...: 9.26955448510129e-06 ph/cm2/s/MeV at 1e-06 TeV - 2017-08-27T08:13:44: Integral flux limit .......: 6.43104918005e-05 ph/cm2/s within [7.5e-08-3e-05] TeV - 2017-08-27T08:13:44: Energy flux limit .........: 1.54555137999306e-10 erg/cm2/s within [7.5e-08-3e-05] TeV + 2019-04-04T15:47:22: +=====================+ + 2019-04-04T15:47:22: | Upper limit results | + 2019-04-04T15:47:22: +=====================+ + 2019-04-04T15:47:22: Differential flux limit ...: 9.26467937144868e-06 ph/cm2/s/MeV at 1e-06 TeV + 2019-04-04T15:47:22: Integral flux limit .......: 6.42766691440729e-05 ph/cm2/s within [7.5e-08-3e-05] TeV + 2019-04-04T15:47:22: Energy flux limit .........: 1.54473853201363e-10 erg/cm2/s within [7.5e-08-3e-05] TeV diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png b/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png old mode 100755 new mode 100644 index fa5416a..80b8ef0 --- a/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png +++ b/doc/source/users/tutorials/howto/fermi/howto_fermi_butterfly.png diff --git a/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.rst b/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.rst index 9b614d3..9be47d1 100755 --- a/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.rst +++ b/doc/source/users/tutorials/howto/fermi/howto_fermi_spectrum.rst @@ -8,7 +8,7 @@ Generate a spectral energy distribution for a source You will learn how to use :ref:`csspec` to derive the spectral energy density of a source from Fermi-LAT data. -To determine a spectral of the Vela pulsar you can use the :ref:`csspec` +To determine a spectrum of the Vela pulsar you can use the :ref:`csspec` script. To do this, type .. code-block:: bash diff --git a/models/Makefile.am b/models/Makefile.am index 6294930..859f199 100755 --- a/models/Makefile.am +++ b/models/Makefile.am @@ -1,7 +1,7 @@ ############################################################################# # Makefile for models # # ------------------------------------------------------------------------- # -# Copyright (C) 2016-2018 Juergen Knoedlseder # +# Copyright (C) 2016-2019 Juergen Knoedlseder # # ------------------------------------------------------------------------- # # # # This program is free software: you can redistribute it and/or modify # @@ -24,6 +24,7 @@ models = 1dc_howto.xml \ 1dc_howto_unbinned.xml \ 1dc_howto_vela.xml \ 1dc_howto_velapsr.xml \ + comptel_howto_gc.xml \ bkg_irf.xml \ bkg_cube.xml \ test_irf.xml \ @@ -31,7 +32,7 @@ models = 1dc_howto.xml \ crab.xml \ crab_beyond.xml \ crab_2bkg.xml \ - crab_nobkg.xml \ + crab_nobkg.xml \ crab_poly.xml \ crab_profile.xml \ disk.xml \ diff --git a/models/comptel_howto_gc.xml b/models/comptel_howto_gc.xml new file mode 100644 index 0000000..4219193 --- /dev/null +++ b/models/comptel_howto_gc.xml @@ -0,0 +1,98 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +