#include "readEvents.C"

void countEvents() {
    readEvents();
    int dataEvents = 0;
    int npi0Events = 0;
    int pi0Events = 0;
    for (struct eventData* ev = gEventData; 
         ev < &gEventData[gNEvents]; 
         ++ev) {
        if (ev->type == 0) ++dataEvents;
        if (ev->type == 1) ++npi0Events;
        if (ev->type == 2) ++pi0Events;
    }
    std::cout << "Data Events:   " << dataEvents << std::endl;
    std::cout << "MC Events:     " << npi0Events + pi0Events << std::endl;
    std::cout << "   Non-pi0:       " << npi0Events << std::endl;
    std::cout << "   pi0:           " << pi0Events << std::endl;
}