else:
value.append(parseerror(*globals()[k](m)))
return value
m = model('traditional/new.par')
value1 = addmodeldata(m)
del m
m = model('1713.sns.par')
value2 = addmodeldata(m)
del m
#m = model('Feb.T2.RN.par')
#value3 = addmodeldata(m)
#del m
#data = [parameter, value1, value2, value3]
data = [parameter, value1, value2]#, value3]
comments=[
#'Numbers in parentheses indicate the uncertainties on the last digit(s). Uncertainties on parameters are estimated from a combination of 23 indepedent MCMC chains randomly generated using the Metropolis algorithm. Each chain run for 100000 points after 10000 burn-in steps.',
#"We used {\it tempo2}'s {\it T2} binary model, which implicitly account for the change of the projected semi-major axis ($\dot{x}$ and allow us to fit for the position angle of ascending node ($\Omega$). Numbers in parentheses indicate the uncertainties on the last digit(s). Uncertainties on parameters are estimated by the {\it tempo2} program using information in the covariance matrix.",
#' The averaged DM value; See Section 3.2 and Figure 2 for the more discussion.',
#' See Figure 2 of \citealt{sns+05} for definition.'
]
#print len(parameter), len(value)
table = deluxetable(Caption=Caption, colsetting='lcc', colnames = colnames, data=data, label="tab:sns05", comments=comments, fontsize='scriptsize')
print table
except:
value.append(SF(globals()[k](m)))
data = [parameter, value]
parameter.append(r"\textit{Derived Parameters}")
value.append("")
for k in Derived:
parameter.append(Derived[k])
try:
value.append(parseerror(*m.__dict__[k]))
except:
value.append(SF(globals()[k](m)))
data = [parameter, value]
table = deluxetable(
Caption=Caption,
colsetting="lc",
colnames=colnames,
data=data,
label="tab:par",
comments=[
"Numbers in parentheses indicate the uncertainties on the last disgit(s). Uncertainties on parameters are estimated using MCMC simulation. ",
"The averaged DM value; See Section \label{sec:dmx} and Figure \label{fig:dmx} for the more discussion.",
"See Figure 2 of \citealt{sns+05} for definition.",
],
)
print table
try:
value.append(parseerror(*m.__dict__[k]))
except:
value.append(SF(globals()[k](m)))
return value
m = model('normaltempo.par')
value1 = addmodeldata(m)
del m
m = model('glstempo.par')
value2 = addmodeldata(m)
del m
m = model('RNtempo.par')
value3 = addmodeldata(m)
del m
data = [parameter, value1, value2, value3]
comments=[
#'Numbers in parentheses indicate the uncertainties on the last disgit(s). Uncertainties on parameters are estimated from a combination of 23 indepedent MCMC chains randomly generated using the Metropolis algorithm. Each chain run for 100000 points after 10000 burn-in steps.',
'Numbers in parentheses indicate the uncertainties on the last disgit(s). Uncertainties on parameters are estimated by the {\it tempo} program using information in the covariance matrix.',
' The averaged DM value; See Section 3.2 and Figure 2 for the more discussion.',
' See Figure 2 of \citealt{sns+05} for definition.']
#print len(parameter), len(value)
table = deluxetable(Caption=Caption, colsetting='lccc', colnames = colnames, data=data, label="tab:par", comments=comments)
print table
for k in Fixed:
parameter.append(Fixed[k])
try:
value.append(parseerror(*m.__dict__[k]))
except:
if type(m.__dict__[k]) == type(''):
value.append(m.__dict__[k])
elif type(m.__dict__[k]) == type(m.__dict__['PEPOCH']):
value.append(m.__dict__[k].quantize(50000))
else:
value.append(SF(globals()[k](m)))
data = [parameter, value]
parameter.append(r'\textit{Derived Parameters}')
value.append('')
for k in Derived:
parameter.append(Derived[k])
try:
value.append(parseerror(*m.__dict__[k]))
except:
value.append(SF(globals()[k](m)))
data = [parameter, value]
table = deluxetable(Caption=Caption, colsetting='lc', colnames = colnames, data=data, label="tab:par", comments=['Numbers in parentheses indicate the uncertainties on the last digit(s). Uncertainties on parameters are estimated from a general least square fit using {\sc tempo}.', 'The averaged DM value; See Section 3.2 and Figure 2 for the more discussion.', 'See Figure 2 of \citealt{sns+05} for definition.'])
print table
maxd = max(mjds)
mind = min(mjds)
return '%s$-$%s' % (MJD_to_datetime(mind).strftime('%Y %b'),
MJD_to_datetime(maxd).strftime('%Y %b'))
for k in keys:
krange = findrange(k)
#print k, krange
DateRange.append(krange)
NoOfEpoch = []
for k in keys:
NoOfEpoch.append(len(m.averes[k]))
#print NoOfEpoch
Caption = '21 year J1713+0747 observations.'
colnames = [
'System', 'Dates', 'Number of', 'Epochs', 'Bandwidth',
'Typical Integration'
]
data = [keys, DateRange, NoOfTOA, NoOfEpoch, BandWidth, Integration]
table = deluxetable(Caption=Caption,
colnames=colnames,
data=data,
label='tab:obs',
colsetting='lccccc')
print table