Hardroc1Reader Class Reference

#include <Hardroc1Reader.hh>

Inheritance diagram for Hardroc1Reader:

[legend]Collaboration diagram for Hardroc1Reader:

[legend]List of all members.

Public Member Functions
	Hardroc1Reader (Run &aRun, FILE *aFile, ui32 lastEventId)
	~Hardroc1Reader ()
int	getNextEvent (Event &eventToFill)
Static Public Member Functions
static const char *	type (const int subtype=0)
Protected Types
typedef std::set< int >	SkippedChipSet_t
Protected Attributes
unsigned short	CRC
int	difId
ui32	difTriggerCounter
ui32	acqTriggerCounter
ui32	globalTriggerCounter
ui32	bcPeriod
i64	bcId_Abs
i64	bcId_AbsOrg
i64	bcId_Dif
i64	bcId_Hit
i64	eventTime
SkippedChipSet_t	skippedChips
int	SkippedBcId
Private Member Functions
void	reset ()
int	getScData ()
int	getAcqData (Event &event)
int	storeEvent (Event &event)
int	newHit (Event &event, const Detector &detector, int data, const int chNum, const int chipId, const int difId)

---

Detailed Description

CLASS DATAREADER //

Definition at line 26 of file Hardroc1Reader.hh.

---

Member Typedef Documentation

typedef std::set<int> Hardroc1Reader::SkippedChipSet_t [protected]

Definition at line 59 of file Hardroc1Reader.hh.

---

Constructor & Destructor Documentation

Hardroc1Reader::Hardroc1Reader	(	Run &	aRun,
		FILE *	aFile,
		ui32	lastEventId
	)

Definition at line 53 of file Hardroc1Reader.cpp.

                                                                       : DifReader(aRun, aFile, lastEventId)
{
  FILE_LOG(logDEBUG1) << "--Hardroc1Reader constructor()"<< endl;
  reset();
}

Hardroc1Reader::~Hardroc1Reader

(

)

Definition at line 60 of file Hardroc1Reader.cpp.

{
  FILE_LOG(logDEBUG1)<< "----Hardroc1Reader destructor"<< endl;
}

---

Member Function Documentation

static const char* Hardroc1Reader::type

(

const int

subtype = 0

)

[inline, static]

Definition at line 30 of file Hardroc1Reader.hh.

Referenced by main().

{ return CROSSDAQHR1; };

int Hardroc1Reader::getNextEvent

(

Event &

eventToFill

)

[virtual]

Implements AcquisitionParser.

Definition at line 471 of file Hardroc1Reader.cpp.

                                             {
  // Start data reconstruction
  Level(0);

        ui16 dataFormat = run.getDataFormat();

  if (dataFormat == 0) {
    // first read of the file
    bcIdDifAbsOrgMap.clear();
    // read file header
    dataFormat = getData(1);
                run.setDataFormat(dataFormat);
    if (fileError == EOF) {
      FILE_LOG(logERROR) << "parsing file[" << inputFile << "]  - error " << errno << ":" << strerror(errno) << endl;
      reset();
      return(LAST_FILE_EVENT);
    }
    display(+1);
    FILE_LOG(logINFO) << "  HARDROC1 Dif VHDL data format : " << dataFormat << endl;

    if (dataFormat != 2 && dataFormat != 3 && dataFormat != 4 && dataFormat != 7 && dataFormat != 8) {
      // unknown format
      FILE_LOG(logERROR) << " HARDROC1 Dif VHDL data format " << dataFormat << " is not implemented" << endl;
      reset();
      return(LAST_FILE_EVENT);
    }

    FILE_LOG(logDEBUG) << run.getDetector() << endl;
    display(-1);
    return NEW_CONFIG;
  }

  // empty the list of chips for which first line has to be skipped
  skippedChips.clear();

  // loops until a non-empty event has been found
  for (unsigned int nbHits = 0; !nbHits; ) {

    if (dataFormat == 2 || dataFormat == 3 || dataFormat == 4) {
      ui32 timestampTime = getData(4);
      ui32 timestampMilliseconds = getData(4);
      globalTime = timestampTime * 1000 + timestampMilliseconds;
    }
    else if (dataFormat == 7) {
      // read in Slow Control data
    }
    else if (dataFormat == 8) {
      ui32 timestampTime = getData(4);
      globalTime = timestampTime * 1000;
    }

    // look for a global data header and process corresponding data
    unsigned int globalHeader = getData(1);
    if (fileError == EOF) {
      FILE_LOG(logINFO) << "  END OF DATA FILE" << endl;
      reset();
      return(LAST_FILE_EVENT);
    }

    if (globalHeader == 0xB1) {
      // slow control data
      getScData();
    }
    else if (globalHeader == 0xB0) {
      // acquisition data
      nbHits = getAcqData(event);
    }
    else {
      display();
      FILE_LOG(logERROR) << "  parsing file error. global header error : " << hex << globalHeader << dec << endl;
    }
  } // for (nbHits)

  return EVENT_CORRECT;
}

void Hardroc1Reader::reset

(

)

[private]

Definition at line 66 of file Hardroc1Reader.cpp.

Referenced by getNextEvent(), and Hardroc1Reader().

                           {
  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;
} // Hardroc1Reader.reset

int Hardroc1Reader::getScData

(

)

[private]

Definition at line 81 of file Hardroc1Reader.cpp.

Referenced by getNextEvent().

                              {
  unsigned int data;

  display(+1);
  FILE_LOG(logDEBUG) << " Slow Control data header B1 found" << endl;
        ui16 dataFormat = run.getDataFormat();

  if (dataFormat == 7) {
    // get time stamp
    i64 timestampTime = getData(4);
    globalTime = timestampTime * 1000;
    display();
    FILE_LOG(logDEBUG) << "   timestampTime = 0x" << hex << timestampTime << " = " << dec << timestampTime << endl;
  }

  int scBoardNumber = getData(1);
  display();
 FILE_LOG(logDEBUG) << "NB chip configuration[" << scBoardNumber << "]" << endl;

  // skip al data
//  for (int dataNum = 0; ; ) {
//    data = getData(1);  ++dataNum;
//    if ((dataNum % 72) == 0)


  for ( int nbChip = 0 ; nbChip < scBoardNumber ; nbChip++)
  {
      unsigned char *buf = initGetBits(72, false);
      if (!buf)
      {
        FILE_LOG(logERROR) << "  malloc error while getting SC buffer" << endl;
        break;
      }
      FILE_LOG(logDEBUG) << "------------------- Config chip [" << nbChip << "]" << endl;
      ui16 en = getBits(4);
      //FILE_LOG(logDEBUG) << "En 1-4[" << BinOut(en,4) <<"]" << endl;
      ui16 chipId = getBits(8);
      FILE_LOG(logDEBUG) << "header chip id: 5-12 [" << BinOut(chipId,8) <<"]=" << dec<< chipId  << "]" <<  endl;
      ui16 bypass_chip = getBits(1);
      //FILE_LOG(logDEBUG) << "bypass_chip 13 [ " << BinOut(bypass_chip,1) <<"]" << endl;

      ui16 en2 = getBits(6);
      //FILE_LOG(logDEBUG) << "En2 14-19 [" << BinOut(en2,6) <<"]" << endl;

      getBits(1);

      ui32 validTrig = getBits(32);
      //FILE_LOG(logDEBUG) << "validTrig 21-52[" << BinOut(validTrig,32) <<"]" << endl;
      validTrig = getBits(32);
     // FILE_LOG(logDEBUG) << "validTrig 53-84[" << BinOut(validTrig,32) <<"]" << endl;

      ui32 bzero = getBits(10);
      FILE_LOG(logDEBUG) << "b0 85-94[" << BinOut(bzero) <<"]"  << dec<< bzero << "]" << endl;


      ui32 bun = getBits(10);
      FILE_LOG(logDEBUG) << "b1 95-104[" << BinOut(bun,20) <<"]" << dec<< bun  << "]" <<  endl;

      ui32 on = getBits(3);
      //FILE_LOG(logDEBUG) << "on 105-107[" << BinOut(on,3) <<"]" << endl;

      ui32 ctest = getBits(32);
      //FILE_LOG(logDEBUG) << "ctest 108-139[" << BinOut(ctest,32) <<"]" << endl;
      ctest = getBits(32);
      //FILE_LOG(logDEBUG) << "ctest 140 - 171[" << BinOut(ctest,32) <<"]" << endl;

      for ( int chIndex = 0 ; chIndex <=63 ; chIndex++ )
      {
          ui32 pream = getBits(6);
                                //FILE_LOG(logDEBUG) << "preamp ch"<< chIndex << " " << 172+(chIndex*6) << "-" << 172+(chIndex*6)+ 6 << "["<< BinOut(pream,6) <<"]" << endl;
                                FILE_LOG(logDEBUG) << "gain channel"<< chIndex <<  "["<< BinOut(pream,6) <<"]=" << dec << pream << endl;
      }

      on = getBits(7);
      //FILE_LOG(logDEBUG) << "on 556 - 562  " << BinOut(on,7) <<"]" << endl;

      ui16 valid_dc = getBits(1);
      //FILE_LOG(logDEBUG) << "valid_dc  563[" << BinOut(valid_dc,1) <<"]" << endl;

      ui16 sw = getBits(4);
      //FILE_LOG(logDEBUG) << "sw 564 - 567 [" << BinOut(sw,4) <<"]" << endl;

      ui16 choix_caisson = getBits(1);
      //FILE_LOG(logDEBUG) << "choix_caisson 568 [" << BinOut(choix_caisson,1) <<"]" << endl;

      sw = getBits(3);
      //FILE_LOG(logDEBUG) << "sw  569 - 571[" << BinOut(sw,3) <<"]" << endl;

      getBits(4);


      finishGetBits();



      // Store config information to Hardroc1Chip
      
      const Detector &detector = run.getDetector();
      try
      {
          Dif& dif = detector.getDifById(difId);
          try {
                  Hardroc1Chip& chip = dynamic_cast<Hardroc1Chip &> (dif.getChipById(chipId));
                  FILE_LOG(logDEBUG1) << " Hardroc1 slow control for chip " << chipId << " et Dif "  << difId << "  configuration set: " <<    endl;
                  chip.setThresholdDac_0(bzero);
                  chip.setThresholdDac_1(bun);
          }
          catch ( ...  )
          {
          }
      }
      catch ( MicroException e )
      {
          FILE_LOG(logINFO) << " Hardroc1 slow control for chip " << chipId << " et Dif "  << difId << "  configuration: " <<   e.getMessage() << endl;

      }

  }

    while ( getData(1) != 0xA1)
    {
      FILE_LOG(logERROR) << "0xA1 trailer not found" << endl;
    }
  data = 0xA1;
  display();
  //  if (data == 0xA1)
 //     break;
  //}


  display(-1);
  if (data != 0xA1)
    FILE_LOG(logERROR) << "  parsing file error : Slow control trailer 0xA1 not found [" << hex << data << dec << "]" << endl;

  return(1);
} // Hardroc1Reader.getScData

int Hardroc1Reader::getAcqData

(

Event &

event

)

[private]

Definition at line 220 of file Hardroc1Reader.cpp.

Referenced by getNextEvent().

                                           {
  unsigned int nbHits = 0;

  ui16 dataFormat = run.getDataFormat();
  display(+1);
  FILE_LOG(logDEBUG) << " global header B0 found" << endl;

  // reset the data for the block
//  bcId_AbsOrg = -1;
  bcId_Abs = 0; // max. 48 bits
  bcId_Dif = 0;
  // the following data is repeated one time for each Dif which has sent data
  difId = getData(1);
  difTriggerCounter = getData(4);
  if (dataFormat == 4 || dataFormat == 7 || dataFormat == 8)
    acqTriggerCounter = getData(4);  // trigger counter version >= 4
  globalTriggerCounter = getData(4);
  if (dataFormat == 3 || dataFormat == 4) {
    bcId_Abs = getData(3);
  }
  else if (dataFormat == 7 || dataFormat == 8) {
    ui64 bcIdHigh = getData(2);
    ui64 bcIdLow = getData(4);
    bcId_Abs = (bcIdHigh << 32) + bcIdLow;
    bcId_Dif = getData(3);
  }

//  if (bcId_AbsOrg == -1) {
//    // set the absolute bcId origin
//    bcId_AbsOrg = bcId_Abs;
//    FILE_LOG(logDEBUG)   << "    set bc origin         " << bcId_AbsOrg << endl;
//  }
  bcId_AbsOrg = getbcIdAbsOrg(difId, bcId_Abs);
  display();

  // the bcTime in milliseconds
  i64 bcTime = ((bcId_Abs - bcId_AbsOrg) * bcPeriod) / 1000;
  // the reading time in milliseconds
  i64 readTime = globalTime + bcTime;

  FILE_LOG(logDEBUG)   << "  START EVENT ------------- time " << readTime << " ms" << endl;
  FILE_LOG(logDEBUG)   << "    bc time               " << bcTime << " ms" << endl;
  FILE_LOG(logDEBUG)   << "    DIF Id                0x" << hex << (int)difId << dec << endl;
  FILE_LOG(logDEBUG)   << "    DIF trigger counter   " << difTriggerCounter << endl;
  if (dataFormat == 4 || dataFormat == 7 || dataFormat == 8)
    FILE_LOG(logDEBUG) << "    ACQ trigger counter   " << acqTriggerCounter << endl;
  FILE_LOG(logDEBUG)   << "    total trigger counter " << globalTriggerCounter << endl;
  if (dataFormat == 3 || dataFormat == 4 || dataFormat == 7 || dataFormat == 8)
    FILE_LOG(logDEBUG) << "    absolute bcId         " << bcId_Abs - bcId_AbsOrg << endl;
  // end of global header

  bool difFound = false;
  Detector &detector = run.getDetector();
  { 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);
      }
      catch (...) {
        // any error case - try next chamber
        continue;
      }

      // the DIF has been found in a chamber. Process all data and continue
      difFound = true;
      break;
    } // for
    if (!difFound) {
      FILE_LOG(logERROR) << "    No dif found with Id 0x" << hex << difId << dec << endl;
    }
  }

  for (bool finish = false; !finish; ) {
    int delimiter = 0;

    display();
    for (; delimiter != 0xA0 && delimiter != 0xB4; ) {
      // seek for first data frame header (0xB4)
      delimiter = getData(1);
      if (fileError == EOF) {
        delimiter = 0;
        break;
      }
      if (delimiter == 0xA3)
        FILE_LOG(logERROR) << "    frame trailer A3 found -- should not" << endl;
    }
    if (delimiter == 0x00) {
      FILE_LOG(logERROR) << "  end of file -- should not" << endl;
      break;
    }
    if (delimiter == 0xA0) {
      unsigned short CRC_check = getData(2);
      display(-1);
      FILE_LOG(logDEBUG) << " global trailer A0 found" << endl;
      break;
    }
    display(+1);
    FILE_LOG(logDEBUG) << "   frame header B4 found" << endl;

    int chipId = -1;
    for (; ; ) {
      chipId = getData(1);
      if (fileError == EOF) {
        finish = true;
        FILE_LOG(logINFO) << "    end of file found -- should not" << endl;
        break;
      }
      if (chipId == 0x00) {
        // dummy null bytes
        continue;
      }
      else if (chipId == 0xA3) {
        display(-1);
        FILE_LOG(logDEBUG) << "   frame trailer A3 found" << endl;
        break;
      }
      else if (chipId == 0xA0) {
        finish = true;
        display(-1);
        FILE_LOG(logERROR) << "    global trailer A0 found -- should not" << endl;
        int CRC_check = getData(1);
        break;
      }
      else if (chipId == 0xC3) {
        display(-1);
        FILE_LOG(logERROR) << "    frame trailer C3 found - error occured during data transfer" << endl;
        break;
      }

      // bcId from hardroc chip is counted as Gray code
      bcId_Hit = grayToBinary(getData(3));
      // event time in milliseconds
      eventTime = readTime + (((bcId_Dif - bcId_Hit) * bcPeriod) / 1000);

      // for each chip, always skip first line of data which is corrupted
      //   <=> if chip not found in the list then this is the first line
      // also skip lines where bcId is < 2 (most of the time, data is corrupted)
      int logicalChipId = (difId << 8) + chipId;
      SkippedChipSet_t::const_iterator chipIt = skippedChips.find(logicalChipId);
      if (chipIt == skippedChips.end()
       || (bcId_Hit <= 1)) {
        if (chipIt == skippedChips.end()) {
          FILE_LOG(logDEBUG) << "      skip data (first line for chip 0x" << hex << logicalChipId << dec << ")" << endl;
          skippedChips.insert(logicalChipId);
        }
        else {
          FILE_LOG(logDEBUG) << "      skip data (bcIdHit = " << (int)bcId_Hit << ")" << endl;
          SkippedBcId++;
        }

        for (int chNum = 48; chNum >= 0; chNum -= 16) {
          unsigned int data4 = getData(4); // lit 4 bytes => 8 char => 16 canaux
        }
      }
      else {
        // not on first line and bcId OK : get and check data
        if (difFound) {
          FILE_LOG(logDEBUG) << "      chip Id               0x" << hex << (int)chipId << dec << endl;
          FILE_LOG(logDEBUG) << "      bcId Hardroc          " << (int)bcId_Hit << endl;
          FILE_LOG(logDEBUG) << "      event time            " << eventTime << endl;
        }

        // trigger data for 64 channels : 1 bit trigger1, 1 bit trigger0 for each channel from 63 downto 0.
        // this makes 2-bit values which should not give bit 1 without bit 0
        for (int chNum = 48; chNum >= 0; chNum -= 16) {
          unsigned int data4 = getData(4); // lit 4 bytes => 8 char => 16 canaux

          if (difFound) {
            for (int index = 0; index < 16; index++) {
              // sur les 16 canaux, recupere les donnees canal par canal
              int trigger = (data4 >> ((15 - index) * 2)) & 0x03; // 2 bits
              // store new event data if any
              if (trigger) {
                nbHits += newHit(event, detector, trigger, index + chNum, chipId, difId);
              }
            } // for (index)
          }
        } // for (chNum)
      }
      display();
    } // for(;;)

  } // for (;;)
  display(-1);

  if (nbHits) {
    // CREATE EVENT for current Run
    event.setTimeStamp(globalTime + (bcId_AbsOrg * bcPeriod) / 1000);

    event.setId(++lastEventId);
    FILE_LOG(logDEBUG1) << "  STORE EVENT " << lastEventId << ", " << nbHits << " hits" << endl;
  }
  else
    FILE_LOG(logDEBUG) << "  EMPTY EVENT" << endl;


  if (SkippedBcId) {
    FILE_LOG(logWARNING) << "    skipped " << SkippedBcId << " data because of bad bcId" << endl;
    SkippedBcId = 0;
  }

  return(nbHits);
} // Hardroc1Reader.getAcqData

int Hardroc1Reader::storeEvent

(

Event &

event

)

[private]

int Hardroc1Reader::newHit	(	Event &	event,
		const Detector &	detector,
		int	data,
		const int	chNum,
		const int	chipId,
		const int	difId
	)			[private]

Definition at line 426 of file Hardroc1Reader.cpp.

Referenced by getAcqData().

                                                                                                                               {
  // cannot be binary pattern 10
  if (data == 0x2) {
    FILE_LOG(logERROR) << "      CHANNEL DATA = 0x02" << endl;
//    return(0);
  }

  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);
} // Hardroc1Reader.newHit

---

Member Data Documentation

unsigned short Hardroc1Reader::CRC [protected]

Definition at line 45 of file Hardroc1Reader.hh.

Referenced by reset().

int Hardroc1Reader::difId [protected]

Definition at line 48 of file Hardroc1Reader.hh.

Referenced by getAcqData(), getScData(), and reset().

ui32 Hardroc1Reader::difTriggerCounter [protected]

Definition at line 49 of file Hardroc1Reader.hh.

Referenced by getAcqData(), and reset().

ui32 Hardroc1Reader::acqTriggerCounter [protected]

Definition at line 50 of file Hardroc1Reader.hh.

Referenced by getAcqData(), and reset().

ui32 Hardroc1Reader::globalTriggerCounter [protected]

Definition at line 51 of file Hardroc1Reader.hh.

Referenced by getAcqData(), and reset().

ui32 Hardroc1Reader::bcPeriod [protected]

Definition at line 52 of file Hardroc1Reader.hh.

Referenced by getAcqData(), and reset().

i64 Hardroc1Reader::bcId_Abs [protected]

Definition at line 53 of file Hardroc1Reader.hh.

Referenced by getAcqData(), newHit(), and reset().

i64 Hardroc1Reader::bcId_AbsOrg [protected]

Definition at line 54 of file Hardroc1Reader.hh.

Referenced by getAcqData(), and reset().

i64 Hardroc1Reader::bcId_Dif [protected]

Definition at line 55 of file Hardroc1Reader.hh.

Referenced by getAcqData(), newHit(), and reset().

i64 Hardroc1Reader::bcId_Hit [protected]

Definition at line 56 of file Hardroc1Reader.hh.

Referenced by getAcqData(), newHit(), and reset().

i64 Hardroc1Reader::eventTime [protected]

Definition at line 57 of file Hardroc1Reader.hh.

Referenced by getAcqData().

SkippedChipSet_t Hardroc1Reader::skippedChips [protected]

Definition at line 60 of file Hardroc1Reader.hh.

Referenced by getAcqData(), and getNextEvent().

int Hardroc1Reader::SkippedBcId [protected]

Definition at line 61 of file Hardroc1Reader.hh.

Referenced by getAcqData(), and reset().

---

The documentation for this class was generated from the following files:

• /data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/include/parser/Hardroc1Reader.hh
• /data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/parser/Hardroc1Reader.cpp

---