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

Go to the documentation of this file.

/* @version $Revision: 1337 $ * @modifiedby $Author: jacquem $ * @lastmodified $Date: 2011-10-06 10:48:29 +0200 (Thu, 06 Oct 2011) $ */

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

#include "event/Event.hh"

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

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

#include <cstdio>
#include <iostream>
#include <iomanip>



#define bigTolittle Toolbox::bigTolittle


using namespace std;


Centaure::Centaure(Run& aRun, FILE *aFile, ui32 lastEventId) : AcquisitionParser(aRun, aFile, lastEventId)
{
  //FILE_LOG(logDEBUG) << "--Centaure constructor()"<< endl;
    newConfig =true;
                globalTriggerCounter = 1;
}

//--------------- Destructeur -------------//
Centaure::~Centaure()
{
  //FILE_LOG(logDEBUG)<< "----Centaure destructeur"<< endl;
}

int  Centaure::getNextEvent(Event& event)
{
  if ( newConfig == true)
  {
      newConfig = false;
      return NEW_CONFIG;
  }
  Detector& aDet = run.getDetector();

// --- totalGassiplex96 = TOTAL Nb of Gassiplex Cards ---
ui32 totalGassiplex96 = 0;


totalGassiplex96 = aDet.getNumberOfGaziplex96();
//FILE_LOG(logDEBUG) << "TOTAL Nb of Gassiplex Cards  " << totalGassiplex96 << endl;


// Data File Variables
//ui16 readValue = 0;    // 16bit equivalent UShort_t? yes
ui16 readValue = 0;          // 16bit equivalent UShort_t? yes
ui16 eventLength = 0;
ui32 centaureEvtId =0;     // 32 bits = 4 octets
ui64  value = 0;        // 64bit? no it seems to be 32 bits as well
//ULong64_t  timestamp = 0;        //this one is 64bit!! ;)
ui64  timestamp = 0;        //this one is 64bit!! ;)

// Tree Variables
i32 eventId = 0;

ui32 nbHiChannels[totalGassiplex96]; // number of channels above VME threshold
i32 nEvents  = lastEventId ;    // number of events  in DATA Files
double GoodRate = 0.0;  // (number of entries in Tree)/(number of events  in DATA Files)
ui32 checklength=0;        // expected event length

// Start data reconstruction

  // LOOP over the DATA Files
bool finish = false;

  checklength=0;
  // get CENTAURE event ID
  fread(&readValue,sizeof(i16),1,inputFile);  // char -> 8 bits
  centaureEvtId = bigTolittle(readValue);//// WATCH OUT EVENTID CAN START AT N != 0  !!!!
  nEvents++;
  // eventId = COUTING event ID
  eventId = nEvents;
  // get event Length
  fread(&readValue,sizeof(i16),1,inputFile);
  eventLength =  bigTolittle(readValue);

  if ( centaureEvtId == 0 && eventLength == 0 ) // Last event code
  {
    return LAST_FILE_EVENT;
  }

  else  // Start a new event
    {
        //FILE_LOG(logDEBUG) << "START EVENT --------------------  "  <<endl;
        //FILE_LOG(logDEBUG) << "  centaureEvtId  " << centaureEvtId <<endl;
        //FILE_LOG(logDEBUG) << "  nEvents    " << nEvents <<endl;
        //FILE_LOG(logDEBUG) << "  eventId     " << eventId <<endl;
        //FILE_LOG(logDEBUG) << "  event length " << eventLength << endl;

//    const map<ui16, Chamber*> &chambers = aDet.getChambers();
//    for (map<ui16, Chamber*>::const_iterator iter=chambers.begin(); iter!=chambers.end(); iter++)
    const ChamberMap_t &chambers = aDet.getChambers();
    for (ChamberMap_t::const_iterator iter=chambers.begin(); iter!=chambers.end(); iter++)
    {
         Chamber& chamber = *(*iter).second;

         // LOOP over number of Gassiplex96 per chamber
         for ( ui32 jj = 0; jj < chamber.getNumberOfGassiplexBoards(); jj++ )
            {
                // Get number of channels above VME threshold for each Gassiplex card :
                fread(&readValue,sizeof(i16),1,inputFile);  // word count
                nbHiChannels[jj] = 1023 & bigTolittle(readValue); // get rid of MSB, 1023=0x3FF
                checklength += nbHiChannels[jj];

                 //FILE_LOG(logINFO)<< "nbHiChannels["<<jj<<"] : " << nbHiChannels[jj] << " HIT for chamber id["<< chamber.getId() << " ,pos: "<< chamber.getZ() <<endl;


                // LOOP over number of Channels above VME threshold for each Gassiplex card
                for ( ui32 j = 0 ; j < nbHiChannels[jj] ; j++  )
                {
                  // Get Channel number and its corresponding channel value
                  //Get  bit 16  to 31
                  fread(&readValue,sizeof(i16),1,inputFile);
                  readValue = bigTolittle(readValue);
                  // fill 32 bit data with  bit 16  to 31
                  value = ( readValue << 16)  ;
                  // Get  bit 0 a 15
                  fread(&readValue,sizeof(i16),1,inputFile);
                  readValue = bigTolittle(readValue);

                  // fill 32 bit data with  bit 0  to 15
                  value = value + readValue;

                  ui16 channelData = value & Toolbox::MASK_CHANNEL_DATA;
                  ui16 channelNumber = ( ( value & Toolbox::MASK_CHANNEL_NUMBER) >> 12 ) & Toolbox::MASK_12BIT;

                //FILE_LOG(logDEBUG) << "ChannelNumber[" << channelNumber << "]type:  ["<<  chamber.getType() << endl;
                if ( chamber.getType() == GASSIPLEXCHAMBER4 )
                {
                  channelNumber = channelNumber + 96*jj ;
                }

            i16 channelValue = channelData;  // to have the negative values as well
              if (channelValue>2000)
                {
                  channelValue=channelValue-4096;
                }
              //chmax = max(chmax,channelValue);
              i16 channelNumber2 = channelNumber ;
//              i32 posX, posY, posZ = 0;
  

                    try
                    {
                        const Channel&  channel = chamber.getChannelById(channelNumber);
 //                       posX = channel.getX();
  //                      posY = channel.getY();
   //                     posZ=chamber.getZ(); //distance between chambers, should be entered by steering file

                        i32 digitalValue = 0;
                                                                                                ChannelHit *chHit =new ChannelHit(channel,digitalValue, channelValue) ;
                        event.insertHit(chamber,chHit);
                    }
                      catch ( MicroException &e)
                    {
                        FILE_LOG(logERROR) << "---------------" << "Chamber->id=["<< chamber.getId() << "], channelNuber[" << channelNumber << "]"  << e.getMessage() << endl;
                    }

                }//end of LOOP over number of Channels above VME threshold
            } // end of LOOP over number of GassiplexChamber1 cards in current chamber
    } // end of loop over chamber


    checklength=checklength*2+3+4+totalGassiplex96;
    FILE_LOG(logDEBUG) <<"  checklength: " << checklength<< endl;
    FILE_LOG(logDEBUG) <<"  eventLength: " << eventLength<< endl;
    FILE_LOG(logDEBUG) <<"  totalGassiplex96: " << totalGassiplex96<< endl;
    //FILE_LOG(logDEBUG) << "END EVENT ------------------  "  <<endl;
    if (checklength != eventLength)
          {
            FILE_LOG(logERROR) << "CRAP EVENT checklength != eventLength: NOTHING CAN BE DONE!" << endl;
            FILE_LOG(logERROR) << "totalGassiplex96[" << totalGassiplex96 << "]" << endl;
          }
          else
          {
          //now read the timestamp and store it via event.SetTimestamp(timestamp)
            fread(&readValue,sizeof(i16),1,inputFile);
            readValue = bigTolittle(readValue);
            timestamp = ((i64) readValue<<48);
            fread(&readValue,sizeof(i16),1,inputFile);
            readValue = bigTolittle(readValue);
            timestamp =timestamp+ ((i64) readValue<<32);
            fread(&readValue,sizeof(i16),1,inputFile);
            readValue = bigTolittle(readValue);
            timestamp =timestamp+ ((i64) readValue<<16);
            fread(&readValue,sizeof(i16),1,inputFile);
            readValue = bigTolittle(readValue);
            timestamp =timestamp+(i64) readValue;
            timestamp = timestamp - TIMEREF_OFFSET  ;  // TIMEREF_OFFSET diff from 1904 to 1970

          }

          // Read end event code  0xffff
          fread(&readValue,sizeof(i16),1,inputFile);
          if ( readValue != 0xffff )
          {
              FILE_LOG(logDEBUG) << "ERROR last code event != 0xffff: "<< readValue << endl;
              return EVENT_ERROR;
          }
          else
          {
              FILE_LOG(logDEBUG) <<   "last code event == 0xffff -> reconstruction successfull"<< endl;
              // CREATE EVENT  for current Run
              event.setTimeStamp(timestamp);
              event.setId(eventId);
                                                        event.setGlobalTriggerCounter(globalTriggerCounter);
                                                        event.setDifTriggerCounter(globalTriggerCounter++);


          }
     } //end of event
     lastEventId = eventId;
     return EVENT_CORRECT;
}

---