示例#1
0
def main():
    pfdfn = sys.argv[1]

    cand = candidate.read_pfd_file(pfdfn)

    print "Loaded %s" % cand.pfdfn
    print "    Best topo period (s): %f" % cand.topo_period
    print "    Best bary period (s): %f" % cand.bary_period
    print "    Best DM (cm^-3/pc): %f" % cand.dm
    print "    RA (J2000 - deg): %f" % cand.raj_deg
    print "    Dec (J2000 - deg): %f" % cand.decj_deg

    #print "-"*10
    #pprint.pprint(cand.__dict__)
    #print "-"*10

    mgauss.add_data(cand)

    print "Added multi-gaussian fit to cand"
    #pprint.pprint(cand.__dict__)
    #print "-"*10

    print cand.multigaussfit

    data = cand.profile.copy()
    data /= np.sqrt(cand.pfd.varprof)
    data -= data.mean()
    cand.multigaussfit.plot_comparison(data, True)
示例#2
0
def main():
    pfdfn = sys.argv[1]

    cand = candidate.read_pfd_file(pfdfn)

    print "Loaded %s" % cand.pfdfn
    print "    Best topo period (s): %f" % cand.topo_period
    print "    Best bary period (s): %f" % cand.bary_period
    print "    Best DM (cm^-3/pc): %f" % cand.dm
    print "    RA (J2000 - deg): %f" % cand.raj_deg
    print "    Dec (J2000 - deg): %f" % cand.decj_deg

    #print "-"*10
    #pprint.pprint(cand.__dict__)
    #print "-"*10

    mgauss.add_data(cand)

    print "Added multi-gaussian fit to cand"
    #pprint.pprint(cand.__dict__)
    #print "-"*10

    print cand.multigaussfit

    data = cand.profile.copy()
    data /= np.sqrt(cand.pfd.varprof)
    data -= data.mean()
    cand.multigaussfit.plot_comparison(data, True)
示例#3
0
def rate_pfd(pfdfn, rater_instances):
    """Given the name of a *.pfd file and a list of Rater instances
        compute the ratings.

        Inputs:
            pfdfn: Name of the *.pfd file.
            rater_instances: A list of Rater instances to compute ratings.

        Outputs:
            cand: The resulting (rated) Candidate object.

        ***NOTE: RatingValues are added directly to the Candidate object.
    """
    cand = candidate.read_pfd_file(pfdfn)
    for rater in rater_instances:
        ratval = rater.rate(cand)
        cand.add_rating(ratval)
    return cand
示例#4
0
def rate_pfd(pfdfn, rater_instances):
    """Given the name of a *.pfd file and a list of Rater instances
        compute the ratings.

        Inputs:
            pfdfn: Name of the *.pfd file.
            rater_instances: A list of Rater instances to compute ratings.

        Outputs:
            cand: The resulting (rated) Candidate object.

        ***NOTE: RatingValues are added directly to the Candidate object.
    """
    cand = candidate.read_pfd_file(pfdfn)
    for rater in rater_instances:
        ratval = rater.rate(cand)
        cand.add_rating(ratval)
    return cand