ctobssim.py
| 1 |
import gammalib |
|---|---|
| 2 |
import ctools |
| 3 |
import cscripts |
| 4 |
|
| 5 |
obssim = ctools.ctobssim() |
| 6 |
obssim['ra'] = 186.151 |
| 7 |
obssim['dec'] = -4.82077 |
| 8 |
obssim['rad'] = 2.16 |
| 9 |
obssim['tmin'] = 0 |
| 10 |
obssim['tmax'] = 3420 |
| 11 |
obssim['emin'] = 0.1 |
| 12 |
obssim['emax'] = 100.0 |
| 13 |
obssim['caldb'] = 'prod3b-v2' |
| 14 |
obssim['irf'] = 'South_z20_0.5h' |
| 15 |
obssim['inmodel'] = './model.xml' |
| 16 |
obssim['outevents'] = 'event1.fits' |
| 17 |
obssim.execute() |
| 18 |
|
| 19 |
# Second pointing
|
| 20 |
obssim = ctools.ctobssim() |
| 21 |
obssim['ra'] = 186.151 |
| 22 |
obssim['dec'] = -7.61578 |
| 23 |
obssim['rad'] = 2.16 |
| 24 |
obssim['tmin'] = 0 |
| 25 |
obssim['tmax'] = 3420 |
| 26 |
obssim['emin'] = 0.1 |
| 27 |
obssim['emax'] = 100.0 |
| 28 |
obssim['caldb'] = 'prod3b-v2' |
| 29 |
obssim['irf'] = 'South_z20_0.5h' |
| 30 |
obssim['inmodel'] = './model.xml' |
| 31 |
obssim['outevents'] = 'event2.fits' |
| 32 |
obssim.execute() |