Commit 87ac426b authored by Marcella Wijngaarden's avatar Marcella Wijngaarden

Add cbcOnly runmode for only full CBC model (intrins + extrinsic)

parent 2176551c
......@@ -298,6 +298,7 @@ struct Data
int signalModelFlag;
int cleanModelFlag;
int cleanOnlyFlag;
int cbcOnlyFlag;
// models available
int cleanFlag;
......
......@@ -864,12 +864,13 @@ void print_run_flags(FILE *fptr, struct Data *data, struct Prior *prior)
fprintf(fptr, " signal model is .......... ");
if(!data->signalFlag) fprintf(fptr, "DISABLED");
else if(!data->glitchFlag) fprintf(fptr, "REQUIRED");
else if(!data->glitchFlag || data->cbcOnlyFlag) fprintf(fptr, "REQUIRED");
else fprintf(fptr, "ENABLED");
fprintf(fptr, "\n");
fprintf(fptr, " cbc model is ............. ");
if(!data->bayesCBCFlag) fprintf(fptr, "DISABLED");
else if(data->cbcOnlyFlag) fprintf(fptr, "REQUIRED");
else fprintf(fptr, "ENABLED");
fprintf(fptr, "\n");
......@@ -986,7 +987,7 @@ void print_chain_files(struct Data *data, struct Chain *chain, struct Model **mo
//print bayescbc parameters
if(data->bayesCBCFlag)
{
print_cbc_model(chain->cbcChainFile[ic], bayescbc, ic, data->Tobs);
print_cbc_model(chain->cbcChainFile[ic], bayescbc, chain->index[ic], data->Tobs);
}
if(data->verboseFlag)fflush(chain->intChainFile[chain->index[ic]]);
......@@ -1415,6 +1416,7 @@ void parse_command_line(struct Data *data, struct Chain *chain, struct Prior *pr
//models off by default
data->fullModelFlag = 0;
data->cleanOnlyFlag = 0;
data->cbcOnlyFlag = 0;
//assume the data are analyzable
data->loudGlitchFlag = 0;
......@@ -1554,6 +1556,17 @@ void parse_command_line(struct Data *data, struct Chain *chain, struct Prior *pr
data->fullModelFlag = 0;
}
ppt = LALInferenceGetProcParamVal(commandLine, "--cbcOnly");
if(ppt)
{
data->noiseModelFlag = 0;
data->glitchModelFlag = 0;
data->signalModelFlag = 1;
data->fullModelFlag = 0;
data->cbcOnlyFlag = 1;
data->bayesCBCFlag = 1;
}
ppt = LALInferenceGetProcParamVal(commandLine, "--fullOnly");
if(ppt)
{
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment