|
|
|
#include <DCCTowerBlock.h>
Public Member Functions | |
| DCCTowerBlock (DCCDataUnpacker *u, EcalElectronicsMapper *m, DCCEventBlock *e, bool unpack, bool forceToKeepFRdata) | |
| void | updateCollectors () |
Protected Member Functions | |
| void | fillEcalElectronicsError (std::auto_ptr< EcalElectronicsIdCollection > *) |
| int | unpackXtalData (unsigned int stripID, unsigned int xtalID) |
Protected Attributes | |
| std::auto_ptr< EBDigiCollection > * | digis_ |
| std::auto_ptr < EBDetIdCollection > * | invalidChIds_ |
| std::auto_ptr < EBDetIdCollection > * | invalidGains_ |
| std::auto_ptr < EBDetIdCollection > * | invalidGainsSwitch_ |
| EBDetId * | pDetId_ |
Definition at line 22 of file DCCTowerBlock.h.
| DCCTowerBlock::DCCTowerBlock | ( | DCCDataUnpacker * | u, |
| EcalElectronicsMapper * | m, | ||
| DCCEventBlock * | e, | ||
| bool | unpack, | ||
| bool | forceToKeepFRdata | ||
| ) |
Definition at line 10 of file DCCTowerBlock.cc.
: DCCFEBlock(u,m,e,unpack,forceToKeepFRdata){}
| void DCCTowerBlock::fillEcalElectronicsError | ( | std::auto_ptr< EcalElectronicsIdCollection > * | errorColection | ) | [protected, virtual] |
Reimplemented from DCCFEBlock.
Definition at line 259 of file DCCTowerBlock.cc.
References DCCDataBlockPrototype::event_, DCCFEBlock::expTowerID_, EcalElectronicsMapper::getActiveSM(), EcalElectronicsMapper::getTTEleIdPointer(), DCCEventBlock::l1A(), DCCDataBlockPrototype::mapper_, NUMB_SM_EB_MIN_MIN, NUMB_SM_EB_PLU_MAX, DCCDataUnpacker::silentMode_, and TCCID_SMID_SHIFT_EB.
Referenced by unpackXtalData().
{
const int activeDCC = mapper_->getActiveSM();
if(NUMB_SM_EB_MIN_MIN<=activeDCC && activeDCC<=NUMB_SM_EB_PLU_MAX){
EcalElectronicsId * eleTp = mapper_->getTTEleIdPointer(activeDCC+TCCID_SMID_SHIFT_EB,expTowerID_);
(*errorColection)->push_back(*eleTp);
}else{
if( ! DCCDataUnpacker::silentMode_ ){
edm::LogWarning("IncorrectBlock")
<<"For event "<<event_->l1A()<<" there's fed: "<< activeDCC
<<" activeDcc: "<<mapper_->getActiveSM()
<<" but that activeDcc is not valid in EB.";
}
}
}
| int DCCTowerBlock::unpackXtalData | ( | unsigned int | stripID, |
| unsigned int | xtalID | ||
| ) | [protected, virtual] |
Reimplemented from DCCFEBlock.
Definition at line 29 of file DCCTowerBlock.cc.
References BLOCK_UNPACKED, DCCEventBlock::bx(), runTheMatrix::data, DCCDataBlockPrototype::data_, DCCDataBlockPrototype::event_, fillEcalElectronicsError(), EcalElectronicsMapper::getActiveDCC(), DCCDataUnpacker::getChannelValue(), EcalElectronicsMapper::getDetIdPointer(), i, DCCFEBlock::invalidZSXtalIds_, DCCEventBlock::l1A(), DCCFEBlock::lastStripId_, DCCFEBlock::lastXtalId_, DCCDataBlockPrototype::mapper_, min, DCCFEBlock::nTSamples_, NUMB_STRIP, NUMB_XTAL, DCCFEBlock::numbDWInXtalBlock_, pDetId_, asciidump::s, DCCDataUnpacker::silentMode_, SKIP_BLOCK_UNPACKING, TOWER_DIGI_MASK, TOWER_STRIPID_MASK, TOWER_XTALID_B, TOWER_XTALID_MASK, DCCFEBlock::towerId_, DCCDataBlockPrototype::unpacker_, DCCFEBlock::xtalGains_, and DCCFEBlock::zs_.
{
bool errorOnXtal(false);
uint16_t * xData_= reinterpret_cast<uint16_t *>(data_);
// Get xtal data ids
unsigned int stripId = (*xData_) & TOWER_STRIPID_MASK;
unsigned int xtalId = ((*xData_)>>TOWER_XTALID_B ) & TOWER_XTALID_MASK;
// check id in case data are not 0suppressed
if( !zs_ && (expStripID != stripId || expXtalID != xtalId)){
if(! DCCDataUnpacker::silentMode_){
edm::LogWarning("IncorrectBlock")
<<"For event L1A: "<<event_->l1A()<<", fed "<<mapper_->getActiveDCC()<<" and tower "<<towerId_
<<"\n The expected strip is "<<expStripID<<" and "<<stripId<<" was found"
<<"\n The expected xtal is "<<expXtalID <<" and "<<xtalId<<" was found";
}
// using expected cry_di to raise warning about xtal_id problem
pDetId_ = (EBDetId*) mapper_->getDetIdPointer(towerId_,expStripID,expXtalID);
(*invalidChIds_)->push_back(*pDetId_);
stripId = expStripID;
xtalId = expXtalID;
errorOnXtal= true;
// return here, so to skip all following checks
lastXtalId_++;
if (lastXtalId_ > NUMB_XTAL) {lastXtalId_=1; lastStripId_++;}
data_ += numbDWInXtalBlock_;
return BLOCK_UNPACKED;
}
// check id in case of 0suppressed data
else if(zs_){
// Check for valid Ids 1) values out of range
if(stripId == 0 || stripId > 5 || xtalId == 0 || xtalId > 5){
if( ! DCCDataUnpacker::silentMode_ ){
edm::LogWarning("IncorrectBlock")
<<"For event L1A: "<<event_->l1A()<<", fed "<<mapper_->getActiveDCC()<<" and tower "<<towerId_
<<"\n Invalid strip : "<<stripId<<" or xtal : "<<xtalId
<<" ids ( last strip was: " << lastStripId_ << " last ch was: " << lastXtalId_ << ")";
}
int st = lastStripId_;
int ch = lastXtalId_;
ch++;
if (ch > NUMB_XTAL) {ch=1; st++;}
if (st > NUMB_STRIP) {ch=1; st=1;}
// adding channel following the last valid
//pDetId_ = (EBDetId*) mapper_->getDetIdPointer(towerId_,st,ch);
//(*invalidChIds_)->push_back(*pDetId_);
fillEcalElectronicsError(invalidZSXtalIds_);
errorOnXtal = true;
lastStripId_ = st;
lastXtalId_ = ch;
// return here, so to skip all following checks
return SKIP_BLOCK_UNPACKING;
}else{
// Check for zs valid Ids 2) if channel-in-strip has increased wrt previous xtal
// Check for zs valid Ids 2) if channel-in-strip has increased wrt previous xtal
// 3) if strip has increased wrt previous xtal
if( ( stripId == lastStripId_ && xtalId <= lastXtalId_ ) ||
(stripId < lastStripId_))
{
if (! DCCDataUnpacker::silentMode_) {
edm::LogWarning("IncorrectBlock")
<< "Xtal id was expected to increase but it didn't - last valid xtal id was " << lastXtalId_ << " while current xtal is " << xtalId
<< " (LV1 " << event_->l1A() << " fed " << mapper_->getActiveDCC() << " tower " << towerId_ << ")";
}
int st = lastStripId_;
int ch = lastXtalId_;
ch++;
if (ch > NUMB_XTAL) {ch=1; st++;}
if (st > NUMB_STRIP) {ch=1; st=1;}
// adding channel following the last valid
//pDetId_ = (EBDetId*) mapper_->getDetIdPointer(towerId_,st,ch);
//(*invalidChIds_)->push_back(*pDetId_);
fillEcalElectronicsError(invalidZSXtalIds_);
errorOnXtal = true;
lastStripId_ = st;
lastXtalId_ = ch;
// return here, so to skip all following checks
return SKIP_BLOCK_UNPACKING;
}
// if channel id not proven wrong, update lastStripId_ and lastXtalId_
lastStripId_ = stripId;
lastXtalId_ = xtalId;
}//end else
}// end if (zs_)
bool addedFrame=false;
// if there is an error on xtal id ignore next error checks
// otherwise, assume channel_id is valid and proceed with making and checking the data frame
if(errorOnXtal) return SKIP_BLOCK_UNPACKING;
pDetId_ = (EBDetId*) mapper_->getDetIdPointer(towerId_,stripId,xtalId);
(*digis_)->push_back(*pDetId_);
EBDataFrame df( (*digis_)->back() );
addedFrame=true;
bool wrongGain(false);
//set samples in the data frame
for(unsigned int i =0; i< nTSamples_ ;i++){ // loop on samples
xData_++;
unsigned int data = (*xData_) & TOWER_DIGI_MASK;
unsigned int gain = data>>12;
xtalGains_[i]=gain;
if(gain == 0){
wrongGain = true;
// although gain==0 found, produce the dataFrame in order to have it, for saturation case
}
df.setSample(i,data);
}// loop on samples
bool isSaturation(true);
if(wrongGain){
// check whether the gain==0 has features of saturation or not
// gain==0 occurs either in case of data corruption or of ADC saturation
// \->reject digi \-> keep digi
// determine where gainId==0 starts
short firstGainZeroSampID(-1); short firstGainZeroSampADC(-1);
for (unsigned int s=0; s<nTSamples_; s++ ) {
if(df.sample(s).gainId()==0 && firstGainZeroSampID==-1)
{
firstGainZeroSampID = s;
firstGainZeroSampADC = df.sample(s).adc();
break;
}
}
// check whether gain==0 and adc() stays constant for (at least) 5 consecutive samples
unsigned int plateauEnd = std::min(nTSamples_,(unsigned int)(firstGainZeroSampID+5));
for (unsigned int s=firstGainZeroSampID; s<plateauEnd; s++)
{
if( df.sample(s).gainId()==0 && df.sample(s).adc()==firstGainZeroSampADC ) {;}
else
{ isSaturation=false; break;} //it's not saturation
}
// get rid of channels which are stuck in gain0
if(firstGainZeroSampID<3) {isSaturation=false; }
if (! DCCDataUnpacker::silentMode_) {
if (unpacker_->getChannelValue(mapper_->getActiveDCC(), towerId_, stripId, xtalId) != 10) {
edm::LogWarning("IncorrectBlock")
<< "Gain zero" << (isSaturation ? " with features of saturation" : "" ) << " was found in Tower Block"
<< " (L1A " << event_->l1A() << " bx " << event_->bx() << " fed " << mapper_->getActiveDCC()
<< " tower " << towerId_ << " strip " << stripId << " xtal " << xtalId << ")";
}
}
if (! isSaturation)
{
(*invalidGains_)->push_back(*pDetId_);
(*digis_)->pop_back();
errorOnXtal = true;
//Point to begin of next xtal Block
data_ += numbDWInXtalBlock_;
//return here, so to skip all the rest
//make special collection for gain0 data frames when due to saturation
return BLOCK_UNPACKED;
}//end isSaturation
else {
data_ += numbDWInXtalBlock_;
return BLOCK_UNPACKED;
}
}//end WrongGain
// from here on, care about gain switches
short firstGainWrong=-1;
short numGainWrong=0;
for (unsigned int i=1; i<nTSamples_; i++ ) {
if (i>0 && xtalGains_[i-1]>xtalGains_[i]) {
numGainWrong++;
if (firstGainWrong == -1) { firstGainWrong=i;}
}
}
if (numGainWrong > 0) {
if (! DCCDataUnpacker::silentMode_) {
edm::LogWarning("IncorrectBlock")
<< "A wrong gain transition switch was found for Tower Block in strip " << stripId << " and xtal " << xtalId
<< " (L1A " << event_->l1A() << " bx " << event_->bx() << " fed " << mapper_->getActiveDCC() << " tower " << towerId_ << ")";
}
(*invalidGainsSwitch_)->push_back(*pDetId_);
errorOnXtal = true;
}
//Add frame to collection only if all data format and gain rules are respected
if (errorOnXtal && addedFrame) {
(*digis_)->pop_back();
}
//Point to begin of next xtal Block
data_ += numbDWInXtalBlock_;
return BLOCK_UNPACKED;
}
| void DCCTowerBlock::updateCollectors | ( | ) | [virtual] |
Reimplemented from DCCFEBlock.
Definition at line 14 of file DCCTowerBlock.cc.
References digis_, DCCDataUnpacker::ebDigisCollection(), invalidChIds_, DCCDataUnpacker::invalidChIdsCollection(), invalidGains_, DCCDataUnpacker::invalidGainsCollection(), invalidGainsSwitch_, DCCDataUnpacker::invalidGainsSwitchCollection(), and DCCDataBlockPrototype::unpacker_.
{
DCCFEBlock::updateCollectors();
// needs to be update for eb/ee
digis_ = unpacker_->ebDigisCollection();
invalidGains_ = unpacker_->invalidGainsCollection();
invalidGainsSwitch_ = unpacker_->invalidGainsSwitchCollection();
invalidChIds_ = unpacker_->invalidChIdsCollection();
}
std::auto_ptr<EBDigiCollection>* DCCTowerBlock::digis_ [protected] |
Definition at line 37 of file DCCTowerBlock.h.
Referenced by updateCollectors().
std::auto_ptr<EBDetIdCollection>* DCCTowerBlock::invalidChIds_ [protected] |
Definition at line 44 of file DCCTowerBlock.h.
Referenced by updateCollectors().
std::auto_ptr<EBDetIdCollection>* DCCTowerBlock::invalidGains_ [protected] |
Definition at line 42 of file DCCTowerBlock.h.
Referenced by updateCollectors().
std::auto_ptr<EBDetIdCollection>* DCCTowerBlock::invalidGainsSwitch_ [protected] |
Definition at line 43 of file DCCTowerBlock.h.
Referenced by updateCollectors().
EBDetId* DCCTowerBlock::pDetId_ [protected] |
Definition at line 39 of file DCCTowerBlock.h.
Referenced by unpackXtalData().
1.7.3