/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/parser/Hardroc2Reader.cpp

Go to the documentation of this file.

/* @version $Revision: 1328 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2011-10-03 17:04:17 +0200 (Mon, 03 Oct 2011) $ */

////////////////////////////////////////////////////////////////////////
// MICROMEGAS
//
//       EVENT RECONSTRUCTION SOFTWARE FOR HARDROC BOARDS
//
////////////////////////////////////////////////////////////////////////
// v7 8.6.2009
//    29.6.2009 : added timestamps with Hardroc data
////////////////////////////////////////////////////////////////////////
// timestamp (v7) :
//   t(hit) = t(lecture) - t2 + t3
//   l(lecture) = t0 + (t1 - t10)
//   avec
//     t0 = timestamp unix (32 bits, avec ou sans millisecondes)
//     t1 = BcId absolu (bcId_Abs, 48 bits)
//     t10 = premiere lecture de t1
//     t2 = BcId pour la DIF (bcId_HRHit, 24 bits)
//     t3 = BCId pour la HR (bcId_Hit, 24 bits en CODE GRAY)
////////////////////////////////////////////////////////////////////////

#include <stdio.h>
#include <errno.h>

#include "parser/Hardroc2Reader.hh"
#include "parser/AcquisitionParser.hh"

#include "event/Event.hh"

#include "geometry/Chamber.hh"
#include "geometry/Dif.hh"
#include "geometry/Hardroc2Chip.hh"
#include "geometry/Detector.hh"

#include "tools/Toolbox.hh"
#include "tools/SteerDesc.hh"
#include "tools/Log.hh"
#include "tools/MicroException.hh"

#include <iostream>
#include <iomanip>

using namespace std;

////////////////////////////////////////////////////////////////////////
Hardroc2Reader::Hardroc2Reader(Run& aRun, FILE *aFile, ui32 lastEventId) : DifReader(aRun, aFile, lastEventId)
{
  FILE_LOG(logDEBUG1) << "--Hardroc2Reader constructor()"<< endl;
  reset();
}

//--------------- Destructeur -------------//
Hardroc2Reader::~Hardroc2Reader()
{
  FILE_LOG(logDEBUG1)<< "----Hardroc2Reader destructor"<< endl;
}

//-------------------------------------------------------------------
void Hardroc2Reader::reset() {
  bcId_AbsOrg = -1; // the bcId origin in milliseconds
  CRC = 0; // current CRC for the whole file
  SkippedBcId = 0;

  // following initializations are useless
  globalTime = 0;
  difId = -1;
  difTriggerCounter = acqTriggerCounter = globalTriggerCounter = 0;
  bcPeriod = 200;
  bcId_Abs = bcId_Dif = bcId_Hit = 0;
} // Hardroc2Reader.reset


//-------------------------------------------------------------------
int Hardroc2Reader::newHit(Event& event, const Detector &detector, int data, const int chNum, const int chipId, const int difId) {
  // cannot be binary pattern 10

  const ChamberMap_t &chambers = detector.getChambers();
  for (ChamberMap_t::const_iterator it = chambers.begin(); it != chambers.end(); ++it) {
    try {
      Chamber& chamber = *(it->second);
      const BoardMap_t& boards = chamber.getBoardsByDifId(difId);
      // look for type of chamber
//      if (chamber.getType() != HARDROC1CHAMBER)
//        continue;
      // store new event data
      try {
        const Channel& channel = chamber.getChannelById(chNum, chipId, difId);
        i32 analogValue = 0;
        ChannelHit *hit = new ChannelHit(channel, data, analogValue);

        hit->setHardrocTime(globalTime, bcId_Abs, bcId_Dif, bcId_Hit,0);
        event.insertHit(chamber, hit);
        FILE_LOG(logDEBUG) << "      new Hit : ch " << chNum << ", chip " << chipId << ", dif " << difId
                           << ", Time=" << globalTime << " abs=" << bcId_Abs << " dif=" << bcId_Dif << " hit=" << bcId_Hit << " Data=" << data << endl;
        return(1);
      }
      catch (MicroException &e) {
        FILE_LOG(logERROR) << "ch " << chNum << ", chip " << chipId << ", dif " << difId << " : " << e.getMessage() << endl;
        return(0);
      }
    }
    catch (...) {
      FILE_LOG(logDEBUG) << "      board list not found - try next" << endl;
      continue;
    }

  } // for
  FILE_LOG(logERROR) << "      ch " << chNum << ", chip " << chipId << ", dif " << difId << " : geometry element not found" << endl;
  return(0);
} // Hardroc2Reader.newHit

---