test_write.py
1 |
from gammalib import * |
---|---|
2 |
from ctools import * |
3 |
import math |
4 |
import os |
5 |
from gfuncs import * |
6 |
from userclass import * |
7 |
datadir = '/home/hfm/chiachun/disk1/gammalib/CubeAnalysis/cta1dc/crab/data/'
|
8 |
run = 23523
|
9 |
|
10 |
# files for this run
|
11 |
f_event = '{datadir}/run_00{run}_std_eventlist.fits'.format(**globals()) |
12 |
f_aeff = '{datadir}/run_00{run}_std_arf.fits'.format(**globals()) |
13 |
f_psf = '{datadir}/run_00{run}_std_psf.fits'.format(**globals()) |
14 |
f_bg = '/home/hfm/chiachun/disk1/gammalib/CubeAnalysis/irf_file.fits'.format(**globals()) |
15 |
|
16 |
# set up a cta observation
|
17 |
cta = GCTAObservation() |
18 |
cta.id('{run}'.format(**globals())) |
19 |
cta.load(f_event) |
20 |
rsp = cta.response() |
21 |
rsp.load_aeff(f_aeff) |
22 |
rsp.load_psf(f_psf) |
23 |
rsp.load_background(f_bg) |
24 |
bg = GCTABackground3D(rsp.background()) |
25 |
bg_table = bg.rsp_table() |
26 |
newhdu = GFitsBinTable() |
27 |
bg_table.write(newhdu) |
28 |
fits = GFits() |
29 |
fits.append(newhdu) |
30 |
fits.saveto('newfile.fits', True) |
31 |
fits.close() |
32 |
|
33 |
# set up a new cta observation
|
34 |
ncta = GCTAObservation() |
35 |
ncta.id('{run}'.format(**globals())) |
36 |
ncta.load(f_event) |
37 |
nrsp = cta.response() |
38 |
nrsp.load_aeff(f_aeff) |
39 |
nrsp.load_psf(f_psf) |
40 |
nrsp.load_background('newfile.fits')
|
41 |
bg = GCTABackground3D(rsp.background()) |
42 |
logE = 0
|
43 |
print( bg(logE, 0 ,0 ) ) |