DCCEBTCCBlock.cc

Go to the documentation of this file.

#include "EventFilter/EcalRawToDigi/interface/DCCEBTCCBlock.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "EventFilter/EcalRawToDigi/interface/EcalElectronicsMapper.h"
#include "EventFilter/EcalRawToDigi/interface/DCCDataUnpacker.h"

DCCEBTCCBlock::DCCEBTCCBlock ( DCCDataUnpacker  * u,  EcalElectronicsMapper * m, DCCEventBlock * e, bool unpack) 
: DCCTCCBlock(u,m,e,unpack)
{
  blockLength_ = mapper_->getEBTCCBlockLength();
  expNumbTTs_  = TCC_EB_NUMBTTS;
}

void DCCEBTCCBlock::updateCollectors(){
  tps_ = unpacker_->ecalTpsCollection();
}

bool DCCEBTCCBlock::checkTccIdAndNumbTTs(){

  expTccId_ = mapper_->getActiveSM()+TCCID_SMID_SHIFT_EB;

  if( tccId_ != expTccId_ ){

    if( ! DCCDataUnpacker::silentMode_ ){
      edm::LogWarning("IncorrectBlock")
        <<"Error on event "<<event_->l1A()<<" with bx "<<event_->bx()<<" in fed "<<mapper_->getActiveDCC()
        <<"\n TCC id is "<<tccId_<<" while expected is "<<expTccId_
        <<"\n TCC Block Skipped ...";  
         //todo : add this to error colection
     }
     return false;
  }
  
  //Check number of TT Flags
  if( nTTs_ != expNumbTTs_ ){
    if( ! DCCDataUnpacker::silentMode_ ){
      edm::LogWarning("IncorrectBlock")
       <<"Unable to unpack TCC block for event "<<event_->l1A()<<" in fed "<<mapper_->getActiveDCC()
       <<"\n Number of TTs "<<nTTs_<<" is different from expected "<<expNumbTTs_
       <<"\n TCC Block Skipped ..."; 
         //todo : add this to error colection
     }
     return false;
  }  
  return true;
}


void DCCEBTCCBlock::addTriggerPrimitivesToCollection(){

  //point to trigger data
  data_++;

  unsigned int towersInPhi = EcalElectronicsMapper::kTowersInPhi;
  
  const uint16_t * tccP_= reinterpret_cast<const uint16_t * >(data_);
 

  for( unsigned int i = 1; i <= expNumbTTs_; i++){

    unsigned int theTT = i;

    if(NUMB_SM_EB_PLU_MIN<= mapper_->getActiveSM() && mapper_->getActiveSM()<=NUMB_SM_EB_PLU_MAX)
    {
        unsigned int u = (i-1)%towersInPhi;
        u      = towersInPhi-u;
        theTT  = ( (i-1)/towersInPhi )*towersInPhi + u;
    }
   
    pTP_ =  mapper_->getTPPointer(tccId_,theTT);
    for(unsigned int ns = 0; ns<nTSamples_;ns++,tccP_++){
      
      pTP_->setSampleValue(ns, (*tccP_));
      (*tps_)->push_back(*pTP_);
                  
    }
  }


}