pSEOBNRv4HM test of GR branch

Description

This is a pipeline to do a test of GR using gravitational wave signals from binary black hole systems, where we treat the (complex) QNM frequencies as free parameters and estimate them directly from the data. The code involves changes to the SEOBNRv4HM model, where the QNM parameters are sampled over as non-GR parameters, in addition to the GR parameters. It also introduces changes to the parameter estimation code to recognise those free parameters, sample over them and plot them. We further make changes to the injection infrastructure to allow us to simulate signals with free QNM parameters.

API Changes and Justification

Backwards Compatible Changes

• This change introduces no API changes
• This change adds new API calls

Backwards Incompatible Changes

• This change modifies an existing API
• This change removes an existing API

If any of the Backwards Incompatible check boxes are ticked please provide a justification why this change is necessary and why it needs to be done in a backwards incompatible way.

Review Status

This pipeline has been reviewed and used to produce results for the pSEOBNRv4HM ringdown test of GR that has gone into the O3a TGR paper. The review wiki is here: https://git.ligo.org/cbc-testinggr/reviews/nohairtesteobhm/-/wikis/home

and specifically, the code review page is here: https://git.ligo.org/cbc-testinggr/reviews/nohairtesteobhm/-/wikis/Code-review_v4

lalapps/src/inspiral/inspinj.c

@@ -543,6 +543,9 @@ LALSimInspiralApplyTaper taperInj = LAL_SIM_INSPIRAL_TAPER_NONE;
 AlignmentType alignInj = notAligned;
 REAL8 redshift;
 
+REAL8 domega220=0.0;
+REAL8 dtau220=0.0;
+
 REAL8 single_IFO_SNR_threshold=0.0;
 char ** ifonames=NULL;
 int numifos=0;
@@ -2022,6 +2025,8 @@ int main( int argc, char *argv[] )
     {"stdev-spin2",             required_argument, 0,                 1004},
     {"mean-spin1",              required_argument, 0,                 1005},
     {"mean-spin2",              required_argument, 0,                 1006},
+    {"domega220",                   required_argument, 0,                 1100},
+    {"dtau220",                   required_argument, 0,                 1101},
     {0, 0, 0, 0}
   };
   int c;
@@ -3186,6 +3191,19 @@ int main( int argc, char *argv[] )
         next_process_param( long_options[option_index].name,
           "float", "%le", meanSpin2 );
         break;
+      case 1100:
+        domega220 = atof( LALoptarg );
+        this_proc_param = this_proc_param->next =
+        next_process_param( long_options[option_index].name,
+                            "float", "%le", domega220 );
+        break;
+      case 1101:
+        dtau220 = atof( LALoptarg );
+        this_proc_param = this_proc_param->next =
+        next_process_param( long_options[option_index].name,
+                            "float", "%le", dtau220 );
+        break;
+
 
 
       default:
@@ -4428,6 +4446,9 @@ int main( int argc, char *argv[] )
     /* populate the bandpass options */
     simTable->bandpass = bandPassInj;
 
+    simTable->domega220=domega220;
+    simTable->dtau220=dtau220;
+
 
     /* increment current time, avoiding roundoff error;
        check if end of loop is reached */