---

HcalUnpacker Class Reference

#include <EventFilter/HcalRawToDigi/interface/HcalUnpacker.h>

List of all members.

Public Member Functions
	HcalUnpacker (int sourceIdOffset)
	For histograms, no begin and end.
	HcalUnpacker (int sourceIdOffset, int beg, int end)
	for normal data
void	unpack (const FEDRawData &raw, const HcalElectronicsMap &emap, std::vector< HFDataFrame > &precision, std::vector< HcalTriggerPrimitiveDigi > &tp)
void	unpack (const FEDRawData &raw, const HcalElectronicsMap &emap, std::vector< HODataFrame > &precision, std::vector< HcalTriggerPrimitiveDigi > &tp)
void	unpack (const FEDRawData &raw, const HcalElectronicsMap &emap, std::vector< HBHEDataFrame > &precision, std::vector< HcalTriggerPrimitiveDigi > &tp)
void	unpack (const FEDRawData &raw, const HcalElectronicsMap &emap, Collections &conts, HcalUnpackerReport &report, bool silent=false)
void	unpack (const FEDRawData &raw, const HcalElectronicsMap &emap, std::vector< HcalHistogramDigi > &histoDigis)
Private Attributes
int	endSample_
	last sample from fed raw data to copy (if present)
int	sourceIdOffset_
	number to subtract from the source id to get the dcc id
int	startSample_
	first sample from fed raw data to copy
std::set< HcalElectronicsId >	unknownIds_
std::set< HcalElectronicsId >	unknownIdsTrig_
	Recorded to limit number of times a log message is generated.
Classes
struct	Collections

---

Detailed Description

Definition at line 18 of file HcalUnpacker.h.

---

Constructor & Destructor Documentation

HcalUnpacker::HcalUnpacker	(	int	sourceIdOffset,
		int	beg,
		int	end
	)			[inline]

for normal data

Definition at line 33 of file HcalUnpacker.h.

: sourceIdOffset_(sourceIdOffset), startSample_(beg), endSample_(end) { }

HcalUnpacker::HcalUnpacker

(

int

sourceIdOffset

)

[inline]

For histograms, no begin and end.

Definition at line 35 of file HcalUnpacker.h.

: sourceIdOffset_(sourceIdOffset), startSample_(-1), endSample_(-1) { }

---

Member Function Documentation

void HcalUnpacker::unpack	(	const FEDRawData &	raw,
		const HcalElectronicsMap &	emap,
		std::vector< HFDataFrame > &	precision,
		std::vector< HcalTriggerPrimitiveDigi > &	tp
	)

Definition at line 318 of file HcalUnpacker.cc.

References c, HcalUnpacker::Collections::hfCont, r, HcalUnpacker::Collections::tpCont, and unpack().

                                                                                                                                                           {
  Collections c;
  c.hfCont=&container;
  c.tpCont=&tp;
  HcalUnpackerReport r;
  unpack(raw,emap,c,r);
}

void HcalUnpacker::unpack	(	const FEDRawData &	raw,
		const HcalElectronicsMap &	emap,
		std::vector< HODataFrame > &	precision,
		std::vector< HcalTriggerPrimitiveDigi > &	tp
	)

Definition at line 310 of file HcalUnpacker.cc.

References c, HcalUnpacker::Collections::hoCont, r, HcalUnpacker::Collections::tpCont, and unpack().

                                                                                                                                                           {
  Collections c;
  c.hoCont=&container;
  c.tpCont=&tp;
  HcalUnpackerReport r;
  unpack(raw,emap,c,r);
}

void HcalUnpacker::unpack	(	const FEDRawData &	raw,
		const HcalElectronicsMap &	emap,
		std::vector< HBHEDataFrame > &	precision,
		std::vector< HcalTriggerPrimitiveDigi > &	tp
	)

Definition at line 302 of file HcalUnpacker.cc.

References c, HcalUnpacker::Collections::hbheCont, r, HcalUnpacker::Collections::tpCont, and unpack().

                                                                                                                                                             {
  Collections c;
  c.hbheCont=&container;
  c.tpCont=&tp;
  HcalUnpackerReport r;
  unpack(raw,emap,c,r);
}

void HcalUnpacker::unpack	(	const FEDRawData &	raw,
		const HcalElectronicsMap &	emap,
		Collections &	conts,
		HcalUnpackerReport &	report,
		bool	silent = false
	)

work through the samples

work through the samples

Definition at line 53 of file HcalUnpacker.cc.

References HcalUnpacker::Collections::calibCont, DetId::Calo, HcalHTRData::check(), HOUnrolledTP::checked, HcalUnpackerReport::countSpigotFormatError(), HcalUnpackerReport::countUnmappedDigi(), HcalUnpackerReport::countUnmappedTPDigi(), FEDRawData::data(), HcalHTRData::dataPointers(), DetId::det(), endSample_, HcalQIESample::fiber(), HcalQIESample::fiberAndChan(), HcalQIESample::fiberChan(), HcalHTRData::getFirmwareFlavor(), HcalHTRData::getFormatVersion(), HcalHTRData::getNPS(), HcalDCCHeader::getSourceId(), HcalDCCHeader::getSpigotCRCError(), HcalDCCHeader::getSpigotData(), HcalDCCHeader::getSpigotPresent(), HcalHTRData::getSubmodule(), HcalUnpacker::Collections::hbheCont, DetId::Hcal, HcalBarrel, HcalCalibration, HcalEmpty, HcalEndcap, HcalForward, HcalOther, HcalOuter, HcalUnpacker::Collections::hfCont, HcalUnpacker::Collections::hoCont, i, id, HOUnrolledTP::ieta, HcalDetId::ieta(), HOUnrolledTP::iphi, HcalDetId::iphi(), HcalHTRData::isHistogramEvent(), isTPGSOI(), LogDebug, HcalElectronicsMap::lookup(), HcalElectronicsMap::lookupTrigger(), DetId::null(), HcalTriggerPrimitiveSample::raw(), HcalQIESample::raw(), HOUnrolledTP::samples, HOUnrolledTP::setbit(), HcalElectronicsId::setHTR(), FEDRawData::size(), HcalTriggerPrimitiveSample::slb(), HcalTriggerPrimitiveSample::slbAndChan(), HcalTriggerPrimitiveSample::slbChan(), HOUnrolledTP::soi, sourceIdOffset_, HcalDCCHeader::SPIGOT_COUNT, startSample_, HcalOtherDetId::subdet(), HcalZDCDetId::SubdetectorId, DetId::subdetId(), HcalUnpacker::Collections::tpCont, HcalUnpacker::Collections::tphoCont, HcalTrigTowerDetId::Undefined, unknownIds_, unknownIdsTrig_, HOUnrolledTP::valid, TrackValidation_HighPurity_cff::valid, and HcalUnpacker::Collections::zdcCont.

                                                                                       {

  if (raw.size()<16) {
    if (!silent) edm::LogWarning("Invalid Data") << "Empty/invalid DCC data, size = " << raw.size();
    return;
  }

  // get the DCC header
  const HcalDCCHeader* dccHeader=(const HcalDCCHeader*)(raw.data());
  int dccid=dccHeader->getSourceId()-sourceIdOffset_;

  // check the summary status
  
  // walk through the HTR data...
  HcalHTRData htr;
  const unsigned short* daq_first, *daq_last, *tp_first, *tp_last;
  const HcalQIESample* qie_begin, *qie_end, *qie_work;
  const HcalTriggerPrimitiveSample *tp_begin, *tp_end, *tp_work; 
  for (int spigot=0; spigot<HcalDCCHeader::SPIGOT_COUNT; spigot++) {
    if (!dccHeader->getSpigotPresent(spigot)) continue;

    int retval=dccHeader->getSpigotData(spigot,htr,raw.size());
    if (retval!=0) {
      if (retval==-1) {
        if (!silent) edm::LogWarning("Invalid Data") << "Invalid HTR data (data beyond payload size) observed on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId();
        report.countSpigotFormatError();
      }
      continue;
    }
    // check
    if (dccHeader->getSpigotCRCError(spigot)) {
      if (!silent) 
        edm::LogWarning("Invalid Data") << "CRC Error on HTR data observed on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId();
      report.countSpigotFormatError();
      continue;
    }  
    if (!htr.check()) {
      if (!silent) 
        edm::LogWarning("Invalid Data") << "Invalid HTR data observed on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId();
      report.countSpigotFormatError();
      continue;
    }  
    if (htr.isHistogramEvent()) {
      if (!silent) edm::LogWarning("Invalid Data") << "Histogram data passed to non-histogram unpacker on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId();
      continue;
    }
    if (htr.getFirmwareFlavor()==0x81) {
      LogDebug("HcalTechTrigProcessor") << "Skipping data on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId() << " which is from the TechTrigProcessor (use separate unpacker!)";
      continue;
    }
    if (htr.getFirmwareFlavor()>=0x80) {
      if (!silent) edm::LogWarning("HcalUnpacker") << "Skipping data on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId() << " which is of unknown flavor " << htr.getFirmwareFlavor();
      continue;
    }

    // calculate "real" number of presamples
    int nps=htr.getNPS()-startSample_;
    
    // get pointers
    htr.dataPointers(&daq_first,&daq_last,&tp_first,&tp_last);
    unsigned int smid=htr.getSubmodule();
    int htr_tb=smid&0x1;
    int htr_slot=(smid>>1)&0x1F;
    int htr_cr=(smid>>6)&0x1F;
    
    tp_begin=(HcalTriggerPrimitiveSample*)tp_first;
    tp_end=(HcalTriggerPrimitiveSample*)(tp_last+1); // one beyond last..
    
    int currFiberChan=0x3F; // invalid fiber+channel...
    int ncurr=0;
    bool valid=false;

    bool tpgSOIbitInUse=htr.getFormatVersion()>=3; // version 3 and later
    bool isHOtpg=htr.getFormatVersion()>=3 && htr.getFirmwareFlavor()==0; // HO is flavor zero
    int npre=0;
    /*
      Unpack the trigger primitives
    */
    if (isHOtpg) {
      HOUnrolledTP unrolled[24];
      for (tp_work=tp_begin; tp_work!=tp_end; tp_work++) {
        if (tp_work->raw()==0xFFFF) continue; // filler word
        int sector=tp_work->slbChan();
        if (sector>2) continue;

        for (int ibit=0; ibit<8; ibit++) {
          int linear=sector*8+ibit; 
          if (!unrolled[linear].checked) {
            unrolled[linear].checked=true;
            int fiber=(linear/3)+1;
            int fc=(linear%3);
            // electronics id (use precision match for HO TP)
            HcalElectronicsId eid(fc,fiber,spigot,dccid);       
            eid.setHTR(htr_cr,htr_slot,htr_tb);
            DetId did=emap.lookup(eid);
            if (!did.null()) {
              if (did.det()==DetId::Hcal && ((HcalSubdetector)did.subdetId())==HcalOuter ) {
                HcalDetId hid(did);
                unrolled[linear].valid=true;
                unrolled[linear].ieta=hid.ieta();
                unrolled[linear].iphi=hid.iphi();
              }
            } else {
              report.countUnmappedTPDigi();
            }
          }
          if (unrolled[linear].valid) {
            if (isTPGSOI(*tp_work)) unrolled[linear].soi=unrolled[linear].samples;
            if (tp_work->raw()&(1<<ibit)) unrolled[linear].setbit(unrolled[linear].samples);
            unrolled[linear].samples++;
          }
        }
      }
      for (int i=0; i<24; i++) {
        if (unrolled[i].valid) 
          colls.tphoCont->push_back(HOTriggerPrimitiveDigi(
                                                           unrolled[i].ieta,
                                                           unrolled[i].iphi,
                                                           unrolled[i].samples,
                                                           unrolled[i].soi,
                                                           unrolled[i].databits));
      }
    } else {
      for (tp_work=tp_begin; tp_work!=tp_end; tp_work++) {
        if (tp_work->raw()==0xFFFF) continue; // filler word
        if (tp_work->slbAndChan()!=currFiberChan) { // start new set
          npre=0;
          currFiberChan=tp_work->slbAndChan();
          // lookup the right channel
          HcalElectronicsId eid(tp_work->slbChan(),tp_work->slb(),spigot,dccid,htr_cr,htr_slot,htr_tb);
          DetId did=emap.lookupTrigger(eid);
          if (did.null()) {
            report.countUnmappedTPDigi();
            if (unknownIdsTrig_.find(eid)==unknownIdsTrig_.end()) {
              if (!silent) edm::LogWarning("HCAL") << "HcalUnpacker: No trigger primitive match found for electronics id :" << eid;
              unknownIdsTrig_.insert(eid);
            }
            valid=false;
            continue;
          } else if (did==HcalTrigTowerDetId::Undefined || 
                     (did.det()==DetId::Hcal && did.subdetId()==0)) {
            // known to be unmapped
            valid=false;
            continue;
          }
          HcalTrigTowerDetId id(did);
          colls.tpCont->push_back(HcalTriggerPrimitiveDigi(id));
          // set the various bits
          if (!tpgSOIbitInUse) colls.tpCont->back().setPresamples(nps);
          // no hits recorded for current
          ncurr=0;
          valid=true;
        }      
        // add the word (if within settings or recent firmware [recent firmware ignores startSample/endSample])
        if (valid && ((tpgSOIbitInUse && ncurr<10) || (ncurr>=startSample_ && ncurr<=endSample_))) {
          colls.tpCont->back().setSample(colls.tpCont->back().size(),*tp_work);
          colls.tpCont->back().setSize(colls.tpCont->back().size()+1);
        }
        // set presamples,if SOI
        if (valid && tpgSOIbitInUse && isTPGSOI(*tp_work)) {
          colls.tpCont->back().setPresamples(ncurr);
        }
        ncurr++;
        npre++;
      }
    }
  

    qie_begin=(HcalQIESample*)daq_first;
    qie_end=(HcalQIESample*)(daq_last+1); // one beyond last..

    currFiberChan=0x3F; // invalid fiber+channel...
    ncurr=0;
    valid=false;
    
    for (qie_work=qie_begin; qie_work!=qie_end; ) {
      if (qie_work->raw()==0xFFFF) {
        qie_work++;
        continue; // filler word
      }
      // always at the beginning ...
      currFiberChan=qie_work->fiberAndChan();

      // lookup the right channel
      HcalElectronicsId eid(qie_work->fiberChan(),qie_work->fiber(),spigot,dccid);
      eid.setHTR(htr_cr,htr_slot,htr_tb);
      DetId did=emap.lookup(eid);

      if (!did.null()) {
        if (did.det()==DetId::Calo && did.subdetId()==HcalZDCDetId::SubdetectorId) {
          colls.zdcCont->push_back(ZDCDataFrame(HcalZDCDetId(did)));
          qie_work=HcalUnpacker_impl::unpack<ZDCDataFrame>(qie_work, qie_end, colls.zdcCont->back(), nps, eid, startSample_, endSample_); 
        } else if (did.det()==DetId::Hcal) {
          switch (((HcalSubdetector)did.subdetId())) {
          case (HcalBarrel):
          case (HcalEndcap): {
            colls.hbheCont->push_back(HBHEDataFrame(HcalDetId(did)));
            qie_work=HcalUnpacker_impl::unpack<HBHEDataFrame>(qie_work, qie_end, colls.hbheCont->back(), nps, eid, startSample_, endSample_);
          } break;
          case (HcalOuter): {
            colls.hoCont->push_back(HODataFrame(HcalDetId(did)));
            qie_work=HcalUnpacker_impl::unpack<HODataFrame>(qie_work, qie_end, colls.hoCont->back(), nps, eid, startSample_, endSample_);
          } break;
          case (HcalForward): {
            colls.hfCont->push_back(HFDataFrame(HcalDetId(did)));
            qie_work=HcalUnpacker_impl::unpack<HFDataFrame>(qie_work, qie_end, colls.hfCont->back(), nps, eid, startSample_, endSample_);
          } break;
          case (HcalOther) : {
            HcalOtherDetId odid(did);
            if (odid.subdet()==HcalCalibration) {
              colls.calibCont->push_back(HcalCalibDataFrame(HcalCalibDetId(did)));
              qie_work=HcalUnpacker_impl::unpack<HcalCalibDataFrame>(qie_work, qie_end, colls.calibCont->back(), nps, eid, startSample_, endSample_); 
            }
          } break;
          case (HcalEmpty): 
          default: {
            for (int fiberC=qie_work->fiberAndChan();
                 qie_work!=qie_end && qie_work->fiberAndChan()==fiberC;
                 qie_work++);
          }
          break;
          }
        }
      } else {
        report.countUnmappedDigi();
        if (unknownIds_.find(eid)==unknownIds_.end()) {
          if (!silent) edm::LogWarning("HCAL") << "HcalUnpacker: No match found for electronics id :" << eid;
          unknownIds_.insert(eid);
        }
        for (int fiberC=qie_work->fiberAndChan();
             qie_work!=qie_end && qie_work->fiberAndChan()==fiberC;
             qie_work++);
      }
    }
  }
}

void HcalUnpacker::unpack	(	const FEDRawData &	raw,
		const HcalElectronicsMap &	emap,
		std::vector< HcalHistogramDigi > &	histoDigis
	)

Definition at line 326 of file HcalUnpacker.cc.

References HcalHTRData::check(), FEDRawData::data(), DetId::det(), f, HcalHistogramDigi::getArray(), HcalHTRData::getHistogramFibers(), HcalDCCHeader::getSourceId(), HcalDCCHeader::getSpigotData(), HcalDCCHeader::getSpigotPresent(), HcalHTRData::getSubmodule(), DetId::Hcal, HcalHTRData::isHistogramEvent(), HcalElectronicsMap::lookup(), DetId::null(), HcalElectronicsId::setHTR(), FEDRawData::size(), sourceIdOffset_, HcalDCCHeader::SPIGOT_COUNT, DetId::subdetId(), unknownIds_, and HcalHTRData::unpackHistogram().

Referenced by HcalRawToDigi::produce(), HcalHistogramRawToDigi::produce(), and unpack().

                                                                                                                         {

  // get the DCC header
  const HcalDCCHeader* dccHeader=(const HcalDCCHeader*)(raw.data());
  int dccid=dccHeader->getSourceId()-sourceIdOffset_;
  
  // check the summary status
  
  // walk through the HTR data...
  HcalHTRData htr;
  for (int spigot=0; spigot<HcalDCCHeader::SPIGOT_COUNT; spigot++) {
    if (!dccHeader->getSpigotPresent(spigot)) continue;
    
    int retval=dccHeader->getSpigotData(spigot,htr,raw.size());
    // check
    if (retval || !htr.check()) {
      edm::LogWarning("Invalid Data") << "Invalid HTR data observed on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId();
      continue;
    }
    if (!htr.isHistogramEvent()) {
      edm::LogWarning("Invalid Data") << "Non-histogram data passed to histogram unpacker on spigot " << spigot << " of DCC with source id " << dccHeader->getSourceId();
      continue;
    }

    unsigned int smid=htr.getSubmodule();
    int htr_tb=smid&0x1;
    int htr_slot=(smid>>1)&0x1F;
    int htr_cr=(smid>>6)&0x1F;
    
    // find out the fibers
    int f[2],fc;
    htr.getHistogramFibers(f[0],f[1]);

    for (int nf=0; nf<2; nf++) {
      if (f[nf]<0 || nf==1 && f[0]==f[1]) continue; // skip if invalid or the same
      for (fc=0; fc<=2; fc++) {
        HcalElectronicsId eid(fc,f[nf],spigot,dccid);     
        eid.setHTR(htr_cr,htr_slot,htr_tb);
        DetId did=emap.lookup(eid);

        if (did.null() || did.det()!=DetId::Hcal || did.subdetId()==0) {
          if (unknownIds_.find(eid)==unknownIds_.end()) {
            edm::LogWarning("HCAL") << "HcalHistogramUnpacker: No match found for electronics id :" << eid;
            unknownIds_.insert(eid);
          }       
          continue;
        }
        histoDigis.push_back(HcalHistogramDigi(HcalDetId(did))); // add it!
        HcalHistogramDigi& digi=histoDigis.back();
        
        // unpack the four capids
        for (int capid=0; capid<4; capid++) 
          htr.unpackHistogram(f[nf],fc,capid,digi.getArray(capid));
        
      }
    }
  }
}      

---

Member Data Documentation

int HcalUnpacker::endSample_ [private]

last sample from fed raw data to copy (if present)

Definition at line 47 of file HcalUnpacker.h.

Referenced by unpack().

int HcalUnpacker::sourceIdOffset_ [private]

number to subtract from the source id to get the dcc id

Definition at line 45 of file HcalUnpacker.h.

Referenced by unpack().

int HcalUnpacker::startSample_ [private]

first sample from fed raw data to copy

Definition at line 46 of file HcalUnpacker.h.

Referenced by unpack().

std::set<HcalElectronicsId> HcalUnpacker::unknownIds_ [private]

Definition at line 48 of file HcalUnpacker.h.

Referenced by unpack().

std::set<HcalElectronicsId> HcalUnpacker::unknownIdsTrig_ [private]

Recorded to limit number of times a log message is generated.

Definition at line 48 of file HcalUnpacker.h.

Referenced by unpack().

---

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

• EventFilter/HcalRawToDigi/interface/HcalUnpacker.h
• EventFilter/HcalRawToDigi/src/HcalUnpacker.cc

---