CalibMicrorocParser Class Reference

#include <CalibMicrorocParser.hh>

Inheritance diagram for CalibMicrorocParser:

[legend]Collaboration diagram for CalibMicrorocParser:

[legend]List of all members.

Public Member Functions
	CalibMicrorocParser (Run &aRun, FILE *aFile, ui32 firstEventId, bool _testMicroroc=false)
virtual	~CalibMicrorocParser ()
int	getNextEvent (Event &eventToFill)
Static Public Member Functions
static const char *	type ()
Protected Member Functions
int	getLine (FILE *file, std::string &line)
Private Member Functions
unsigned int	getData (FILE *inputFile, const int nBytes, bool computeCRC=true)
unsigned int	getDecimalData (FILE *inputFile, const int nBytes)
unsigned int	checkData (FILE *inputFile, const int nBytes)
int	getDacValue (std::string line)
void	setMotifValue (std::string line)
void	setChargeValue (std::string line)
unsigned int	grayToBinary (unsigned int gray)
int	newHit (Event &event, const Detector &detector, int data, const int chNum, const int chipId, const int difId)
std::string	binaire (unsigned int nombre)
void	reset ()
void	initParsing (FILE *inputFile)
bool	checkNextLine (FILE *inputFile, std::string target)
virtual std::string	getCtestString (void)
Private Attributes
MicrorocSCReader *	scReader
bool	init
bool	motif [MOTIFNB]
float	charge
bool	newCalibration
bool	newDif
bool	newHR
int	difId
int	dac [3]
char	shaper [2]
bool	testMicroroc

---

Detailed Description

CLASS DATAREADER //

Definition at line 32 of file CalibMicrorocParser.hh.

---

Constructor & Destructor Documentation

CalibMicrorocParser::CalibMicrorocParser	(	Run &	aRun,
		FILE *	aFile,
		ui32	firstEventId,
		bool	_testMicroroc = false
	)

Definition at line 55 of file CalibMicrorocParser.cpp.

: AcquisitionParser(aRun, aFile,  _firstEventId), scReader(NULL) 
{

  testMicroroc = _testMicroroc;
  run.setCalibrationRun(true);
  run.setName("Microroc calibration Run");

  reset();
}

CalibMicrorocParser::~CalibMicrorocParser

(

)

[virtual]

Definition at line 67 of file CalibMicrorocParser.cpp.

{
  FILE_LOG(logDEBUG1)<< "----CalibMicrorocParser destructeur"<< endl;
}

---

Member Function Documentation

static const char* CalibMicrorocParser::type

(

)

[inline, static]

Reimplemented in TestMicrorocParser.

Definition at line 38 of file CalibMicrorocParser.hh.

Referenced by main().

{ return CALIBMICRORC; }

int CalibMicrorocParser::getNextEvent

(

Event &

eventToFill

)

[virtual]

Implements AcquisitionParser.

Reimplemented in TestMicrorocParser.

Definition at line 252 of file CalibMicrorocParser.cpp.

Referenced by TestMicrorocParser::getNextEvent().

{
    unsigned int nbHits = 0;
    // The first line does not contains needed informations
    // init fill the motifCtest and place streaming just before the first event
    if ( ! init)
    {
        initParsing(inputFile);
        init = true;
    }

    string line;

    if ( newCalibration)
    { // get calibration parameters ( DAC_2, DAC_1, D AC_0 values )
        FILE_LOG(logDEBUG1) <<  "New calibration !" << endl;
        size_t found;

        do
        {
            int result = getLine ( inputFile, line);
             if (result == EOF)
             {
                 FILE_LOG(logDEBUG1) <<  "End of file !" << endl;
                 reset();
                 FILE_LOG(logDEBUG1) <<  "return LAST_FILE_EVENT:["  << LAST_FILE_EVENT <<  endl;
                 return LAST_FILE_EVENT;
             }

             // Get DAC_2
             found=line.find("DAC_2");
        }
        while (found==string::npos);

        dac[0] = -1;
        dac[1] = -1;
        dac[2] = getDacValue(line);

        // Get DAC_1
        if ( getLine ( inputFile, line) != EOF )
        {
            if ( line.find("DAC_1")!=string::npos) { dac[1] = getDacValue(line); }
            else { return EVENT_ERROR;  }  // DAC_1 value not found
        }

        // Get DAC_0
        if ( getLine ( inputFile, line) != EOF )
        {
            if ( line.find("DAC_0")!=string::npos) { dac[0] = getDacValue(line); }
            else { return EVENT_ERROR;  }  // DAC_0 value not found
        }
        FILE_LOG(logDEBUG) <<  "dac_2[" << dac[2] << "]" << endl;
        FILE_LOG(logDEBUG) <<  "dac_1[" << dac[1] << "]" << endl;
        FILE_LOG(logDEBUG) <<  "dac_0[" << dac[0] << "]" << endl;
        // Set chip configuration
        const Detector &detector = run.getDetector();
        try {
          ChamberMap_t chambers = detector.getChambers();
          for (ChamberMap_t::iterator chamberIt = chambers.begin(); chamberIt != chambers.end(); chamberIt++) {
            Chamber& chamber = *chamberIt->second;
            DifMap_t difs = chamber.getDifs();
            for (DifMap_t::iterator difIt = difs.begin(); difIt != difs.end(); difIt++) {
              Dif& dif = *difIt->second;
              BoardMap_t boards = dif.getBoards();
              for (BoardMap_t::iterator boardIt = boards.begin(); boardIt != boards.end(); boardIt++) {
                Board& b = *boardIt->second;
                ChipMap_t chips = b.getChips();
                for (ChipMap_t::iterator chipIt = chips.begin(); chipIt != chips.end(); chipIt++) {
                  MicrorocChip &chip = dynamic_cast<MicrorocChip &> (*chipIt->second);
                  chip.setThresholdDac_0(dac[0]);
                  chip.setThresholdDac_1(dac[1]);
                  chip.setThresholdDac_2(dac[2]);
                  //for (int chNum = 0; chNum < 64; chNum++) {
                  //  Channel& ch = chip.getChannelById(chNum);
                  //  ch.setEnable(1);
                  //}
                }
              }
            }
          }
        }
        catch (MicroException &e) { cout << e.getMessage() << endl; }

        newCalibration = false; // nany event could be attach to these calibration parametre ( DAC_2 ... )
        return NEW_CONFIG;
     } // end if new calibratio
    if ( newDif)
    {
      FILE_LOG(logDEBUG1) <<  "New Dif !" << endl;
      // Place file streamer after B4
      int space = 0;  // nb bytes beetween B4 and B0 tags
      int b0Pos = 6;  // B0 tag position
      bool find = false;
      do
      {
        b0Pos = ftell(inputFile);

        FILE_LOG(logDEBUG) <<  "inputFile [" << inputFile << "]" << endl;
        FILE_LOG(logDEBUG) <<  "b0Pos [" << b0Pos << "]" << endl;
        if (getLine ( inputFile, line) == EOF) {
          return(LAST_FILE_EVENT);
        }
      //FILE_LOG(logINFO) << "line[" << line << "]" << endl;
        // On verifie que l'on ne retombe pas sur le mot DAC_2, ce qui ce produit si aucune Dif n'a repondu ( du a un BUG )

        if (checkNextLine(inputFile,"DAC_") )
        {
          newCalibration = true;
          FILE_LOG(logWARNING) <<  "No data for DAC values:[" << dac[0] << "," << dac[1] << "," << dac[2] << "]" << endl;
          return getNextEvent( event);
        }

        int b0find = line.find("B0",18);
        FILE_LOG(logDEBUG) <<  "b0find [" << b0find << "]" << endl;
        int b4pos = 0;
        do
        {
          b4pos = line.find("B4",b4pos+1);
          FILE_LOG(logDEBUG) <<  "b4pos [" << b4pos << "]" << endl;
          FILE_LOG(logDEBUG) <<  "b4pos - b0find [" << b4pos - b0find << "]" << endl;

          ui16 gap = 48;
          
          if ( testMicroroc == true ) { gap = 46 ;}
    
          if ( b4pos - b0find == gap )
          {
            find = true;
          }
        }
        while  (b4pos != string::npos && find == false );
      }
      while ( !find );

      b0Pos = b0Pos + 18 + 2 ;  // position du B0
      FILE_LOG(logDEBUG) <<  " b0Pos = b0Pos + 18 + 2 [" << b0Pos << "]" << endl;
      fseek(inputFile, b0Pos, SEEK_SET);

      // File after position "B0"
      // Here we are sure to be at the correct position in the file to start

      this->setComputeCrcFlag(true);
      
      FILE_LOG(logDEBUG) << "crc->compute(0xB0)[ "<< crc <<"]"   << endl;
      crc->compute(0xB0);

      // New DIF
      difId = getData(inputFile, 1);
      if ( scReader   ) { scReader->parseSC(difId); } 
      FILE_LOG(logDEBUG) << "!!!  difId[0x"<< hex  << difId << "]" <<  endl;

      //fseek(inputFile, 42, SEEK_CUR);
      //fseek(inputFile, 44, SEEK_CUR);
      ui16 skip = 22;
      if ( testMicroroc == true ) { skip = 21 ;}
      for ( int i = 0 ; i < skip ; i++) {
        int foo = getData(inputFile,1);
//      FILE_LOG(logINFO) << (hex) << "  foo[" << foo << "]" <<  endl;
      }
      // Fichier avant la position "B4"
      // New B4 Tag

      int dataHeader = getData(inputFile, 1);
      FILE_LOG(logDEBUG1) << (hex) << "  dataHeader[" << dataHeader << "]" <<  endl;

      if (dataHeader == 0xB4)
      {
        FILE_LOG(logDEBUG1) << "  data header B4 found" << endl;
      }
      else
      {
         FILE_LOG(logERROR) << "  parsing file error. data header error : " << hex << dataHeader << dec << endl;
         return EVENT_ERROR;
      }

      newDif = false;
    }



    if ( newHR)
    {
          FILE_LOG(logDEBUG) <<  "New MR !" << endl;

       // VERIF TAG A3 pas de skip. Dans le cas ou aucun des chips ne repond
      //    if ( checkData(inputFile, 1) == 0xA3)
       //   {
        //     FILE_LOG(logERROR) <<  "ERROR TAG A3 just after newHR for DAC values:[" << dac[0] << "," << dac[1] << "," << dac[2] << "]" << endl;
         //    return EVENT_ERROR;
         // }
       // On skip les insultes du microroc
       //fseek ( inputFile, 40, SEEK_CUR);
       newHR = false;
    }


    // On Verifie en hexadecimal que l'on a pas de tag A0
    // Pour une Dif, lit ChipId + BCId + donnee
    // get Chip Id. Warning ! Decimal

    int chipId = -1;
    // On doit lire 2 fois car le tag A3 est en Hexedecimal alors que le chipId est en decimal
    int foo = getData(inputFile , 1);
    if ( foo == 0xA3)
    { // Fin de donnees chip(HR)
      FILE_LOG(logDEBUG) << " -------- TAG A3  force newHR to true " <<  endl;
      newHR = true;
      foo = getData(inputFile , 1);
      if ( foo == 0xA0 )
      {
        newDif = true;
        FILE_LOG(logDEBUG1) << " -------- Tag A0 derriere A3 " <<  endl;
        // Stop CRC computation 

        this->setComputeCrcFlag(false);
        unsigned short CRC_check = getData(inputFile , 2 ); // READ CRC
        FILE_LOG(logDEBUG1) << hex << " Compute CRC[" << crc->getCrc() << "]" <<  endl;
        FILE_LOG(logDEBUG1) << hex << " Received CRC[" << CRC_check << "]" <<  endl;
        if ( CRC_check != crc->getCrc() )
        {
          FILE_LOG(logERROR) << " CRC Error " << endl;
          exit(-1);
        }
        
        try {
          foo = getData(inputFile, 1 );
        }
        catch ( MicroException &e) 
        {
          foo = 0;
        }
        if ( foo != 0xB0)
        {
          // On verifie si apres la ligne separatrice  "-----------------"c'est une nouvelle acquisition on des nouvelles valeur de DAC
          if (  getLine ( inputFile, line) == EOF ) { return LAST_FILE_EVENT ;}  // ligne separatrice ---------------
          if (checkNextLine(inputFile,"acquisition") )
          {
            newDif = true;
          }
          else 
          {
            newCalibration = true;
            FILE_LOG(logDEBUG1) << " -------- Pas de B0 derrire le CRC donc newCalibration = true " <<  endl;
          }
        }
        else 
        {
          FILE_LOG(logDEBUG1) << " -------- B0 derrire le CRC donc Dif  " <<  endl;
        }
      }
      return getNextEvent( event);
    }
    else
    { // Pas Fin de donnees chip(HR)
      fseek ( inputFile, -2, SEEK_CUR);
      chipId = getData(inputFile, 1, false);


    // get bcid
    int bcId = getData(inputFile, 3);
    FILE_LOG(logDEBUG1) << (hex) << "  bcid[" << bcId << "]" <<  endl;


    // the following data is repeated one time for each Dif which has sent data

    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;
        return EVENT_ERROR;
      }
    }

   // 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

        int canal[16]={0};
        for (int chNum = 48; chNum >= 0; chNum -= 16) {
                unsigned int data4;
          try {
              data4 = getData(inputFile, 4); // on lit 4 Bytes => 8 cara => 16 cannaux
          }
          catch ( MicroException &e) 
          {
              continue;
          }


          int index = 0;
          for ( index = 0; index <= 15 ; index++)
          {
            // sur les 16 canaux, on recupere les donnees canal par canal
            canal[index] = (data4 >> ((15 - index )* 2 )) & 0x03; // 2 Most Significant Bits of the digi

// MODIF BY RENAUD 9 march 2012 to deal with 2 testboard without changing Labview
//            if ( motif[index+chNum])
            {
//              if (canal[index]!=0){
              {
                if ( newHit(event, detector, canal[index], index + chNum    , chipId, difId) != 1 )
                                        {
// Format de donées non valid, on cherche une balise A3 pour tenter de recuperer les prochain evemements
                  try 
                  {
                        foo =  getData(inputFile, 1) ;
                        while (  foo != 0xA3 )
                                                      {
                          foo =  getData(inputFile, 1) ;
                                                      }
                    }
                                                  catch ( MicroException &e) 
                    {
                      FILE_LOG(logINFO) << "End of line 0xA3 not Found !!!" << endl; 
                                                          newDif = true;
//                    event.setId(++lastEventId);
                    //FILE_LOG(logDEBUG1) << "  STORE EVENT " << lastEventId << ", " << nbHits << " hits" << endl;
                    //return EVENT_CORRECT;
                    }
                                        }
                                        else 
                                        {//, dac);
                        nbHits++;
                                }}
               }
          }
      }

      event.setId(++lastEventId);
      FILE_LOG(logDEBUG1) << "  STORE EVENT " << lastEventId << ", " << nbHits << " hits" << endl;
      return EVENT_CORRECT;
    }

}

unsigned int CalibMicrorocParser::getData	(	FILE *	inputFile,
		const int	nBytes,
		bool	computeCRC = true
	)			[private]

Definition at line 632 of file CalibMicrorocParser.cpp.

Referenced by checkData(), and getNextEvent().

                                                                                            {
  unsigned int result = 0;
  for (int byte = 0; byte < nBytes; ++byte) {
    int data = 0;
    if (fscanf(inputFile, "%2x", &data) != 1) {
            throw MicroException("get data: fcanf error");

//      return(0);
    }
    result = (result << 8) + data;
    if ( computeCrc && computeCRC) {
      crc->compute(data) ;
    }

  }
  return(result);
} 

unsigned int CalibMicrorocParser::getDecimalData	(	FILE *	inputFile,
		const int	nBytes
	)			[private]

Definition at line 653 of file CalibMicrorocParser.cpp.

                                                                                  {
  unsigned int result = 0;
  for (int byte = 0; byte < nBytes; ++byte) {
    int data = 0;
    if (fscanf(inputFile, "%2d", &data) != 1) {
      return(0);
    }
    result = (result << 8) + data;
  }
  return(result);
} 

unsigned int CalibMicrorocParser::checkData	(	FILE *	inputFile,
		const int	nBytes
	)			[private]

Definition at line 622 of file CalibMicrorocParser.cpp.

                                                                             {

    int curPos = ftell(inputFile);
    unsigned int result = getData(inputFile, nBytes);
    fseek(inputFile,curPos,SEEK_SET);
    return result;
}

int CalibMicrorocParser::getDacValue

(

std::string

line

)

[private]

Referenced by getNextEvent().

void CalibMicrorocParser::setMotifValue

(

std::string

line

)

[private]

Referenced by initParsing().

void CalibMicrorocParser::setChargeValue

(

std::string

line

)

[private]

Referenced by initParsing().

unsigned int CalibMicrorocParser::grayToBinary

(

unsigned int

gray

)

[private]

Reimplemented from AcquisitionParser.

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

Definition at line 668 of file CalibMicrorocParser.cpp.

Referenced by getNextEvent().

{

  if ( data != 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);
          i16 analogValue = 0;
          //ChannelHit *hit = new ChannelHit(channel, data, analogValue, dac);
          ChannelHit *hit = new ChannelHit(channel, data, analogValue);
          event.insertHit(chamber, hit);
          FILE_LOG(logDEBUG1) << "      new Hit : ch " << chNum << ", chip " << chipId << ", dif " << difId
                            << " Data=" << data << " , x:" << channel.getX() << ", y:" << channel.getY() << endl;
          return(1);
        }
        catch (MicroException &e) {
          FILE_LOG(logERROR) << "ch " << chNum << ", chip " << chipId << ", dif " << difId << " : " << e.getMessage() << endl;
          return(0);
        }
      }
      catch (...) {
        continue;
      }

    } // for
  
    return(0);
  }
  else
  { 
    return 1;
  }
} 

string CalibMicrorocParser::binaire

(

unsigned int

nombre

)

[private]

Definition at line 715 of file CalibMicrorocParser.cpp.

{
string res = "";
unsigned int val = nombre;

while (val != 0) {
/* Invariant de boucle : code_binaire(val) + res = code_binaire(nombre) */

   if (val & 1 == 1) {
res = string("1") + res;
}
    else {
/* val & 1 == 0 */
res = string("0") + res;
}

val = val >> 1;
}
return res;
}

void CalibMicrorocParser::reset

(

)

[private]

Definition at line 73 of file CalibMicrorocParser.cpp.

Referenced by CalibMicrorocParser(), and getNextEvent().

{
init = false;
newCalibration = true;
newDif = true;
newHR = true;
FILE_LOG(logDEBUG) << " -------- reset  force newHR to true " <<  endl;
difId = -1;
dac[0] = -1;
dac[1] = -1;
dac[2] = -1;
shaper[0] = -1;
shaper[1] = -1;
difId = -1;
charge = 0.;

    for ( int i = 0; i < MOTIFNB; i++)
    {
        motif[i] = 0;
        }


}

void CalibMicrorocParser::initParsing

(

FILE *

inputFile

)

[private]

Definition at line 169 of file CalibMicrorocParser.cpp.

Referenced by getNextEvent().

{

    string line;
    size_t found;
    // search charge (fC)
    do
    {
      if ( getLine ( inputFile, line) == EOF ) { return ;}
      found = line.find("charge (fC)" );
      if  (found!=string::npos)
      { 
        setChargeValue(line);
      }
      line.clear();
    }
    while  (found==string::npos);
    do
    {
       if ( getLine ( inputFile, line) == EOF ) { return ;}
       found = line.find("Chips en calibration" );
       line.clear();
    }
    while  (found==string::npos);

    // get motifCTest
    if ( getLine ( inputFile, line) == EOF ) { return ;}
    found = line.find(getCtestString() );
    if  (found!=string::npos)
    {
      setMotifValue(line);
    }
        else {
            found = line.find(CalibMicrorocParser::getCtestString() );
        if  (found!=string::npos)
        {
          setMotifValue(line);
        }
        }


// Find slow control data
   do
    {
       if ( getLine ( inputFile, line) == EOF ) { return ;}
       found = line.find("donnees Slow-Control" );
       line.clear();
    }
    while  (found==string::npos);
     // "donnees Slow-Control"  found get value in following line
    if ( getLine ( inputFile, line) == EOF ) { return ;}
        // supression du retour chariot
     line = line.substr(0, line.size());
//difId=80;
  //run.getDetector().getChambers().begin().first.getDifs().begin().fisrt.getId();
  scReader = new MicrorocSCReader(run,line,difId); 

}

bool CalibMicrorocParser::checkNextLine	(	FILE *	inputFile,
		std::string	target
	)			[private]

Referenced by getNextEvent().

virtual std::string CalibMicrorocParser::getCtestString

(

void

)

[inline, private, virtual]

Reimplemented in TestMicrorocParser.

Definition at line 58 of file CalibMicrorocParser.hh.

Referenced by initParsing().

{return "motif Ctest";} ;

int CalibMicrorocParser::getLine	(	FILE *	file,
		std::string &	line
	)			[protected]

Referenced by TestMicrorocParser::getNextEvent(), getNextEvent(), and initParsing().

---

Member Data Documentation

MicrorocSCReader* CalibMicrorocParser::scReader [private]

Definition at line 58 of file CalibMicrorocParser.hh.

Referenced by getNextEvent(), and initParsing().

bool CalibMicrorocParser::init [private]

Definition at line 63 of file CalibMicrorocParser.hh.

Referenced by getNextEvent(), and reset().

bool CalibMicrorocParser::motif[MOTIFNB] [private]

Definition at line 64 of file CalibMicrorocParser.hh.

Referenced by reset().

float CalibMicrorocParser::charge [private]

Definition at line 65 of file CalibMicrorocParser.hh.

Referenced by reset().

bool CalibMicrorocParser::newCalibration [private]

Definition at line 66 of file CalibMicrorocParser.hh.

Referenced by getNextEvent(), and reset().

bool CalibMicrorocParser::newDif [private]

Definition at line 67 of file CalibMicrorocParser.hh.

Referenced by getNextEvent(), and reset().

bool CalibMicrorocParser::newHR [private]

Definition at line 68 of file CalibMicrorocParser.hh.

Referenced by getNextEvent(), and reset().

int CalibMicrorocParser::difId [private]

Definition at line 69 of file CalibMicrorocParser.hh.

Referenced by getNextEvent(), initParsing(), and reset().

int CalibMicrorocParser::dac[3] [private]

Definition at line 70 of file CalibMicrorocParser.hh.

Referenced by getNextEvent(), and reset().

char CalibMicrorocParser::shaper[2] [private]

Definition at line 71 of file CalibMicrorocParser.hh.

Referenced by reset().

bool CalibMicrorocParser::testMicroroc [private]

Definition at line 72 of file CalibMicrorocParser.hh.

Referenced by CalibMicrorocParser(), and getNextEvent().

---

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

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

---