#include <cstdio>

#include "TTree.h"

#include "base/EventProc.h"
#include "base/Event.h"
#include "hadaq/TdcSubEvent.h"

class SecondProc : public base::EventProc {
   protected:

      std::string fTdcId;      ///< tdc id where channels will be selected

      double      fHits[8];    ///< 8 channel, abstract hits

      base::H1handle  hNumHits; ///< histogram with hits number

   public:
      SecondProc(const char* procname, const char* _tdcid) :
         base::EventProc(procname),
         fTdcId(_tdcid),
         hNumHits(0)
      {
         printf("Create %s for %s\n", GetName(), fTdcId.c_str());

         hNumHits = MakeH1("NumHits","Number of hits", 32, 0, 32, "number");

         // enable storing already in constructor
         SetStoreEnabled();
      }

      virtual void CreateBranch(TTree* t)
      {
         // only called when tree is created in first.C
         // one can ignore
         t->Branch(GetName(), fHits, "hits[8]/D");
      }

      virtual bool Process(base::Event* ev)
      {
         for (unsigned n=0;n<8;n++) fHits[n] = 0.;

         hadaq::TdcSubEvent* sub =
               dynamic_cast<hadaq::TdcSubEvent*> (ev->GetSubEvent(fTdcId));

         // printf("%s process sub %p %s\n", GetName(), sub, fTdcId.c_str());

         if (sub==0) return false;

         double num(0), ch0tm(0);

         for (unsigned cnt=0;cnt<sub->Size();cnt++) {
            const hadaq::TdcMessageExt& ext = sub->msg(cnt);

            unsigned chid = ext.msg().getHitChannel();
            if (chid==0) { ch0tm = ext.GetGlobalTime(); continue; }

            // full time
            double tm = ext.GetGlobalTime() + ch0tm;

            // printf("  ch:%3d tm:%f\n", chid, tm);
            num+=1;
         }

         FillH1(hNumHits, num);

         return true;
      }
};


void second()
{
   new SecondProc("Second1", "TDC_1000");
}