diff --git a/doc/source/users/notebooks.rst b/doc/source/users/notebooks.rst index f932a44..0b62da1 100644 --- a/doc/source/users/notebooks.rst +++ b/doc/source/users/notebooks.rst @@ -1,7 +1,7 @@ .. _sec_notebooks: Jupyter notebooks -============ +================= Some of the tutorials are provided as Jupyter notebooks. If you see this icon diff --git a/doc/source/users/tutorials/quickstart/butterfly.rst b/doc/source/users/tutorials/quickstart/butterfly.rst index 3a5a50a..35512d3 100644 --- a/doc/source/users/tutorials/quickstart/butterfly.rst +++ b/doc/source/users/tutorials/quickstart/butterfly.rst @@ -38,10 +38,10 @@ default is set to 100. Below an excerpt of the first lines of the .. code-block:: none - 103514.216667934 8.08876441182872e-15 7.86364405017512e-15 8.31388477348232e-15 - 110917.48152624 6.81370736945431e-15 6.62965024227895e-15 6.99776449662967e-15 - 118850.222743702 5.73964152654309e-15 5.58918876927557e-15 5.89009428381061e-15 - 127350.308101666 4.83488401643174e-15 4.71191629605204e-15 4.95785173681143e-15 + 103514.216667934 8.08876441182872e-15 7.86364405017512e-15 8.31388477348232e-15 + 110917.48152624 6.81370736945431e-15 6.62965024227895e-15 6.99776449662967e-15 + 118850.222743702 5.73964152654309e-15 5.58918876927557e-15 5.89009428381061e-15 + 127350.308101666 4.83488401643174e-15 4.71191629605204e-15 4.95785173681143e-15 ... diff --git a/doc/source/users/tutorials/quickstart/fitting.rst b/doc/source/users/tutorials/quickstart/fitting.rst index 4db4aea..3e926de 100644 --- a/doc/source/users/tutorials/quickstart/fitting.rst +++ b/doc/source/users/tutorials/quickstart/fitting.rst @@ -49,27 +49,27 @@ that was produced by the run: .. code-block:: xml - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + In this example, the ``Prefactor`` and ``Index`` of the spectral model for the Crab nebula as well as the ``Prefactor`` and ``Index`` of the background spectral diff --git a/doc/source/users/tutorials/quickstart/selecting.rst b/doc/source/users/tutorials/quickstart/selecting.rst index bffc9b8..699f321 100644 --- a/doc/source/users/tutorials/quickstart/selecting.rst +++ b/doc/source/users/tutorials/quickstart/selecting.rst @@ -35,43 +35,41 @@ event selection: .. code-block:: none - 2017-11-28T14:21:57: +===================+ - 2017-11-28T14:21:57: | Input observation | - 2017-11-28T14:21:57: +===================+ - 2017-11-28T14:21:57: === GObservations === - 2017-11-28T14:21:57: Number of observations ....: 1 - 2017-11-28T14:21:57: Number of models ..........: 0 - 2017-11-28T14:21:57: Number of observed events .: 200833 - 2017-11-28T14:21:57: Number of predicted events : 0 - 2017-11-28T14:21:57: - 2017-11-28T14:21:57: +=================+ - 2017-11-28T14:21:57: | Event selection | - 2017-11-28T14:21:57: +=================+ - 2017-11-28T14:21:57: === CTA observation === - 2017-11-28T14:21:57: Input filename ............: events.fits - 2017-11-28T14:21:57: Event extension name ......: EVENTS - 2017-11-28T14:21:57: GTI extension name ........: GTI - 2017-11-28T14:21:58: === Event selection === - 2017-11-28T14:21:58: Time range (MJD) ..........: 58849.0077451852 - 58849.0285785185 days - 2017-11-28T14:21:58: Time range (UTC) ..........: 2020-01-01T00:10:00 - 2020-01-01T00:40:00 - 2017-11-28T14:21:58: Time range (MET) ..........: 631109469.184 - 631111269.184 seconds - 2017-11-28T14:21:58: Selected energy range .....: 0.1 - 100 TeV - 2017-11-28T14:21:58: Requested RoI .............: Centre(RA,DEC)=(83.63, 22.01) deg, Radius=3 deg - 2017-11-28T14:21:58: RoI of data ...............: Centre(RA,DEC)=(83.5, 22.8) deg, Radius=5 deg - 2017-11-28T14:21:58: Selected RoI ..............: Centre(RA,DEC)=(83.63, 22.01) deg, Radius=3 deg - 2017-11-28T14:21:58: cfitsio selection .........: TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000 - 2017-11-28T14:21:58: FITS filename .............: /var/tmp/tmp.0.Nc6oWk[EVENTS][TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000] - 2017-11-28T14:21:58: - 2017-11-28T14:21:58: +====================+ - 2017-11-28T14:21:58: | Output observation | - 2017-11-28T14:21:58: +====================+ - 2017-11-28T14:21:58: === GObservations === - 2017-11-28T14:21:58: Number of observations ....: 1 - 2017-11-28T14:21:58: Number of models ..........: 0 - 2017-11-28T14:21:58: Number of observed events .: 20877 - 2017-11-28T14:21:58: Number of predicted events : 0 - 2017-11-28T14:21:58: - 2017-11-28T14:21:58: +=================+ + 2017-11-28T14:21:57: +===================+ + 2017-11-28T14:21:57: | Input observation | + 2017-11-28T14:21:57: +===================+ + 2017-11-28T14:21:57: === GObservations === + 2017-11-28T14:21:57: Number of observations ....: 1 + 2017-11-28T14:21:57: Number of models ..........: 0 + 2017-11-28T14:21:57: Number of observed events .: 200833 + 2017-11-28T14:21:57: Number of predicted events : 0 + 2017-11-28T14:21:57: + 2017-11-28T14:21:57: +=================+ + 2017-11-28T14:21:57: | Event selection | + 2017-11-28T14:21:57: +=================+ + 2017-11-28T14:21:57: === CTA observation === + 2017-11-28T14:21:57: Input filename ............: events.fits + 2017-11-28T14:21:57: Event extension name ......: EVENTS + 2017-11-28T14:21:57: GTI extension name ........: GTI + 2017-11-28T14:21:58: === Event selection === + 2017-11-28T14:21:58: Time range (MJD) ..........: 58849.0077451852 - 58849.0285785185 days + 2017-11-28T14:21:58: Time range (UTC) ..........: 2020-01-01T00:10:00 - 2020-01-01T00:40:00 + 2017-11-28T14:21:58: Time range (MET) ..........: 631109469.184 - 631111269.184 seconds + 2017-11-28T14:21:58: Selected energy range .....: 0.1 - 100 TeV + 2017-11-28T14:21:58: Requested RoI .............: Centre(RA,DEC)=(83.63, 22.01) deg, Radius=3 deg + 2017-11-28T14:21:58: RoI of data ...............: Centre(RA,DEC)=(83.5, 22.8) deg, Radius=5 deg + 2017-11-28T14:21:58: Selected RoI ..............: Centre(RA,DEC)=(83.63, 22.01) deg, Radius=3 deg + 2017-11-28T14:21:58: cfitsio selection .........: TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000 + 2017-11-28T14:21:58: FITS filename .............: /var/tmp/tmp.0.Nc6oWk[EVENTS][TIME >= 631109469.18399978 && TIME <= 631111269.18400002 && ENERGY >= 0.10000000 && ENERGY <= 100.00000000 && ANGSEP(83.630000,22.010000,RA,DEC) <= 3.000000] + 2017-11-28T14:21:58: + 2017-11-28T14:21:58: +====================+ + 2017-11-28T14:21:58: | Output observation | + 2017-11-28T14:21:58: +====================+ + 2017-11-28T14:21:58: === GObservations === + 2017-11-28T14:21:58: Number of observations ....: 1 + 2017-11-28T14:21:58: Number of models ..........: 0 + 2017-11-28T14:21:58: Number of observed events .: 20877 + 2017-11-28T14:21:58: Number of predicted events : 0 From the 200833 events in the input event list 20877 events were selected. diff --git a/doc/source/users/tutorials/quickstart/simulating.rst b/doc/source/users/tutorials/quickstart/simulating.rst index 0207e62..80920e6 100644 --- a/doc/source/users/tutorials/quickstart/simulating.rst +++ b/doc/source/users/tutorials/quickstart/simulating.rst @@ -24,18 +24,18 @@ characteristics of the simulation: .. code-block:: bash - $ ctobssim - RA of pointing (degrees) (0-360) [83.63] 83.5 - Dec of pointing (degrees) (-90-90) [22.01] 22.8 - Radius of FOV (degrees) (0-180) [5.0] - Start time (UTC string, JD, MJD or MET in seconds) [2020-01-01T00:00:00] - Stop time (UTC string, JD, MJD or MET in seconds) [2020-01-01T01:00:00] 2020-01-01T01:00:00 - Lower energy limit (TeV) [0.03] - Upper energy limit (TeV) [200.] - Calibration database [prod2] - Instrument response function [South_0.5h] - Input model definition XML file [$CTOOLS/share/models/crab.xml] - Output event data file or observation definition XML file [events.fits] + $ ctobssim + RA of pointing (degrees) (0-360) [83.63] 83.5 + Dec of pointing (degrees) (-90-90) [22.01] 22.8 + Radius of FOV (degrees) (0-180) [5.0] + Start time (UTC string, JD, MJD or MET in seconds) [2020-01-01T00:00:00] + Stop time (UTC string, JD, MJD or MET in seconds) [2020-01-01T01:00:00] 2020-01-01T01:00:00 + Lower energy limit (TeV) [0.03] + Upper energy limit (TeV) [200.] + Calibration database [prod2] + Instrument response function [South_0.5h] + Input model definition XML file [$CTOOLS/share/models/crab.xml] + Output event data file or observation definition XML file [events.fits] Each line represents a query for one parameter value. The line starts with a short description of the parameter, followed by @@ -211,29 +211,29 @@ from this file are shown here: .. code-block:: none - 2017-11-28T14:13:40: === CTA observation === - 2017-11-28T14:13:40: Simulation cone ...........: RA=83.5 deg, Dec=22.8 deg, radius=5.5 deg - 2017-11-28T14:13:40: Time interval .............: 6.31109e+08 - 6.31112e+08 s - 2017-11-28T14:13:40: Photon energy range .......: 30 GeV - 72.3622611060088 GeV - 2017-11-28T14:13:40: Event energy range ........: 30 GeV - 72.3622611060088 GeV - 2017-11-28T14:13:40: Simulation area ..........: 1.97769e+09 cm2 - 2017-11-28T14:13:40: Use model ................: Crab - 2017-11-28T14:13:40: Normalization ............: 1 [Crab] - 2017-11-28T14:13:40: Flux .....................: 2.5413e-09 [Crab] photons/cm2/s - 2017-11-28T14:13:40: Normalized flux ..........: 2.5413e-09 [Crab] photons/cm2/s - 2017-11-28T14:13:40: Photon rate ..............: 5.0259 photons/s [Crab] - 2017-11-28T14:13:40: MC source photons ........: 18186 [Crab] - 2017-11-28T14:13:40: MC source events .........: 3544 [Crab] - 2017-11-28T14:13:40: MC source events .........: 3544 (all source models) - 2017-11-28T14:13:40: Photon energy range .......: 72.3622611060088 GeV - 174.543227745807 GeV - ... - 2017-11-28T14:13:40: MC source photons .........: 47219 [Crab] - 2017-11-28T14:13:40: MC source events ..........: 11356 [Crab] - 2017-11-28T14:13:49: MC events outside ROI .....: 0 - 2017-11-28T14:13:49: MC background events ......: 189477 - 2017-11-28T14:13:49: MC identifier 1 ...........: Crab - 2017-11-28T14:13:49: MC identifier 2 ...........: CTABackgroundModel - 2017-11-28T14:13:49: MC events .................: 200833 (all models) + 2017-11-28T14:13:40: === CTA observation === + 2017-11-28T14:13:40: Simulation cone ...........: RA=83.5 deg, Dec=22.8 deg, radius=5.5 deg + 2017-11-28T14:13:40: Time interval .............: 6.31109e+08 - 6.31112e+08 s + 2017-11-28T14:13:40: Photon energy range .......: 30 GeV - 72.3622611060088 GeV + 2017-11-28T14:13:40: Event energy range ........: 30 GeV - 72.3622611060088 GeV + 2017-11-28T14:13:40: Simulation area ..........: 1.97769e+09 cm2 + 2017-11-28T14:13:40: Use model ................: Crab + 2017-11-28T14:13:40: Normalization ............: 1 [Crab] + 2017-11-28T14:13:40: Flux .....................: 2.5413e-09 [Crab] photons/cm2/s + 2017-11-28T14:13:40: Normalized flux ..........: 2.5413e-09 [Crab] photons/cm2/s + 2017-11-28T14:13:40: Photon rate ..............: 5.0259 photons/s [Crab] + 2017-11-28T14:13:40: MC source photons ........: 18186 [Crab] + 2017-11-28T14:13:40: MC source events .........: 3544 [Crab] + 2017-11-28T14:13:40: MC source events .........: 3544 (all source models) + 2017-11-28T14:13:40: Photon energy range .......: 72.3622611060088 GeV - 174.543227745807 GeV + ... + 2017-11-28T14:13:40: MC source photons .........: 47219 [Crab] + 2017-11-28T14:13:40: MC source events ..........: 11356 [Crab] + 2017-11-28T14:13:49: MC events outside ROI .....: 0 + 2017-11-28T14:13:49: MC background events ......: 189477 + 2017-11-28T14:13:49: MC identifier 1 ...........: Crab + 2017-11-28T14:13:49: MC identifier 2 ...........: CTABackgroundModel + 2017-11-28T14:13:49: MC events .................: 200833 (all models) Each line starts with the UTC time at which the line has been written. In this run, 47219 Crab photons have been thrown. 11356 of these photons have been diff --git a/doc/source/users/tutorials/quickstart/spectrum.rst b/doc/source/users/tutorials/quickstart/spectrum.rst index aeec703..657bb81 100644 --- a/doc/source/users/tutorials/quickstart/spectrum.rst +++ b/doc/source/users/tutorials/quickstart/spectrum.rst @@ -26,19 +26,19 @@ bins between 100 GeV and 100 TeV: .. code-block:: bash - $ csspec - Input event list, counts cube, or observation definition XML file [events.fits] cntcube.fits - Input exposure cube file [NONE] expcube.fits - Input PSF cube file [NONE] psfcube.fits - Input background cube file [NONE] bkgcube.fits - Input model definition XML file [$CTOOLS/share/models/crab.xml] crab_results.xml - Source name [Crab] - Spectrum generation method (SLICE|NODES|AUTO) [AUTO] - Binning algorithm (LIN|LOG|FILE) [LOG] - Lower energy limit (TeV) [0.1] - Upper energy limit (TeV) [100.0] - Number of energy bins [20] 10 - Output spectrum file [spectrum.fits] + $ csspec + Input event list, counts cube, or observation definition XML file [events.fits] cntcube.fits + Input exposure cube file [NONE] expcube.fits + Input PSF cube file [NONE] psfcube.fits + Input background cube file [NONE] bkgcube.fits + Input model definition XML file [$CTOOLS/share/models/crab.xml] crab_results.xml + Source name [Crab] + Spectrum generation method (SLICE|NODES|AUTO) [AUTO] + Binning algorithm (LIN|LOG|FILE) [LOG] + Lower energy limit (TeV) [0.1] + Upper energy limit (TeV) [100.0] + Number of energy bins [20] 10 + Output spectrum file [spectrum.fits] There are different methods to create a binned spectrum: diff --git a/models/Makefile.am b/models/Makefile.am index 67dde84..c46f343 100644 --- a/models/Makefile.am +++ b/models/Makefile.am @@ -21,7 +21,7 @@ # Models to install models = crab.xml \ - crab_beyond.xml \ + crab_beyond.xml \ crab_poly.xml \ crab_profile.xml \ disk.xml \