Difference between revisions of "LST1Analysis MVA"
(→Ctools) |
|||
Line 76: | Line 76: | ||
ctskymap | ctskymap | ||
ds9 skymap.fits | ds9 skymap.fits | ||
− | + | ctbin | |
+ | fv cntcube.fits | ||
+ | ctexpcube | ||
+ | ctpsfcube | ||
+ | ctbkgcube | ||
</pre> | </pre> |
Revision as of 09:07, 19 August 2021
Contents
Location of LST1 data at PIC
- Directory: /pnfs/pic.es/data/cta/LST/Data/ZFITS/
- Web: LST1 ZFITS files on the web
Same user/password as in CTA portal
Location of LST1 at La Palma
- Directory: /fefs/onsite/data
Use cp01/cp02 to copy out data, never the login machine
Run summary info/data-check
- Run-summary web in GAE
LST1 account
CTApipe
git clone https://github.com/cta-observatory/ctapipe.git cd ctapipe CTAPIPE_VER=0.10.5 wget https://raw.githubusercontent.com/cta-observatory/ctapipe/v$CTAPIPE_VER/environment.yml conda env create -n cta -f environment.yml conda activate cta conda install -c conda-forge ctapipe=$CTAPIPE_VER
Analysis of LST1 data
Donwload LST analysis software cta-lstchain (wrapper that uses centrally supported ctapipe routines)
git clone https://github.com/cta-observatory/cta-lstchain.git cd cta-lstchain LSTCHAIN_VER=0.7.2 wget https://raw.githubusercontent.com/cta-observatory/cta-lstchain/v$LSTCHAIN_VER/environment.yml conda env create -n lst -f environment.yml conda activate lst pip install lstchain==$LSTCHAIN_VER rm environment.yml In case of problems with conda env try mamba: mamba env create -f environment.yml
Files needed to analyze golden Run 442 (20190527):
- Pedestal run: pedestal_file_run446_0000.fits
- Low level DRS4 calibration: calibration.hdf5
Run directly jupyter notebook
cd cta-lstchain/notebooks jupyter notebook -> select Analyze_real_muon_data.ipynb
Convert example jupyter notebook to python
cd cta-lstchain/notebooks jupyter nbconvert --to python Analyze_real_muon_data.ipynb
Fix path of data, pedestal and calibration and run!
Lines to fix: source = event_source(input_url="../data/LST-1.4.Run00442.0001.fits.fz",max_events=None) pedestal_path="../pedestal/pedestal_file_run446_0000.fits", with HDF5TableReader('../calibration/calibration.hdf5') as h5_table: python Analyze_real_muon_data.py
Ctools
ctobssim edisp=yes fv events.fits ctskymap ds9 skymap.fits ctbin fv cntcube.fits ctexpcube ctpsfcube ctbkgcube