#include "cwb.hh"#include "config.hh"#include "wavearray.hh"#include "TString.h"#include "TObjArray.h"#include "TObjString.h"#include "TRandom.h"#include "TComplex.h"#include "TGraphAsymmErrors.h"#include "TPaveStats.h"#include "TMath.h"#include "mdc.hh"#include "cwb2G.hh"#include "watplot.hh"#include "gwavearray.hh"#include "network.hh"#include "gnetwork.hh"#include "constants.hh"#include <fstream>#include <vector>#include "TKDE.h"#include "TSpline.h"Macros | |
| #define | _USE_LAL |
| #define | CCUT_BCHIRP 1.35 |
| #define | CCUT_FMAX 512.0 |
| #define | CCUT_LTIME 0.5 |
| #define | CCUT_RTIME 0.0 |
| #define | CCUT_UCHIRP 1.65 |
| #define | CCUT_WDM_FRES 16 |
| #define | LOGL_ARTHR 0.001 |
| #define | LOGL_ENABLED false |
| #define | LOGL_FHIGH 8192 |
| #define | LOGL_FLOW 0 |
| #define | LOGL_IFO_MASK 0x1111111 |
| #define | LOGL_RESAMPLE 1 |
| #define | MAX_ECOV 1.0 |
| #define | MAX_TRIALS 30 |
| #define | PE_CCUT "" |
| #define | PE_COVERAGE "local" |
| #define | PE_DATA_LABEL "wave_" |
| #define | PE_DUMP_CR true |
| #define | PE_EFRACTION 0.01 |
| #define | PE_GPS_ERROR 0.01 |
| #define | PE_GPS_EVENT 0 |
| #define | PE_GW_NAME "GWXXYYZZ" |
| #define | PE_IFAR_THR 0 |
| #define | PE_INJ_ONSRC false |
| #define | PE_INJ_TIME_STEP 150 |
| #define | PE_LOGL "" |
| #define | PE_MAKE_FF true |
| #define | PE_MAKE_NE false |
| #define | PE_MAKE_OF false |
| #define | PE_MAKE_ONSRC false |
| #define | PE_MAKE_RE true |
| #define | PE_MAX_EVENTS 15000 |
| #define | PE_NCYCLES 800 |
| #define | PE_NIFO 2 |
| #define | PE_OFF_EVENT 0 |
| #define | PE_PLOT_99PERC false |
| #define | PE_PLOT_CR true |
| #define | PE_PLOT_EFRACTION 0.999 |
| #define | PE_PLOT_FNAME "pewave" |
| #define | PE_PLOT_ODIR "" |
| #define | PE_PLOT_ONSRC false |
| #define | PE_PLOT_TYPE "time" |
| #define | PE_PLOT_ZOOM 1.0 |
| #define | PE_RCUT "" |
| #define | PE_RESIDUALS false |
| #define | PE_RSTAT "" |
| #define | PE_SENSITIVITY false |
| #define | PE_SIDE "full" |
| #define | PE_STATISTICS false |
| #define | PE_SYNC_PHASE false |
| #define | PE_SYNC_TIME false |
| #define | PE_TCOA_SHIFT 0. |
| #define | PE_USE_ONSRC "" |
| #define | PE_WAVE_DIR "output" |
| #define | RCUT_THR 0.95 |
| #define | RCUT_WDM_FRES 8 |
| #define | RSTAT_ENABLED false |
| #define | RSTAT_JTRIALS 0 |
| #define | RSTAT_RTRIALS 0 |
| #define | RSTAT_TYPE "median" |
| #define | XIFO 4 |
Functions | |
| void | ComputeAmpSensitivity (TString match, double lower, double upper, double &dmatch, double &damp) |
| int | ComputeCCutStatistics (int nIFO) |
| void | ComputePhaseSensitivity (TString match, double lower, double upper, double &dmatch, double &dphase) |
| int | ComputeStatistics (int nIFO) |
| void | ComputeTimeSensitivity (TString match, double lower, double upper, double &dmatch, double &dtime) |
| int | ComputeWaveformCR (int nIFO) |
| void | ComputeWaveformCR (int nIFO, int selected, double *cov50, double *cov90, double *cov99) |
| int | ComputeWaveformFCR (int nIFO) |
| void | cwb_pereport (TString options) |
| void | DumpStatistics (TString ofname, TString params, bool app) |
| void | DumpWaveformCR (int nIFO) |
| double | GetCCut (wavearray< double > *ts, double tcoa, double m1, double m2, double s1z, double s2z) |
| void | GetCCutParms (TString options) |
| double | GetDelay (network *NET, TString ifo, double theta, double phi) |
| void | GetFullCoverage (int nIFO, int selected, double *fcov50, double *fcov90, double *fcov99) |
| double | GetInjTcoa (double geocentric_tcoa, network *NET, TString ifo, double theta, double phi) |
| void | GetLoglParms (TString options) |
| void | GetOffSourceStatistic (int nIFO, int recId, vector< int > injId, vector< double > &oFF, vector< double > &oOF, vector< double > &oRE) |
| int | GetOnSourceStatistic (int nIFO, double &mFF, double &mOF, double &mRE) |
| TString | GetParameterFromInjLog (TString log, TString param) |
| wavearray< double > | GetRCut (wavearray< double > *ts, wavearray< double > *aux) |
| void | GetRCutParms (TString options) |
| void | GetRstatParms (TString options) |
| double | GetSyncTcoa (wavearray< double > wf, double sync_time, double sync_phase, double tcoa) |
| void | LoadOutputRootFiles (int nIFO, TString data_label, TString output_dir, int max_events=PE_MAX_EVENTS, double ifar_thr=0) |
| void | MakeDistributions (TString type="FF") |
| void | MakePvalueDistribution (TString type="FF") |
| void | PlotWaveformAsymmErrors (TString ofname, TString title, wavearray< double > wrec, wavearray< double > wmed, wavearray< double > wl50, wavearray< double > wu50, wavearray< double > wl90, wavearray< double > wu90, wavearray< double > wl99, wavearray< double > wu99, wavearray< double > wref, TString pdir, double P, double Q, double tref, double tcoa, int ifoid) |
| void | PlotWaveforms (int nIFO, int id, bool residual=false) |
| void | PlotWaveformsCR (int nIFO) |
| void | PrintUserOptions (TString ofname="") |
| void | ReadUserOptions (TString options) |
| void | ResetUserOptions () |
| int | SyncWaveforms (int nIFO, TString type) |
| float | TestStatistic (int nIFO, int id, float ref) |
Variables | |
| int | gEVENTS |
| static TString | gIFO [3] = {"L1","H1","V1"} |
| uoptions | gOPT |
| WDM< double > * | gWDM |
| std::vector< wavearray< double > > | pFRQ |
| std::vector< wavearray< double > > | pSNR |
| std::vector< wavearray< double > > | pTIM |
| std::vector< wavearray< double > > | vAUX |
| float | vDOF |
| double | vFACTOR |
| std::vector< double > | vFF |
| std::vector< wavearray< double > > | vINJ |
| std::vector< wavearray< double > > | vL50 |
| std::vector< wavearray< double > > | vL90 |
| std::vector< wavearray< double > > | vL99 |
| float | vM1 |
| float | vM2 |
| std::vector< wavearray< double > > | vMED |
| std::vector< double > | vNE |
| std::vector< wavearray< double > > | vNUL |
| std::vector< double > | vOF |
| float | vPHI |
| double | vPSYNC [NIFO_MAX] |
| std::vector< double > | vRE |
| std::vector< wavearray< double > > | vREC |
| float | vS1 |
| float | vS2 |
| float | vSNR [NIFO_MAX] |
| double | vTCOA |
| float | vTHETA |
| std::vector< double > | vTREC |
| double | vTSYNC [NIFO_MAX] |
| std::vector< wavearray< double > > | vU50 |
| std::vector< wavearray< double > > | vU90 |
| std::vector< wavearray< double > > | vU99 |
| std::vector< wavearray< double > > | wAUX [PE_MAX_EVENTS] |
| float | wDOF [PE_MAX_EVENTS] |
| double | wFACTOR [PE_MAX_EVENTS] |
| std::vector< wavearray< double > > | wINJ [PE_MAX_EVENTS] |
| float | wM1 [PE_MAX_EVENTS] |
| float | wM2 [PE_MAX_EVENTS] |
| std::vector< wavearray< double > > | wNUL [PE_MAX_EVENTS] |
| float | wPHI [PE_MAX_EVENTS] |
| std::vector< wavearray< double > > | wREC [PE_MAX_EVENTS] |
| float | wS1 [PE_MAX_EVENTS] |
| float | wS2 [PE_MAX_EVENTS] |
| float | wSNR [PE_MAX_EVENTS][NIFO_MAX] |
| double | wTCOA [PE_MAX_EVENTS] |
| float | wTHETA [PE_MAX_EVENTS] |
| std::vector< double > | xFF |
| std::vector< double > | xNE |
| std::vector< double > | xOF |
| std::vector< double > | xRE |
| #define _USE_LAL |
Definition at line 27 of file cwb_pereport.C.
| #define CCUT_BCHIRP 1.35 |
Definition at line 258 of file cwb_pereport.C.
| #define CCUT_FMAX 512.0 |
Definition at line 262 of file cwb_pereport.C.
| #define CCUT_LTIME 0.5 |
Definition at line 260 of file cwb_pereport.C.
| #define CCUT_RTIME 0.0 |
Definition at line 261 of file cwb_pereport.C.
| #define CCUT_UCHIRP 1.65 |
Definition at line 259 of file cwb_pereport.C.
| #define CCUT_WDM_FRES 16 |
Definition at line 257 of file cwb_pereport.C.
| #define LOGL_ARTHR 0.001 |
Definition at line 272 of file cwb_pereport.C.
| #define LOGL_ENABLED false |
Definition at line 269 of file cwb_pereport.C.
| #define LOGL_FHIGH 8192 |
Definition at line 271 of file cwb_pereport.C.
| #define LOGL_FLOW 0 |
Definition at line 270 of file cwb_pereport.C.
| #define LOGL_IFO_MASK 0x1111111 |
Definition at line 273 of file cwb_pereport.C.
| #define LOGL_RESAMPLE 1 |
Definition at line 274 of file cwb_pereport.C.
| #define MAX_ECOV 1.0 |
| #define MAX_TRIALS 30 |
| #define PE_CCUT "" |
Definition at line 101 of file cwb_pereport.C.
| #define PE_COVERAGE "local" |
Definition at line 122 of file cwb_pereport.C.
| #define PE_DATA_LABEL "wave_" |
Definition at line 89 of file cwb_pereport.C.
| #define PE_DUMP_CR true |
Definition at line 108 of file cwb_pereport.C.
| #define PE_EFRACTION 0.01 |
Definition at line 119 of file cwb_pereport.C.
| #define PE_GPS_ERROR 0.01 |
Definition at line 96 of file cwb_pereport.C.
| #define PE_GPS_EVENT 0 |
Definition at line 94 of file cwb_pereport.C.
| #define PE_GW_NAME "GWXXYYZZ" |
Definition at line 88 of file cwb_pereport.C.
| #define PE_IFAR_THR 0 |
Definition at line 93 of file cwb_pereport.C.
| #define PE_INJ_ONSRC false |
Definition at line 99 of file cwb_pereport.C.
| #define PE_INJ_TIME_STEP 150 |
Definition at line 98 of file cwb_pereport.C.
| #define PE_LOGL "" |
Definition at line 104 of file cwb_pereport.C.
| #define PE_MAKE_FF true |
Definition at line 132 of file cwb_pereport.C.
| #define PE_MAKE_NE false |
Definition at line 135 of file cwb_pereport.C.
| #define PE_MAKE_OF false |
Definition at line 133 of file cwb_pereport.C.
| #define PE_MAKE_ONSRC false |
Definition at line 136 of file cwb_pereport.C.
| #define PE_MAKE_RE true |
Definition at line 134 of file cwb_pereport.C.
| #define PE_MAX_EVENTS 15000 |
Definition at line 86 of file cwb_pereport.C.
| #define PE_NCYCLES 800 |
Definition at line 92 of file cwb_pereport.C.
| #define PE_NIFO 2 |
Definition at line 91 of file cwb_pereport.C.
| #define PE_OFF_EVENT 0 |
Definition at line 97 of file cwb_pereport.C.
| #define PE_PLOT_99PERC false |
Definition at line 130 of file cwb_pereport.C.
| #define PE_PLOT_CR true |
Definition at line 124 of file cwb_pereport.C.
| #define PE_PLOT_EFRACTION 0.999 |
Definition at line 128 of file cwb_pereport.C.
| #define PE_PLOT_FNAME "pewave" |
Definition at line 126 of file cwb_pereport.C.
| #define PE_PLOT_ODIR "" |
Definition at line 127 of file cwb_pereport.C.
| #define PE_PLOT_ONSRC false |
Definition at line 111 of file cwb_pereport.C.
| #define PE_PLOT_TYPE "time" |
Definition at line 125 of file cwb_pereport.C.
| #define PE_PLOT_ZOOM 1.0 |
Definition at line 129 of file cwb_pereport.C.
| #define PE_RCUT "" |
Definition at line 102 of file cwb_pereport.C.
| #define PE_RESIDUALS false |
Definition at line 117 of file cwb_pereport.C.
| #define PE_RSTAT "" |
Definition at line 106 of file cwb_pereport.C.
| #define PE_SENSITIVITY false |
Definition at line 116 of file cwb_pereport.C.
| #define PE_SIDE "full" |
Definition at line 121 of file cwb_pereport.C.
| #define PE_STATISTICS false |
Definition at line 118 of file cwb_pereport.C.
| #define PE_SYNC_PHASE false |
Definition at line 113 of file cwb_pereport.C.
| #define PE_SYNC_TIME false |
Definition at line 114 of file cwb_pereport.C.
| #define PE_TCOA_SHIFT 0. |
Definition at line 95 of file cwb_pereport.C.
| #define PE_USE_ONSRC "" |
Definition at line 109 of file cwb_pereport.C.
| #define PE_WAVE_DIR "output" |
Definition at line 90 of file cwb_pereport.C.
| #define RCUT_THR 0.95 |
Definition at line 265 of file cwb_pereport.C.
| #define RCUT_WDM_FRES 8 |
Definition at line 264 of file cwb_pereport.C.
| #define RSTAT_ENABLED false |
Definition at line 281 of file cwb_pereport.C.
| #define RSTAT_JTRIALS 0 |
Definition at line 284 of file cwb_pereport.C.
| #define RSTAT_RTRIALS 0 |
Definition at line 283 of file cwb_pereport.C.
| #define RSTAT_TYPE "median" |
Definition at line 282 of file cwb_pereport.C.
| #define XIFO 4 |
Definition at line 23 of file cwb_pereport.C.
| void ComputeAmpSensitivity | ( | TString | match, |
| double | lower, | ||
| double | upper, | ||
| double & | dmatch, | ||
| double & | damp | ||
| ) |
Definition at line 2783 of file cwb_pereport.C.
Definition at line 786 of file cwb_pereport.C.
| void ComputePhaseSensitivity | ( | TString | match, |
| double | lower, | ||
| double | upper, | ||
| double & | dmatch, | ||
| double & | dphase | ||
| ) |
Definition at line 2718 of file cwb_pereport.C.
Definition at line 496 of file cwb_pereport.C.
| void ComputeTimeSensitivity | ( | TString | match, |
| double | lower, | ||
| double | upper, | ||
| double & | dmatch, | ||
| double & | dtime | ||
| ) |
Definition at line 2749 of file cwb_pereport.C.
Definition at line 946 of file cwb_pereport.C.
Definition at line 4169 of file cwb_pereport.C.
Definition at line 4091 of file cwb_pereport.C.
| void cwb_pereport | ( | TString | options | ) |
Definition at line 402 of file cwb_pereport.C.
Definition at line 2707 of file cwb_pereport.C.
| void DumpWaveformCR | ( | int | nIFO | ) |
Definition at line 2182 of file cwb_pereport.C.
| wavearray< double > GetCCut | ( | wavearray< double > * | ts, |
| double | tcoa, | ||
| double | m1, | ||
| double | m2, | ||
| double | s1z, | ||
| double | s2z | ||
| ) |
Definition at line 3551 of file cwb_pereport.C.
| void GetCCutParms | ( | TString | options | ) |
Definition at line 3778 of file cwb_pereport.C.
Definition at line 2697 of file cwb_pereport.C.
Definition at line 4264 of file cwb_pereport.C.
Definition at line 2685 of file cwb_pereport.C.
| void GetLoglParms | ( | TString | options | ) |
Definition at line 3955 of file cwb_pereport.C.
| void GetOffSourceStatistic | ( | int | nIFO, |
| int | recId, | ||
| vector< int > | injId, | ||
| vector< double > & | oFF, | ||
| vector< double > & | oOF, | ||
| vector< double > & | oRE | ||
| ) |
Definition at line 3089 of file cwb_pereport.C.
Definition at line 2945 of file cwb_pereport.C.
Definition at line 3526 of file cwb_pereport.C.
Definition at line 3855 of file cwb_pereport.C.
| void GetRCutParms | ( | TString | options | ) |
Definition at line 3909 of file cwb_pereport.C.
| void GetRstatParms | ( | TString | options | ) |
Definition at line 4026 of file cwb_pereport.C.
| double GetSyncTcoa | ( | wavearray< double > | wf, |
| double | sync_time, | ||
| double | sync_phase, | ||
| double | tcoa | ||
| ) |
Definition at line 2854 of file cwb_pereport.C.
| void LoadOutputRootFiles | ( | int | nIFO, |
| TString | data_label, | ||
| TString | output_dir, | ||
| int | max_events = PE_MAX_EVENTS, |
||
| double | ifar_thr = 0 |
||
| ) |
Definition at line 1200 of file cwb_pereport.C.
| void MakeDistributions | ( | TString | type = "FF" | ) |
Definition at line 2281 of file cwb_pereport.C.
| void MakePvalueDistribution | ( | TString | type = "FF" | ) |
Definition at line 3207 of file cwb_pereport.C.
| void PlotWaveformAsymmErrors | ( | TString | ofname, |
| TString | title, | ||
| wavearray< double > | wrec, | ||
| wavearray< double > | wmed, | ||
| wavearray< double > | wl50, | ||
| wavearray< double > | wu50, | ||
| wavearray< double > | wl90, | ||
| wavearray< double > | wu90, | ||
| wavearray< double > | wl99, | ||
| wavearray< double > | wu99, | ||
| wavearray< double > | wref, | ||
| TString | pdir, | ||
| double | P, | ||
| double | Q, | ||
| double | tref, | ||
| double | tcoa, | ||
| int | ifoid | ||
| ) |
Definition at line 1557 of file cwb_pereport.C.
Definition at line 834 of file cwb_pereport.C.
| void PlotWaveformsCR | ( | int | nIFO | ) |
Definition at line 1498 of file cwb_pereport.C.
| void PrintUserOptions | ( | TString | ofname = "" | ) |
Definition at line 2062 of file cwb_pereport.C.
| void ReadUserOptions | ( | TString | options | ) |
Definition at line 1743 of file cwb_pereport.C.
| void ResetUserOptions | ( | ) |
Definition at line 2019 of file cwb_pereport.C.
Definition at line 3348 of file cwb_pereport.C.
Definition at line 2876 of file cwb_pereport.C.
| int gEVENTS |
Definition at line 251 of file cwb_pereport.C.
Definition at line 252 of file cwb_pereport.C.
| uoptions gOPT |
Definition at line 250 of file cwb_pereport.C.
| WDM<double>* gWDM |
Definition at line 253 of file cwb_pereport.C.
| std::vector<wavearray<double> > pFRQ |
Definition at line 303 of file cwb_pereport.C.
| std::vector<wavearray<double> > pSNR |
Definition at line 302 of file cwb_pereport.C.
| std::vector<wavearray<double> > pTIM |
Definition at line 304 of file cwb_pereport.C.
| std::vector<wavearray<double> > vAUX |
Definition at line 293 of file cwb_pereport.C.
| float vDOF |
Definition at line 316 of file cwb_pereport.C.
| double vFACTOR |
Definition at line 309 of file cwb_pereport.C.
| std::vector<double> vFF |
Definition at line 338 of file cwb_pereport.C.
| std::vector<wavearray<double> > vINJ |
Definition at line 292 of file cwb_pereport.C.
| std::vector<wavearray<double> > vL50 |
Definition at line 295 of file cwb_pereport.C.
| std::vector<wavearray<double> > vL90 |
Definition at line 297 of file cwb_pereport.C.
| std::vector<wavearray<double> > vL99 |
Definition at line 299 of file cwb_pereport.C.
| float vM1 |
Definition at line 312 of file cwb_pereport.C.
| float vM2 |
Definition at line 313 of file cwb_pereport.C.
| std::vector<wavearray<double> > vMED |
Definition at line 294 of file cwb_pereport.C.
| std::vector<double> vNE |
Definition at line 338 of file cwb_pereport.C.
| std::vector<wavearray<double> > vNUL |
Definition at line 290 of file cwb_pereport.C.
| std::vector<double> vOF |
Definition at line 338 of file cwb_pereport.C.
| float vPHI |
Definition at line 311 of file cwb_pereport.C.
| double vPSYNC[NIFO_MAX] |
Definition at line 320 of file cwb_pereport.C.
| std::vector<double> vRE |
Definition at line 338 of file cwb_pereport.C.
| std::vector<wavearray<double> > vREC |
Definition at line 291 of file cwb_pereport.C.
| float vS1 |
Definition at line 314 of file cwb_pereport.C.
| float vS2 |
Definition at line 315 of file cwb_pereport.C.
| float vSNR[NIFO_MAX] |
Definition at line 317 of file cwb_pereport.C.
| double vTCOA |
Definition at line 308 of file cwb_pereport.C.
| float vTHETA |
Definition at line 310 of file cwb_pereport.C.
| std::vector<double> vTREC |
Definition at line 306 of file cwb_pereport.C.
| double vTSYNC[NIFO_MAX] |
Definition at line 319 of file cwb_pereport.C.
| std::vector<wavearray<double> > vU50 |
Definition at line 296 of file cwb_pereport.C.
| std::vector<wavearray<double> > vU90 |
Definition at line 298 of file cwb_pereport.C.
| std::vector<wavearray<double> > vU99 |
Definition at line 300 of file cwb_pereport.C.
| std::vector<wavearray<double> > wAUX[PE_MAX_EVENTS] |
Definition at line 325 of file cwb_pereport.C.
| float wDOF[PE_MAX_EVENTS] |
Definition at line 335 of file cwb_pereport.C.
| double wFACTOR[PE_MAX_EVENTS] |
Definition at line 328 of file cwb_pereport.C.
| std::vector<wavearray<double> > wINJ[PE_MAX_EVENTS] |
Definition at line 324 of file cwb_pereport.C.
| float wM1[PE_MAX_EVENTS] |
Definition at line 331 of file cwb_pereport.C.
| float wM2[PE_MAX_EVENTS] |
Definition at line 332 of file cwb_pereport.C.
| std::vector<wavearray<double> > wNUL[PE_MAX_EVENTS] |
Definition at line 323 of file cwb_pereport.C.
| float wPHI[PE_MAX_EVENTS] |
Definition at line 330 of file cwb_pereport.C.
| std::vector<wavearray<double> > wREC[PE_MAX_EVENTS] |
Definition at line 322 of file cwb_pereport.C.
| float wS1[PE_MAX_EVENTS] |
Definition at line 333 of file cwb_pereport.C.
| float wS2[PE_MAX_EVENTS] |
Definition at line 334 of file cwb_pereport.C.
| float wSNR[PE_MAX_EVENTS][NIFO_MAX] |
Definition at line 336 of file cwb_pereport.C.
| double wTCOA[PE_MAX_EVENTS] |
Definition at line 327 of file cwb_pereport.C.
| float wTHETA[PE_MAX_EVENTS] |
Definition at line 329 of file cwb_pereport.C.
| std::vector<double> xFF |
Definition at line 342 of file cwb_pereport.C.
| std::vector<double> xNE |
Definition at line 342 of file cwb_pereport.C.
| std::vector<double> xOF |
Definition at line 342 of file cwb_pereport.C.
| std::vector<double> xRE |
Definition at line 342 of file cwb_pereport.C.