GalacticBinary.h

Go to the documentation of this file.

/*
 *  Copyright (C) 2019 Tyson B. Littenberg (MSFC-ST12), Neil J. Cornish
 *
 *  This program is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation; either version 2 of the License, or
 *  (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with with program; see the file COPYING. If not, write to the
 *  Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston,
 *  MA  02111-1307  USA
 */
 
 
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
 
#define MAXSTRINGSIZE 1024 
 
struct Data
{
    int N;        
    int Nmax;     
    int NT;       
    int Nchannel; 
    int DMAX;     
    double logN;  
 
    long cseed; 
    long nseed; 
    long iseed; 
 
    double T; 
    double sqT; 
    double *t0;   
    double *tgap; 
 
    
    int qmin; 
    int qmax; 
    int qpad; 
    
    double fmin; 
    double fmax; 
    double sine_f_on_fstar; 
 
    //some manipulations of f,fmin for likelihood calculation
    double sum_log_f; 
    double logfmin; 
    double logfmax; 
 
    struct TDI **tdi; 
    struct TDI **raw; 
    struct Noise **noise; 
    char format[16];
    char dataDir[MAXSTRINGSIZE]; 
 
    //Spectrum proposal
    double *p; 
    double pmax; 
    double SNR2; 
 
    int Nwave; 
    int downsample; 
 
    double ****h_rec; 
    double ****h_res; 
    double ****r_pow; 
    double ****h_pow; 
    double ****S_pow; 
    char fileName[MAXSTRINGSIZE]; 
 
    int NP; 
    struct Source *inj; 
 
    
    int Ncache; 
    char **cache; 
    struct Catalog *catalog; 
};
 
struct Flags
{
    int verbose;    
    int quiet;      
    int NMCMC;      
    int NBURN;      
    int NINJ;       
    int NDATA;      
    int NT;         
    int NVB;        
    int DMAX;       
    int simNoise;   
    int fixSky;     
    int fixFdot;
    int fixFreq;    
    int galaxyPrior;
    int snrPrior;   
    int emPrior;    
    int knownSource;
    int detached;   
    int strainData; 
    int hdf5Data;   
    int orbit;      
    int prior;      
    int debug;      
    int cheat;      
    int burnin;     
    int maximize;   
    int update;     
    int updateCov;  
    int match;      
    int rj;         
    int gap;        
    int calibration;
    int confNoise;  
    int resume;     
    int catalog;    
    int threads;    
    int psd;        
 
    
    char runDir[MAXSTRINGSIZE];       
    char vbFile[MAXSTRINGSIZE];       
    char **injFile;                   
    char noiseFile[MAXSTRINGSIZE];    
    char cdfFile[MAXSTRINGSIZE];      
    char gmmFile[MAXSTRINGSIZE];      
    char covFile[MAXSTRINGSIZE];      
    char matchInfile1[MAXSTRINGSIZE]; 
    char matchInfile2[MAXSTRINGSIZE]; 
    char pdfFile[MAXSTRINGSIZE];      
    char psdFile[MAXSTRINGSIZE];      
    char catalogFile[MAXSTRINGSIZE];  
};
 
struct Chain
{
    int NC;
    
    int NP;
 
    int *index;
    
    int **dimension;
    
    double *acceptance;
    
    double *temperature;
    
    double *avgLogL;
    
    double annealing;
    
    double logLmax;
    
    const gsl_rng_type **T;
    
    gsl_rng **r;
    
    
    FILE **noiseFile;
    FILE **foregroundFile;
    
    FILE **chainFile;
    
    FILE **calibrationFile;
    
    FILE **dimensionFile;
    
    FILE **parameterFile;
    
    FILE *likelihoodFile;
    
    FILE *temperatureFile;
    
    char chainDir[MAXSTRINGSIZE]; 
    char chkptDir[MAXSTRINGSIZE]; 
 
};
 
struct Source
{
    int NP;
 
    double *params;
 
    struct TDI *tdi;
    
 
    double f0;    
    double dfdt;  
    double d2fdt2;
    double amp;   
    
    double psi;      
    double cosi;     
    double phi0;     
    double phi;      
    double costheta; 
    
    double m1;  
    double m2;  
    double Mc;  
    double D;   
    
    int BW;   
    int qmin; 
    int qmax; 
    int imin; 
    int imax; 
 
    double **fisher_matrix; 
    double **fisher_evectr; 
    double *fisher_evalue;  
    int fisher_update_flag; 
};
 
struct Noise
{
    int N;
    
    int Nchannel;
    
    double *eta;
    
    double *f;
    
    double ***C;    
    double ***invC; 
    double *detC;   
    
    double *transfer;
};
 
struct Calibration
{
    double dampA;
    double dampE;
    double dampX;
    double dampY;
    double dampZ;
    
    double dphiA;
    double dphiE;
    double dphiX;
    double dphiY;
    double dphiZ;
 
    double real_dphiA;
    double real_dphiE;
    double real_dphiX;
    double real_dphiY;
    double real_dphiZ;
    double imag_dphiA;
    double imag_dphiE;
    double imag_dphiX;
    double imag_dphiY;
    double imag_dphiZ;
 
};
 
struct Model
{
    int NT;     
    int NP;     
    int Nmax;   
    int Nlive;  
    struct Source **source; 
    
    struct Noise **noise;
    
    struct Calibration **calibration;
    
    struct TDI **tdi; 
    struct TDI **residual; 
    
    double *t0; 
    double *t0_min; 
    double *t0_max; 
    
    double **prior; 
    double *logPriorVolume; 
    
    double logL; 
    double logLnorm; 
};