/data3/calcul/jacquem/working_dir/Micromegas/micromegasFrameWork/src/parser/DiracLabview.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) $ */

#include "parser/DiracLabview.hh"
#include "parser/DifReader.hh"

#include "event/Event.hh"

#include "geometry/Detector.hh"
#include "geometry/Chamber.hh"
#include "geometry/Dif.hh"
#include "geometry/Chip.hh"
#include "geometry/DiracChip.hh"

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

#include <iostream>
#include <cstdio>



#define bigTolittle Toolbox::bigTolittle


using namespace std;


DiracLabview::DiracLabview(Run& aRun, FILE* aFile, ui32 _firstEventId) :  DifReader(aRun, aFile, _firstEventId)
{
  //FILE_LOG(logDEBUG) << "--DiracLabview constructor()"<< endl;
}

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


int  DiracLabview::getNextEvent(Event& event)
{
        ui16  dataFormat = run.getDataFormat();
  if (dataFormat == 0)
  {
      dataFormat = 0;
                        run.setDataFormat(dataFormat);
      return  NEW_CONFIG;
  }
        char strReadLine[256];
        static short iNbEventInMemory = 0;
        static unsigned char fChannelValues[NB_EVENT_DUMP][NB_CHANNELS], fStimChannel;
        static short fDacLoA,fDacMidA,fDacHiA,fDacLoB,fDacMidB,fDacHiB;
        static int fNbDump, fBcid[NB_EVENT_DUMP];
        static double fCharge;
//      Int_t fEventId=0;
        unsigned short fNbChannelHit;
        bool bFlagRun = true;
        bool bFlagReadEventHeader = true;
        Detector& aDet = run.getDetector();

        if (iNbEventInMemory==0){
                while((feof(inputFile)==0) && bFlagRun){
                        fgets(strReadLine,256,inputFile) ;
                        if (!strncmp(strReadLine,"voie  0 = ",10)){
                                bFlagReadEventHeader=0;
                        }
                        if (!strncmp(strReadLine,"voie 63 = ",10)){
                                bFlagRun=0;
                        }
                        // Read Chip Identifier
                        if (bFlagReadEventHeader){
                                if (!strncmp(strReadLine,"Fin de run",10)){
                                        return LAST_FILE_EVENT;
                                }
                                if (!strncmp(strReadLine,"Channel = ",10)){
                                        sscanf(&strReadLine[10],"%hhd",&fStimChannel);
                                }
                                if (!strncmp(strReadLine,"DAC_A (L-M-H) = ",16)){
                                        sscanf(&strReadLine[16],"%hd %hd %hd",&fDacLoA,&fDacMidA,&fDacHiA);
                                }
                                if (!strncmp(strReadLine,"DAC_B (L-M-H) = ",16)){
                                        sscanf(&strReadLine[16],"%hd %hd %hd",&fDacLoB,&fDacMidB,&fDacHiB);
                                }
                                if (!strncmp(strReadLine,"acquisition = ",14)){
                                        sscanf(&strReadLine[14],"%d",&fNbDump);
                                }
                                if (!strncmp(strReadLine,"Charge = ",9)){
                                        sscanf(&strReadLine[9],"%lE",&fCharge);
                                }
                                if (!strncmp(strReadLine,"BCID = ",7)){
                                        sscanf(&strReadLine[7],"%d %d %d %d %d %d %d %d",&fBcid[0],&fBcid[1],&fBcid[2],&fBcid[3],&fBcid[4],&fBcid[5],&fBcid[6],&fBcid[7]);
                                }
                        }
                        if (!bFlagReadEventHeader){
                                //if (!strncmp(strReadLine,"voie ",5)){
                                        short iChannel;
                                        sscanf(&strReadLine[5],"%hd",&iChannel);
                                        for (unsigned short iIndex=0;iIndex < NB_EVENT_DUMP;iIndex++){
                                                fChannelValues[iIndex][iChannel]=0;
                                                if (strReadLine[10+iIndex]=='1') fChannelValues[iIndex][iChannel]+=2;
                                                if (strReadLine[19+iIndex]=='1') fChannelValues[iIndex][iChannel]+=1;
                                        }
                                //}
                        }
                }
                iNbEventInMemory=8;
        }

        if (iNbEventInMemory!=0){
//              if (lastEventId<17){
                        FILE_LOG(logDEBUG1) << endl;
                        FILE_LOG(logDEBUG1) << "_________________________________________________________" << endl;
                        FILE_LOG(logDEBUG1) << "Event ID : " << lastEventId << " - N° acq. : " << fNbDump << " - Stim' Channel : " << static_cast<int>(fStimChannel) << endl;
                        FILE_LOG(logDEBUG1) << "Chip BCID : " << fBcid[8-iNbEventInMemory] << " - Dif BCID : " << "??" << endl;
                        FILE_LOG(logDEBUG1) << "Hits :";
                        for(unsigned char iIndex=0;iIndex<NB_CHANNELS;iIndex++){
                                if (fChannelValues[8-iNbEventInMemory][iIndex]!=0){
                                        FILE_LOG(logDEBUG1) << " Ch"<<static_cast<int>(iIndex)<<"=>'"<<static_cast<int>(fChannelValues[8-iNbEventInMemory][iIndex])<<"'" << endl;
                                        fNbChannelHit++;
                                }
                        }

      run.setInjectedCharge(fCharge);
                        FILE_LOG(logDEBUG1) << "fCharge[" <<  fCharge << "]"  << endl;
      try {
          Dif& dif = aDet.getDifById(1);
          DiracChip& diracChip = dynamic_cast<DiracChip &> (dif.getChipById(1));

          FILE_LOG(logDEBUG1) << " Dirac slow control for chip configuration set: " <<    endl;
          diracChip.setBLoThreshold(fDacLoB);
          diracChip.setBMidThreshold(fDacMidB);
          diracChip.setBHiThreshold(fDacHiB);

          diracChip.setALoThreshold(fDacLoA);
          diracChip.setAMidThreshold(fDacMidA);
          diracChip.setAHiThreshold(fDacHiA);
          diracChip.setConfigured(true);

                            FILE_LOG(logDEBUG1) << endl;
          FILE_LOG(logDEBUG1) << "Dac A : " << fDacLoA << " / " << fDacMidA << " / " << fDacHiA << endl;
          FILE_LOG(logDEBUG1) << "Dac B : " << fDacLoB << " / " << fDacMidB << " / " << fDacHiB << endl;
          FILE_LOG(logDEBUG1) << "Charge: " << fCharge << endl;
      }
      catch(...) {
             FILE_LOG(logDEBUG1) << "Ouh la la tu va avoir des gros problemes toi.: " << endl;
      }

                //event.Build(fEventId,1,123,123,123);
                //event.setQin(fCharge);
                //event.SetNbDump(fNbDump);
                //event.SetChipBcid(fBcid[8-iNbEventInMemory]);
                //event.SetDifBcid(1234);                                       // !!!!!

                fNbChannelHit=0;
                //event.SetNchannel(fNbChannelHit);
                const ChamberMap_t &chambers = aDet.getChambers();
    for (ChamberMap_t::const_iterator iter=chambers.begin(); iter!=chambers.end(); iter++)
    {
      Chamber& chamber = *(*iter).second;
      FILE_LOG(logDEBUG1)<< "Chamber->id="<< chamber.getId() << " , descrition:" << chamber.getDescription() << endl;
                        //if (chamber.getId()==0){
                                for(unsigned char iIndex=0;iIndex<NB_CHANNELS;iIndex++){
                                        if (fChannelValues[8-iNbEventInMemory][iIndex]!=0){
                                        try{

                                                const Channel&  channel = chamber.getChannelById(iIndex,1,1);
            i16 analogValue = 0;
            FILE_LOG(logDEBUG1)<< "insertHit" << endl;
                                         event.insertHit(chamber,new ChannelHit(channel,fChannelValues[8-iNbEventInMemory][iIndex],analogValue));
                                                }
                                                catch ( MicroException &e){
                                                        FILE_LOG(logERROR) << e.getMessage() << endl;
                                                }
                                                //event.AddChannel(iIndex+1,fChannelValues[8-iNbEventInMemory][iIndex]);
                                                fNbChannelHit++;
                                        }
                                }
                        //}
                }
                //event.SetNchannel(fNbChannelHit);
                event.setId(++lastEventId);
                event.setTimeStamp(fBcid[8-iNbEventInMemory]);
                iNbEventInMemory--;
        }
  FILE_LOG(logDEBUG1)<< "event correct" << endl;
        return EVENT_CORRECT;
}

---