diff options
Diffstat (limited to 'src/indexamajig.c')
| -rw-r--r-- | src/indexamajig.c | 13 |
1 files changed, 11 insertions, 2 deletions
diff --git a/src/indexamajig.c b/src/indexamajig.c index b8d3df88..b118593e 100644 --- a/src/indexamajig.c +++ b/src/indexamajig.c @@ -61,6 +61,7 @@ struct process_args int config_sanity; int config_satcorr; int config_sa; + float threshold; struct detector *det; IndexingMethod indm; IndexingPrivate *ipriv; @@ -145,6 +146,7 @@ static void show_help(const char *s) " included in the HDF5 file.\n" " --no-sa Don't correct for the differing solid angles of\n" " the pixels.\n" +" --threshold=<n> Only accept peaks above <n> ADU. Default: 800.\n" "\n" "\nOptions for greater performance or verbosity:\n\n" " --verbose Be verbose about indexing.\n" @@ -325,7 +327,7 @@ static struct process_result process_image(struct process_args *pargs) } /* Perform 'fine' peak search */ - search_peaks(&image); + search_peaks(&image, pargs->threshold); /* Get rid of noise-filtered version at this point * - it was strictly for the purposes of peak detection. */ @@ -468,6 +470,7 @@ int main(int argc, char *argv[]) int config_sanity = 0; int config_satcorr = 0; int config_sa = 1; + float threshold = 800.0; struct detector *det; char *geometry = NULL; IndexingMethod indm; @@ -512,11 +515,12 @@ int main(int argc, char *argv[]) {"check-sanity", 0, &config_sanity, 1}, {"sat-corr", 0, &config_satcorr, 1}, {"no-sa", 0, &config_sa, 0}, + {"threshold", 1, NULL, 't'}, {0, 0, NULL, 0} }; /* Short options */ - while ((c = getopt_long(argc, argv, "hi:wp:j:x:g:", + while ((c = getopt_long(argc, argv, "hi:wp:j:x:g:t:", longopts, NULL)) != -1) { switch (c) { @@ -552,6 +556,10 @@ int main(int argc, char *argv[]) geometry = strdup(optarg); break; + case 't' : + threshold = strtof(optarg, NULL); + break; + case 0 : break; @@ -707,6 +715,7 @@ int main(int argc, char *argv[]) pargs->indm = indm; pargs->intensities = intensities; pargs->gctx = gctx; + pargs->threshold = threshold; pargs->id = i; pthread_mutex_lock(&pargs->control_mutex); pargs->done = 0; |