DiracReader Class Reference

#include <DiracReader.hh>

Inheritance diagram for DiracReader:

[legend]Collaboration diagram for DiracReader:

[legend]List of all members.

Public Member Functions
	DiracReader (Run &aRun, FILE *aFile, ui32 firstEventId)
	~DiracReader ()
int	getNextEvent (Event &eventToFill)
Static Public Member Functions
static const char *	type (const int subtype=0)
Protected Attributes
unsigned short	CRC
int	nbDifs
std::vector< DiracDifConfig >	difConfigs
i64	bcId_AbsOrg
int	acqDifId
DiracTriggerInfo	trigInfos [8]
Private Member Functions
void	reset ()
int	getConfigData ()
int	getAcqData ()
int	storeHits (Event &event)
ChannelHit *	newHit (Event &event, const DiracDifConfig &difConfig, const DiracTriggerInfo &trigInfo, const int chipNum, const int chNum)
void	storeDifInfo (const int difId)

---

Detailed Description

CLASS DATAREADER //

Definition at line 108 of file DiracReader.hh.

---

Constructor & Destructor Documentation

DiracReader::DiracReader	(	Run &	aRun,
		FILE *	aFile,
		ui32	firstEventId
	)

Definition at line 228 of file DiracReader.cpp.

                                                                   : DifReader(aRun, aFile, _firstEventId) {
    reset();
}

DiracReader::~DiracReader

(

)

Definition at line 234 of file DiracReader.cpp.

                          {
} // ~DiracReader

---

Member Function Documentation

int DiracReader::getNextEvent

(

Event &

eventToFill

)

[virtual]

Implements AcquisitionParser.

Definition at line 524 of file DiracReader.cpp.

                                          {
    // Start data reconstruction
    Level(0);

                i16 dataFormat = run.getDataFormat();
    if (dataFormat == -1) {
        // first read of the file
        try
        {
            bcIdDifAbsOrgMap.clear();
        }

        catch(exception & e) {
            FILE_LOG(logERROR) << "bcIdDifAbsOrgMap.clear():" << e.what() << endl;
        }
        // read file header
        dataFormat = getBinData(1);
                                run.setDataFormat(dataFormat);
        if (fileError == EOF) {
            FILE_LOG(logDEBUG) << "end of file found" << endl;
            reset();
            return (LAST_FILE_EVENT);
        }

        if (dataFormat != 1) {
            // unknown format
            FILE_LOG(logERROR) << "DIRAC data format " << dataFormat << " is not implemented" << endl;
            reset();
            return (LAST_FILE_EVENT);
        }

        nbDifs = getBinData(1);
        display(+1);
        FILE_LOG(logDEBUG) << "DIRAC data format : " << dataFormat << ", nb Dif = " << nbDifs << endl;

        for (int difNum = 0; difNum < nbDifs; difNum++) {
            // read DIF configuration
            DiracDifConfig difConfig;
            unsigned char buf[21];
            getBinData(20, buf);
            display();

            difConfig.read(buf);
            difConfigs.push_back(difConfig);

            // store data into Run info
            storeDifInfo(difConfig.id);

            // store info for the whole run
            run.setInjectedCharge(difConfigs[0].charge);
            run.setWidthCtest(difConfigs[0].ctestWidth);
            run.setPeriodCtest(difConfigs[0].ctestPeriod);
            run.setDelayCtest(difConfigs[0].ctestDelay);
            run.setDelayTrigger(difConfigs[0].triggerDelay);
        }
        display(-1);

        // get unix timestamp for the whole file (seconds since jan 1st, 1970)
        globalTime = getBinData(4);
        display();
        FILE_LOG(logDEBUG1) << "global time : " << globalTime << endl;
        return NEW_CONFIG;
    }

    for (unsigned int nbHits = 0; !nbHits;) {
        unsigned char header = getBinData(1);
        if (fileError == EOF) {
            FILE_LOG(logDEBUG) << "  END OF DATA FILE" << endl;
            reset();
            return (LAST_FILE_EVENT);
        }

        if (header == 0xB1) {
            // store all configuration data
            display(+1);
            getConfigData();
            if (getBinData(1) != 0xA1) {
                display(-1);
                FILE_LOG(logERROR) << "    frame trailer A1 not found" << endl;
                break;
            }
            display(-1);
        } else if (header == 0xB0) {
            // store all acquisition data
            display(+1);
            getAcqData();
            if (getBinData(1) != 0xA0) {
                display(-1);
                FILE_LOG(logERROR) << "    frame trailer A0 not found" << endl;
                break;
            }
            display(-1);
            CRC = getBinData(2);
            display();

            // then store all data
            nbHits = storeHits(event);
        } else {
            display();
            FILE_LOG(logERROR) << "    expected header not found" << endl;
        }

        if (fileError == EOF)
            break;
    }

    return EVENT_CORRECT;
} // DiracReader.getNextEvent

static const char* DiracReader::type

(

const int

subtype = 0

)

[inline, static]

Definition at line 115 of file DiracReader.hh.

Referenced by main().

{ return LABVIEW_DIRAC_READER; };

void DiracReader::reset

(

)

[private]

Definition at line 239 of file DiracReader.cpp.

Referenced by DiracReader(), and getNextEvent().

                        {
    nbDifs = 0;
    CRC = 0;
    bcId_AbsOrg = -1; // the bcId origin
} // DiracReader.reset

int DiracReader::getConfigData

(

)

[private]

Definition at line 248 of file DiracReader.cpp.

Referenced by getNextEvent().

                               {
    int difId = getBinData(1); // read 8 bits

    for (int difNum = 0; difNum < difConfigs.size(); difNum++) {
        DiracDifConfig& difConfig = difConfigs[difNum];
        if (difConfig.id == difId) {
            difConfig.nbChips = getBinData(1); // read 8 bits

            display();
            FILE_LOG(logDEBUG1) << " DIF #" << difConfig.id << " : " << difConfig.nbChips << " chips" << endl;

            for (int chipNum = 0; chipNum < difConfig.nbChips; chipNum++) {
                DiracChipConfig chipConfig;

                unsigned char buf[24];
                getBinData(23, buf);
                display();
                chipConfig.read(buf);
                // InitGetBits(23);
                // chipconfg.read();   <--- GetBits(n);
                // FinishGetBits();

                difConfig.chipConfigs.push_back(chipConfig);

                const Detector &detector = run.getDetector();
                try
                {
                    Dif& dif = detector.getDifById(difId);
                    BoardMap_t boards = dif.getBoards();
                    for (BoardMap_t::iterator boardIt = boards.begin(); boardIt != boards.end(); boardIt++) {
                        Board& b = *boardIt->second;
                        try
                        {
                            DiracChip& chip = dynamic_cast<DiracChip&> (b.getChipById(chipNum));
                            chip.setBLoThreshold(chipConfig.bLoThreshold);
                            chip.setBMidThreshold(chipConfig.bMidThreshold);
                            chip.setBHiThreshold(chipConfig.bHiThreshold);
                            chip.setALoThreshold(chipConfig.aLoThreshold);
                            chip.setAMidThreshold(chipConfig.aMidThreshold);
                            chip.setAHiThreshold(chipConfig.aHiThreshold);
                            chip.setGain(chipConfig.gain);
                            for (int chNum = 0; chNum < 64; chNum++) {
                                Channel& ch = chip.getChannelById(chNum);
                                ch.setEnable(chipConfig.chConfig[chNum].enable);
                                ch.setStimulate(chipConfig.chConfig[chNum].test);
                            }
                            break;
                        }

                        catch(...) {
                        }
                    }
                }

                catch(...) {
                }
            }

            if (difConfigs.size() > nbDifs) {
                FILE_LOG(logERROR) << "DIF config error" << endl;
                return (0);
            }
            return (1);
        } // if
    } // for

    FILE_LOG(logERROR) << "DIF #" << difId << " not found" << endl;
    return (0);
} // DiracReader.getConfigData

int DiracReader::getAcqData

(

)

[private]

Definition at line 320 of file DiracReader.cpp.

Referenced by getNextEvent().

                            {
    acqDifId = getBinData(1); // read 8 bits
    unsigned int nbTotalTrigger = getBinData(4);
    unsigned int nbTrigger = getBinData(4);
    display();

    FILE_LOG(logDEBUG) << "  data for DIF #" << acqDifId << ", " << nbTotalTrigger << " total, " << nbTrigger << " triggers" << endl;

    // look for current DIF configuration
    for (int difNum = 0; difNum < difConfigs.size(); difNum++) {
        const DiracDifConfig& difConfig = difConfigs[difNum];
        if (difConfig.id == acqDifId) {
            // read 8 trigger informations (10 bytes each)
            for (int trigInfoNum = 0; trigInfoNum < 8; trigInfoNum++) {
                unsigned char buf[11];
                getBinData(10, buf);
                display();
                trigInfos[trigInfoNum].read(buf);
            }

            // set the origin for absolute bcId
            //      if (bcId_AbsOrg == -1) {
            //        // set the absolute bcId origin
            //        bcId_AbsOrg = trigInfos[0].bcIdAbs;
            //        FILE_LOG(logDEBUG)   << "    set bc origin         " << bcId_AbsOrg << endl;
            //      }
            bcId_AbsOrg = getbcIdAbsOrg(acqDifId, trigInfos[0].bcIdAbs);
            FILE_LOG(logDEBUG) << " getbcIdAbsOrg(" << acqDifId << " , trigInfos[0].bcIdAbs)  " << getbcIdAbsOrg(acqDifId, trigInfos[0].bcIdAbs) << endl;

            // the following is sent for every chip for the given DIF
            for (int chipNum = 0; chipNum < difConfig.nbChips; chipNum++) {
                unsigned int header = getBinData(1);
                if (header != 0xB4) {
                    FILE_LOG(logERROR) << "parsing file error : header 0xB4 not found" << endl;
                    return (0);
                }
                display(+1);

                // get 8 bcIds (96 bits = 12 bytes) <-> 8 triggers * 12-bit bcId
                for (int bitNum = 11; bitNum >= 0; bitNum--) {
                    const unsigned int data = getBinData(1);
                    unsigned int mask = 0x80;
                    for (int trigNum = 0; trigNum < 8; trigNum++, mask >>= 1) {
                        DiracTriggerInfo& trigInfo = trigInfos[trigNum];
                        unsigned int bit = ((data & mask) != 0) ? 1 : 0;
                        trigInfo.bcId |= bit << bitNum;
                    }
                }
                display();

                // get 8 channel data (1024 bits = 128 bytes) <-> 8 triggers * 64 channels * 2-bit data
                for (int chNum = 63; chNum >= 0; chNum--) {
                    for (int bitNum = 0; bitNum < 2; bitNum++) {
                        const unsigned int data = getBinData(1);
                        unsigned int mask = 0x80;
                        for (int trigNum = 0; trigNum < 8; trigNum++, mask >>= 1) {
                            DiracTriggerInfo& trigInfo = trigInfos[trigNum];
                            unsigned int bit = ((data & mask) != 0) ? 1 : 0;
                            trigInfo.data[chNum] |= bit << bitNum;
                        }
                    }
                    if (chNum == 32)
                        display();
                }
                display();

                header = getBinData(1);
                display(-1);
                if (header != 0xA3) {
                    FILE_LOG(logERROR) << "parsing file error : trailer 0xA3 not found" << endl;
                    return (0);
                }
            } // for chipNum

            return (1);
        } // if
    } // for difNum

    return (0);
} // DiracReader.getAcqData

int DiracReader::storeHits

(

Event &

event

)

[private]

Definition at line 403 of file DiracReader.cpp.

Referenced by getNextEvent().

                                       {
    int nbHits = 0;

    for (int difNum = 0; difNum < difConfigs.size(); difNum++) {
        const DiracDifConfig& difConfig = difConfigs[difNum];
        if (difConfig.id == acqDifId) {
            for (int chipNum = 0; chipNum < difConfig.nbChips; chipNum++) {

                for (int trigNum = 0; trigNum < 8; trigNum++) {
                    const DiracTriggerInfo& trigInfo = trigInfos[trigNum];

                    for (int chNum = 0; chNum < 64; chNum++) {
                        if (trigInfo.data[chNum]) {
                            if (newHit(event, difConfig, trigInfo, chipNum, chNum))
                                nbHits++;
                        }
                    }
                }
            } // for chipNum

            if (nbHits) {
                // CREATE EVENT for current Run
                event.setId(++lastEventId);
                FILE_LOG(logDEBUG1) << "  STORE EVENT " << lastEventId << ", " << nbHits << " hits" << endl;
            } else
                FILE_LOG(logDEBUG) << "  EMPTY EVENT" << endl;
            return (nbHits);
        } // if
    } // for difNum

    return (0);
} // DiracReader.storeHits

ChannelHit * DiracReader::newHit	(	Event &	event,
		const DiracDifConfig &	difConfig,
		const DiracTriggerInfo &	trigInfo,
		const int	chipNum,
		const int	chNum
	)			[private]

Definition at line 438 of file DiracReader.cpp.

Referenced by storeHits().

                                                                                                                                                   {
    const Detector &detector = run.getDetector();
    const ChamberMap_t &chambers = detector.getChambers();
    for (ChamberMap_t::const_iterator it = chambers.begin(); it != chambers.end(); ++it) {
        try
        {
            DiracChamber1& chamber = dynamic_cast<DiracChamber1&> (*(it->second));
            const BoardMap_t& boards = chamber.getBoardsByDifId(difConfig.id);
            // store new event data
            try
            {
                const Channel& channel = chamber.getChannelById(chNum, chipNum, difConfig.id);
                const int data = trigInfo.data[chNum];
                ChannelHit *hit = new ChannelHit(channel, data, 0);
                // all in bcPeriods - this depends on the acquisition mode
                i64 bcId_Abs; // absolute bc Id for current DIF ALWAYS @5MHz
                i64 bcId_Dif; // current bcId for current DIF
                i64 bcId_Hit; // bcId of current hit
                i64 bcId_Hit_Fine = trigInfos[0].trigger;

                //        if (difConfig.acqMode == 0x0A) {
                bcId_Abs = trigInfos[0].bcIdAbs;
                bcId_Dif = trigInfos[0].bcIdFpga;
                //        }
                //        else {
                //          bcId_Abs = trigInfo.bcIdAbs;
                //          bcId_Dif = trigInfo.bcIdFpga;
                //        }

                bcId_Hit = trigInfo.bcId;

                hit->setHardrocTime(globalTime, bcId_Abs, bcId_Dif, bcId_Hit,bcId_Hit_Fine);
                event.insertHit(chamber, hit);
                i64 timeStamp = (globalTime * 1000) // all in milliseconds
                        + ((trigInfo.bcIdAbs - bcId_AbsOrg) / 5000000) // ALWAYS @5MHz
                        + (((trigInfo.bcId - trigInfo.bcIdFpga) * difConfig.bcPeriod) / 1000);
                event.setTimeStamp(timeStamp);
                FILE_LOG(logDEBUG) << "bcId_Abs[" << bcId_Abs << "] - bcId_AbsOrg[" << bcId_AbsOrg << "] = [" << bcId_Abs - bcId_AbsOrg << "]" << endl;
                FILE_LOG(logDEBUG) << "  new Hit : ch " << chNum << ", chip " << chipNum << ", dif " << difConfig.id
                        << ", Time=" << timeStamp << "ms, bcId:abs=" << (bcId_Abs - bcId_AbsOrg) << ",dif=" << bcId_Dif << ",hit=" << bcId_Hit << " Data=" << data << endl;
                return (hit);
            }

            catch(MicroException & e) {
                FILE_LOG(logERROR) << "ch " << chNum << ", chip " << chipNum << ", dif " << difConfig.id << " : " << e.getMessage() << endl;
                return (0);
            }
        }

        catch(...) {
            continue;
        }

    } // for
    FILE_LOG(logERROR) << "ch " << chNum << ", chip " << chipNum << ", dif " << difConfig.id << " : geometry element not found" << endl;
    return (0);
} // DiracReader.newHit

void DiracReader::storeDifInfo

(

const int

difId

)

[private]

Definition at line 498 of file DiracReader.cpp.

Referenced by getNextEvent().

                                              {
    const Detector &detector = run.getDetector();
    try
    {
        Dif& dif = detector.getDifById(difId);
        for (int difNum = 0; difNum < difConfigs.size(); difNum++) {
            DiracDifConfig& difConfig = difConfigs[difNum];
            if (difConfig.id == difId) {
                dif.setAcqMode(difConfig.acqMode);
                dif.setLowRegister(difConfig.lowState);
                dif.setHighRegister(difConfig.highState);
                dif.setAfterRegister(difConfig.afterState);
                dif.setBeforeRegister(difConfig.beforeState);
                return;
            }
        }
    }

    catch(...) {
    };
}; // DiracReader.storeDifInfo

---

Member Data Documentation

unsigned short DiracReader::CRC [protected]

Definition at line 128 of file DiracReader.hh.

Referenced by getNextEvent(), and reset().

int DiracReader::nbDifs [protected]

Definition at line 131 of file DiracReader.hh.

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

std::vector<DiracDifConfig> DiracReader::difConfigs [protected]

Definition at line 132 of file DiracReader.hh.

Referenced by getAcqData(), getConfigData(), getNextEvent(), storeDifInfo(), and storeHits().

i64 DiracReader::bcId_AbsOrg [protected]

Definition at line 135 of file DiracReader.hh.

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

int DiracReader::acqDifId [protected]

Definition at line 136 of file DiracReader.hh.

Referenced by getAcqData(), and storeHits().

DiracTriggerInfo DiracReader::trigInfos[8] [protected]

Definition at line 137 of file DiracReader.hh.

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

---

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

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

---