io_online.hh File Reference

#include "km3net-dataformat/offline/Evt.hh"
#include "km3net-dataformat/offline/Hit.hh"
#include "km3net-dataformat/offline/Trk.hh"
#include "km3net-dataformat/offline/Exception.hh"
#include "km3net-dataformat/online/JDAQEvent.hh"
#include "km3net-dataformat/online/JDAQKeyHit.hh"
#include "km3net-dataformat/online/JDAQTimeslice.hh"
#include "km3net-dataformat/online/JDAQSummaryslice.hh"
#include "TStreamerInfo.h"
#include "TFile.h"
#include "TTree.h"
#include "km3net-dataformat/definitions/root.hh"
#include <map>
#include <vector>

Go to the source code of this file.

Functions
void	read (Hit &hit, const JDAQHit &daqhit)
	Read a hit from a DAQ hit. More...
void	read (Hit &hit, const JDAQKeyHit &daqhit)
	Read a hit from a DAQ key hit. More...
void	read (Evt &evt, const JDAQEvent &de)
	Read an event from a DAQ event. More...
void	read (Evt &evt, const JDAQTimeslice &ts)
	Read an event from a DAQ time slice. More...
JDAQSummaryslice *	get_summary_slice (TFile *f, int frame_index)
	Get summary slice from given file with given frame index. More...

Function Documentation

void read ( Hit &  hit, const JDAQHit &  daqhit  )

inline

Read a hit from a DAQ hit.

Parameters

hit	hit
daqhit	DAQ hit

Definition at line 31 of file io_online.hh.

{
  hit.channel_id = daqhit.getPMT();
  hit.tot        = daqhit.getToT();
  hit.tdc        = daqhit.getT(); // GetT() just return the bare TDC
}

void read ( Hit &  hit, const JDAQKeyHit &  daqhit  )

inline

Read a hit from a DAQ key hit.

Parameters

hit	hit
daqhit	DAQ key hit

Definition at line 44 of file io_online.hh.

{
  hit.id = hit.pmt_id = 0;

  hit.dom_id     = daqhit.getModuleID();
  hit.channel_id = daqhit.getPMT();
  hit.tot        = daqhit.getToT();
  hit.tdc        = daqhit.getT(); // GetT() just return the bare TDC
}

void read ( Evt &  evt, const JDAQEvent &  de  )

inline

Read an event from a DAQ event.

Parameters

evt	evt
de	DAQ event

Definition at line 60 of file io_online.hh.

{
  evt.run_id          = de.getRunNumber();
  evt.det_id          = de.getDetectorID();
  evt.frame_index     = de.getFrameIndex();
  evt.trigger_counter = de.getCounter();
  evt.overlays        = de.getOverlays();
  evt.trigger_mask    = de.getTriggerMask();
  evt.t.SetSec(     de.getTimesliceStart().getUTCseconds() );
  evt.t.SetNanoSec( de.getTimesliceStart().getUTC16nanosecondcycles() * 16 );

  // The only way to the hits that are also in the triggeredhits collection
  // is by dom and channel id and time.

  Hit h;
  std::map<int, std::map< int, std::map < unsigned int, Hit*> > > M;

  const std::vector<JDAQSnapshotHit>&   snapshotHits = de.getHits<JDAQSnapshotHit>();
  const std::vector<JDAQTriggeredHit>& triggeredHits = de.getHits<JDAQTriggeredHit>();

  // http://stackoverflow.com/questions/10735135/reallocation-in-stdvector-after-stdvector-reserve
  evt.hits.clear();
  evt.hits.reserve(snapshotHits.size());

  for (auto& daqhit : snapshotHits ) // JDAQSnapshotHit
  {
    read( h, daqhit );
    h.trig = 0;
    evt.hits.push_back( h );
    M[ h.dom_id ][ h.channel_id ][ h.tdc ] = &(*evt.hits.rbegin());
  }

  for (auto& daqtrighit : triggeredHits)
  {
    Hit* g = M[daqtrighit.getModuleID()][daqtrighit.getPMT()][daqtrighit.getT()];

    if (g)
      g->trig = daqtrighit.getTriggerMask ();
    else
      THROW(Exception, "Failed to flag snaphot hit " << (int) daqtrighit.getModuleID() << "." << (int) daqtrighit.getPMT() << " " << daqtrighit.getT());
  }
}

void read ( Evt &  evt, const JDAQTimeslice &  ts  )

inline

Read an event from a DAQ time slice.

Parameters

evt	evt
ts	DAQ time slice

Definition at line 110 of file io_online.hh.

{
  evt.run_id          = ts.getRunNumber();
  evt.det_id          = ts.getDetectorID();
  evt.frame_index     = ts.getFrameIndex();
  evt.t.SetSec( ts.getTimesliceStart().getUTCseconds() );
  evt.t.SetNanoSec( ts.getTimesliceStart().getUTC16nanosecondcycles() * 16 );
  evt.hits.clear();

  // a timeslice is a vector of JDAQSuperFrame's, which is a JDAQFrame, which
  // is an stl::vector-like object (supporting stl-like iteration.)

  Hit h; h.id = 0; h.pmt_id = 0;

  for (auto& sf : ts )
  {
    for (auto& daqhit : sf )
    {
      read( h, daqhit);
      h.dom_id = sf.getModuleID();
      evt.hits.push_back(h);
    }
  }
}

JDAQSummaryslice* get_summary_slice ( TFile *  f, int  frame_index  )

inline

Get summary slice from given file with given frame index.

Parameters

f	pointer to ROOT file
frame_index	frame index

Returns

pointer to summary slice

Definition at line 143 of file io_online.hh.

{
  static TFile* _f = 0;
  static TTree*  S = 0;
  static TBranch* BS = 0;
  JDAQSummaryslice *r = 0 ;

  if (_f != f) // setup for tree reading and build tree index
  {
    _f = f;

    // first we have to deal with the following....
    // The streamer of JDAQSummaryFrame needs to know what to do since
    // this information is not written to root file.

    const char* name = JDAQSummaryslice::Class()->GetName();

    JDAQSummaryFrame::ROOT_IO_VERSION = ((TStreamerInfo*)_f -> GetStreamerInfoList()->FindObject(name))->GetClassVersion();

    //print ("JDAQSummaryFrame::ROOT_IO_VERSION set to", JDAQSummaryFrame::ROOT_IO_VERSION );

    S = (TTree*) _f->Get( TTREE_ONLINE_TIMESLICE );

    if (!S)
    {
      THROW(Exception, "Failed to get timeslice TTree : " <<  TBRANCH_ONLINE_TIMESLICE );
    }
    BS = S->GetBranch( TBRANCH_ONLINE_TIMESLICE );
    if (!BS)
    {
      THROW(Exception, "Failed to get brach :" << TBRANCH_ONLINE_TIMESLICE );
    }

    std::cout << "building index to lookup timeslices..." << std::endl;
    int n = S->BuildIndex("frame_index");
    (void) n;
    BS->SetAddress( &r );
  }

  int nbytes = S->GetEntryWithIndex( frame_index );
  if ( nbytes == 0 )
  {
    THROW(Exception, "Failed to find entry with frame_index " << frame_index);
  }
  return r;

}