aemspillmipp.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// $Id: aemspillmipp.cc,v 1.1 2005/11/21 20:35:30 hmeyer Exp $
//
// Simple utility to get time and # of events for spills in a run
// Based on rawdump.cc from MippIo/test
//
// hmeyer@fnal.gov
////////////////////////////////////////////////////////////////////////
#include <iostream>
#include "MippIo/MippIoFile.h"
#include "MippIo/MippIoBlock.h"
#include "MippIo/MippIoFileHeader.h"
#include "MippIo/MippIoEvent.h"
#include "MippIo/MippIoEOFBlock.h"
#include "MippIo/MippIoEventHeader.h"
#include "MippIo/MippIoMDDBlock.h"
#include "MippIo/MippIoStatsBlock.h"

std::ostream& os = std::cout;

static void gsDumpFileHeader(const MippIoFileHeader* h) 
{
  os <<
    "************************************"
    "************************************"       << std::endl;
  os << "* Run:         " << h->Data().fRunNumber     << 
    "."              << h->Data().fSubRunNumber  << std::endl <<
    "* Run type:        "   << MippIo::GetRunTypeName(h->GetRunType());
  time_t t = h->Data().fCreationDate;
  os << "\n* Date opened:   "     << ctime(&t) <<
    "* Beam momentum:   " << h->GetMomentum() <<
    "\n* Enabled triggers:";
  for (int i = 0; i < 16; ++i) 
    if(h->IsTrigOn(i)) std::cout << " " << i;
  os << "\n* File size:       " << h->Data().fFileSize      <<
    "\n* Number of events: "  << h->Data().fNevents       << std::endl <<
    "* Run Comment: "         << h->Data().fRunComment    << std::endl <<
    "* End of Run Comment: "  << h->Data().fEndRunComment <<
    "\n*";

  for (int i = 0; i < MippIoFileHeader::fsNTrigBit; ++i) {
    os << "\n* Trigger " << i << ": " << h->GetTrigCfg(i).fName 
       << ", prescale: " << h->GetTrigCfg(i).fPrescale;
  }
  os << "\n*";
  for (int i = 0; i < MippIoFileHeader::fsNTrigInput; ++i) {
    os << "\n* Trigger input " << i << ": " << h->GetTrigInputName(i);
  }

  os << "\n*****" << std::endl;
}

//......................................................................

static void gsDumpSpills(const MippIoEvent* evt) 
{
//======================================================================
// Dump Spill data from MDD block
//======================================================================

  static int fNTrigInSpill = 0;

  //Skip inter-spill calibration events
  const MippIoStatsBlock* stBlk =
    evt->Block<MippIoStatsBlock>(MippIoBlock::kStatsId);
  bool calibEvent = false;
  if (stBlk) {
    const MippIoStat104* st =
      dynamic_cast<const MippIoStat104*>(stBlk->GetStat());
    if (st) {
      if (!st->fLatch.empty()) {
        if (st->fLatch[0] & 0x2000)  calibEvent = true;
      }
    }
  }
  if (calibEvent) return;
  
  //Count In-spill triggers for each spill
  const MippIoMDDBlock* mddb =
    evt->Block<MippIoMDDBlock>(MippIoBlock::kBeamPPCId);
  if (mddb) {
    const MippIoMDDBlock::Data_t& mddData = mddb->Data();

    if (mddData.fIRQ == 3) {
      os << fNTrigInSpill << " events_at " << mddData.fSec << std::endl;
      fNTrigInSpill = 0;
    }
    else {
      fNTrigInSpill++;
    }

  } // valid mdd block
}

//......................................................................


static void gsDumpEOF(const MippIoEOFBlock* eof) 
{  
  os << 
    "*" <<
    "** File closed" <<
    " on "        << eof->Data().fDateClosed <<
    " size = "    << eof->Data().fFileSize   << " bytes" <<
    " nEvents = " << eof->Data().fNevents    << std::endl << 
    "************************************" <<
    "************************************" << std::endl;
}

//......................................................................

static void gsDumpEvent(const MippIoEvent* evt) 
{
  gsDumpSpills(evt);
}

//......................................................................

static void gsRawDump(const char* filename) 
{
  MippIoFile f(MippIoFile::kReadable);
  
  f.Open(filename);

  const MippIoBlock* b   = 0;
  bool               eof = false;
  while (1) {
    b = f.ReadBlock();
    if (b<=0) {
      std::cerr << "rawdump: ReadBlock failed..." << std::endl;
      break;
    }
    switch (b->Prefix().fId) {
    case MippIoBlock::kFileHeaderId: 
      gsDumpFileHeader((const MippIoFileHeader*)b); 
      break;
    case MippIoBlock::kEventId:
      gsDumpEvent((const MippIoEvent*)b);
      break;
    case MippIoBlock::kEOFId:
      gsDumpEOF((const MippIoEOFBlock*)b); 
      eof=true; 
      break;
    default:
      std::cerr << "rawdump: Bad block. prefix="
        << (void*)b->Prefix().fId        << " "
        << (void*)b->Prefix().fVersion   << " "
        << (void*)b->Prefix().fBlockSize << " " 
        << std::endl;
    }
    if (eof) break;
  }
  f.Close();
}

//......................................................................

int main(int argc, char** argv) 
{
  for (int i=1; i<argc; ++i) {
    gsRawDump(argv[i]); 
  }
  return 1;
}

////////////////////////////////////////////////////////////////////////

---