![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
Go to the documentation of this file.
11 #include "GaudiKernel/MsgStream.h"
17 #include "GaudiKernel/Bootstrap.h"
18 #include "GaudiKernel/ISvcLocator.h"
23 #define MYLEVEL (MSG::FATAL)
24 #define LARBSDBG(text) logstr<<MYLEVEL<<text<<endmsg
26 #define LARBSDBG(text)
41 m_onlineHelper(nullptr)
54 ISvcLocator* svcLoc = Gaudi::svcLocator( );
57 std::cout <<
"Unable to locate DetectorStore" << std::endl;
62 std::cout <<
"Could not get LArOnlineID helper !" << std::endl;
132 int off1_tmp, dim1_tmp;
133 int off2_tmp, dim2_tmp;
134 int off3_tmp, dim3_tmp;
138 off1_tmp = dim1_tmp = 0;
139 off2_tmp = dim2_tmp = 0;
140 off3_tmp = dim3_tmp = 0;
151 dim1_tmp = 83+(
ns+1)/2+
n1;
154 off2_tmp = off1_tmp+dim1_tmp;
155 dim2_tmp = (n2*
ns+1)/2;
156 off3_tmp = off2_tmp+dim2_tmp;
158 if(dim2_tmp==0) off2_tmp = 0;
159 if(dim3_tmp==0) off3_tmp = 0;
162 if(off1 != off1_tmp) problem=1;
163 if(dim1 != dim1_tmp) problem=2;
164 if(off2 != off2_tmp) problem=3;
165 if(dim2 != dim2_tmp) problem=4;
166 if(off3 != off3_tmp) problem=5;
167 if(dim3 != dim3_tmp) problem=6;
168 if(
n1 != n1_tmp) problem=7;
169 if(n2 != n2_tmp) problem=8;
177 std::cout <<
"LArByteStreamProblem " << problem << std::endl;
178 std::cout <<
"NSamples = " << std::dec <<
ns << std::endl;
179 std::cout <<
"getHeader32(NGains) = " << std::hex <<
getHeader32(
NGains) << std::endl;
180 std::cout <<
"NWTot: " << std::hex <<
getNumberOfWords() <<
" n1=" <<
n1 <<
" (" << n1_tmp <<
") n2=" << n2 <<
" (" << n2_tmp <<
")" << std::endl;
181 std::cout <<
"Found 1: " << off1 <<
" " << dim1 << std::endl;
182 std::cout <<
"Found 2: " << off2 <<
" " << dim2 << std::endl;
183 std::cout <<
"Found 3: " << off3 <<
" " << dim3 << std::dec << std::endl;
185 if(
n1==n1_tmp && n2==n2_tmp) {
192 std::cout <<
"Recomputed 1: " << std::hex << off1 <<
" " << dim1 << std::endl;
193 std::cout <<
"Recomputed 2: " << off2 <<
" " << dim2 << std::endl;
194 std::cout <<
"Recomputed 3: " << off3 <<
" " << dim3 << std::dec << std::endl;
224 if(off1< off2 && off1 + dim1 > off2) problem = 1;
225 if(off1< off3 && off1 + dim1 > off3) problem = 2;
226 if(off2< off1 && off2 + dim2 > off1) problem = 3;
227 if(off2< off3 && off2 + dim2 > off3) problem = 4;
228 if(off3< off1 && off3 + dim3 > off1) problem = 5;
229 if(off3< off2 && off3 + dim3 > off2) problem = 6;
233 std::cout <<
"LArByteStreamProblem " << problem << std::endl;
234 std::cout <<
"Unrecoverable problem" << std::endl;
300 #ifdef LARBSDBGOUTPUT
303 logstr <<
MYLEVEL <<
"Let s go in getNextRawData..." <<
endmsg;
310 #ifdef LARBSDBGOUTPUT
311 logstr <<
MYLEVEL <<
"Maximum number of channels reached" <<
endmsg;
323 channelNumber=((rodChannelNumber&0xe)<<2) + ((rodChannelNumber&0x1)<<6) + (rodChannelNumber>>4);
329 #ifdef LARBSDBGOUTPUT
331 logstr <<
MYLEVEL <<
"This FEB has " << ngains <<
" gains" <<
endmsg;
334 if(ngains==0 ||
nsamples==0)
return 0;
381 #ifdef LARBSDBGOUTPUT
383 logstr <<
MYLEVEL <<
" ===> FEB Channel = " << channelNumber <<
endmsg;
394 unsigned rearrangeFirstSample=0;
400 if (rearrangeFirstSample && rearrangeFirstSample<
samples.size())
403 for (
unsigned i=1;
i<=rearrangeFirstSample;
i++)
405 samples[rearrangeFirstSample]=movedSample;
407 #ifdef LARBSDBGOUTPUT
416 #ifdef LARBSDBGOUTPUT
419 logstr <<
MYLEVEL <<
"Let s go in getNextDigits..." <<
endmsg;
427 #ifdef LARBSDBGOUTPUT
428 logstr <<
MYLEVEL <<
"Maximum number of channels reached" <<
endmsg;
434 #ifdef LARBSDBGOUTPUT
435 logstr <<
MYLEVEL <<
"No Digits Block in this FEB" <<
endmsg;
440 #ifdef LARBSDBGOUTPUT
441 logstr <<
MYLEVEL <<
"No Mask Digits Block in this FEB" <<
endmsg;
451 while(hasDigits==0) {
454 #ifdef LARBSDBGOUTPUT
455 logstr <<
MYLEVEL <<
"Maximum number of channels reached" <<
endmsg;
464 channelNumber=((rodChannelNumber&0xe)<<2) + ((rodChannelNumber&0x1)<<6) + (rodChannelNumber>>4);
472 }
else gain=0xffffffff;
474 #ifdef LARBSDBGOUTPUT
516 #ifdef LARBSDBGOUTPUT
518 logstr <<
MYLEVEL <<
" ===> FEB Channel = " << channelNumber <<
endmsg;
526 unsigned rearrangeFirstSample=0;
532 if (rearrangeFirstSample && rearrangeFirstSample<
samples.size())
535 for (
unsigned i=1;
i<=rearrangeFirstSample;
i++)
537 samples[rearrangeFirstSample]=movedSample;
539 #ifdef LARBSDBGOUTPUT
562 for(
int j=0;j<32;j++)
572 for(
int j=0;j<32;j++)
614 if(
adc>=8)
return x>>8;
679 unsigned int sizeRead=0;
681 while (sizeRead<fragment.size()) {
683 FebIter=fragment.begin()+sizeRead;
687 if (FebIter+currFebSize>fragment.end()) {
692 m_mFebBlocks[currFEBid].assign(FebIter,FebIter+currFebSize);
759 thequality = (
int16_t) quality;
766 if ((abse>8192)&&(abse<65536))
768 EncodedE=((abse>>3)|0x4000);
770 else if ((abse>65535)&&(abse<524288))
772 EncodedE=((abse>>6)|0x8000);
774 else if ((abse>524288))
776 EncodedE=((abse>>9)|0xc000);
781 if (
sign<0) EncodedE |= 0x2000;
890 for(
unsigned int i=0;
i<
n;
i++){
902 for (
unsigned i=0;
i<
n;
i++){
915 for (
unsigned i=0;
i<
n;
i++){
926 radd_nANC = (radd_nANC<<16);
947 for(
unsigned int i=0;
i<
n/2;
i++) {
961 unsigned int imax =
n/2;
962 for(
unsigned int i=0;
i<
imax;
i++){
976 for(
unsigned int i=0;
i<
imax;
i++){
989 for(
size_t ii=0;ii<
endtag/2;ii++){
990 sum+=((*m_vFragment)[ii]);
1006 FEBMAPTYPE::const_iterator feb_it_b=
m_mFebBlocks.begin();
1007 FEBMAPTYPE::const_iterator feb_it_e=
m_mFebBlocks.end();
1008 FEBMAPTYPE::const_iterator feb_it;
1009 for (feb_it=feb_it_b;feb_it!=feb_it_e;++feb_it) {
1010 if (feb_it!=feb_it_b)
1015 feb_it->second.begin(), feb_it->second.end());
1023 template<
class RAWDATA>
1024 bool LArRodBlockPhysicsV5::operator ()
1025 (
const RAWDATA* ch1,
const RAWDATA* ch2)
const
1034 int cId1 = m_onlineHelper->channel(id1);
1035 int cId2 = m_onlineHelper->channel(
id2);
1036 return FebToRodChannel(cId1) < FebToRodChannel(cId2);
1043 #ifdef LARBSDBGOUTPUT
1044 #undef LARBSDBGOUTPUT
JetConstituentVector::iterator iterator
virtual uint16_t getRawDataSize() const
virtual uint16_t getNbSweetCells2() const
std::vector< uint32_t > m_SumBlkBlockE2
std::vector< uint32_t > * m_vFragment
void setRawData(const int, const std::vector< short > &, const uint32_t)
virtual uint32_t getNumberOfSamples() const
singleton-like access to IMessageSvc via open function and helper
const LArOnlineID * m_onlineHelper
uint32_t getNumberOfWords() const
const uint32_t * m_GainPointer
const uint16_t * m_DigitsPointer
uint16_t m_EnergyThreshold2
virtual uint16_t getCtrl2(uint32_t adc) const
const uint32_t * m_FebBlock
void setHeader16(const unsigned n, const uint16_t w)
virtual int getNextRawData(int &channelNumber, std::vector< short > &samples, uint32_t &gain)
virtual uint32_t getNumberOfGains() const
setScaleOne setStatusOne setSaturated int16_t
std::vector< uint32_t > * m_pRODblock
IMessageSvc * getMessageSvc(bool quiet=false)
uint32_t OfflineToRawGain(const uint32_t gain) const
void setNextEnergy(const int channel, const int32_t energy, const int32_t time, const int32_t quality, const uint32_t gain)
std::vector< uint16_t > m_DigitsEncode
const uint16_t * m_TimeQualityPointer
void initializeFragment(std::vector< uint32_t > &fragment)
unsigned short m_iHeadBlockSize
unsigned int m_rearrangeFirstSample
const int32_t * m_SumPointer
uint32_t getHeader32(const unsigned n) const
The Athena Transient Store API.
virtual uint16_t getResults1Size() const
uint32_t RawToOfflineGain(const uint32_t gain) const
std::vector< uint32_t > m_GainBlock
virtual void resetPointers()
::StatusCode StatusCode
StatusCode definition for legacy code.
int getNextDigits(int &channelNumber, std::vector< short > &samples, uint32_t &gain)
const uint16_t * m_EnergyPointer
unsigned short m_requiredNSamples
const uint32_t * m_MaskDigitsPointer
const uint32_t * m_MaskTimeQualityPointer
uint16_t m_numberHotCellOffTime
std::vector< uint16_t > m_TimeQualityBlock
uint16_t getNbSweetCells2FromMask() const
const uint16_t * m_RawDataPointer
uint16_t getFirstSampleIndex() const
virtual int getNextEnergy(int &channelNumber, int32_t &energy, int32_t &time, int32_t &quality, uint32_t &gain)
virtual uint32_t getRadd(uint32_t adc, uint32_t sample) const
virtual uint32_t getStatus() const
virtual uint16_t getCtrl1(uint32_t adc) const
virtual uint16_t getCtrl3(uint32_t adc) const
int32_t m_MiddleHeaderSize
def time(flags, cells_name, *args, **kw)
const uint16_t * m_RaddPointer
std::vector< uint16_t > m_EnergyBlockEncode
void initializeFEB(const uint32_t id)
uint16_t m_EnergyThreshold1
virtual uint16_t getNbSweetCells1() const
std::vector< uint32_t > m_SumBlkBlockE1
uint16_t getNbSweetCells1FromMask() const
virtual int FebToRodChannel(int ch) const
uint16_t getHeader16(const unsigned n) const
virtual bool setPointers()
void setHeader32(const unsigned n, const uint32_t w)
virtual uint16_t getResults2Size() const