Research Journal: Jamil Batshoun

CTA simulations with gammapy

June 28, 2021

  1. What was done:
    1. Met with Massimo to copy cta files over to my directory and had ds9 set up to view fits files.
      1. Massimo sent me an email with a pdf of the CTA DC1 to read.
    2. Nevis VERITAS/CTA meeting
  2. Things left to do:
  3. Notes:
Week 4
  1. What was done:
    1. Completed tutorials:
      1. Overview overview.ipynb
        CTA with gammapy cta.ipynb
        3D map simulation simulate_3d.ipynb
        1D spectrum simulation spectrum_simulation.ipynb
        Point source sensitivity cta_sensitivity.ipynb
    2. Went through some exercises in the tutorials. Included plots I added as attatchments to this page.
    3. Attended VERITAS summer meeting. Went to blazar science and ACG talks.
  2. Things left to do:
    1. Complete exercises from the tutorials i've completed so far to practice what they covered.
    2. Will go through irf section of cta.ipynb again to take notes on the plots the code produces.
    3. Need to complete taking notes on the papers I have read/need to read so far.
June 24, 2021
  1. What was done today:
    1. Attended 1st half of today's session of the VERITAS summer meeting.
    2. Attended REU talk
    3. Took notes on cta.ipynb. Worked on exercises in overview.ipynb but struggled getting through the first exercise.
      1. Made a sky map of hess observations!
  2. Things left to do:
    1. Finish notes on cta.ipynb
  3. Notes
    1. overview.ipynb exercise notes: my best guess is that i am supposed to use SphericalCircle from astropy.visualization.wcsaxes.
    2. ^^^Figured it out
June 16, 2021
  1. What was done today:
    1. read about fits, found the standard paper for fits
    2. read about the icrs, dec, and ra
    3. looked for a formal definition of the dl3-dr1 data from hess. did not find much info about what a data level is
    4. ran into same issue with the path associated with the GAMMAPY_DATA environment variable
      1. I fixed this by adding a line to my .profile on the nevis computers. I don't see this becoming a problem again.
    5. completed tutorial 2
  2. What is left to do:
    1. Read up on background modeling.
    2. Understand what the classes of gammapy do
    3. get a better understanding of how the analysis (both high level and lower level gammapy api)
June 17-18, 2021
  1. What was done today:
    1. Went over again the analysis steps found in tutorial "analysis 1"
  2. Things left to do:
    1. started working through tutorial "overview"
  3. Notes:
    1. How higher-level analysis works in gammapy:
      1. tutorial intro intro definition:
        1. (data reduction): observations reduced to a specified geometry (1D/3D). Combined with IRFS, this makes a dataset
        2. (Analysis): set a physical model and fit it to the datasets
      2. definition based on code:
        1. create config object for the settings used in your analysis
        2. set config object's type, stack, geom, background, fit, and flux_point settings settings
        3. (optional) write config object to a YAML file to use later.
        4. select observations to use
        5. data reduction
        6. (optional) save dataset to disk
        7. create model_config object with model settings
        8. set model_config object to analysis object
        9. use run_fit on analysis object
        10. create model best fit yaml file.
        11. estimate flux points
    2. questions:
      1. what is the model-best-fit.yaml file for? what about the model-best-fit_covariance.dat file?
      2. is data really reduced?
June 9, 2021
  1. What was done today:
    1. Added gammapy resources to .myprofile file
    2. started working on gammapy tutorials
      1. code to start tutorials:
        1. ssh -L 9888:localhost:9888

        2. module load gammapy/0.18
        3. gammapy_start
        4. jupyter notebook --no-browser --port=9888
        5. copy & paste resulting link to chrome and use jupyter notebook normally
      2. COMPLETED: 1st gammapy tutorial analysis_1.ipynb
  2. Things left to do:
    1. create a .myprofile and copy contents of .profile into it.
    2. Remove gammapy line in .profile so that it is the same as when I first started (This might be unnecessary).
  3. Notes:
    1. Useful commands
      1. To check where the "juypter notebook" command is being used from (i.e. gammapy vs python): which jupyter-notebook
    2. Made an error when executing "export $GAMMAPY_DATA". executed in jbatshoun directory and not datasets directory. Noticed in tutorial what should have been the real path and made corrections.
June 9, 2021
  1. What was done today:
  2. Things Left to do:
  3. Notes:
    1. aa5351-06: Is the IRF a good measure of telescope performance? what about the psf?
June 8, 2021
  1. Read through gammapy and cta explanation of IRFs
  2. Read gammapy overview

June 7, 2021 (Meeting Massimo, Jamil)

  1. Introduction to project
  2. To-do list
    2. Read overview carefully
    3. Installation: no need to worry about it, we have shared tools on the nevis cluster
      1. log into the nevis cluster
      2. module load gammapy/0.18
      3. gammapy_start
    4. download tutorial datasets (gammapy download tutorials --release 0.18.2 , export GAMMAPY_DATA=<>)
    5. go through tutorials
  3. Weekly meeting with Massimo: Thursdays, 11 AM

-- Massimo Capasso - 2021-06-27

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng cta_gps_highestenergyphotoninfo.png r1 manage 97.9 K 2021-06-27 - 00:30 JamilBatshoun  
PNGpng hess_surveymap.png r1 manage 75.6 K 2021-06-27 - 00:25 JamilBatshoun  
PNGpng markerandcirclegalacticcenter.png r1 manage 185.1 K 2021-06-27 - 00:24 JamilBatshoun  
Edit | Attach | Watch | Print version | History: r22 | r16 < r15 < r14 < r13 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r14 - 2021-06-28 - JamilBatshoun
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2022 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback