DCCMemBlock Class Reference

#include <DCCMemBlock.h>

Inheritance diagram for DCCMemBlock:

Public Member Functions
	DCCMemBlock (DCCDataUnpacker *u, EcalElectronicsMapper *m, DCCEventBlock *e)
void	display (std::ostream &o) override
virtual int	unpack (const uint64_t **data, unsigned int *dwToEnd)
int	unpack (const uint64_t **data, unsigned int *dwToEnd, unsigned int expectedTowerID)
void	updateCollectors () override
	~DCCMemBlock () override
Public Member Functions inherited from DCCDataBlockPrototype
	DCCDataBlockPrototype (DCCDataUnpacker *unpacker, EcalElectronicsMapper *mapper, DCCEventBlock *event, bool unpack=true)
void	enableSyncChecks ()
virtual unsigned int	getLength ()
virtual int	unpack (const uint64_t **data, unsigned int *dwToEnd)
virtual void	updateEventPointers ()
virtual	~DCCDataBlockPrototype ()

Protected Member Functions
void	fillPnDiodeDigisCollection ()
void	unpackMemTowerData ()

Protected Attributes
unsigned int	bx_
unsigned int	expTowerID_
unsigned int	expXtalTSamples_
std::unique_ptr< EcalElectronicsIdCollection > *	invalidMemBlockSizes_
std::unique_ptr< EcalElectronicsIdCollection > *	invalidMemChIds_
std::unique_ptr< EcalElectronicsIdCollection > *	invalidMemGains_
std::unique_ptr< EcalElectronicsIdCollection > *	invalidMemTtIds_
unsigned int	kSamplesPerPn_
unsigned int	l1_
unsigned int	lastStripId_
unsigned int	lastTowerBeforeMem_
unsigned int	lastXtalId_
unsigned int	nTSamples_
unsigned int	numbDWInXtalBlock_
std::vector< short >	pn_
std::unique_ptr< EcalPnDiodeDigiCollection > *	pnDiodeDigis_
unsigned int	towerId_
unsigned int	unfilteredTowerBlockLength_
unsigned int	xtalBlockSize_
Protected Attributes inherited from DCCDataBlockPrototype
unsigned int	blockLength_
const uint64_t *	data_
const uint64_t **	datap_
unsigned int *	dwToEnd_
bool	error_
DCCEventBlock *	event_
EcalElectronicsMapper *	mapper_
bool	sync_
DCCDataUnpacker *	unpacker_
bool	unpackInternalData_

Detailed Description

Definition at line 34 of file DCCMemBlock.h.

Constructor & Destructor Documentation

DCCMemBlock()

DCCMemBlock::DCCMemBlock	(	DCCDataUnpacker *	u,
		EcalElectronicsMapper *	m,
		DCCEventBlock *	e
	)

Definition at line 7 of file DCCMemBlock.cc.

References expXtalTSamples_, EcalElectronicsMapper::getUnfilteredTowerBlockLength(), mps_fire::i, kSamplesPerPn_, DCCDataBlockPrototype::mapper_, numbDWInXtalBlock_, EcalElectronicsMapper::numbXtalTSamples(), pn_, unfilteredTowerBlockLength_, and xtalBlockSize_.

    : DCCDataBlockPrototype(u, m, e) {
  unfilteredTowerBlockLength_ = mapper_->getUnfilteredTowerBlockLength();
  expXtalTSamples_ = mapper_->numbXtalTSamples();

  numbDWInXtalBlock_ = (expXtalTSamples_ - 2) / 4 + 1;
  xtalBlockSize_ = numbDWInXtalBlock_ * 8;
  kSamplesPerPn_ = expXtalTSamples_ * 5;

  unsigned int numbOfXtalBlocks = (unfilteredTowerBlockLength_ - 1) / numbDWInXtalBlock_;
  unsigned int numbOfPnBlocks = numbOfXtalBlocks / 5;  //change 5 by a variable
  unsigned int vectorSize = numbOfPnBlocks * 10 * expXtalTSamples_;

  //Build pnDiodevector
  for (unsigned int i = 0; i < vectorSize; i++) {
    pn_.push_back(-1);
  }
}

~DCCMemBlock()

DCCMemBlock::~DCCMemBlock ( )

inlineoverride

Definition at line 38 of file DCCMemBlock.h.

{}

Member Function Documentation

display()

void DCCMemBlock::display ( std::ostream &  o )

overridevirtual

Reimplemented from DCCDataBlockPrototype.

Definition at line 315 of file DCCMemBlock.cc.

References DCCDataBlockPrototype::blockLength_, bx_, l1_, nTSamples_, EcalTangentSkim_cfg::o, and towerId_.

                                       {
  o << "\n Unpacked Info for DCC MEM Block"
    << "\n DW1 ============================="
    << "\n Mem Tower Block Id " << towerId_ << "\n Numb Samp " << nTSamples_ << "\n Bx " << bx_ << "\n L1 " << l1_
    << "\n blockLength " << blockLength_;
}

fillPnDiodeDigisCollection()

void DCCMemBlock::fillPnDiodeDigisCollection ( )

protected

Definition at line 260 of file DCCMemBlock.cc.

References EcalBarrel, EcalEndcap, EcalElectronicsMapper::getActiveSM(), kSamplesPerPn_, lastTowerBeforeMem_, DCCDataBlockPrototype::mapper_, NUMB_SM_EB_MIN_MIN, NUMB_SM_EB_PLU_MAX, NUMB_SM_EE_MIN_MAX, NUMB_SM_EE_MIN_MIN, NUMB_SM_EE_PLU_MAX, NUMB_SM_EE_PLU_MIN, pn_, EcalPnDiodeDigi::setSample(), EcalPnDiodeDigi::setSize(), DCCDataUnpacker::silentMode_, and towerId_.

Referenced by unpack().

                                             {
  //todo change pnId max
  for (int pnId = 1; pnId <= 5; pnId++) {
    bool errorOnPn(false);
    unsigned int realPnId = pnId;

    if (towerId_ == 70) {
      realPnId += 5;
    }

    // Note : we are assuming always 5 VFE channels enabled
    // This means we all have 5 pns per tower

    // solution before sending creation of PnDigi's in mapper as done with crystals
    //     mapper_->getActiveSM()  : this is the 'dccid'
    //     number ranging internally in ECAL from 1 to 54, according convention specified here:
    //     http://indico.cern.ch/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=11621

    //     mapper_->getActiveDCC() : this is the FED_id (601 - 654 for ECAL at CMS)

    const int activeSM = mapper_->getActiveSM();
    int subdet(0);
    if (NUMB_SM_EB_MIN_MIN <= activeSM && activeSM <= NUMB_SM_EB_PLU_MAX) {
      subdet = EcalBarrel;
    } else if ((NUMB_SM_EE_MIN_MIN <= activeSM && activeSM <= NUMB_SM_EE_MIN_MAX) ||
               (NUMB_SM_EE_PLU_MIN <= activeSM && activeSM <= NUMB_SM_EE_PLU_MAX)) {
      subdet = EcalEndcap;
    } else {
      if (!DCCDataUnpacker::silentMode_) {
        edm::LogWarning("IncorrectMapping") << "\n mapper points to non existing dccid: " << activeSM;
      }
    }

    EcalPnDiodeDetId PnId(subdet, activeSM, realPnId);

    EcalPnDiodeDigi thePnDigi(PnId);
    thePnDigi.setSize(kSamplesPerPn_);

    for (unsigned int ts = 0; ts < kSamplesPerPn_; ts++) {
      short pnDiodeData = pn_[(towerId_ - lastTowerBeforeMem_) * 250 + (pnId - 1) * kSamplesPerPn_ + ts];
      if (pnDiodeData == -1) {
        errorOnPn = true;
        break;
      }

      EcalFEMSample thePnSample(pnDiodeData);
      thePnDigi.setSample(ts, thePnSample);
    }

    if (!errorOnPn) {
      (*pnDiodeDigis_)->push_back(thePnDigi);
    }
  }
}

unpack() [1/2]

virtual int DCCDataBlockPrototype::unpack

inline

Definition at line 40 of file DCCDataBlockPrototype.h.

{ return BLOCK_UNPACKED; }

unpack() [2/2]

int DCCMemBlock::unpack	(	const uint64_t **	data,
		unsigned int *	dwToEnd,
		unsigned int	expectedTowerID
	)

Definition at line 34 of file DCCMemBlock.cc.

References BLOCK_UNPACKED, DCCDataBlockPrototype::blockLength_, DCCEventBlock::bx(), bx_, data, DCCDataBlockPrototype::data_, DCCDataBlockPrototype::datap_, DCCDataBlockPrototype::dwToEnd_, DCCDataBlockPrototype::error_, DCCDataBlockPrototype::event_, expTowerID_, expXtalTSamples_, FE_MEM, fillPnDiodeDigisCollection(), DCCEventBlock::fov(), EcalElectronicsMapper::getActiveDCC(), EcalElectronicsMapper::getActiveSM(), H_FOV_MASK, l1ctLayer2EG_cff::id, isSynced(), l1_, DCCEventBlock::l1A(), lastStripId_, lastXtalId_, DCCDataBlockPrototype::mapper_, nTSamples_, DCCEventBlock::setFESyncNumbers(), DCCDataUnpacker::silentMode_, SKIP_BLOCK_UNPACKING, STOP_EVENT_UNPACKING, DCCDataBlockPrototype::sync_, TOWER_BX_B, TOWER_BX_MASK, TOWER_ID_MASK, TOWER_L1_B, TOWER_L1_MASK, TOWER_LENGTH_B, TOWER_LENGTH_MASK, TOWER_NSAMP_B, TOWER_NSAMP_MASK, towerId_, unfilteredTowerBlockLength_, unpackMemTowerData(), and DCCDataBlockPrototype::updateEventPointers().

Referenced by DCCEEEventBlock::unpack(), and DCCEBEventBlock::unpack().

                                                                                                  {
  error_ = false;
  datap_ = data;
  data_ = *data;
  dwToEnd_ = dwToEnd;

  if ((*dwToEnd_) < 1) {
    if (!DCCDataUnpacker::silentMode_) {
      edm::LogWarning("IncorrectEvent") << "\nUnable to unpack MEM block for event " << event_->l1A() << " in fed "
                                        << mapper_->getActiveDCC() << "\nThe end of event was reached !";
    }
    return STOP_EVENT_UNPACKING;
  }

  lastStripId_ = 0;
  lastXtalId_ = 0;
  expTowerID_ = expectedTowerID;

  //Point to begin of block
  data_++;

  towerId_ = (*data_) & TOWER_ID_MASK;
  nTSamples_ = (*data_ >> TOWER_NSAMP_B) & TOWER_NSAMP_MASK;
  bx_ = (*data_ >> TOWER_BX_B) & TOWER_BX_MASK;
  l1_ = (*data_ >> TOWER_L1_B) & TOWER_L1_MASK;
  blockLength_ = (*data_ >> TOWER_LENGTH_B) & TOWER_LENGTH_MASK;

  event_->setFESyncNumbers(l1_, bx_, short(expectedTowerID - 1));

  //debugging
  //display(cout);

  // Block Length Check (1)
  if (unfilteredTowerBlockLength_ != blockLength_) {
    // chosing channel 1 as representative of a dummy...
    EcalElectronicsId id(mapper_->getActiveSM(), expTowerID_, 1, 1);
    (*invalidMemBlockSizes_)->push_back(id);
    if (!DCCDataUnpacker::silentMode_) {
      edm::LogWarning("IncorrectEvent") << "\nFor event " << event_->l1A() << ", fed " << mapper_->getActiveDCC()
                                        << " and tower block " << towerId_ << "\nExpected mem block size is "
                                        << (unfilteredTowerBlockLength_ * 8) << " bytes while " << (blockLength_ * 8)
                                        << " was found";
    }
    return STOP_EVENT_UNPACKING;
  }

  // Block Length Check (2)
  if ((*dwToEnd_) < blockLength_) {
    if (!DCCDataUnpacker::silentMode_) {
      edm::LogWarning("IncorrectEvent") << "\nUnable to unpack MEM block for event " << event_->l1A() << " in fed "
                                        << mapper_->getActiveDCC() << "\n Only " << ((*dwToEnd_) * 8)
                                        << " bytes are available while " << (blockLength_ * 8) << " are needed!";
      // chosing channel 1 as representative of a dummy...
    }
    EcalElectronicsId id(mapper_->getActiveSM(), expTowerID_, 1, 1);
    (*invalidMemBlockSizes_)->push_back(id);
    return STOP_EVENT_UNPACKING;
  }

  // Synchronization Check
  if (sync_) {
    const unsigned int dccBx = (event_->bx()) & TOWER_BX_MASK;
    const unsigned int dccL1 = (event_->l1A()) & TOWER_L1_MASK;
    const unsigned int fov = (event_->fov()) & H_FOV_MASK;

    if (!isSynced(dccBx, bx_, dccL1, l1_, FE_MEM, fov)) {
      if (!DCCDataUnpacker::silentMode_) {
        edm::LogWarning("IncorrectEvent")
            << "Synchronization error for Mem block"
            << " (L1A " << event_->l1A() << " bx " << event_->bx() << " fed " << mapper_->getActiveDCC() << ")\n"
            << "  dccBx = " << dccBx << " bx_ = " << bx_ << " dccL1 = " << dccL1 << " l1_ = " << l1_ << "\n"
            << "  => Stop event unpacking";
      }
      //Note : add to error collection ?
      // need of a new collection
      return STOP_EVENT_UNPACKING;
    }
  }

  // Number Of Samples Check
  if (nTSamples_ != expXtalTSamples_) {
    if (!DCCDataUnpacker::silentMode_) {
      edm::LogWarning("IncorrectEvent") << "\nUnable to unpack MEM block for event " << event_->l1A() << " in fed "
                                        << mapper_->getActiveDCC() << "\nNumber of time samples " << nTSamples_
                                        << " is not the same as expected (" << expXtalTSamples_ << ")";
    }
    //Note : add to error collection ?
    return STOP_EVENT_UNPACKING;
  }

  //Channel Id Check
  if (expTowerID_ != towerId_) {
    // chosing channel 1 as representative as a dummy...
    EcalElectronicsId id(mapper_->getActiveSM(), expTowerID_, 1, 1);
    (*invalidMemTtIds_)->push_back(id);
    if (!DCCDataUnpacker::silentMode_) {
      edm::LogWarning("IncorrectBlock") << "For event " << event_->l1A() << " and fed " << mapper_->getActiveDCC()
                                        << " and sm: " << mapper_->getActiveSM() << "\nExpected mem tower block is "
                                        << expTowerID_ << " while " << towerId_ << " was found ";
    }

    towerId_ = expTowerID_;

    // todo : go to the next mem
    error_ = true;

    updateEventPointers();
    return SKIP_BLOCK_UNPACKING;
  }

  //point to xtal data
  data_++;

  unpackMemTowerData();

  if (!error_) {
    fillPnDiodeDigisCollection();
  }

  updateEventPointers();

  return BLOCK_UNPACKED;
}

unpackMemTowerData()

void DCCMemBlock::unpackMemTowerData ( )

protected

Definition at line 158 of file DCCMemBlock.cc.

References DCCDataBlockPrototype::data_, DCCDataBlockPrototype::error_, DCCDataBlockPrototype::event_, PedestalClient_cfi::gain, EcalElectronicsMapper::getActiveDCC(), EcalElectronicsMapper::getActiveSM(), mps_fire::i, cuy::ib, l1ctLayer2EG_cff::id, DCCEventBlock::l1A(), lastTowerBeforeMem_, DCCDataBlockPrototype::mapper_, nTSamples_, numbDWInXtalBlock_, pn_, ecalGpuTask_cfi::sample, DCCDataUnpacker::silentMode_, groupFilesInBlocks::temp, TOWER_DIGI_MASK, TOWER_STRIPID_MASK, TOWER_XTALID_B, TOWER_XTALID_MASK, and towerId_.

Referenced by unpack().

                                     {
  //todo: move EcalPnDiodeDetId to electronics mapper

  lastTowerBeforeMem_ = 0;
  // differentiating the barrel and the endcap case
  if (9 < mapper_->getActiveSM() || mapper_->getActiveSM() < 46) {
    lastTowerBeforeMem_ = 69;
  } else {
    lastTowerBeforeMem_ = 69;
  }

  for (unsigned int expStripId = 1; expStripId <= 5; expStripId++) {
    for (unsigned int expXtalId = 1; expXtalId <= 5; expXtalId++) {
      const uint16_t* xData_ = reinterpret_cast<const uint16_t*>(data_);

      // Get xtal data ids
      unsigned int stripId = (*xData_) & TOWER_STRIPID_MASK;
      unsigned int xtalId = ((*xData_) >> TOWER_XTALID_B) & TOWER_XTALID_MASK;

      bool errorOnDecoding(false);

      if (expStripId != stripId || expXtalId != xtalId) {
        // chosing channel and strip as EcalElectronicsId
        EcalElectronicsId id(mapper_->getActiveSM(), towerId_, expStripId, expXtalId);
        (*invalidMemChIds_)->push_back(id);

        if (!DCCDataUnpacker::silentMode_) {
          edm::LogWarning("IncorrectBlock")
              << "For event " << event_->l1A() << ", fed " << mapper_->getActiveDCC() << " and tower mem block "
              << towerId_ << "\nThe expected strip is " << expStripId << " and " << stripId << " was found"
              << "\nThe expected xtal  is " << expXtalId << " and " << xtalId << " was found";
        }

        stripId = expStripId;
        xtalId = expXtalId;

        errorOnDecoding = true;

        //Note : move to the next ...
      }

      unsigned int ipn, index;

      if ((stripId - 1) % 2 == 0) {
        ipn = (towerId_ - lastTowerBeforeMem_) * 5 + xtalId - 1;
      } else {
        ipn = (towerId_ - lastTowerBeforeMem_) * 5 + 5 - xtalId;
      }

      //Cooking samples
      for (unsigned int i = 0; i < nTSamples_; i++) {
        xData_++;

        index = ipn * 50 + (stripId - 1) * nTSamples_ + i;

        //edm::LogDebug("EcalRawToDigiMemChId")<<"\n Strip id "<<std::dec<<stripId<<" Xtal id "<<xtalId
        //  <<" tsamp = "<<i<<" 16b = 0x "<<std::hex<<(*xData_)<<dec;

        unsigned int temp = (*xData_) & TOWER_DIGI_MASK;

        short sample(0);

        if ((stripId - 1) % 2) {
          // If strip number is even, 14 bits are reversed in order
          for (int ib = 0; ib < 14; ib++) {
            sample <<= 1;
            sample |= (temp & 1);
            temp >>= 1;
          }

        } else {
          sample = temp;
        }

        sample ^= 0x800;
        unsigned int gain = sample >> 12;

        if (gain >= 2) {
          EcalElectronicsId id(mapper_->getActiveSM(), towerId_, stripId, xtalId);
          (*invalidMemGains_)->push_back(id);

          if (!DCCDataUnpacker::silentMode_) {
            edm::LogWarning("IncorrectGain")
                << "For event " << event_->l1A() << ", fed " << mapper_->getActiveDCC() << " , mem tower block "
                << towerId_ << "\nIn strip " << stripId << " xtal " << xtalId << " the gain is " << gain
                << " in sample " << (i + 1);
          }

          errorOnDecoding = true;
        }

        if (!errorOnDecoding && !error_) {
          pn_[index] = sample;
        }  //Note : move to the next versus flag...

      }  // loop over samples ended

      data_ += numbDWInXtalBlock_;
    }  //loop over xtals
  }    // loop over strips
}

updateCollectors()

void DCCMemBlock::updateCollectors ( )

overridevirtual

Reimplemented from DCCDataBlockPrototype.

Definition at line 26 of file DCCMemBlock.cc.

References invalidMemBlockSizes_, DCCDataUnpacker::invalidMemBlockSizesCollection(), invalidMemChIds_, DCCDataUnpacker::invalidMemChIdsCollection(), invalidMemGains_, DCCDataUnpacker::invalidMemGainsCollection(), invalidMemTtIds_, DCCDataUnpacker::invalidMemTtIdsCollection(), pnDiodeDigis_, DCCDataUnpacker::pnDiodeDigisCollection(), and DCCDataBlockPrototype::unpacker_.

Referenced by DCCEventBlock::updateCollectors().

                                   {
  invalidMemChIds_ = unpacker_->invalidMemChIdsCollection();
  invalidMemBlockSizes_ = unpacker_->invalidMemBlockSizesCollection();
  invalidMemTtIds_ = unpacker_->invalidMemTtIdsCollection();
  invalidMemGains_ = unpacker_->invalidMemGainsCollection();
  pnDiodeDigis_ = unpacker_->pnDiodeDigisCollection();
}

Member Data Documentation

bx_

unsigned int DCCMemBlock::bx_

protected

Definition at line 66 of file DCCMemBlock.h.

Referenced by display(), and unpack().

expTowerID_

unsigned int DCCMemBlock::expTowerID_

protected

Definition at line 52 of file DCCMemBlock.h.

Referenced by unpack().

expXtalTSamples_

unsigned int DCCMemBlock::expXtalTSamples_

protected

Definition at line 53 of file DCCMemBlock.h.

Referenced by DCCMemBlock(), and unpack().

invalidMemBlockSizes_

std::unique_ptr<EcalElectronicsIdCollection>* DCCMemBlock::invalidMemBlockSizes_

protected

Definition at line 70 of file DCCMemBlock.h.

Referenced by updateCollectors().

invalidMemChIds_

std::unique_ptr<EcalElectronicsIdCollection>* DCCMemBlock::invalidMemChIds_

protected

Definition at line 69 of file DCCMemBlock.h.

Referenced by updateCollectors().

invalidMemGains_

std::unique_ptr<EcalElectronicsIdCollection>* DCCMemBlock::invalidMemGains_

protected

Definition at line 72 of file DCCMemBlock.h.

Referenced by updateCollectors().

invalidMemTtIds_

std::unique_ptr<EcalElectronicsIdCollection>* DCCMemBlock::invalidMemTtIds_

protected

Definition at line 71 of file DCCMemBlock.h.

Referenced by updateCollectors().

kSamplesPerPn_

unsigned int DCCMemBlock::kSamplesPerPn_

protected

Definition at line 54 of file DCCMemBlock.h.

Referenced by DCCMemBlock(), and fillPnDiodeDigisCollection().

l1_

unsigned int DCCMemBlock::l1_

protected

Definition at line 67 of file DCCMemBlock.h.

Referenced by display(), and unpack().

lastStripId_

unsigned int DCCMemBlock::lastStripId_

protected

Definition at line 56 of file DCCMemBlock.h.

Referenced by unpack().

lastTowerBeforeMem_

unsigned int DCCMemBlock::lastTowerBeforeMem_

protected

Definition at line 58 of file DCCMemBlock.h.

Referenced by fillPnDiodeDigisCollection(), and unpackMemTowerData().

lastXtalId_

unsigned int DCCMemBlock::lastXtalId_

protected

Definition at line 57 of file DCCMemBlock.h.

Referenced by unpack().

nTSamples_

unsigned int DCCMemBlock::nTSamples_

protected

Definition at line 63 of file DCCMemBlock.h.

Referenced by display(), unpack(), and unpackMemTowerData().

numbDWInXtalBlock_

unsigned int DCCMemBlock::numbDWInXtalBlock_

protected

Definition at line 61 of file DCCMemBlock.h.

Referenced by DCCMemBlock(), and unpackMemTowerData().

pn_

std::vector<short> DCCMemBlock::pn_

protected

Definition at line 50 of file DCCMemBlock.h.

Referenced by DCCMemBlock(), fillPnDiodeDigisCollection(), and unpackMemTowerData().

pnDiodeDigis_

std::unique_ptr<EcalPnDiodeDigiCollection>* DCCMemBlock::pnDiodeDigis_

protected

Definition at line 73 of file DCCMemBlock.h.

Referenced by updateCollectors().

towerId_

unsigned int DCCMemBlock::towerId_

protected

Definition at line 60 of file DCCMemBlock.h.

Referenced by display(), fillPnDiodeDigisCollection(), unpack(), and unpackMemTowerData().

unfilteredTowerBlockLength_

unsigned int DCCMemBlock::unfilteredTowerBlockLength_

protected

Definition at line 64 of file DCCMemBlock.h.

Referenced by DCCMemBlock(), and unpack().

xtalBlockSize_

unsigned int DCCMemBlock::xtalBlockSize_

protected

Definition at line 62 of file DCCMemBlock.h.

Referenced by DCCMemBlock().