Konstancja Satalecka / ctools

Commit 4e61d90f99342be9beda259a8eb3b3c193f718b6

Authored by Jürgen Knödlseder
1 parent 890ca618

Longitude in skymaps and counts maps now increases towards the left. Only releva… 

…nt projections are proposed in parameter file.
Inline Side-by-side
Showing 9 changed files with 19 additions and 12 deletions
ChangeLog
  View file @4e61d90
1 -2013-11-20 Juergen Knoedlseder <jurgen.knodlseder@irap.omp.eu> 1 +2013-11-29 Juergen Knoedlseder <jurgen.knodlseder@irap.omp.eu>
2 2
3 * Version 0.7.0 released 3 * Version 0.7.0 released
4 ======================== 4 ========================
@@ -7,6 +7,8 @@ @@ -7,6 +7,8 @@
7 ctobssim now only simulates events within ROI and supports 7 ctobssim now only simulates events within ROI and supports
8 all kinds of data space models 8 all kinds of data space models
9 ctbin now implements LIN and FILE binning algorithms 9 ctbin now implements LIN and FILE binning algorithms
  10 + ctbin, ctskymap and ctmodel now generate maps with
  11 + increasing longitude to the left
10 12
11 2012-12-21 Juergen Knoedlseder <jurgen.knodlseder@irap.omp.eu> 13 2012-12-21 Juergen Knoedlseder <jurgen.knodlseder@irap.omp.eu>
12 14
@@ -2,7 +2,7 @@ New Features and Important Changes in ctools 00-07-00 @@ -2,7 +2,7 @@ New Features and Important Changes in ctools 00-07-00
2 2
3 Author(s) : Juergen Knoedlseder 3 Author(s) : Juergen Knoedlseder
4 4
5 -20 November 2013 5 +29 November 2013
6 6
7 7
8 1. Introduction 8 1. Introduction
@@ -52,6 +52,9 @@ ENERGYBINS is expected. @@ -52,6 +52,9 @@ ENERGYBINS is expected.
52 52
53 The ebinalg and ebinfile parameters are now automatic. 53 The ebinalg and ebinfile parameters are now automatic.
54 54
  55 +Counts maps generated by ctbin now have the longitude increasing
  56 +towards the left.
  57 +
55 58
56 7. ctlike - CTA maximum likelihood fitting 59 7. ctlike - CTA maximum likelihood fitting
57 ------------------------------------------ 60 ------------------------------------------
@@ -60,12 +63,14 @@ None @@ -60,12 +63,14 @@ None
60 63
61 8. ctskymap - CTA sky mapping 64 8. ctskymap - CTA sky mapping
62 ----------------------------- 65 -----------------------------
63 -None 66 +Skymaps generated by ctskymap now have the longitude increasing
  67 +towards the left.
64 68
65 69
66 9. ctmodel - CTA counts map modelling 70 9. ctmodel - CTA counts map modelling
67 ------------------------------------- 71 -------------------------------------
68 -None 72 +Model maps generated by ctmodel now have the longitude increasing
  73 +towards the left.
69 74
70 75
71 10. Examples 76 10. Examples
1 ctools information 1 ctools information
2 ================== 2 ==================
3 -Version: 00-06-00 (December 21th, 2012) 3 +Version: 00-07-00 (November 29th, 2013)
4 Author: Juergen Knodlseder (jurgen.knodlseder@irap.omp.eu) 4 Author: Juergen Knodlseder (jurgen.knodlseder@irap.omp.eu)
5 -GammaLib dependency: 00-07-00 5 +GammaLib dependency: 00-08-00
6 6
7 7
8 License information 8 License information
src/ctbin/ctbin.cpp
  View file @4e61d90
@@ -590,7 +590,7 @@ void ctbin::bin_events(GCTAObservation* obs) @@ -590,7 +590,7 @@ void ctbin::bin_events(GCTAObservation* obs)
590 590
591 // Create skymap 591 // Create skymap
592 GSkymap map = GSkymap(m_proj, m_coordsys, 592 GSkymap map = GSkymap(m_proj, m_coordsys,
593 - xref, yref, m_binsz, m_binsz, 593 + xref, yref, -m_binsz, m_binsz,
594 m_nxpix, m_nypix, m_enumbins); 594 m_nxpix, m_nypix, m_enumbins);
595 595
596 // Initialise binning statistics 596 // Initialise binning statistics
src/ctbin/ctbin.par
  View file @4e61d90
@@ -39,7 +39,7 @@ coordsys, s, a, &quot;CEL&quot;, CEL|GAL,,&quot;Coordinate system (CEL - celestial, GAL - galac @@ -39,7 +39,7 @@ coordsys, s, a, &quot;CEL&quot;, CEL|GAL,,&quot;Coordinate system (CEL - celestial, GAL - galac
39 xref, r, a, 83.63,,, "First coordinate of image center in degrees (RA or galactic l)" 39 xref, r, a, 83.63,,, "First coordinate of image center in degrees (RA or galactic l)"
40 yref, r, a, 22.01,,, "Second coordinate of image center in degrees (DEC or galactic b)" 40 yref, r, a, 22.01,,, "Second coordinate of image center in degrees (DEC or galactic b)"
41 axisrot, r, h, 0.0,,, "Rotation angle of image axis, in degrees" 41 axisrot, r, h, 0.0,,, "Rotation angle of image axis, in degrees"
42 -proj, s, a, "CAR",,, "Projection method e.g. AIT|ARC|CAR|GLS|MER|NCP|SIN|STG|TAN" 42 +proj, s, a, "CAR", AZP|CAR|MER|STG|TAN,, "Projection method e.g. AZP|CAR|MER|STG|TAN"
43 43
44 # 44 #
45 # Standard parameters 45 # Standard parameters
src/ctmodel/ctmodel.cpp
  View file @4e61d90
@@ -495,7 +495,7 @@ void ctmodel::setup_obs(void) @@ -495,7 +495,7 @@ void ctmodel::setup_obs(void)
495 495
496 // Setup skymap 496 // Setup skymap
497 GSkymap map = GSkymap(m_proj, m_coordsys, 497 GSkymap map = GSkymap(m_proj, m_coordsys,
498 - m_xref, m_yref, m_binsz, m_binsz, 498 + m_xref, m_yref, -m_binsz, m_binsz,
499 m_nxpix, m_nypix, m_enumbins); 499 m_nxpix, m_nypix, m_enumbins);
500 500
501 // Create model cube from sky map 501 // Create model cube from sky map
src/ctmodel/ctmodel.par
  View file @4e61d90
@@ -52,7 +52,7 @@ coordsys, s, a, &quot;CEL&quot;, CEL|GAL,,&quot;Coordinate system (CEL - celestial, GAL - galac @@ -52,7 +52,7 @@ coordsys, s, a, &quot;CEL&quot;, CEL|GAL,,&quot;Coordinate system (CEL - celestial, GAL - galac
52 xref, r, a, 83.63,,, "First coordinate of image center in degrees (RA or galactic l)" 52 xref, r, a, 83.63,,, "First coordinate of image center in degrees (RA or galactic l)"
53 yref, r, a, 22.01,,, "Second coordinate of image center in degrees (DEC or galactic b)" 53 yref, r, a, 22.01,,, "Second coordinate of image center in degrees (DEC or galactic b)"
54 axisrot, r, h, 0.0,,, "Rotation angle of image axis, in degrees" 54 axisrot, r, h, 0.0,,, "Rotation angle of image axis, in degrees"
55 -proj, s, a, "CAR",,, "Projection method e.g. AIT|ARC|CAR|GLS|MER|NCP|SIN|STG|TAN" 55 +proj, s, a, "CAR", AZP|CAR|MER|STG|TAN,, "Projection method e.g. AZP|CAR|MER|STG|TAN"
56 # 56 #
57 # Standard parameters 57 # Standard parameters
58 #==================== 58 #====================
src/ctskymap/ctskymap.cpp
  View file @4e61d90
@@ -437,7 +437,7 @@ void ctskymap::init_map(GCTAObservation* obs) @@ -437,7 +437,7 @@ void ctskymap::init_map(GCTAObservation* obs)
437 } // endelse: map centre set to pointing 437 } // endelse: map centre set to pointing
438 438
439 // Create skymap 439 // Create skymap
440 - m_skymap = GSkymap(m_proj, m_coordsys, xref, yref, m_binsz, m_binsz, 440 + m_skymap = GSkymap(m_proj, m_coordsys, xref, yref, -m_binsz, m_binsz,
441 m_nxpix, m_nypix, 1); 441 m_nxpix, m_nypix, 1);
442 442
443 // Return 443 // Return
src/ctskymap/ctskymap.par
  View file @4e61d90
@@ -34,7 +34,7 @@ coordsys, s, a, &quot;CEL&quot;, CEL|GAL,,&quot;Coordinate system (CEL - celestial, GAL - galac @@ -34,7 +34,7 @@ coordsys, s, a, &quot;CEL&quot;, CEL|GAL,,&quot;Coordinate system (CEL - celestial, GAL - galac
34 xref, r, a, 83.63,,, "First coordinate of image center in degrees (RA or galactic l)" 34 xref, r, a, 83.63,,, "First coordinate of image center in degrees (RA or galactic l)"
35 yref, r, a, 22.01,,, "Second coordinate of image center in degrees (DEC or galactic b)" 35 yref, r, a, 22.01,,, "Second coordinate of image center in degrees (DEC or galactic b)"
36 axisrot, r, h, 0.0,,, "Rotation angle of image axis, in degrees" 36 axisrot, r, h, 0.0,,, "Rotation angle of image axis, in degrees"
37 -proj, s, a, "CAR",,, "Projection method e.g. AIT|ARC|CAR|GLS|MER|NCP|SIN|STG|TAN" 37 +proj, s, a, "CAR", AZP|CAR|MER|STG|TAN,, "Projection method e.g. AZP|CAR|MER|STG|TAN"
38 38
39 # 39 #
40 # Standard parameters 40 # Standard parameters