def phaseogram(self, weights=None, bins=100, rotate=0.0, size=5, alpha=0.25, file=False): """ Make a nice 2-panel phaseogram for the current model """ mjds = self.toas.table['tdbld'].astype(np.float64) phss = self.get_event_phases() fermi.phaseogram(mjds, phss, weights=self.weights, bins=bins, rotate=rotate, size=size, alpha=alpha, file=file)
def phaseogram(self, bins=100, rotate=0.0, size=5, alpha=0.25, file=False): """ Make a nice 2-panel phaseogram for the current model """ mjds = self.toas.table['tdbld'].astype(np.float64) phss = self.get_event_phases() fermi.phaseogram(mjds, phss, weights=self.weights, bins=bins, rotate=rotate, size=size, alpha=alpha, file=file)
if tt.mjd < maxT: tlnew.append(tt) tl = tlnew print("post len : ", len(tlnew)) # Now convert to TOAs object and compute TDBs and posvels ts = toa.TOAs(toalist=tl) ts.filename = args.eventfile ts.compute_TDBs() ts.compute_posvels(ephem=args.ephem, planets=args.planets) print(ts.get_summary()) mjds = ts.get_mjds() print(mjds.min(), mjds.max()) # Read in model modelin = pint.models.get_model(args.parfile) # Remove the dispersion delay as it is unnecessary modelin.delay_funcs['L1'].remove(modelin.dispersion_delay) # Compute model phase for each TOA phss = modelin.phase(ts.table)[1] # ensure all postive phases = np.where(phss < 0.0, phss + 1.0, phss) mjds = ts.get_mjds() weights = np.array([w['weight'] for w in ts.table['flags']]) h = float(hmw(phases, weights)) print("Htest : {0:.2f} ({1:.2f} sigma)".format(h, h2sig(h))) phaseogram(mjds, phases, weights, bins=100, file=args.outfile)
if tt.mjd < maxT: tlnew.append(tt) tl=tlnew print("post len : ",len(tlnew)) # Now convert to TOAs object and compute TDBs and posvels ts = toa.TOAs(toalist=tl) ts.filename = args.eventfile ts.compute_TDBs() ts.compute_posvels(ephem=args.ephem,planets=args.planets) print(ts.get_summary()) mjds = ts.get_mjds() print(mjds.min(),mjds.max()) # Read in model modelin = pint.models.get_model(args.parfile) # Remove the dispersion delay as it is unnecessary modelin.delay_funcs['L1'].remove(modelin.dispersion_delay) # Compute model phase for each TOA phss = modelin.phase(ts.table)[1] # ensure all postive phases = np.where(phss < 0.0, phss + 1.0, phss) mjds = ts.get_mjds() weights = np.array([w['weight'] for w in ts.table['flags']]) h = float(hmw(phases,weights)) print("Htest : {0:.2f} ({1:.2f} sigma)".format(h,h2sig(h))) phaseogram(mjds,phases,weights,bins=100,file = args.outfile)