Skip to content

oLIB and MLy VOEvents do not have central_frequency and duration values

The content of VO alert voevents are missing p_central_freq p_duration for MLy, and p_duration for oLIB. Fir oLIB is still present the "fluency" that the collaboration has decided not ti distribute.

Related issue: emfollow/gwcelery#594 (closed)

The code should be changed (my guess to):

        elif isinstance(event, LalInferenceBurstEvent):
            p_freq = vp.Param(
                "frequency",
                value=float(event.frequency_mean),
                ucd="gw.frequency",
                unit="Hz",
                ac=True,
            )
            p_freq.Description = "Mean frequency of GW burst signal"
            v.What.append(p_freq)

            # Calculate the fluence. 
            # From Min-A Cho: fluence = pi*(c**3)*(freq**2)*(hrss_max**2)*(10**3)/(4*G)
            # Note that hrss here actually has units of s^(-1/2)
            # XXX obviously need to refactor here.
            # try:
            #     fluence = pi * pow(c,3) * pow(event.frequency,2) 
            #     fluence = fluence * pow(event.hrss,2)
            #     fluence = fluence / (4.0*G)
            # 
            #     p_fluence = vp.Param(
            #         "Fluence",
            #         value=fluence,
            #         ucd="gw.fluence",
            #        unit="erg/cm^2",
            #         ac=True
            #     )
            #     p_fluence.Description = "Estimated fluence of GW burst signal"
            #     v.What.append(p_fluence)
            ### Duration
            p_duration = vp.Param(
                "Duration",
                value=float(event.duration),
                unit="s",
                ucd="time.duration",
                ac=True,
            )
            p_duration.Description = "Measured duration of GW burst signal"
            v.What.append(p_duration)
        elif isinstance(event, MLyBurstEvent):
            p_central_freq = vp.Param(
                "CentralFreq",
                value=float(event.central_freq),
                ucd="gw.frequency",
                unit="Hz",
                ac=True,
            )
            p_central_freq.Description = \
                "Central frequency of GW burst signal"
            v.What.append(p_central_freq)

            ### Duration
            duration = event.quality_mean / (2 * np.pi * event.frequency_mean)
            p_duration = vp.Param(
                "Duration",
                value=float(duration),
                unit="s",
                ucd="time.duration",
                ac=True,
            )
            p_duration.Description = "Measured duration of GW burst signal"
            v.What.append(p_duration)
        except Exception as e:
            logger.exception(e)
Edited by Roberto DePietri
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information