/*
 * cdccount - an example program for accessing the contents
 *            of events stored in a hddm file.
 *
 * Richard Jones
 * GlueX collaboration
 * January 10, 2005
 *
 */

#include <stdio.h>
#include <stdlib.h>
#include "hddm_s.h"

int process_event(s_HDDM_t *event);

int main(int argc, char **argv)
{
   s_HDDM_t *thisInputEvent = 0;
   s_iostream_t *thisInputFile = 0;
   int events_with_hits = 0;
   int events = 0;
   int input;
   for (input=1; input<argc; input++) {
      if (! (thisInputFile = open_s_HDDM(argv[input]))) {
         fprintf(stderr,"Error - could not open input file %s\n",
                 argv[input]);
         exit(1);
      }
      while (thisInputEvent = read_s_HDDM(thisInputFile)) {
         events_with_hits += process_event(thisInputEvent);
         flush_s_HDDM(thisInputEvent,0);
         ++events;
      }
      close_s_HDDM(thisInputFile);
   }
   printf("Total events seen %d, with cdc hits %d.\n",events,events_with_hits);
}

int process_event(s_HDDM_t *event)
{
   s_HitView_t *hits;
   s_CdcTruthPoints_t *points;
   hits = event->physicsEvents->in[0].hitView;
   if (hits == HDDM_NULL ||
       hits->centralDC == HDDM_NULL ||
       hits->centralDC->cdcTruthPoints == HDDM_NULL) {
      return 0;
   }
   //printf("New event number %d,",event->physicsEvents->in[0].eventNo);
   //printf(" run number %d\n",event->physicsEvents->in[0].runNo);
   points = hits->centralDC->cdcTruthPoints;
   //printf(" found %d cdcTruthPoints!\n",points->mult);
   int ipoint;
   int count=0;
   for (ipoint=0; ipoint<points->mult; ipoint++) {
		s_CdcTruthPoint_t *point = &points->in[ipoint];
      if (fabs(point->dradius-19.5) < 0.5e5) {
         //printf("  dradius=%f,",point->dradius);
         //printf("  phi=%f,",point->phi);
         //printf("  primary=%s,",point->primary ? "true":"false");
         //printf("  r=%f,",point->r);
         //printf("  track=%f,",point->track);
         //printf("  z=%f,",point->z);
         //printf("  dE/dx=%f\n",point->dEdx * 1e6);
      }
      ++count;
   }
   return count;
}