diff --git a/runspecS6H1feedforward b/runspecS6H1feedforward index 3256f85b0892771d8f40da86425a49112a8f9bbe..0f140b44cd4f2f8c72e004d7e0512716409a2dfd 100755 --- a/runspecS6H1feedforward +++ b/runspecS6H1feedforward @@ -1,3 +1,3 @@ #/bin/csh -/home/keithr/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000feedforward +/home/gmeadors/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/hann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000feedforward diff --git a/runspecS6H1test b/runspecS6H1test index ef658a8739d922bcd156b16dfd3c01f4a4319faa..94a701e678cbf4062789781e9f0acd4492dd3675 100755 --- a/runspecS6H1test +++ b/runspecS6H1test @@ -1,4 +1,4 @@ #/bin/csh -/home/keithr/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000test +/home/gmeadors/master/opt/lscsoft/lalapps/bin/lalapps_spec_avg --SFTs /archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-931/\*.sft\;/archive/frames/S6/pulsar/sfts/feedforward/beforeHann/H-1_H1_1800SFT_AMPSC02L2-932/\*.sft --startGPS 931000000 --endGPS 932810000 --fMin 40. --fMax 1999.999 --freqRes 0.10 --timeBaseline 1800 -v 3 --IFO H1 --outputBname S6H1_Hann_40_2000test diff --git a/spec_avg.c b/spec_avg.c index 64410b9738146eeb37b39de1a8bea730ff96fc05..17354febedbbe805e2ac21bf387c5c006cbb217a 100644 --- a/spec_avg.c +++ b/spec_avg.c @@ -253,19 +253,22 @@ int main(int argc, char **argv) } } thisavepower /= count; - if (i >= 1457995 && i <= 1458005) + // Print out the calibrated ASD bin value over about 5/16 Hz, + // centered on 850 Hz. + if (i >= 1457719 && i <= 1458281) { f = f0 + ((REAL4)i)*deltaF; - printf("SFT number=%d\n",j); - printf("Frequency=%f\n",f); - printf("Bin=%d\n",i); - printf("Real component=%f\n", sft_vect->data[0].data->data[i].re); - printf("Imaginary component=%f\n", sft_vect->data[0].data->data[i].im); - printf("Bin power=%f\n", thispower); - printf("Bin-neighbor average power=%f\n", thisavepower); - printf("Scale-factor=%f\n", scalefactor); - printf("Time baseline=%f\n", timebaseline); - fprintf(fp5,"%d %16.8f %g %g %g %g\n",j, f, thispower, thisavepower, scalefactor, timebaseline); + //printf("SFT number=%d\n",j); + //printf("Frequency=%f\n",f); + //printf("Bin=%d\n",i); + //printf("Real component=%f\n", sft_vect->data[0].data->data[i].re); + //printf("Imaginary component=%f\n", sft_vect->data[0].data->data[i].im); + //printf("Bin power=%f\n", thispower); + //printf("Bin-neighbor average power=%f\n", thisavepower); + //printf("Scale-factor=%f\n", scalefactor); + //printf("Time baseline=%f\n", timebaseline); + printf("%f\n", pow(2*thispower/scalefactor/scalefactor/timebaseline, 0.5)); + fprintf(fp5,"%d %16.8f %g %g %g %g\n",j, f, thispower, pow(2*thispower/scalefactor/scalefactor/timebaseline,0.5), scalefactor, timebaseline); } weight = 1./thisavepower; // weight = 1.;