DCCEETCCBlock Class Reference

#include <DCCEETCCBlock.h>

Inheritance diagram for DCCEETCCBlock:

Public Member Functions
void	addTriggerPrimitivesToCollection () override
	DCCEETCCBlock (DCCDataUnpacker *u, EcalElectronicsMapper *m, DCCEventBlock *e, bool unpacking)
unsigned int	getLength () override
void	updateCollectors () override
Public Member Functions inherited from DCCTCCBlock
	DCCTCCBlock (DCCDataUnpacker *u, EcalElectronicsMapper *m, DCCEventBlock *e, bool unpack)
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, short tccChId=0)
Public Member Functions inherited from DCCDataBlockPrototype
	DCCDataBlockPrototype (DCCDataUnpacker *unpacker, EcalElectronicsMapper *mapper, DCCEventBlock *event, bool unpack=true)
void	enableSyncChecks ()
virtual int	unpack (const uint64_t **data, unsigned int *dwToEnd)
virtual void	updateEventPointers ()
virtual	~DCCDataBlockPrototype ()

Protected Member Functions
bool	checkTccIdAndNumbTTs () override

Additional Inherited Members
Protected Attributes inherited from DCCTCCBlock
unsigned int	bx_
unsigned int	expNumbTTs_
unsigned int	expTccId_
unsigned int	l1_
unsigned int	nTSamples_
unsigned int	nTTs_
EcalPseudoStripInputDigi *	pPS_
unsigned int	ps_
std::unique_ptr< EcalPSInputDigiCollection > *	pss_
EcalTriggerPrimitiveDigi *	pTP_
unsigned int	tccId_
std::unique_ptr< EcalTrigPrimDigiCollection > *	tps_
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 29 of file DCCEETCCBlock.h.

Constructor & Destructor Documentation

DCCEETCCBlock()

DCCEETCCBlock::DCCEETCCBlock	(	DCCDataUnpacker *	u,
		EcalElectronicsMapper *	m,
		DCCEventBlock *	e,
		bool	unpacking
	)

Class constructor

Definition at line 6 of file DCCEETCCBlock.cc.

References DCCDataBlockPrototype::blockLength_.

    : DCCTCCBlock(u, m, e, unpack) {
  blockLength_ = 0;
}

Member Function Documentation

addTriggerPrimitivesToCollection()

void DCCEETCCBlock::addTriggerPrimitivesToCollection ( )

overridevirtual

Reimplemented from DCCTCCBlock.

Definition at line 16 of file DCCEETCCBlock.cc.

References DCCEventBlock::bx(), DCCDataBlockPrototype::data_, dcc_FOV_0, dcc_FOV_1, dcc_FOV_2, DCCDataBlockPrototype::event_, DCCEventBlock::fov(), EcalElectronicsMapper::getActiveDCC(), EcalElectronicsMapper::getPSInputDigiPointer(), EcalElectronicsMapper::getTPPointer(), mps_fire::i, EcalElectronicsMapper::isTCCExternal(), DCCEventBlock::l1A(), DCCDataBlockPrototype::mapper_, DCCTCCBlock::nTTs_, NUMB_PSEUDOSTRIPS, NUMB_TTS_TPG1, NUMB_TTS_TPG2, NUMB_TTS_TPG2_DUPL, HLT_IsoTrack_cff::offset, DCCTCCBlock::pPS_, DCCTCCBlock::ps_, DCCTCCBlock::pTP_, EcalTriggerPrimitiveDigi::setSample(), EcalPseudoStripInputDigi::setSampleValue(), DCCDataUnpacker::silentMode_, and DCCTCCBlock::tccId_.

                                                     {
  //point to trigger data
  data_++;

  bool processTPG2(true);
  unsigned int psInputCounter(0);

  const uint16_t* tccP_ = reinterpret_cast<const uint16_t*>(data_);

  int dccFOV = event_->fov();
  if (!(dccFOV == dcc_FOV_0 || dccFOV == dcc_FOV_1 || dccFOV == dcc_FOV_2)) {
    if (!DCCDataUnpacker::silentMode_) {
      edm::LogWarning("IncorrectEvent") << "\n FOV value in data is: " << dccFOV << "At event: " << event_->l1A()
                                        << " with bx " << event_->bx() << " in fed <<" << mapper_->getActiveDCC()
                                        << "\n TCC id " << tccId_ << " FOV " << dccFOV
                                        << " which is not a foreseen value. Setting it to: " << dcc_FOV_2;
    }
    dccFOV = dcc_FOV_2;
  }

  // MC raw data based on CMS NOTE 2005/021
  // (and raw data when FOV was unassigned, earlier than mid 2008)
  if (dccFOV == dcc_FOV_0) {
    // Unpack TPG1 pseudostrip input block
    if (ps_) {
      for (unsigned int i = 1; i <= NUMB_PSEUDOSTRIPS; i++, tccP_++, psInputCounter++) {
        pPS_ = mapper_->getPSInputDigiPointer(tccId_, psInputCounter);
        if (!pPS_)
          continue;
        pPS_->setSampleValue(0, *tccP_);
        (*pss_)->push_back(*pPS_);
      }
    }

    // Unpack TPG1 trigger primitive block
    // loop over tp_counter=i and navigate forward in raw (tccP_)
    for (unsigned int i = 1; i <= NUMB_TTS_TPG1; i++, tccP_++) {  //16

      if (i <= nTTs_) {
        pTP_ = mapper_->getTPPointer(tccId_, i);  // pointer to tp digi container
        if (pTP_)
          pTP_->setSample(0, *tccP_);  // fill it
      } else {
        processTPG2 = false;
        break;
      }
      // adding trigger primitive digi to the collection
      if (pTP_)
        (*tps_)->push_back(*pTP_);
      else
        edm::LogError("IncorrectBlock") << "trigger primitive digi was not aquired";
    }

    if (processTPG2) {
      // Unpack TPG2 pseudostrip input block
      if (ps_) {
        for (unsigned int i = 1; i <= NUMB_PSEUDOSTRIPS; i++, tccP_++, psInputCounter++) {
          if (i <= NUMB_TTS_TPG2_DUPL)
            continue;
          //fill pseudostrip container
          pPS_ = mapper_->getPSInputDigiPointer(tccId_, psInputCounter);
          if (!pPS_)
            continue;
          pPS_->setSampleValue(0, *tccP_);
          (*pss_)->push_back(*pPS_);
        }
      }

      // Unpack TPG2 trigger primitive block
      for (unsigned int i = 1; i <= NUMB_TTS_TPG2; i++, tccP_++) {  //12
        unsigned int tt = i + NUMB_TTS_TPG1;

        if (tt <= nTTs_) {
          pTP_ = mapper_->getTPPointer(tccId_, tt);
          if (pTP_)
            pTP_->setSample(0, *tccP_);
        } else
          break;
        // adding trigger primitive digi to the collection
        (*tps_)->push_back(*pTP_);
      }

    }  // end if(processTPG2)

  }  // end FOV==0

  // real data since ever FOV was initialized; only 2 used >= June 09
  else if (dccFOV == dcc_FOV_1 || dccFOV == dcc_FOV_2) {
    // Unpack TPG1 pseudostrip input block
    if (ps_) {
      for (unsigned int i = 1; i <= NUMB_PSEUDOSTRIPS; i++, tccP_++, psInputCounter++) {
        pPS_ = mapper_->getPSInputDigiPointer(tccId_, psInputCounter);
        if (!pPS_)
          continue;
        pPS_->setSampleValue(0, *tccP_);
        (*pss_)->push_back(*pPS_);
      }
    }

    int offset(0);
    // Unpack TPG1 trigger primitive block
    // loop over tp_counter=i and navigate forward in raw (tccP_)
    for (unsigned int i = 1; i <= NUMB_TTS_TPG1; i++, tccP_++) {  //16

      if (i <= nTTs_) {
        if (mapper_->isTCCExternal(tccId_)) {
          if (i > 8 && i <= 16)
            continue;  // skip blank tp's [9,16]
          if (i > 8)
            offset = 8;  //
          if (i > 24)
            continue;  // skip blank tp's [25, 28]
        }

        pTP_ = mapper_->getTPPointer(tccId_, i - offset);  // pointer to tp digi container
        if (pTP_)
          pTP_->setSample(0, *tccP_);  // fill it

      } else {
        processTPG2 = false;
        break;
      }
      // adding trigger primitive digi to the collection
      if (pTP_)
        (*tps_)->push_back(*pTP_);
      else
        edm::LogError("IncorrectBlock") << "trigger primitive digi was not aquired";
    }

    if (processTPG2) {
      // Unpack TPG2 pseudostrip input block
      if (ps_) {
        for (unsigned int i = 1; i <= NUMB_PSEUDOSTRIPS; i++, tccP_++, psInputCounter++) {
          if (i <= NUMB_TTS_TPG2_DUPL)
            continue;
          //fill pseudostrip container
          pPS_ = mapper_->getPSInputDigiPointer(tccId_, psInputCounter);
          if (!pPS_)
            continue;
          pPS_->setSampleValue(0, *tccP_);
          (*pss_)->push_back(*pPS_);
        }
      }

      // Unpack TPG2 trigger primitive block
      for (unsigned int i = 1; i <= NUMB_TTS_TPG2; i++, tccP_++) {  //12
        unsigned int tt = i + NUMB_TTS_TPG1;

        if (i <= nTTs_) {
          if (mapper_->isTCCExternal(tccId_)) {
            if (tt > 8 && tt <= 16)
              continue;  // skip blank tp's [9,16]
            if (tt > 8)
              offset = 8;  //
            if (tt > 24)
              continue;  // skip blank tp's [25, 28]
          }

          pTP_ = mapper_->getTPPointer(tccId_, tt - offset);
          if (pTP_)
            pTP_->setSample(0, *tccP_);
        } else
          break;
        // adding trigger primitive digi to the collection
        (*tps_)->push_back(*pTP_);
      }
    }  // end if(processTPG2)
  }  // end FOV==1 or 2
}

checkTccIdAndNumbTTs()

bool DCCEETCCBlock::checkTccIdAndNumbTTs ( )

overrideprotectedvirtual

Reimplemented from DCCTCCBlock.

Definition at line 188 of file DCCEETCCBlock.cc.

References DCCEventBlock::bx(), DCCDataBlockPrototype::event_, EcalElectronicsMapper::getActiveDCC(), EcalElectronicsMapper::getActiveSM(), EcalElectronicsMapper::getTccs(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, DCCEventBlock::l1A(), visualization-live-secondInstance_cfg::m, DCCDataBlockPrototype::mapper_, DCCTCCBlock::nTTs_, DCCDataUnpacker::silentMode_, and DCCTCCBlock::tccId_.

                                         {
  bool tccFound(false);
  bool errorOnNumbOfTTs(false);
  const int activeDCC = mapper_->getActiveSM();
  std::vector<unsigned int>* m = mapper_->getTccs(activeDCC);
  std::vector<unsigned int>::iterator it;
  for (it = m->begin(); it != m->end(); it++) {
    if ((*it) == tccId_) {
      tccFound = true;

      /*
            expNumbTTs_= 28; //separate from inner and outer tcc 
            For inner TCCs you expect 28 TTs (4 in phi, 7 in eta).
            For outer TCCs you expect 16 TTs (4 in phi, 4 in eta).
            to implement : map tccid-> number of tts

            Pascal P.: For outer board, only 16 TPs are meaningful but 28 are sent.
            - [1,8] are meaningful 
            - [9,16] are empty ie should be skipped by unpacker 
            - [17, 24] are meaningful 
            - [25, 28] are empty ie should be skipped by unpacker 
          */

      if (nTTs_ != 28 && nTTs_ != 16) {
        if (!DCCDataUnpacker::silentMode_) {
          edm::LogWarning("IncorrectBlock")
              << "Error on event " << event_->l1A() << " with bx " << event_->bx() << " in fed <<"
              << mapper_->getActiveDCC() << "\n TCC id " << tccId_ << " has " << nTTs_
              << " Trigger Towers (only 28 or 16 are the expected values in EE)"
              << "\n => Skipping to next fed block...";
          errorOnNumbOfTTs = true;
          //todo : add to error collection
        }
      }
    }
  }

  if (!tccFound) {
    if (!DCCDataUnpacker::silentMode_) {
      edm::LogWarning("IncorrectBlock") << "Error on event " << event_->l1A() << " with bx " << event_->bx()
                                        << " in fed <<" << mapper_->getActiveDCC() << "\n TCC id " << tccId_
                                        << " is not valid for this dcc "
                                        << "\n => Skipping to next fed block...";
      //todo : add to error collection
    }
  }

  return (tccFound || errorOnNumbOfTTs);
}

getLength()

unsigned int DCCEETCCBlock::getLength ( )

overridevirtual

Reimplemented from DCCDataBlockPrototype.

Definition at line 238 of file DCCEETCCBlock.cc.

References B_MASK, DCCDataBlockPrototype::data_, NUMB_PSEUDOSTRIPS, NUMB_TTS_TPG1, NUMB_TTS_TPG2, DCCTCCBlock::ps_, TCC_PS_B, and groupFilesInBlocks::temp.

                                      {
  const uint64_t* temp = data_;
  temp++;

  ps_ = (*temp >> TCC_PS_B) & B_MASK;

  unsigned int numbTps = (NUMB_TTS_TPG1 + NUMB_TTS_TPG2);
  if (ps_) {
    numbTps += 2 * NUMB_PSEUDOSTRIPS;
  }

  unsigned int length = numbTps / 4 + 2;  //header and trailer
  if (numbTps % 4)
    length++;

  return length;
}

updateCollectors()

void DCCEETCCBlock::updateCollectors ( )

overridevirtual

Reimplemented from DCCDataBlockPrototype.

Definition at line 11 of file DCCEETCCBlock.cc.

References DCCDataUnpacker::ecalPSsCollection(), DCCDataUnpacker::ecalTpsCollection(), DCCTCCBlock::pss_, DCCTCCBlock::tps_, and DCCDataBlockPrototype::unpacker_.

                                     {
  tps_ = unpacker_->ecalTpsCollection();
  pss_ = unpacker_->ecalPSsCollection();
}