diff options
| -rw-r--r-- | src/fit.c | 11 |
1 files changed, 10 insertions, 1 deletions
@@ -5766,6 +5766,7 @@ void usage(void) fprintf(stderr," --max-time maximum time in seconds per fit (default: 3600)\n"); fprintf(stderr," --skip-second-event only fit the first event after a MAST bank\n"); fprintf(stderr," --max-particles maximum number of particles to fit for (default: 1)\n"); + fprintf(stderr," --min-nhit minimum nhit to fit an event (default: 100)\n"); fprintf(stderr," -h display this help message\n"); exit(1); } @@ -5915,6 +5916,7 @@ int main(int argc, char **argv) size_t combos[100]; size_t len; char tmp[256]; + size_t nhit, min_nhit = 100; for (i = 1; i < argc; i++) { if (strlen(argv[i]) >= 2 && !strncmp(argv[i], "--", 2)) { @@ -5927,6 +5929,9 @@ int main(int argc, char **argv) } else if (!strcmp(argv[i]+2,"max-particles")) { max_particles = atoi(argv[++i]); continue; + } else if (!strcmp(argv[i]+2,"min-nhit")) { + min_nhit = atoi(argv[++i]); + continue; } } else if (argv[i][0] == '-') { switch (argv[i][1]) { @@ -6166,9 +6171,13 @@ skip_mc: rv = get_event(f,&ev,&b); + nhit = get_nhit(&ev); + + if (nhit < min_nhit) break; + if (fout) { fprintf(fout, " - gtid: %i\n", ev.gtid); - fprintf(fout, " nhit: %zu\n", get_nhit(&ev)); + fprintf(fout, " nhit: %zu\n", nhit); fprintf(fout, " fit:\n"); } |