Commit ddc6fb07 authored by Cristina Valeria Torres's avatar Cristina Valeria Torres
Browse files

Merge branch 'master' of git+ssh://cristina.torres@ligo-vcs.phys.uwm.edu/usr/local/git/lalsuite

Original: d6cbb87a273b58359828205be9dfdeee7a9f222e
parents 082d4a7a 36d285b3
......@@ -79,7 +79,7 @@ lalapps_ComputeFStatistic_v2_cuda_LDADD = $(LDADD) -lcudart
lalapps_PredictFStat_SOURCES = PredictFStat.c ../transientCW_utils.c ../transientCW_utils.h ../ProbabilityDensity.c ../ProbabilityDensity.h
# force linking in C++ mode for E@H
eah_PredictFStat_SOURCES = PredictFStat.c ../hough/src2/EinsteinAtHome/win_lib.cpp
eah_PredictFStat_SOURCES = $(lalapps_PredictFStat_SOURCES) ../hough/src2/EinsteinAtHome/win_lib.cpp
lalapps_ComputeFStatAnalyticMonteCarloUpperLimit_SOURCES = ComputeFStatAnalyticMonteCarloUpperLimit.c
......
......@@ -18,7 +18,7 @@ EXTRA_PROGRAMS = psinject eah_Makefakedata_v4 $(frameprogs)
lalapps_Makefakedata_v4_SOURCES = makefakedata_v4.c sft_extra.c sft_extra.h ../transientCW_utils.c ../transientCW_utils.h ../ProbabilityDensity.c ../ProbabilityDensity.h
# force linking in C++ mode for E@H
eah_Makefakedata_v4_SOURCES = makefakedata_v4.c sft_extra.c sft_extra.h ../hough/src2/EinsteinAtHome/win_lib.cpp
eah_Makefakedata_v4_SOURCES = $(lalapps_Makefakedata_v4_SOURCES) ../hough/src2/EinsteinAtHome/win_lib.cpp
lalapps_makefakedata_v2_SOURCES = makefakedata_v2.c sft_extra.c sft_extra.h
......
......@@ -170,8 +170,8 @@ for i; do
echo " --check-app=<app> only test the app specified, not necessarily the one just built"
echo " --release use some dark magic to make the App most compatible and add remote debugging."
echo " Implies --static and --rebuild and even more dirty hacks on Linux to work on Woody"
echo " --appname set an application name (only used in --release builds, defaults to einstein_S5GC1HF)"
echo " --appversion set an application version (only used in --release builds, defaults to 0.00)"
echo " --appname=<name> set an application name (only used in --release builds, defaults to einstein_S5GC1HF)"
echo " --appversion=N.NN set an application version (only used in --release builds, defaults to 0.00)"
echo " --norebuild disables --rebuild on --release. DANGEROUS! Use only for testing the build script"
echo " --help show this message and exit"
exit ;;
......@@ -252,8 +252,8 @@ fi
# if --pather and not --altivec, make sure the binary runs on G3
test ."$MACOSX_DEPLOYMENT_TARGET" = ."10.3" -a ."$acc" = ."" &&
CFLAGS="-mcpu=G3 $CFLAGS" &&
CXXFLAGS="-mcpu=G3 $CXXFLAGS"
CFLAGS="-mcpu=G3 $CFLAGS" &&
CXXFLAGS="-mcpu=G3 $CXXFLAGS"
if [ ".$cuda" = ".true" -a ."$build_win32" = ."true" ]; then
export CFLAGS="-g0 $CFLAGS"
......
......@@ -696,6 +696,18 @@ static void worker (void) {
res = HIERARCHICALSEARCH_EFILE;
}
/* boinc_resolve and unzip segment list */
else if (MATCH_START("--segmentList=",argv[arg],l)) {
rargv[rarg] = (char*)calloc(MAX_PATH_LEN,sizeof(char));
if(!rargv[rarg]){
LogPrintf(LOG_CRITICAL, "Out of memory\n");
boinc_finish(HIERARCHICALSEARCH_EMEM);
}
strncpy(rargv[rarg],argv[arg],l);
if (resolve_and_unzip(argv[arg]+l, rargv[rarg]+l, MAX_PATH_LEN-l) < 0)
res = HIERARCHICALSEARCH_EFILE;
}
/* boinc_resolve and unzip ephermeris files */
else if (MATCH_START("--ephemE=",argv[arg],l)) {
rargv[rarg] = (char*)calloc(MAX_PATH_LEN,sizeof(char));
......
......@@ -29,19 +29,24 @@ static struct coh_PTF_params *coh_PTF_get_params( int argc, char **argv )
int main( int argc, char **argv )
{
/* Declarations of parameters */
INT4 i,j,k;
UINT4 ui,uj,sp;
/* process structures */
struct coh_PTF_params *params = NULL;
ProcessParamsTable *procpar = NULL;
/* sky position structures */
UINT4 numSkyPoints;
struct coh_PTF_skyPoints *skyPoints = NULL;
/* FFT structures */
REAL4FFTPlan *fwdplan = NULL;
REAL4FFTPlan *revplan = NULL;
COMPLEX8FFTPlan *invPlan = NULL;
/* input data and spectrum storage */
UINT4 numDetectors = 0;
UINT4 singleDetector = 0;
......@@ -49,17 +54,21 @@ int main( int argc, char **argv )
REAL4FrequencySeries *invspec[LAL_NUM_IFO+1];
RingDataSegments *segments[LAL_NUM_IFO+1];
INT4 numSegments = 0;
/* template counters */
INT4 numTmplts = 0;
INT4 numSpinTmplts = 0;
INT4 numNoSpinTmplts = 0;
/* template indexes */
INT4 startTemplate = -1;
INT4 stopTemplate = -1;
/* bank structures */
UINT4 UNUSED spinBank = 0;
char spinFileName[256];
char noSpinFileName[256];
/* template and findchirp data structures */
InspiralTemplate *PTFSpinTemplate = NULL;
InspiralTemplate *PTFNoSpinTemplate = NULL;
......@@ -75,6 +84,7 @@ int main( int argc, char **argv )
REAL8Array *PTFM[LAL_NUM_IFO+1];
REAL8Array *PTFN[LAL_NUM_IFO+1];
COMPLEX8VectorSequence *PTFqVec[LAL_NUM_IFO+1];
/* triggered sky position and sensitivity structures */
LIGOTimeGPS segStartTime;
struct timeval startTime;
......@@ -85,6 +95,7 @@ int main( int argc, char **argv )
REAL8 *Fplustrig;
REAL8 *Fcrosstrig;
REAL8 detLoc[3];
/* coherent statistic structures */
REAL4TimeSeries *cohSNR = NULL;
REAL4TimeSeries *pValues[10];
......@@ -92,10 +103,12 @@ int main( int argc, char **argv )
REAL4TimeSeries *gammaBeta[2];
REAL4TimeSeries *nullSNR = NULL;
REAL4TimeSeries *traceSNR = NULL;
/* consistency test structures */
REAL4TimeSeries *bankVeto = NULL;
REAL4TimeSeries *autoVeto = NULL;
REAL4TimeSeries *chiSquare = NULL;
/* output event structures */
MultiInspiralTable *eventList = NULL;
MultiInspiralTable *thisEvent = NULL;
......@@ -732,7 +745,7 @@ int main( int argc, char **argv )
}
verbose( "Made filters for ifo %d,segment %d, template %d at %ld \n",
ifoNumber,j,i,timeval_subtract(&startTime));
ifoNumber, j, i, timeval_subtract(&startTime) );
}
}
/* If necessary calculate the null stream filters */
......@@ -757,8 +770,6 @@ int main( int argc, char **argv )
j,i,timeval_subtract(&startTime));
}
/* FIXME: Sky location looping should go here.
Clustering may be required if we are doing this! */
switch(params->skyLooping)
{
case SINGLE_SKY_POINT:
......@@ -768,8 +779,10 @@ int main( int argc, char **argv )
/* set 'segStartTime' to trigger time */
segStartTime = params->trigTime;
verbose( "Begin loop over sky points at %ld \n", timeval_subtract(&startTime));
verbose( "Begin loop over sky points at %ld \n",
timeval_subtract(&startTime) );
/* loop over sky points */
for ( sp = 0; sp < numSkyPoints ; sp++ )
{
......@@ -810,7 +823,7 @@ int main( int argc, char **argv )
verbose( "Made coherent statistic for segment %d, template %d, sky point %d at %ld \n", j, i, sp, timeval_subtract(&startTime) );
// This function adds any loud events to the list of triggers
/* This function construct triggers from loud events */
eventId = coh_PTF_add_triggers( params, &eventList, &thisEvent,
cohSNR, *PTFtemplate, eventId,
spinTemplate, singleDetector,
......@@ -819,14 +832,24 @@ int main( int argc, char **argv )
autoVeto, chiSquare, PTFM );
verbose( "Generated triggers for segment %d, template %d, sky point %d at %ld \n", j, i, sp, timeval_subtract(&startTime) );
// Clustering can happen here. The clustering routine needs refining
// coh_PTF_cluster_triggers(params,&eventList,&thisEvent);
// verbose( "Clustered triggers for segment %d, template %d at %ld \n",
// j,i,time(NULL)-startTime);
}
break;
}
// Then we get a bunch of memory freeing statements
case ALL_SKY:
case TWO_DET_ALL_SKY:
error( "All sky code not implemented yet\n" );
break;
default:
error( "Oops. No sky type set, Ian done broke the code.\n" );
break;
}
/* cluster triggers */
/* XLALClusterMultiInspiralTable( ); */
/* verbose( "Clustered triggers for segment %d, template %d at %ld \n",
j, i, time(NULL)-startTime ); */
/* Then we get a bunch of memory freeing statements */
for ( k = 0 ; k < 10 ; k++ )
{
if (pValues[k])
......@@ -980,21 +1003,24 @@ void coh_PTF_statistic(
)
{
verbose( "Entering the statistic loop at %ld \n",timeval_subtract(&startTime));
// This function generates the SNR for every point in time and, where
// appropriate calculates the desired signal based vetoes.
verbose( "Entering the statistic loop at %ld \n",
timeval_subtract(&startTime) );
/* This function generates the SNR for every point in time and, where
* appropriate calculates the desired signal based vetoes. */
UINT4 i, j, k, m, vecLength, vecLengthTwo, UNUSED vecLengthSquare, UNUSED vecLengthTwoSquare;
INT4 l;
INT4 timeOffsetPoints[LAL_NUM_IFO];
REAL4 deltaT = cohSNR->deltaT;
UINT4 numPoints = floor( params->segmentDuration * params->sampleRate + 0.5 );
UINT4 i, j, k, m, vecLength, vecLengthTwo;
INT4 l;
INT4 timeOffsetPoints[LAL_NUM_IFO];
REAL4 deltaT = cohSNR->deltaT;
UINT4 numPoints = floor( params->segmentDuration * params->sampleRate\
+ 0.5 );
struct bankDataOverlaps *chisqOverlaps = *chisqOverlapsP;
REAL4 *frequencyRangesPlus = *frequencyRangesPlusP;
REAL4 *frequencyRangesCross = *frequencyRangesCrossP;
REAL4 *frequencyRangesPlus = *frequencyRangesPlusP;
REAL4 *frequencyRangesCross = *frequencyRangesCrossP;
// Code works slightly differently if spin/non spin and single/coherent
/* Code works slightly differently if spin/non spin and single/coherent */
if (spinTemplate)
vecLength = 5;
else
......@@ -1003,36 +1029,38 @@ void coh_PTF_statistic(
vecLengthTwo = vecLength;
else
vecLengthTwo = 2* vecLength;
vecLengthSquare = vecLength*vecLength;
vecLengthTwoSquare = vecLengthTwo * vecLengthTwo;
// These arrays are used to store the maximized quantities
// For non spin these are the 4 F-stat parameters (only 2 for one detector)
// For spin these are the P values
/* These arrays are used to store the maximized quantities
* For non spin these are the 4 F-stat parameters (only 2 for one detector)
* For spin these are the P values */
for ( i = 0 ; i < vecLengthTwo ; i++ )
{
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length );
}
if (! spinTemplate)
{
for ( i = vecLengthTwo ; i < 2*vecLengthTwo ; i++ )
{
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length );
}
}
if (spinTemplate)
{
/* These store a amplitude and phase information for PTF search */
gammaBeta[0] = XLALCreateREAL4TimeSeries( "Gamma",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
gammaBeta[1] = XLALCreateREAL4TimeSeries( "Beta",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
gammaBeta[0] = XLALCreateREAL4TimeSeries( "Gamma", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length );
gammaBeta[1] = XLALCreateREAL4TimeSeries( "Beta", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length);
}
/* FIXME: All the time series should be outputtable */
......
......@@ -29,19 +29,24 @@ static struct coh_PTF_params *coh_PTF_get_params( int argc, char **argv )
int main( int argc, char **argv )
{
/* Declarations of parameters */
INT4 i,j,k;
UINT4 ui,uj,sp;
/* process structures */
struct coh_PTF_params *params = NULL;
ProcessParamsTable *procpar = NULL;
/* sky position structures */
UINT4 numSkyPoints;
struct coh_PTF_skyPoints *skyPoints = NULL;
/* FFT structures */
REAL4FFTPlan *fwdplan = NULL;
REAL4FFTPlan *revplan = NULL;
COMPLEX8FFTPlan *invPlan = NULL;
/* input data and spectrum storage */
UINT4 numDetectors = 0;
UINT4 singleDetector = 0;
......@@ -49,17 +54,21 @@ int main( int argc, char **argv )
REAL4FrequencySeries *invspec[LAL_NUM_IFO+1];
RingDataSegments *segments[LAL_NUM_IFO+1];
INT4 numSegments = 0;
/* template counters */
INT4 numTmplts = 0;
INT4 numSpinTmplts = 0;
INT4 numNoSpinTmplts = 0;
/* template indexes */
INT4 startTemplate = -1;
INT4 stopTemplate = -1;
/* bank structures */
UINT4 UNUSED spinBank = 0;
char spinFileName[256];
char noSpinFileName[256];
/* template and findchirp data structures */
InspiralTemplate *PTFSpinTemplate = NULL;
InspiralTemplate *PTFNoSpinTemplate = NULL;
......@@ -75,6 +84,7 @@ int main( int argc, char **argv )
REAL8Array *PTFM[LAL_NUM_IFO+1];
REAL8Array *PTFN[LAL_NUM_IFO+1];
COMPLEX8VectorSequence *PTFqVec[LAL_NUM_IFO+1];
/* triggered sky position and sensitivity structures */
LIGOTimeGPS segStartTime;
struct timeval startTime;
......@@ -85,6 +95,7 @@ int main( int argc, char **argv )
REAL8 *Fplustrig;
REAL8 *Fcrosstrig;
REAL8 detLoc[3];
/* coherent statistic structures */
REAL4TimeSeries *cohSNR = NULL;
REAL4TimeSeries *pValues[10];
......@@ -92,10 +103,12 @@ int main( int argc, char **argv )
REAL4TimeSeries *gammaBeta[2];
REAL4TimeSeries *nullSNR = NULL;
REAL4TimeSeries *traceSNR = NULL;
/* consistency test structures */
REAL4TimeSeries *bankVeto = NULL;
REAL4TimeSeries *autoVeto = NULL;
REAL4TimeSeries *chiSquare = NULL;
/* output event structures */
MultiInspiralTable *eventList = NULL;
MultiInspiralTable *thisEvent = NULL;
......@@ -732,7 +745,7 @@ int main( int argc, char **argv )
}
verbose( "Made filters for ifo %d,segment %d, template %d at %ld \n",
ifoNumber,j,i,timeval_subtract(&startTime));
ifoNumber, j, i, timeval_subtract(&startTime) );
}
}
/* If necessary calculate the null stream filters */
......@@ -757,8 +770,6 @@ int main( int argc, char **argv )
j,i,timeval_subtract(&startTime));
}
/* FIXME: Sky location looping should go here.
Clustering may be required if we are doing this! */
switch(params->skyLooping)
{
case SINGLE_SKY_POINT:
......@@ -768,8 +779,10 @@ int main( int argc, char **argv )
/* set 'segStartTime' to trigger time */
segStartTime = params->trigTime;
verbose( "Begin loop over sky points at %ld \n", timeval_subtract(&startTime));
verbose( "Begin loop over sky points at %ld \n",
timeval_subtract(&startTime) );
/* loop over sky points */
for ( sp = 0; sp < numSkyPoints ; sp++ )
{
......@@ -810,7 +823,7 @@ int main( int argc, char **argv )
verbose( "Made coherent statistic for segment %d, template %d, sky point %d at %ld \n", j, i, sp, timeval_subtract(&startTime) );
// This function adds any loud events to the list of triggers
/* This function construct triggers from loud events */
eventId = coh_PTF_add_triggers( params, &eventList, &thisEvent,
cohSNR, *PTFtemplate, eventId,
spinTemplate, singleDetector,
......@@ -819,14 +832,24 @@ int main( int argc, char **argv )
autoVeto, chiSquare, PTFM );
verbose( "Generated triggers for segment %d, template %d, sky point %d at %ld \n", j, i, sp, timeval_subtract(&startTime) );
// Clustering can happen here. The clustering routine needs refining
// coh_PTF_cluster_triggers(params,&eventList,&thisEvent);
// verbose( "Clustered triggers for segment %d, template %d at %ld \n",
// j,i,time(NULL)-startTime);
}
break;
}
// Then we get a bunch of memory freeing statements
case ALL_SKY:
case TWO_DET_ALL_SKY:
error( "All sky code not implemented yet\n" );
break;
default:
error( "Oops. No sky type set, Ian done broke the code.\n" );
break;
}
/* cluster triggers */
/* XLALClusterMultiInspiralTable( ); */
/* verbose( "Clustered triggers for segment %d, template %d at %ld \n",
j, i, time(NULL)-startTime ); */
/* Then we get a bunch of memory freeing statements */
for ( k = 0 ; k < 10 ; k++ )
{
if (pValues[k])
......@@ -980,21 +1003,24 @@ void coh_PTF_statistic(
)
{
verbose( "Entering the statistic loop at %ld \n",timeval_subtract(&startTime));
// This function generates the SNR for every point in time and, where
// appropriate calculates the desired signal based vetoes.
verbose( "Entering the statistic loop at %ld \n",
timeval_subtract(&startTime) );
/* This function generates the SNR for every point in time and, where
* appropriate calculates the desired signal based vetoes. */
UINT4 i, j, k, m, vecLength, vecLengthTwo, UNUSED vecLengthSquare, UNUSED vecLengthTwoSquare;
INT4 l;
INT4 timeOffsetPoints[LAL_NUM_IFO];
REAL4 deltaT = cohSNR->deltaT;
UINT4 numPoints = floor( params->segmentDuration * params->sampleRate + 0.5 );
UINT4 i, j, k, m, vecLength, vecLengthTwo;
INT4 l;
INT4 timeOffsetPoints[LAL_NUM_IFO];
REAL4 deltaT = cohSNR->deltaT;
UINT4 numPoints = floor( params->segmentDuration * params->sampleRate\
+ 0.5 );
struct bankDataOverlaps *chisqOverlaps = *chisqOverlapsP;
REAL4 *frequencyRangesPlus = *frequencyRangesPlusP;
REAL4 *frequencyRangesCross = *frequencyRangesCrossP;
REAL4 *frequencyRangesPlus = *frequencyRangesPlusP;
REAL4 *frequencyRangesCross = *frequencyRangesCrossP;
// Code works slightly differently if spin/non spin and single/coherent
/* Code works slightly differently if spin/non spin and single/coherent */
if (spinTemplate)
vecLength = 5;
else
......@@ -1003,36 +1029,38 @@ void coh_PTF_statistic(
vecLengthTwo = vecLength;
else
vecLengthTwo = 2* vecLength;
vecLengthSquare = vecLength*vecLength;
vecLengthTwoSquare = vecLengthTwo * vecLengthTwo;
// These arrays are used to store the maximized quantities
// For non spin these are the 4 F-stat parameters (only 2 for one detector)
// For spin these are the P values
/* These arrays are used to store the maximized quantities
* For non spin these are the 4 F-stat parameters (only 2 for one detector)
* For spin these are the P values */
for ( i = 0 ; i < vecLengthTwo ; i++ )
{
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length );
}
if (! spinTemplate)
{
for ( i = vecLengthTwo ; i < 2*vecLengthTwo ; i++ )
{
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
pValues[i] = XLALCreateREAL4TimeSeries( "Pvalue", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length );
}
}
if (spinTemplate)
{
/* These store a amplitude and phase information for PTF search */
gammaBeta[0] = XLALCreateREAL4TimeSeries( "Gamma",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
gammaBeta[1] = XLALCreateREAL4TimeSeries( "Beta",
&cohSNR->epoch,cohSNR->f0,cohSNR->deltaT,
&lalDimensionlessUnit,cohSNR->data->length);
gammaBeta[0] = XLALCreateREAL4TimeSeries( "Gamma", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length );
gammaBeta[1] = XLALCreateREAL4TimeSeries( "Beta", &cohSNR->epoch,
cohSNR->f0, cohSNR->deltaT,
&lalDimensionlessUnit,
cohSNR->data->length);
}
/* FIXME: All the time series should be outputtable */
......
......@@ -718,7 +718,7 @@ void coh_PTF_sky_grid(
REAL4 alpha,detalpha;
alpha = 0;
REAL4 angularResolution; /* angular resolution of grid in radians */
double baseline,distance;
double baseline,lightTravelTime;
REAL4 theta,phi; /* sky parameters */
REAL4 raNp = 0.; /* north */
REAL4 decNp = LAL_PI / 2.; /* pole */
......@@ -753,18 +753,12 @@ void coh_PTF_sky_grid(
params->declination,
&params->trigTime );
distance = sqrt( pow( detectors[i]->location[0] -
detectors[j]->location[0], 2) +
pow( detectors[i]->location[1] -
detectors[j]->location[1], 2) +
pow( detectors[i]->location[2] -
detectors[j]->location[2], 2 ) );
/* convert to time */
distance = distance / LAL_C_SI;
/* get light travel time */
lightTravelTime = XLALLightTravelTime( detectors[i], detectors[j] );
lightTravelTime *= 1e-9;
/* calculate opening angle */
angle = acos( baseline/distance );
angle = acos( baseline/lightTravelTime );
/* generate angular window with sky error */
lambdamin = angle-params->skyError;
......@@ -783,7 +777,7 @@ void coh_PTF_sky_grid(
lambda = lambdamax;
/* calculate alpha */
detalpha = distance * sin( lambda );
detalpha = lightTravelTime * sin( lambda );
if ( detalpha > alpha )
alpha = detalpha;
......@@ -867,9 +861,6 @@ void coh_PTF_sky_grid(
skyPoints->rightAscension[i] = phi;
skyPoints->declination[i] = LAL_PI / 2. - theta;
if ( skyPoints->declination[i] < 0. )
skyPoints->declination[i] += LAL_PI / 2.;
}
for( ifoNumber = 0; ifoNumber < LAL_NUM_IFO; ifoNumber++)
......
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