/////////////////////////////////////////////////////////
// $Id: tutRunSimulation.cxx,v 1.17 2011/09/06 18:58:35 mcgrew Exp $

#include <cmath>
#include <algorithm>
#include <iostream>

#include <tut.h>
#include <TGeoManager.h>

#include "tutTestHits.hxx"

#include "HEPUnits.hxx"
#include "ICOMETEvent.hxx"
#include "IG4HitMonitor.hxx"
#include "IG4Trajectory.hxx"

TestHits::HitChecker::const_iterator& TestHits::HitChecker::const_iterator::operator++(){
    ++fHitIt;
    if (fHitIt == (*fHitConIt)->Get<COMET::IG4HitContainer>(".")->end()) {
        ++fHitConIt;
        if (fHitConIt == (*fEventIt)->Get<COMET::IDataVector>("truth/g4Hits")->end()) {
            ++fEventIt;
            // Exit if we've reached the end
            if (fEventIt == fEventEnd){
                fCurrent = *fHitIt;
                return *this;
            }
            // Set the hit data vector
            fHitConIt = (*fEventIt)->Get<COMET::IDataVector>("truth/g4Hits")->begin();
        }
        // Set thte hit container
        fHitIt = (*fHitConIt)->Get<COMET::IG4HitContainer>(".")->begin();
    }
    fCurrent = *fHitIt;
    // Return the iterator
    return *this;
}