/
gammapy-pfspec
executable file
·54 lines (46 loc) · 2.57 KB
/
gammapy-pfspec
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
#!/usr/bin/env python
# Licensed under a 3-clause BSD style license - see LICENSE.rst
"""Create spectra from VHE event lists in FITS format.
prog [options] FILE [ARF RMF]
FILE can either be an indiviual .fits/.fits.gz file or a batch file.
In case it is a individual file, the ARF and RMF must also be specified.
The bankfile must contain three columns: data file, ARF, and RMF.
TODO: document
"""
# Parse command line arguments
from gammapy.utils.scripts import argparse, GammapyFormatter
parser = argparse.ArgumentParser(description=__doc__,
formatter_class=GammapyFormatter)
# TODO: fix next argument
parser.add_argument('input_file_names', type=str,
help='Input file names.')
parser.add_argument('-p', '--analysis_position', type=str, default=None,
help='Center of the skymap in RA and Dec (J2000) in degrees. '
'Format: \'(RA, Dec)\', including the quotation marks. '
'If no center is given, the source position from the first input file is used.')
parser.add_argument('-r', '--analysis_radius', type=float, default=0.125,
help='Aperture for the analysis in degrees.')
parser.add_argument('-m', '--match_rmf', type=str, default=None,
help='RMF filename to which the average PHA file binning is matched.')
parser.add_argument('-d', '--data_dir', type=str, default='',
help='Directory in which the data is located. Will be added as prefix to the entries in the bankfile.')
parser.add_argument('-w', '--write_output', action='store_true', default=False,
help='Write results to FITS files in current directory.')
parser.add_argument('--no_graphical_output', action='store_false', default=True,
help='Switch off graphical output.')
parser.add_argument('-l', '--loglevel', type=str, default='INFO',
help='Amount of logging e.g. DEBUG, INFO, WARNING, ERROR.')
args = parser.parse_args()
args = vars(args)
# Execute script
import logging
logging.basicConfig(level=logging.DEBUG, format='%(levelname)s - %(message)s')
from gammapy.utils.pyfact import create_spectrum
create_spectrum(input_file_names=args['input_file_names'],
analysis_position=args['analysis_position'],
analysis_radius=args['analysis_radius'],
match_rmf=args['match_rmf'],
datadir=args['datadir'],
write_output_files=args['write_output'],
do_graphical_output=args['no_graphical_output'],
loglevel=args['loglevel'])