The data challenge module requires the MEGAlib code, available here. Among other things, MEGAlib simulates the emission from any (MeV) gamma-ray source, simulates the instrument response, performs the event reconstruction, and performs the high-level data analysis. See the above link for more details regarding the MEGAlib package.
For any help/problems with running the data challenge module please contact Chris Karwin at: [email protected].
- March 2023: The first data challenge is now available! It can be found here.
The main purpose of this repository is to simulate the all-sky data that will be observed by COSI. The primary code is run_dc.py, which can be called with run_sims.py, with the main input parameters passed via inputs.yaml. Additionally, parallel simulations with multiple time bins can be ran using run_parallel_sims.py, which distributes the time bins to seperate compute nodes. The pipeline also supports the use of mcosima with numerous cores per compute node. The modules can be ran directly from the command line, or submitted to a batch system, which allows them to be easily employed for generating multiple/long simulations.
See Source_Library for available sources. Let us know if you want any specific source added!
1. Download cosi-data-challenges directory:
$ git clone https://github.com/ckarwin/cosi-data-challenges.git
2. Install with pip:
$ cd cosi-data-challenges
$ pip install -e .
3. Start a new analysis directory, and enter the commmand-line prompt:
$ make_sim
4. Specify inputs in inputs.yaml
5. Run setup script:
$ python run_sim_setup.py
- This will setup the source directory and copy all needed files for running the code.
6. To run the code:
- Uncomment the functions inside run_sims.py that you want to run.
- The code can be ran directly from the terminal or submitted to a batch system.
- To run from the terminal use python run_sims.py.
- To run parallel jobs in cosima with numerous time bins use python run_parallel_sims.py.
- To submit a single job use python submit_jobs.py.
7. If running parallel jobs:
- In run_sims.py uncomment all functions except mimrec.
- Run: python run_parallel_sims.py.
- After all the jobs finish, uncomment just the mimrec function in run_sims.py, then run: python submit_jobs.py.
8. Note that the batch submission commands in run_parallel_sims.py and submit_jobs.py may need to be modified based on the user's specific batch system.
- The number of iterations in ExtractImage.cxx needs to be changed manually if using a different value than the default (20). Specifically, this is at lines 8 and 10. For x iterations: 20 --> x and 22 --> x+2. This will be automated soon.