Diwa Bhusal - Summer 2023 Research Journal

Weekly Slides:


Slides for my Talks:

REU Lunch Slides: https://docs.google.com/presentation/d/1Ztm_2ItfUjJZ_amRI14e0q_5fLjAqNwGjOg22ILRsQA/edit?usp=sharing

X-ray-TeV Symposium Slides: https://docs.google.com/presentation/d/1WjLEOxRtDnl_bznG56CeCZzGKiz1P6tTCA7lBuYVJdc/edit?usp=sharing

REU Final Talk Slides:

REU Final Report:

SRI Lida Orzeck Session Poster:


Week 1:

  • From 05/30/2023 - 06/02/2023 -- A Backlog:
    • I attended the NEVIS Labs Orientation and Intro to the Standard Model lecture
    • ROOT/PyROOT Computing Tutorial (led by William Seligman)
    • On Friday, I attended the VERITAS Intro lecture session, and then the Blazar/AGN Lecture by Professor Reshmi Mukherjee. And, I attended Colin Adam's VEGAS Intro and tutorial.
    • I worked on setting up VEGAS on my computer.

Week 2:

  • Monday 06/05/2023:
    • I went through the VEGAS Tutorial (I am waiting for verification from the UCLA group to access the data).
    • I started the Crab analysis but am waiting for all the files to download.
    • I also accessed the subthreshold candidates Zennodo data release. I also read/started reading a few papers today.
  • Tuesday 06/06/2023:
    • I attended the weekly GECo group meeting this morning.
    • Since the Crab analysis files were downloaded today, I finished up stages 1-5.
    • I also read a couple of papers and also worked with a Jupyter notebook file teaching me how to work with the Zennodo data release.
    • I also attended the Neutrino physics lecture in the Science Center.
  • Wednesday 06/07/2023:
    • I finished the Stage 6 Crab Analysis and produced all my plots and got my results from the analysis.
    • I read a couple of papers.
    • I attended the weekly VERITAS meeting and presented my project overview slides.
    • I met with Colin and got to look at his code.
  • Thursday 06/08/2023:
    • I worked on setting up a conda environment and setting up a Jupyter Notebook.
    • I spent some time reviewing Colin's code.
    • I also read some of a VEGAS documentation paper.
  • Friday 06/09/2023:
    • I met with Colin and Professor Zsusza Marka in the morning.
    • I continued to read some of the VEGAS Documentation papers.
    • I also worked on setting up my conda environment.

Week 3:

  • Monday 06/12/2023:
    • I set up my conda environment.
    • I read some GW papers sent to me by Prof. Zsusza Marka.
    • I also worked on my REU Lunch Presentation.
  • Tuesday 06/13/2023:
    • I attended the weekly GECo meeting in the morning.
    • I presented at the REU Lunch.
    • I attended the RARAF Lecture in the Science Center and the subsequent lab tour.
  • Wednesday 06/14/2023:
    • I presented some slides at the weekly VERITAS group meeting.
    • I also reviewed Colin's code and started running the O3 data for coincidents.
  • Thursday 06/15/2023:
    • I edited Colin's code to make it work for the O3 data and started running it to look for coincidences. It took most of the day to do so.
    • I read some GW papers while waiting for it to run.
  • Friday 06/16/2023:
    • I met with Professor Mukherjee, Prof. Zsusza Marka and Colin at 10 AM to discuss some updates and next steps.
    • I also attended the REU Lecture on Binary Neutron Stars.
    • I started working on the .csv file and met with Colin to understand what all I should be including.

Week 4:

  • Monday 06/19/2023:
    • Juneteenth Work Holiday
  • Tuesday 06/20/2023:
    • I attended the weekly GECo Meeting.
    • I worked on the .csv file for most of the day.
    • Then I attended the REU Lecture and subsequent lab tour of ATLAS.
  • Wednesday 06/21/2023:
    • I continued to work on the csv file and have integrated all the probabilities into them.
    • I'm now working on getting time stamps and seeing if the times of veritas observations overlap with the LIGO time of event.
    • I also attended the weekly VERITAS meeting.
  • Thursday 06/22/2023:
    • I continued work on the csv file.
    • I have been able to integrate the usable duration of the VERITAS observations and whether or not it overlaps with t_0.
    • Then, I got started on defining a value function for the time intervals.
  • Friday 06/23/2023:
    • I continued work on the csv file.
    • I was able to define a value function and I worked today mostly on double checking any errors that could've been made.
    • I also met with Professor Mukherjee and Professor Zsusza Marka about the project and some future steps, like including distance measurements incorporated.
    • I also attended the VERITAS REU Lecture by Professor Mukherjee and VERITAS Journal Club Weekly Meeting. I also read some GW Papers.

Week 5:

  • Monday 06/26/2023:
    • I continued work on the csv file and fixing up any errors and double checking all the values in it.
    • I read some GW Papers.
  • Tuesday 06/27/2023:
    • I attended the weekly GECo Meeting
    • I continued to fix errors in the csv file and incorporaitng pipeline information, LIGO Event Distance and Stdev (which I also converted to a redshift value and also accounted for the fact that the cWB pipeline did not have a distance measurement) and also sorting my chart and looking only at the events from the pipelines that has the highest pastro values.
    • I read some GW Documentation Papers.
  • Wednesday 06/28/2023:
    • I continued fixing errors in the .csv file in order to create and tripe check a final version for all errors and discrepancies.
    • I created a histogram and worked on my slides.
    • I attended the weekly VERITAS Meeting
    • I read some GW Papers and started a VERITAS Paper to read.
  • Thursday 06/29/2023:
    • I finished up the table of coincidences! Now some future directions are to meet and discuss which VERITAS Runs we should run an analysis for and also discuss if incorportating LIGO Run O2 Data inside. I also want to figure out how to possibly incorporate a relative distance measurement to weight my values depending on how far they are from us.
    • I attended the NanoGrav Announcement and the subsequent IceCube Announcement!
    • I read more GW Background Papers and I read the Journal Club Paper (the JHolder Thesis, Ch. 5 - A new technique for Image Analysis)
    • I reran Colin's code for -10 to 10,000 because I noticed a small error in one of the images that were produced.
  • Friday 06/30/2023:
    • I met with Colin this morning to discuss some issues about the data and discussed some future directions.
    • I explored the table and looking through the coincidences. I am trying to figure out a way to weight the different parameters depending on just how important each variable is to us.
    • I attended Weekly Journal Club, Svanik Covered the JHolder Thesis Ch.5
    • I met with Colin and Zsusza to discuss the table I created and how important each parameter may be. We discussed "VERITAS Observational Capability" which sort of is determined by VERITAS Event Overlap, overlaping with t_0 and whether or not it is a BNS Event. We also discussed whether or not we should pursue a time dependent analysis. We have decided to read up on it first and maybe decide from there.
    • I have finished up a preliminary runlist of interesting VTS Runs that may be relevant/interesting to pursue! I created four separate lists for now, each row sorted from lowest to highest redshift value and the ones with much larger redshift values, I have also noted for reference:
      1. Highest pAstro and VTS Overlap Prob: 95846
      2. Overlaps with t_0 and has pBNS Prob.: 95279, 95280, 95396, 95528, 95529, 95667, 95668
      3. Has pBNS Prob: 95397, 95829, 95718, 96002 & 96001 (redshift 0.125+/-0.054)
      4. Overlaps with t_0?: 95674, 95175, 95553 (redshift 0.351+/-0.115), 95931 (redshift 0.550+/-0.242)

Week 6:

  • Monday 07/03/2023:
    • Work Holiday
  • Tuesday 07/04/2023:
    • Fourth of July Work Holiday
  • Wednesday 07/05/2023:
  • Thursday 07/06/2023:
    • I incorporated the masses (lower mass, ratios, and chirp masses) into my data table. From this, I was able to identify 4 candidates and 1 outlier that will be interesting to analyze.
    • I read through more information about going through to do a time-dependent analysis.
  • Friday 07/07/2023:
    • I met with Prof. Zsuszanna Marka in the morning to discuss some future steps and directions. We mentioned and discussed looking through the 4 candidates and 1 outlier (as a start) and then coming back to it to do a time-dependent analysis later. However, we were concerned about possible collaboration policies in case we were to come back and reanalyzed. I discussed this later with Colin and we found that there would be a trials penalty if we re-analyzed using the time-dependent analysis.
    • Some of the 4 runs I looked at also had a second run that was spatially coincident but wasn't t_0 overlap, but instead was the VERITAS Run right after the event (Ex. 95279 and 95280) and I discussed whether or not we should also consider these runs. The conclusion that was made was that if it is a wobbled observation of the same target, then the second runs should be included in a time-dependent analysis but not included if it is a time-independent analysis.
    • I met with Professor Mukherjee for an update meeting and to discuss some possible next steps. We discussed the possibility of working with just 1 candidate for now (whichever one we find least valuable), in order to avoid a trials penalty in the future. We also discussed the time-dependent analysis and how it is a "flag" on VEGAS we can just turn on.
    • I discussed this time-dependent analysis with Colin and we discussed that it was never designed to deal with observations taken before or including t_0 and since the t^-a function blows up at t=0, it could be an issue. We would have to customize the code to include a time-dep that doesn't explode when t_0 is overlapped.
    • Colin outlined some things I should aim to get an understanding of:
      • The empirical time dependence for the prompt gamma-ray emission of previously detected short GRBs. What sort of time dependence do they have? Are there any trusted models we can look at?
        • Look at the afterglow emission paper and we could do a time-dependent search for afterglow emission and then later doing a prompt emission search, with different time-dependent models.
      • Understand the 1.7 s delay from t_0 from FermiGBM detection for BNS GW Event 170817. Is this a constant delay or is it dependent on other factors we can empirically calculate?

Week 7:

  • Monday 07/10/2023:
    • Read a couple of papers on BNS Mergers. I am trying to understand the GRB Jet that exits from a BNS Merger and how the dynamics of it work. So far, I have found:
      • BNS Mergers, when they first merge, will turn into a Hyper/Super Massive Neutron Star for a bit before, eventually, they core collapse into a Blackhole. This communicates to me, that there is a time delay in the merger event and when the Super/Hyper Massive NS Eventually Collapses. I am still trying to understand how this occurs.
  • Tuesday 07/11/2023:
    • I worked from home today because of some 2nd degree burns on my foot from cooking the night before that made it difficult for me to walk and required me to rebandage and clean often.
    • I attended the weekly GECo meeting at 10 AM.
    • I continued to read more papers, specficially a couple of BNS Review papers to understand the time delays.
  • Wednesday 07/12/2023:
    • I contacted Colin and we discussed dropping the time dependent analysis for now as the system is broken.
    • I started a lot of work on my REU Final Report and REU Final Talk.
  • Thursday 07/13/2023:
    • I continued working on my final report.
    • I ran VEGAS Stages 1-4 on my 5 candidates, now we are trying to figure out which time cuts to make in order to continue to Stage 5.
    • I attended the VERITAS Collaboration Meeting's VEGAS Tutorial.
  • Friday 07/14/2023:
    • I continued work on my final report.
    • I tried to read a couple of papers to understand the time cuts we should be using for the most statistically signifigant results.
    • I attended the REU Lecture at 1 PM on Biophysics and the subsequent lab tour at 2 PM.
    • I attended the Weekly Journal Club, this week led by Claudia & Ana, meeting at 3 PM.

Week 8:

  • Monday 07/17/2023: Back to Nevis!
    • I worked most of today on my final report and I also researched the time-cuts, for now I have gotten that around 0 to 200 is the most signifigant but the issue is whether a span of that much would be helpful or harmful and would it wash out any signifigance from the first couple of seconds.
    • I met with Colin and Professor Mukherjee about this question. Here is what we sort of settled on:
      • Meeting with other VERITAS Collaborators (Specifically David Williams) to discuss next steps and what he suggests we should do. Question is, what is the lowest time range we could analyze on VEGAS...?
      • For now, run stages 1-4 and maybe stage 5/6 as well. However, do not look at the results as to not add bias!
  • Tuesday 07/18/2023:
    • I worked most of today on the final report and awaited response from David Williams.
    • Colin and I discussed this morning about his idea to find what is the lowest time that would still be statistically signifigant, awaiting response and further discussion about this. Colin basically said we can use a run from a dark patch of sky where we do not expect to find anything and see which are the most statistically sound time cuts, testing for a minute, 10 mins, etc.
    • I attended the weekly GECo meeting at 11 AM (rescheduled)
    • I attended the 2 PM VERITAS Lab Tour.
    • I attended the 3 PM Neutrinos Lab Tour.
  • Wednesday 07/19/2023: Brookhaven National Lab Tour!!!
    • Colin and I discussed moving forward with his idea to find the shortest time we could look at.
      • Looking at source: Draco (Run #: 77335) for 1 min, 2.5 min, 5 min, 10 min, 15 min
      • We do not expect to see anything for this source. We will run stages 1-4 and then run stages 5/6 with the different time cuts and then look at the signifigances histograms and see the lowest analysis duration that will give us a gaussian distrubution.
  • Thursday 07/20/2023:
    • I continued working on my paper!
    • I started analyzing the blank sky run. It didn't get very far before I ran into some issues and after some troubleshooting assistance with Colin, we realized I was having problems querrying through to find my runs.
    • I contacted Takeru Hyashi at UCLA to help fix my errors, and I am awaiting response.
    • For now, Colin has uploaded the BlankSky run into the system so I am able to work on it. I ran stages 1-4 on it and thought about how I should integrate the timecuts properly for running stages 5&6 on them.
  • Friday 07/21/2023:
    • No more accessing data issues! It's been fixed with a new Public_Key id.
    • For my 1 minute time cut (60 second) I was able to plot the signifigances and get the histogram using Colin's code!
    • I tested running stage5/6 in the command line instead of on the VEGAS Process Script and giving it manually separate places to go to. And also, figuring out if I can run them in batch!
    • I have submitted runs for my candidates so I have them processing stages 1-4 of VEGAS right now for me to access on Monday.
    • I attended the REU Lecture on Atomic, Molecular and Optical Physics by Professor Ana Asenjo-Garcia at 2 PM and the subsequent lab tour at 3 PM.
    • I attended the weekly Journal Club meeting, this week by Jack and Anna!
    • I have started running all of the Stage5s in the background. Hopefully, it will run successfully without overwriting my files. From what I see, it didn't overwrite anything when I tested it earlier in the day!

Week 9:

  • Over the weekend:
    • The Stages 1-4 for our candidates finished running!
    • I ran all of the stage5s for the Draco timecuts and then moved onto running all the stage6s.
  • Monday 07/24/2023: Nevis BBQ Day!
    • I ran into an issue with the fits files and they were getting overwritten when I was running them in batch.
  • Tuesday 07/25/2023:
    • I attended the weekly GECo Meeting!
    • Finished up slides for Symposium Talk and practiced my talk!
    • Cleared up the fits file issues and had the stages running.
    • We discussed 20 minutes was the timespan we should be analyzing for based off of the results we were getting from the Draco blanksky analysis.
      • However, most of the candidates did not have 20 minutes of time in the same run to be able to analyze so we had to choose one of the sources that had a wobbled observation right after we could pull into.
    • I ran stages 5 and 6 on my candidate and produced my plots.
      • Some scaling and issues with the stage6 plot because theres (I think) a bin or something that is odd so I am trying to fix it. However, the distributions on the SigDist plot look mostly the same.
    • I attended the REU Lecture by IAXO/GAPs and subsequent lab tour.
  • Wednesday 07/26/2023: X-Ray-TeV Symposium!
    • I presented a talk at the Symposium (slides linked above)
    • I drafted up my SRI Poster!
    • I reran stage6 and set S6A _Batch to false so that I could get my plots.
  • Thursday 07/27/2023:
    • I attended the REU Lunch.
    • I worked on my SRI Poster and the REU Paper for most of the day.
  • Friday 07/28/2023:
    • I finished my SRI Poster and submitted it
    • I met with Professor Zsusza Marka in the morning
    • I attended the Alumni Panel!

Week 9:

  • Monday 07/31/2023:
    • I worked on the REU Paper and Presentation!
  • Tuesday 08/01/2023:
    • I gave a practice talk to the VERITAS group
    • I continued working on my Paper and Presentation
  • Wednesday 08/02/2023: SRI Poster Session!
  • Thursday 08/03/2023-Friday 08/04/2023: REU Talks and Presentation!


Edit | Attach | Watch | Print version | History: r23 < r22 < r21 < r20 < r19 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r23 - 2023-08-01 - DiwaBhusal
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback