Verified Commit c740c1d3 authored by Daniel Williams's avatar Daniel Williams
Browse files

Made a change to the production of the projected RA in the log file so that it...

Made a change to the production of the projected RA in the log file so that it is actually correct...
parent 437a0583
Loading
Loading
Loading
Loading
+26 −2
Original line number Diff line number Diff line
"""
88b           d88  88               88                    
888b         d888  ""               88                    
88`8b       d8'88                   88                    
88 `8b     d8' 88  88  8b,dPPYba,   88   ,d8   ,adPPYba,  
88  `8b   d8'  88  88  88P'   `"8a  88 ,a8"   a8P_____88  
88   `8b d8'   88  88  88       88  8888[     8PP"""""""  
88    `888'    88  88  88       88  88`"Yba,  "8b,   ,aa  
88     `8'     88  88  88       88  88   `Y8a  `"Ybbd8"'  
                                                          
--------------------------------------------------------

This file is a part of Minke, a tool for generating simulated
gravitational wave signals, used for characterising and training
search algorithms.

Minke was created by Daniel Williams, based on work started by Chris
Pankow and others, and is built around the LALSimulation library.



"""
from glue.ligolw import ligolw, utils, lsctables
lsctables.use_in(ligolw.LIGOLWContentHandler);
import numpy
@@ -484,6 +506,7 @@ class MDCSet():
        energy = self.egw[row]
        row = self.waveforms[row]
        output = []
        
        output.append(self.name)                  # GravEn_SimID
        output.append(strains[0])                 # SimHrss
        output.append(energy)                     # SimEgwR2
@@ -491,7 +514,8 @@ class MDCSet():
        output.append(np.cos(row.incl))           # Internal_x the cosine of the angle the LOS makes with axis of angular momentum
        output.append(row.phi)                    # Intenal_phi angle between source x-axis and the LOS
        output.append(np.cos(np.pi/2.0 - row.dec)) # cos(External_x) # this needs to be the co-declination
        output.append(np.pi - row.ra)             # External_phi # This is the RA projected onto an Earth-based coordinate system
        output.append(row.ra if row.ra < np.pi else row.ra - 2*np.pi)
        # ^ External_phi # This is the RA projected onto an Earth-based coordinate system
        output.append(row.psi)                    # External_psi # source's polarisation angle
        output.append(frame.start)                # FrameGPS
        output.append(row.time_geocent_gps)       # EarthCtrGPS