Pablo Drake Research Journal (Summer '22)



  • Plotted different graphs containing blazar population data both from TeVCat and 4LAC Fermi-LAT Catalogue. In particular, I completed two skymaps of blazars, two pie charts with blazar classifications and a scatterplot showing Fermi blazars' photon index vs peak frequency, with their classification as well. All the data was retrieved from the catalogues and the plotting was carried out in Jupyter notebooks, using matplotlib. This work is to be included in some of Professor Mukherjee's slides.
  • Read Observing the Energetic Universe at Very High Energies with the VERITAS Gamma Ray Observatory by Prof. Mukherjee ( Wrote a set of reference notes.
  • Read Cherenkov Telescope Array: CTA CERN Courier article. Wrote a set of reference notes.
  • Attended Daniela VEGAS tutorial
  • Started the Crab Analysis tutorial with VEGAS
  • Continued with the Crab Analysis (stage 6)
  • Read Exploring the high-energy gamma-ray spectra of TeV blazars by Doctor Qi Feng ( Wrote a set of reference notes.
  • Centered all my work around the analysis of the May 2022 flare of the FSRQ 4FGL J1350.8+3033. The structure of this analysis will be the following: carry out a fermipy ROI analysis, plot the lightcurve of the source, plot the lightcurve of the flare, carry out an spectral analysis both during flare buildup, maximum, and decrease (and potentially also during quiescent state), carry out VEGAS spectral analysis of the source.
  • Completed the VEGAS tutorial
  • Started working again with fermipy. An error popped up when attempting to carry out the ROI analysis of FSRQ J1350.8+3033.
Fatal Python error: init_fs_encoding: failed to get the Python codec of the filesystem encoding Python

runtime state: core initialized

ModuleNotFoundError: No module named 'encodings'


  • Started writing a log/manual fermipy document, noting down the different errors that appear.
  • Found out the error was caused when simultaneously calling through my .myprofile doc VEGAS and fermipy code, possibly because of a conflict between subpackages. Another error appeared:
Traceback (most recent call last):

File "./", line 348, in <module> free_sources=pipeline_config.get('free_sources'))

File "./", line 50, in run_analysis gta.setup()

File "/a/share/ged/src/miniconda3/envs/fermipy-1.0.1/lib/python3.7/site-packages/fermipy/", line 1082, in setup c.setup(overwrite=overwrite)

File "/a/share/ged/src/miniconda3/envs/fermipy-1.0.1/lib/python3.7/site-packages/fermipy/", line 5121, in setup self._ltc = LTCube.create(self.files['ltcube'])

File "/a/share/ged/src/miniconda3/envs/fermipy-1.0.1/lib/python3.7/site-packages/fermipy/", line 192, in create ltc = cls.create_from_fits(files[0])

IndexError: list index out of range

  • Read EGRET (GeV) Blazars by Prof. Mukherjee. Wrote a set of reference notes.



  • Worked on adjusting the legends and the colors for the skymaps and pie chart plots that I coded last week. Also adjusted the sizes in order for them to look better in the slides. Finally, researched and included a paragraph of information on general blazar populations at the present time.
  • Reinstalled fermipy using the files provided by Massimo ( The ROI analysis worked without any errors.
  • As per Prof. Mukherjee's suggestion, I started working in both a matplotlib/data plotting manual as continued updating the fermipy guide.
  • Completed ROI analysis. Redid the analysis in order for the graphs to be created (sometimes the plots aren't created automatically when the finishes).
  • Read The future of gamma-ray astronomy by Jürgen Knödlseder ( Wrote a set of notes on it.
  • Started working on creating a full dataset lightcurve.





  • As per Deivid's suggestion, I attempted carrying out lightcurve analysis freeing lower TS sources and a larger radius. Carried out this modification for the whole baseline analysis, the 1-week timeframe around the flare, and the source J0658. The fit still won't converge, and by looking at the residmap, it seems like there is no abnormal source, nothing unaccounted by the model.
  • I had a Zoom chat with Janeth Valverde to consult on the lightcurve problem. Janeth pointed out that the problem might have to do with some of fermitools not running properly for a source this faint. Instead of using Ari's pipeline, she recommended doing a general check on that by running fermitools directly on the command line. She also recommended focusing solely in the last 6 months, and doing an unbinned analysis.
  • Janeth also discussed trying out an adaptive binning strategy (in parallel to the other analyses). She shared this example as a demonstration of the usefulness of adaptive binning.
  • Contacted Ari Brill

-- Massimo Capasso - 2022-06-16


Edit | Attach | Watch | Print version | History: r13 | r7 < r6 < r5 < r4 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r5 - 2022-07-07 - PabloDrake
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