diff options
| author | Andrew Martin <amartin@cfelsgi.desy.de> | 2011-03-23 15:00:13 +0100 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:26 +0100 |
| commit | f69a74040715b52de96e493b0a9d23fbb391548a (patch) | |
| tree | 96e23ceb12ec59fdd6b523d0f44916131e1edcd6 /tests/integration_check.c | |
| parent | 4517782297a62a0a9ad4cce553d5bec6a3c71225 (diff) | |
Added background subtraction and sigma
Diffstat (limited to 'tests/integration_check.c')
| -rw-r--r-- | tests/integration_check.c | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/tests/integration_check.c b/tests/integration_check.c index f80181ee..1e6e64db 100644 --- a/tests/integration_check.c +++ b/tests/integration_check.c @@ -27,7 +27,7 @@ int main(int argc, char *argv[]) { struct image image; - double fsp, ssp, intensity, bg, max; + double fsp, ssp, intensity, bg, max, sigma; int fs, ss; image.data = malloc(128*128*sizeof(float)); @@ -63,7 +63,7 @@ int main(int argc, char *argv[]) image.beam->adu_per_photon = 100.0; /* First check: no intensity -> zero intensity and bg */ - integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, 0, 1); + integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, &sigma, 0, 1); STATUS(" First check: intensity = %.2f bg = %.2f max = %.2f\n", intensity, bg, max); if ( intensity != 0.0 ) { @@ -82,7 +82,7 @@ int main(int argc, char *argv[]) image.data[fs+image.width*ss] = 1000.0; } } - integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, 0, 1); + integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, &sigma, 0, 1); STATUS("Second check: intensity = %.2f, bg = %.2f, max = %.2f\n", intensity, bg, max); if ( fabs(intensity - M_PI*9.0*9.0*1000.0) > 4000.0 ) { @@ -104,7 +104,7 @@ int main(int argc, char *argv[]) image.data[fs+image.width*ss] = poisson_noise(10.0) - 10.0; } } - integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, 0, 1); + integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, &sigma, 0, 1); STATUS(" Third check: intensity = %.2f, bg = %.2f, max = %.2f\n", intensity, bg, max); if ( fabs(intensity) > 100.0 ) { @@ -122,7 +122,7 @@ int main(int argc, char *argv[]) image.data[fs+image.width*ss] = 1000.0; } } - integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, 0, 1); + integrate_peak(&image, 64, 64, &fsp, &ssp, &intensity, &bg, &max, &sigma, 0, 1); STATUS("Fourth check: intensity = %.2f, bg = %.2f, max = %.2f\n", intensity, bg, max); if ( fabs(intensity - M_PI*9.0*9.0*1000.0) > 4000.0 ) { |