-
Notifications
You must be signed in to change notification settings - Fork 0
/
pulsar.py
136 lines (115 loc) · 4.55 KB
/
pulsar.py
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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
"""A virtual pulsar class"""
# Just to be forward-compatible with Python 3
from __future__ import division, print_function, unicode_literals
import astropy.coordinates as C
import astropy.table as T
import astropy.units as U
import numpy as np
# TODO: define your Pulsar class here
# Testing routines -------------------------------------------------------------
import random
import testhelper as TH
@TH.register_test
def test_class():
"""Check if there is a Pulsar class"""
if isinstance(Pulsar, type):
print("OK, Pulsar is a class.")
else:
print("Something went wrong, the type of Pulsar is {!r}".format(type(Pulsar)))
@TH.register_test
def test_instance():
"""Try to create a Pulsar instance"""
# Check if a constructor already exists
if '__init__' in Pulsar.__dict__:
print("Not performing this test because you already defined a constructor")
return
# TODO: write your code here
@TH.register_test
def test_attribute():
"""Try to add a 'name' attribute to a Pulsar instance"""
# Check if a constructor already exists
if '__init__' in Pulsar.__dict__:
print("Not performing this test because you already defined a constructor")
return
if 'name' in Pulsar.__dict__:
print("Not performing this test because you already defined a 'name' attribute")
return
# TODO: create a Pulsar object here
# psr = ...
print("OK, we created a Pulsar instance: {0!r}. Its name is {0.name}".format(psr))
@TH.register_test
def test_class_attribute():
"""Check that the Pulsar class now has a 'name' attribute"""
if hasattr(Pulsar, "name"):
print("OK, Pulsar now has a 'name' attribute. Its default value is {!r}".format(Pulsar.name))
else:
print("Something went wrong, the Pulsar class does not have a 'name' attribute!")
@TH.register_test
def test_prettyprint():
"""Try to call a method"""
# Check if a constructor already exists, don't bother for now
if '__init__' in Pulsar.__dict__:
print("Not performing this test because you already defined a constructor")
return
psr = Pulsar()
psr.name = "Vela"
# Set the attributes by hand
psr.coords = C.SkyCoord.from_name("Vela")
psr.period = 0.089308556629 * U.s
# Pretty-print the Pulsar object
psr.prettyprint()
def make_crab():
""" Create a representation of the Crab pulsar """
name = "Crab"
coords = C.SkyCoord.from_name("Crab pulsar")
period = 0.0333924123 * U.s
return Pulsar(name, coords, period)
@TH.register_test
def test_constructor():
"""Try to use the Pulsar constructor"""
# Make a representation of the Crab pulsar
psr = make_crab()
# Print it
psr.prettyprint()
@TH.register_test
def test_repr():
"""Try to show a textual version of a Pulsar object"""
psr = make_crab()
# Print the textual version
print(psr)
@TH.register_test
def test_property():
"""Try out the properties of the Pulsar class"""
# Make a representation of the Crab pulsar without the period derivative
psr = make_crab()
print("Without period derivative")
psr.prettyprint()
# Set the frequency derivative by hand
psr.frequency_derivative = -3.77535E-10 * U.Hz / U.s
print("With period derivative")
psr.prettyprint()
@TH.register_test
def test_classmethod():
"""Try the from_catalog class method"""
# Fetch a pulsar from the ATNF catalog
psr = Pulsar.from_catalog("J1944+2236")
psr.prettyprint()
@TH.register_test
def test_full():
"""Perform a full test"""
# Print a random selection of 20 pulsars from the ATNF catalog
Pulsar.read_catalog()
print("Name".center(12), "RA".center(11), "DEC".center(12), "Freq".center(12), "Fdot".center(12), sep=" | ")
print("".center(12), "h:m:s".center(11), "d:m:s".center(12), "Hz".center(12), "Hz/s".center(12), sep=" | ")
print("-" * 12, "-" * 11, "-" * 12, "-" * 12, "-" * 12, sep="-+-")
for prow in random.sample(Pulsar.catalog, 20): # Choose 20 pulsars to show
psr = Pulsar.from_catalog_row(prow)
name = psr.name
ra = psr.coords.ra.to_string(unit=U.hourangle, sep=':', precision=2, pad=True)
dec = psr.coords.dec.to_string(unit=U.degree, sep=':', precision=2, alwayssign=True, pad=True)
freq = psr.frequency.to(U.Hz).value
fdot = psr.frequency_derivative.to(U.Hz/U.s).value
print("{:12s} | {:11s} | {:12s} | {:12.5e} | {:+12.5e}".format(name, ra, dec, freq, fdot))
# Main entry point: list available tests and ask which one to run --------------
if __name__ == '__main__':
TH.test_main()