ATLAS Offline Software
Loading...
Searching...
No Matches
LArRodBlockPhysicsV5 Class Reference

#include <LArRodBlockPhysicsV5.h>

Inheritance diagram for LArRodBlockPhysicsV5:
Collaboration diagram for LArRodBlockPhysicsV5:

Public Types

enum  {
  NWTot , NWTot_h , FEBID , FEBID_h ,
  FEB_SN , FEB_SN_h , ResultsOff1 , ResultsDim1 ,
  ResultsOff2 , ResultsDim2 , RawDataBlkOff , RawDataBlkDim ,
  EventStatus , EventStatus_h , NGains , NSamples ,
  FirstSampleIndex , FebConfig , InFPGAFormat , InFPGAFormat_h ,
  endtag
}

Public Member Functions

 LArRodBlockPhysicsV5 (IMessageSvc *msgSvc)
virtual int getNextEnergy (int &channelNumber, int32_t &energy, int32_t &time, int32_t &quality, uint32_t &gain)
virtual int getNextRawData (int &channelNumber, std::vector< short > &samples, uint32_t &gain)
int getNextDigits (int &channelNumber, std::vector< short > &samples, uint32_t &gain)
virtual uint16_t getResults1Size () const
virtual uint16_t getResults2Size () const
virtual uint16_t getRawDataSize () const
virtual uint16_t getNbSweetCells1 () const
virtual uint16_t getNbSweetCells2 () const
virtual uint32_t getNumberOfSamples () const
virtual uint32_t getNumberOfGains () const
virtual uint32_t getRadd (uint32_t adc, uint32_t sample) const
virtual uint16_t getCtrl1 (uint32_t adc) const
virtual uint16_t getCtrl2 (uint32_t adc) const
virtual uint16_t getCtrl3 (uint32_t adc) const
virtual uint32_t getStatus () const
virtual uint32_t hasCalibBlock () const
virtual uint32_t hasPhysicsBlock () const
virtual uint32_t hasRawDataBlock () const
virtual uint32_t hasControlWords () const
virtual int FebToRodChannel (int ch) const
virtual int32_t getEx () const
virtual int32_t getEy () const
virtual int32_t getEz () const
virtual uint32_t getVROBFebId ()
virtual int32_t getVROBEx () const
virtual int32_t getVROBEy () const
virtual int32_t getVROBEz () const
virtual int setGain (const int)
void initializeFragment (std::vector< uint32_t > &fragment)
void initializeFEB (const uint32_t id)
void setNextEnergy (const int channel, const int32_t energy, const int32_t time, const int32_t quality, const uint32_t gain)
void setRawData (const int, const std::vector< short > &, const uint32_t)
void finalizeFEB ()
void concatinateFEBs ()
virtual void setEx (double)
virtual void setEy (double)
virtual void setEz (double)
void setRequiredNSamples (unsigned short ns)
template<class RAWDATA>
bool operator() (const RAWDATA *ch1, const RAWDATA *ch2) const
void sortDataVector (std::vector< const LArDigit * > &vDigit)
uint16_t getFebConfig () const
uint16_t getFirstSampleIndex () const
virtual bool canSetCalibration ()
bool canSetEnergy ()
bool canSetRawData ()
bool canIncludeRawData ()
virtual void sortDataVector (std::vector< const LArRawChannel * > &)
virtual void sortDataVector (std::vector< const LArCalibDigit * > &)
virtual void sortDataVector (std::vector< const LArAccumulatedCalibDigit * > &)
virtual void sortDataVector (std::vector< const LArAccumulatedDigit * > &)
virtual void setNumberOfSamples (const uint8_t n)
virtual void setNumberOfGains (const uint8_t n)
virtual void setTDCPhase (const uint8_t n)
virtual void setRawDataFixed (const int channel, const std::vector< short > &samples, const uint32_t gain)
virtual void setEtQ (const int channel, const int32_t energy, const int32_t time, const int32_t quality, const uint32_t gain)
virtual void setDAC (const uint16_t DACValue)
virtual void setDelay (const uint16_t DelayValue)
virtual void setPulsed (const unsigned channelNumber)
virtual void setNTrigger (const uint16_t NTrigger)
virtual void setSumE (double)
virtual bool canSetRawDataFixed ()
virtual bool canSetNTrigger ()
bool setFragment (const uint32_t *p, uint32_t n)
bool nextFEB ()
uint32_t getNumberOfWords () const
uint32_t getFEBID () const
uint32_t getFEBSN () const
virtual int32_t getSumE () const
virtual int32_t getVROBSumE () const
virtual uint32_t hasAccumBlock () const
virtual int getNextAccumulatedCalibDigit (int &channelNumber, std::vector< uint64_t > &SamplesSum, std::vector< uint64_t > &Samples2Sum, uint32_t &nStepTriggers, uint32_t &gain)
virtual int getNextAccumulatedDigit (int &channelNumber, std::vector< uint64_t > &SamplesSum, std::vector< uint64_t > &corr2Sum, uint32_t &gain)
virtual bool getPulsed (unsigned channelNumber) const
virtual uint16_t getDAC () const
virtual uint16_t getDelay () const
virtual uint16_t getNTrigger () const
virtual uint16_t getStepIndex () const
virtual uint16_t getNStep () const
virtual uint8_t getTDCPhase () const
virtual uint32_t getDspCodeVersion () const
virtual int32_t getDspEventCounter () const
virtual uint32_t onlineCheckSum () const
virtual uint32_t offlineCheckSum () const
virtual int setFragmentVirtualROB (const uint32_t *p, uint32_t n)
virtual void dumpFragment ()
void setFirstSample (const int rearrangeFirstSample)
uint32_t RawToOfflineGain (const uint32_t gain) const
uint32_t OfflineToRawGain (const uint32_t gain) const
bool report_error (void) const

Static Public Member Functions

static std::string BlockType ()

Protected Types

enum  {
  NWTot , NWTot_h , FEBID , FEBID_h ,
  FEBSN , FEBSN_h , endtag
}
typedef std::map< uint32_t, std::vector< uint32_t > > FEBMAPTYPE

Protected Member Functions

void setHeader16 (const unsigned n, const uint16_t w)
void setHeader32 (const unsigned n, const uint32_t w)
uint16_t getVectorHeader16 (const unsigned n) const
uint32_t getVectorHeader32 (const unsigned n) const
uint16_t getHeader16 (const unsigned n) const
uint32_t getHeader32 (const unsigned n) const
uint16_t LE_getHeader16 (const unsigned n) const
void LE_setHeader16 (const unsigned n, const uint16_t w)
uint16_t LE_getVectorHeader16 (const unsigned n) const
void setBit (uint32_t *const p, const unsigned chan)
int getBit (const uint32_t *const p, const unsigned chan) const

Protected Attributes

int32_t m_Ex
int32_t m_Ey
int32_t m_Ez
int32_t m_SumE
uint32_t m_numberOfEx
uint32_t m_numberOfEy
uint32_t m_numberOfEz
uint32_t m_numberOfSumE
unsigned short m_iHeadBlockSize
int m_channelsPerFEB
const uint32_t * m_FebBlock
const uint32_t * m_RodBlock
int32_t m_FebBlockSize
int32_t m_RodBlockSize
std::vector< uint32_t > * m_vFragment
std::vector< uint32_t > * m_pRODblock
FEBMAPTYPE m_mFebBlocks
int32_t m_MiddleHeaderSize
unsigned int m_rearrangeFirstSample
bool m_error_next_feb
uint32_t m_virtualROBJump
int32_t m_ROB_to_decode
const uint32_t * m_virtualROBPointer
const uint32_t * m_virtualROBPointerLocal
MsgStream m_logstr

Static Protected Attributes

static const uint32_t m_RawToOfflineGainMap [4] ={0, 2, 1,0}
static const uint32_t m_OfflineToRawGainMap [3] ={3,2,1}

Private Member Functions

void setE (unsigned index, double E)
virtual void resetPointers ()
virtual bool setPointers ()
void setNextEnergy (const uint16_t energy, const int16_t time, const int16_t quality, const uint32_t gain)
uint16_t getNbSweetCells1FromMask () const
uint16_t getNbSweetCells2FromMask () const

Private Attributes

std::vector< uint32_t > m_SumBlkBlockE1
std::vector< uint32_t > m_SumBlkBlockE2
std::vector< uint32_t > m_GainBlock
std::vector< uint32_t > m_FebInfoBlock
std::vector< uint16_t > m_TimeQualityBlock
std::vector< uint32_t > m_RawDataBlock
std::vector< uint16_t > m_EnergyBlockEncode
std::vector< uint16_t > m_DigitsEncode
int m_EnergyIndex = 0
int m_TimeQualityIndex = 0
int m_DigitsIndex = 0
int m_DigitsChannel = 0
int m_RawDataIndex = 0
const uint32_t * m_GainPointer = nullptr
const uint32_t * m_MaskTimeQualityPointer = nullptr
const uint32_t * m_MaskDigitsPointer = nullptr
const uint16_t * m_RaddPointer = nullptr
const uint16_t * m_EnergyPointer = nullptr
const int32_t * m_SumPointer = nullptr
const uint16_t * m_TimeQualityPointer = nullptr
const uint16_t * m_DigitsPointer = nullptr
const uint16_t * m_RawDataPointer = nullptr
int m_fixedGain
uint16_t m_numberHotCell = 0U
uint16_t m_numberHotCellOffTime
uint16_t m_EnergyThreshold1
uint16_t m_EnergyThreshold2
int16_t m_OffTimeCut
unsigned short m_requiredNSamples
const LArOnlineIDm_onlineHelper

Detailed Description

Definition at line 30 of file LArRodBlockPhysicsV5.h.

Member Typedef Documentation

◆ FEBMAPTYPE

typedef std::map<uint32_t, std::vector<uint32_t> > LArRodBlockStructure::FEBMAPTYPE
protectedinherited

Definition at line 227 of file LArRodBlockStructure.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
Enumerator
NWTot 
NWTot_h 
FEBID 
FEBID_h 
FEB_SN 
FEB_SN_h 
ResultsOff1 
ResultsDim1 
ResultsOff2 
ResultsDim2 
RawDataBlkOff 
RawDataBlkDim 
EventStatus 
EventStatus_h 
NGains 
NSamples 
FirstSampleIndex 
FebConfig 
InFPGAFormat 
InFPGAFormat_h 
endtag 

Definition at line 34 of file LArRodBlockPhysicsV5.h.

34 {
35 NWTot, // First words: DSP event header
36 NWTot_h,
37 FEBID,
38 FEBID_h,
39 FEB_SN, // FEB serial number
40 FEB_SN_h, // FEB serial number
41 ResultsOff1, // Size of results (Physics averages in DSP)
42 ResultsDim1, // Offset to results
43 ResultsOff2, // Size of times (in physics)
44 ResultsDim2, // Offset to times (in physics)
46 RawDataBlkDim, // Raw FEB event offset
47 EventStatus, // Bits describing the event
49 NGains,
50 NSamples,
53 InFPGAFormat, // added 08.09.2005 - wrong 28.09.2005?
55 endtag //This tag needs to be an odd number, see *) for constructor
56 };

◆ anonymous enum

anonymous enum
protectedinherited
Enumerator
NWTot 
NWTot_h 
FEBID 
FEBID_h 
FEBSN 
FEBSN_h 
endtag 

Definition at line 50 of file LArRodBlockStructure.h.

50 {
51 NWTot, // Number of words in DSP block
52 NWTot_h,
53 FEBID, // FEB ID
54 FEBID_h,
55 FEBSN, // FEB Serial Number
56 FEBSN_h,
57 endtag // This tag needs to be an odd number, see *) for constructor
58 };

Constructor & Destructor Documentation

◆ LArRodBlockPhysicsV5()

LArRodBlockPhysicsV5::LArRodBlockPhysicsV5 ( IMessageSvc * msgSvc)

Definition at line 34 of file LArRodBlockPhysicsV5.cxx.

36 m_onlineHelper(nullptr)
37{
38 m_iHeadBlockSize=endtag/2; // The implicit cast rounds down to the right size
44 m_OffTimeCut=0; //FIXME: Nowhere set to a sensible value ???
46 // retrieve onlineHelper
47 SmartIF<StoreGateSvc> detStore{Gaudi::svcLocator()->service("DetectorStore")};
48 if (!detStore) {
49 std::cout << "Unable to locate DetectorStore" << std::endl;
50 std::abort();
51 }
52 StatusCode sc = detStore->retrieve(m_onlineHelper, "LArOnlineID");
53 if (sc.isFailure()) {
54 std::cout << "Could not get LArOnlineID helper !" << std::endl;
55 std::abort();
56 }
57}
static Double_t sc
const LArOnlineID * m_onlineHelper
unsigned short m_requiredNSamples
static std::string BlockType()
LArRodBlockStructure(IMessageSvc *msgSvc, const std::string &blockType)
@ LARNGAIN
Definition CaloGain.h:19
::StatusCode StatusCode
StatusCode definition for legacy code.

Member Function Documentation

◆ BlockType()

std::string LArRodBlockPhysicsV5::BlockType ( )
inlinestatic

Definition at line 61 of file LArRodBlockPhysicsV5.h.

61{ return std::string("RodBlockPhysicsV5");}

◆ canIncludeRawData()

bool LArRodBlockPhysicsV5::canIncludeRawData ( )
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 173 of file LArRodBlockPhysicsV5.h.

173{ return true;}

◆ canSetCalibration()

virtual bool LArRodBlockPhysicsV5::canSetCalibration ( )
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 170 of file LArRodBlockPhysicsV5.h.

170{return false;}

◆ canSetEnergy()

bool LArRodBlockPhysicsV5::canSetEnergy ( )
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 171 of file LArRodBlockPhysicsV5.h.

171{ return true;}

◆ canSetNTrigger()

virtual bool LArRodBlockStructure::canSetNTrigger ( )
inlinevirtualinherited

Definition at line 99 of file LArRodBlockStructure.h.

99{return false;}

◆ canSetRawData()

bool LArRodBlockPhysicsV5::canSetRawData ( )
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 172 of file LArRodBlockPhysicsV5.h.

172{ return false;}

◆ canSetRawDataFixed()

virtual bool LArRodBlockStructure::canSetRawDataFixed ( )
inlinevirtualinherited

Reimplemented in LArRodBlockTransparentV0< DSPHEADER >.

Definition at line 97 of file LArRodBlockStructure.h.

97{return false;}

◆ concatinateFEBs()

void LArRodBlockPhysicsV5::concatinateFEBs ( )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 855 of file LArRodBlockPhysicsV5.cxx.

856{
857 FEBMAPTYPE::const_iterator feb_it_b=m_mFebBlocks.begin();
858 FEBMAPTYPE::const_iterator feb_it_e=m_mFebBlocks.end();
859 FEBMAPTYPE::const_iterator feb_it;
860 for (feb_it=feb_it_b;feb_it!=feb_it_e;++feb_it) {
861 if (feb_it!=feb_it_b) //Not first Feb
863
864 //Add feb data to rod data block
865 m_pRODblock->insert (m_pRODblock->end(),
866 feb_it->second.begin(), feb_it->second.end());
867 } //end for feb_it
868
869 m_mFebBlocks.clear();
870 return;
871}
std::vector< uint32_t > * m_pRODblock

◆ dumpFragment()

void LArRodBlockStructure::dumpFragment ( )
virtualinherited

Definition at line 173 of file LArRodBlockStructure.cxx.

174{
175 std::cout << "Error: Function dumpFragment() not implemented in this instance of LArRodBlockStructure!\n";
176 return;
177}

◆ FebToRodChannel()

int LArRodBlockPhysicsV5::FebToRodChannel ( int ch) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 316 of file LArRodBlockPhysicsV5.h.

317{
318return ( (ch&0x7) << 4) | ( (ch&0x38) >>2 ) | ((ch&0x40)>>6);
319
320}

◆ finalizeFEB()

void LArRodBlockPhysicsV5::finalizeFEB ( )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 724 of file LArRodBlockPhysicsV5.cxx.

725{
726//Complete non-complete Energy block
728 setNextEnergy((uint16_t)0,(int16_t)32767,(int32_t)-32767,(uint32_t)0);//E=0,t=32767,q=-32767,G=0
729
730uint16_t n;
732// checkSum value
733uint32_t sum=0;
734
735// Will Hardcode here for the moment FIXME. Minimal 1 sample
737setHeader16(NSamples,nsamples);
738// These will never be used form MC. Nice to put in here thought
739setHeader16(FEB_SN,0xfefe);
740setHeader16(FEB_SN_h,0xdede);
743
744// Gain block...
745n = m_GainBlock.size();
746//Check if Gain-Block exists and is not yet part of the fragment
747if (n)
748 {
749 for(unsigned int i=0;i<n;i++){
750 m_vFragment->push_back(m_GainBlock[i]);
751 sum+=m_GainBlock[i];
752 }
753 }
754
755 // Cells above energy threshold E1
756 n = m_SumBlkBlockE1.size();
757 //Check if Summary Block exists and is not yet part of the fragment
758 if (n)
759 {
760 for (unsigned i=0;i<n;i++){
761 m_vFragment->push_back(m_SumBlkBlockE1[i]);
763 }
764 }
765
766 // Cells above energy threshold E2 (not included so far)
767 n = m_SumBlkBlockE2.size();
768 //Check if Summary Block exists and is not yet part of the fragment
769 if (n)
770 {
771 for (unsigned i=0;i<n;i++){
772 m_vFragment->push_back(m_SumBlkBlockE2[i]);
774 }
775 }
776
777 // fill info from counters
778 // for moment just include 1 fake words (32 bits) to put radd
779 uint32_t radd_nANC=0x0;
780 // Second threshold missing (FIXME)
781 radd_nANC = ((m_numberHotCell<<8))+(m_DigitsEncode.size()/nsamples);
782 radd_nANC = (radd_nANC<<16);
783 m_vFragment->push_back(radd_nANC);
784 sum+=radd_nANC;
785 // Need to include radd nsamples-1
786 // No need to include in sum's for now
787 for( int i=0; i < (nsamples-1)/2; i++)
788 m_vFragment->push_back(0x0);
789
790
791 // Energy block...
792 n = 128 ; // Fixed size m_EnergyBlock.size();
793 // Block also include time, whenever necessary
794 int size_of_block=80+(nsamples+1)/2+(m_TimeQualityBlock.size())/2;
795 //Check if Energy-Block exists and is not yet part of the fragment
796 if (n)
797 {
799 setHeader16(ResultsDim1,size_of_block);
800 for(unsigned int i=0;i<n/2;i++) {
801 // WARNING which one should be >>16 2*i or 2*i+1? To be tested
802 uint32_t Encode = m_EnergyBlockEncode[2*i]+(m_EnergyBlockEncode[2*i+1]<<16);
803 m_vFragment->push_back(Encode);
804 sum+=Encode;
805 }
806 }
807
808 // Magic numbers (3 or 6) for Ex, Ey and Ez
809 n = m_TimeQualityBlock.size();
810 //Check if Time and Quality Block exists and is not yet part of the fragment
811 if (n)
812 {
813 unsigned int imax = n/2;
814 for(unsigned int i=0;i<imax;i++){
815 ShortLong to_push{};
816 to_push.s[0] = m_TimeQualityBlock[i*2];
817 to_push.s[1] = m_TimeQualityBlock[i*2+1];
818 m_vFragment->push_back(to_push.l);
819 sum+=to_push.l;
820 }
821 }
822 // Now include digits
823 n = m_DigitsEncode.size();
824 if ( n ) {
825 // First make sure it is not and odd number to store
826 if ( m_DigitsEncode.size() & 0x1 ) m_DigitsEncode.push_back(0x0);
827 unsigned int imax=m_DigitsEncode.size()/2;
828 for(unsigned int i=0;i<imax;i++){
829 // Better by-swap
830 ShortLong to_push{};
831 to_push.s[1]=m_DigitsEncode[i*2];
832 to_push.s[0]=m_DigitsEncode[i*2+1];
833 m_vFragment->push_back(to_push.l);
834 sum+=to_push.l;
835 }
837 setHeader16(ResultsOff2,18+size_of_block);
838 } // End of check for format
839
840 // Need to add header to check sum
841 for(size_t ii=0;ii<endtag/2;ii++){
842 sum+=((*m_vFragment)[ii]);
843 }
844 // Three final magic words
845 m_vFragment->push_back(0x0); // For the moment
846 m_vFragment->push_back(0x12345678); // For the moment
847 m_vFragment->push_back(sum& 0x7fffffff);
848
850 return;
851
852}
int imax(int i, int j)
std::vector< uint16_t > m_EnergyBlockEncode
std::vector< uint32_t > m_SumBlkBlockE1
std::vector< uint32_t > m_SumBlkBlockE2
std::vector< uint32_t > m_GainBlock
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
std::vector< uint16_t > m_TimeQualityBlock
void setHeader32(const unsigned n, const uint32_t w)
void setHeader16(const unsigned n, const uint16_t w)
std::vector< uint32_t > * m_vFragment
setWord1 uint16_t
setEventNumber uint32_t

◆ getBit()

int LArRodBlockStructure::getBit ( const uint32_t *const p,
const unsigned chan ) const
inlineprotectedinherited

Definition at line 416 of file LArRodBlockStructure.h.

417{// chan = (0,127)
418 // int a = chan/32;
419 // int r = chan%32;
420 int a = chan>>5;
421 int r = chan&0x1f;
422 // a = (0,3), r = ( 0, 31 )
423 return (*(p+a)>>r) & 0x1;
424}
static Double_t a
int r
Definition globals.cxx:22

◆ getCtrl1()

uint16_t LArRodBlockPhysicsV5::getCtrl1 ( uint32_t adc) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 523 of file LArRodBlockPhysicsV5.cxx.

524{
525 if(!m_RawDataPointer) return 0;
526 int index=5;
528 return x;
529}
#define x
const uint16_t * m_RawDataPointer
str index
Definition DeMoScan.py:362

◆ getCtrl2()

uint16_t LArRodBlockPhysicsV5::getCtrl2 ( uint32_t adc) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 531 of file LArRodBlockPhysicsV5.cxx.

532{
533 if(!m_RawDataPointer) return 0;
534 int index=4;
536 return x;
537}

◆ getCtrl3()

uint16_t LArRodBlockPhysicsV5::getCtrl3 ( uint32_t adc) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 539 of file LArRodBlockPhysicsV5.cxx.

540{
541 if(!m_RawDataPointer) return 0;
542 int index=7;
544 return x;
545}

◆ getDAC()

uint16_t LArRodBlockStructure::getDAC ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockCalibrationV0< DSPHEADER >, LArRodBlockCalibrationV1, LArRodBlockCalibrationV2, and LArRodBlockCalibrationV3.

Definition at line 451 of file LArRodBlockStructure.h.

452{
453 std::cout << "Error: Function getDAC not implemented in this instance of LArRodBlockStructure!\n";
454 return 0;
455}

◆ getDelay()

uint16_t LArRodBlockStructure::getDelay ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockCalibrationV0< DSPHEADER >, LArRodBlockCalibrationV1, LArRodBlockCalibrationV2, and LArRodBlockCalibrationV3.

Definition at line 457 of file LArRodBlockStructure.h.

458{
459 std::cout << "Error: Function getDelay not implemented in this instance of LArRodBlockStructure!\n";
460 return 0;
461}

◆ getDspCodeVersion()

uint32_t LArRodBlockStructure::getDspCodeVersion ( ) const
virtualinherited

Definition at line 256 of file LArRodBlockStructure.cxx.

257{
258 int i=getNumberOfWords()-2;
259 return m_FebBlock[i];
260}
uint32_t getNumberOfWords() const

◆ getDspEventCounter()

int32_t LArRodBlockStructure::getDspEventCounter ( ) const
virtualinherited

Definition at line 262 of file LArRodBlockStructure.cxx.

263{
264 int i=getNumberOfWords()-1;
265 return m_FebBlock[i];
266}

◆ getEx()

int32_t LArRodBlockPhysicsV5::getEx ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 266 of file LArRodBlockPhysicsV5.h.

267{
268 if(m_SumPointer) return (m_SumPointer[0]>>9);
269 return 0;
270}

◆ getEy()

int32_t LArRodBlockPhysicsV5::getEy ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 272 of file LArRodBlockPhysicsV5.h.

273{
274 if(m_SumPointer) return (m_SumPointer[1]>>9);
275 return 0;
276}

◆ getEz()

int32_t LArRodBlockPhysicsV5::getEz ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 278 of file LArRodBlockPhysicsV5.h.

279{
280 if(m_SumPointer) return (m_SumPointer[2]>>9);
281 return 0;
282}

◆ getFebConfig()

uint16_t LArRodBlockPhysicsV5::getFebConfig ( ) const
inline

Definition at line 176 of file LArRodBlockPhysicsV5.h.

177{
178 return getHeader16(FebConfig);
179}
uint16_t getHeader16(const unsigned n) const

◆ getFEBID()

uint32_t LArRodBlockStructure::getFEBID ( ) const
inlineinherited

Definition at line 292 of file LArRodBlockStructure.h.

293{return getHeader32(FEBID);}
uint32_t getHeader32(const unsigned n) const

◆ getFEBSN()

uint32_t LArRodBlockStructure::getFEBSN ( ) const
inlineinherited

Definition at line 295 of file LArRodBlockStructure.h.

296{return getHeader32(FEBSN);}

◆ getFirstSampleIndex()

uint16_t LArRodBlockPhysicsV5::getFirstSampleIndex ( ) const
inline

Definition at line 181 of file LArRodBlockPhysicsV5.h.

182{
184}

◆ getHeader16()

uint16_t LArRodBlockStructure::getHeader16 ( const unsigned n) const
inlineprotectedinherited

Definition at line 340 of file LArRodBlockStructure.h.

341{
342 if (n&0x1) //n is a odd number
343 return m_FebBlock[n>>1] & 0xffff; //1,3,5... are fetched from lower bits
344 else //n is a even number
345 return m_FebBlock[n>>1] >> 16; //0,2,4... are fetched from higher bits
346}

◆ getHeader32()

uint32_t LArRodBlockStructure::getHeader32 ( const unsigned n) const
inlineprotectedinherited

Definition at line 348 of file LArRodBlockStructure.h.

349{return m_FebBlock[n>>1];}

◆ getNbSweetCells1()

uint16_t LArRodBlockPhysicsV5::getNbSweetCells1 ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 450 of file LArRodBlockPhysicsV5.cxx.

451{
452 if(!m_RaddPointer) return 0;
453 return m_RaddPointer[1]>>8;
454}
const uint16_t * m_RaddPointer

◆ getNbSweetCells1FromMask()

uint16_t LArRodBlockPhysicsV5::getNbSweetCells1FromMask ( ) const
private

Definition at line 462 of file LArRodBlockPhysicsV5.cxx.

463{
464 if(!m_MaskTimeQualityPointer) return 0;
465 int n=0;
466 for(int i=0;i<4;i++)
467 for(int j=0;j<32;j++)
468 if((m_MaskTimeQualityPointer[i] >> j) &0x1) n++;
469 return n;
470}
const uint32_t * m_MaskTimeQualityPointer
float j(const xAOD::IParticle &, const xAOD::TrackMeasurementValidation &hit, const Eigen::Matrix3d &jab_inv)

◆ getNbSweetCells2()

uint16_t LArRodBlockPhysicsV5::getNbSweetCells2 ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 456 of file LArRodBlockPhysicsV5.cxx.

457{
458 if(!m_RaddPointer) return 0;
459 return m_RaddPointer[1] & 0xff;
460}

◆ getNbSweetCells2FromMask()

uint16_t LArRodBlockPhysicsV5::getNbSweetCells2FromMask ( ) const
private

Definition at line 472 of file LArRodBlockPhysicsV5.cxx.

473{
474 if(!m_MaskDigitsPointer) return 0;
475 int n=0;
476 for(int i=0;i<4;i++)
477 for(int j=0;j<32;j++)
478 if((m_MaskDigitsPointer[i] >> j) &0x1) n++;
479 return n;
480}
const uint32_t * m_MaskDigitsPointer

◆ getNextAccumulatedCalibDigit()

int LArRodBlockStructure::getNextAccumulatedCalibDigit ( int & channelNumber,
std::vector< uint64_t > & SamplesSum,
std::vector< uint64_t > & Samples2Sum,
uint32_t & nStepTriggers,
uint32_t & gain )
virtualinherited

Reimplemented in LArRodBlockCalibrationV1, LArRodBlockCalibrationV2, and LArRodBlockCalibrationV3.

Definition at line 92 of file LArRodBlockStructure.cxx.

93{
94 std::cout << "Error: Function getNextAccumulatedCalibDigit not implemented in this instance of LArRodBlockStructure!\n";
95 return 0;
96}

◆ getNextAccumulatedDigit()

int LArRodBlockStructure::getNextAccumulatedDigit ( int & channelNumber,
std::vector< uint64_t > & SamplesSum,
std::vector< uint64_t > & corr2Sum,
uint32_t & gain )
virtualinherited

Reimplemented in LArRodBlockAccumulatedV3.

Definition at line 98 of file LArRodBlockStructure.cxx.

99{
100 std::cout << "Error: Function getNextAccumulatedDigit not implemented in this instance of LArRodBlockStructure!\n";
101 return 0;
102}

◆ getNextDigits()

int LArRodBlockPhysicsV5::getNextDigits ( int & channelNumber,
std::vector< short > & samples,
uint32_t & gain )

Definition at line 335 of file LArRodBlockPhysicsV5.cxx.

336{
337#ifdef LARBSDBGOUTPUT
338 MsgStream logstr(Athena::getMessageSvc(), BlockType());
339 //Debug output
340 logstr << MYLEVEL << "Let s go in getNextDigits..." << endmsg;
341 logstr << MYLEVEL << "GetNextDigits for FEB 0x" << MSG::hex << (uint32_t)getHeader32(FEBID) << MSG::dec << endmsg;
342 logstr << MYLEVEL << "m_DigitsPointer=" << m_DigitsPointer << " m_DigitsIndex="<< m_DigitsIndex
343 << " m_DigitsChannel="<< m_DigitsChannel
344 << " m_channelsPerFEB=" << m_channelsPerFEB << endmsg;
345#endif
346
347 if (m_DigitsChannel>=m_channelsPerFEB) { //Already beyond maximal number of channels
348#ifdef LARBSDBGOUTPUT
349 logstr << MYLEVEL << "Maximum number of channels reached" << endmsg;
350#endif
351 return 0;
352 }
353 if (!m_DigitsPointer) { //Block does not exist
354#ifdef LARBSDBGOUTPUT
355 logstr << MYLEVEL << "No Digits Block in this FEB" << endmsg;
356#endif
357 return 0;
358 }
359 if (!m_MaskDigitsPointer) { //Block does not exist
360#ifdef LARBSDBGOUTPUT
361 logstr << MYLEVEL << "No Mask Digits Block in this FEB" << endmsg;
362#endif
363 return 0;
364 }
365
366 // Get Digits if the information is present according to summary block
367 uint32_t hasDigits;
368
369 hasDigits = (uint32_t) ((m_MaskDigitsPointer[m_DigitsChannel>>5] >> (m_DigitsChannel&0x1f)) &0x1);
370 // Increment channel number until digits are found
371 while(hasDigits==0) {
373 if (m_DigitsChannel>=m_channelsPerFEB) { //Already beyond maximal number of channels
374#ifdef LARBSDBGOUTPUT
375 logstr << MYLEVEL << "Maximum number of channels reached" << endmsg;
376#endif
377 return 0;
378 }
379 hasDigits = (uint32_t) ((m_MaskDigitsPointer[m_DigitsChannel>>5] >> (m_DigitsChannel&0x1f)) &0x1);
380 }
381
382 // Get next channel
383 unsigned rodChannelNumber=m_DigitsChannel; // Index of Channel in ROD-Block
384 channelNumber=((rodChannelNumber&0xe)<<2) + ((rodChannelNumber&0x1)<<6) + (rodChannelNumber>>4); //channel number of the FEB
385 const unsigned int nsamples = getHeader16(NSamples) & 0xff;
386
387 // gain in 2 bits of a 32 bits word
388 if(m_GainPointer) {
389 gain = (uint32_t) ((m_GainPointer[m_DigitsChannel>>4] >> (m_DigitsChannel&0xf)*2) & 0x3);
391 } else gain=0xffffffff;
392
393#ifdef LARBSDBGOUTPUT
394 logstr << MYLEVEL << "This FEB has " << nsamples << " samples" << endmsg;
395#endif
396
397 if(nsamples==0) return 0;
398 int s_size = nsamples;
399 int index;
400 index = s_size*m_DigitsIndex;
401 if(m_DigitsIndex&0x1) {
402 samples.push_back(m_DigitsPointer[index-1]>>2);
403 samples.push_back(m_DigitsPointer[index+2]>>2);
404 samples.push_back(m_DigitsPointer[index+1]>>2);
405 samples.push_back(m_DigitsPointer[index+4]>>2);
406 samples.push_back(m_DigitsPointer[index+3]>>2);
407 if(nsamples==7) {
408 samples.push_back(m_DigitsPointer[index+6]>>2);
409 samples.push_back(m_DigitsPointer[index+5]>>2);
410 }
411 } else {
412 samples.push_back(m_DigitsPointer[index+1]>>2);
413 samples.push_back(m_DigitsPointer[index+0]>>2);
414 samples.push_back(m_DigitsPointer[index+3]>>2);
415 samples.push_back(m_DigitsPointer[index+2]>>2);
416 samples.push_back(m_DigitsPointer[index+5]>>2);
417 if(nsamples==7) {
418 samples.push_back(m_DigitsPointer[index+4]>>2);
419 samples.push_back(m_DigitsPointer[index+7]>>2);
420 }
421 }
422
423#ifdef LARBSDBGOUTPUT
424 logstr << MYLEVEL << " ===> ROD Channel = " << m_DigitsChannel << endmsg;
425 logstr << MYLEVEL << " ===> FEB Channel = " << channelNumber << endmsg;
426 logstr << MYLEVEL << " ===> Gain = " << gain << endmsg;
427 for(int i=0;i<nsamples;i++)
428 logstr << MYLEVEL << " ===> sample " << i << " = " << samples[i] << endmsg;
429#endif
432 unsigned rearrangeFirstSample=0;
434 rearrangeFirstSample=m_rearrangeFirstSample; //Overwrite by jobOptions
435 else
436 rearrangeFirstSample=getFirstSampleIndex();
437 if (rearrangeFirstSample && rearrangeFirstSample<samples.size()) //FIXME: Very ugly hack! See explanation in LArRodDecoder.h file
438 {//Change e.g. 3 0 1 2 4 to 0 1 2 3 4
439 short movedSample=samples[0];
440 for (unsigned i=1;i<=rearrangeFirstSample;i++)
441 samples[i-1]=samples[i];
442 samples[rearrangeFirstSample]=movedSample;
443 }
444#ifdef LARBSDBGOUTPUT
445 logstr << MYLEVEL << "GetNextDigits for FEB finished 0x" << MSG::hex << (uint32_t)getHeader32(FEBID) << MSG::dec << endmsg;
446#endif
447 return 1;
448}
#define endmsg
#define MYLEVEL
const uint32_t * m_GainPointer
const uint16_t * m_DigitsPointer
uint16_t getFirstSampleIndex() const
uint32_t RawToOfflineGain(const uint32_t gain) const
IMessageSvc * getMessageSvc(bool quiet=false)

◆ getNextEnergy()

int LArRodBlockPhysicsV5::getNextEnergy ( int & channelNumber,
int32_t & energy,
int32_t & time,
int32_t & quality,
uint32_t & gain )
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 186 of file LArRodBlockPhysicsV5.h.

187{
188 #ifdef LARBSDBGOUTPUT
189 MsgStream logstr(Athena::getMessageSvc(), BlockType());
190 #endif
191
192 LARBSDBG("in LArRodBlockPhysicsV5::getNextEnergy.");
193 LARBSDBG("m_channelsPerFEB=" << m_channelsPerFEB);
194 if (m_EnergyIndex>=m_channelsPerFEB) // Already beyond maximal number of channels
195 return 0;
196 if (!m_EnergyPointer) // No data block present
197 return 0;
198
199 unsigned rodChannelNumber=m_EnergyIndex; // Index of Channel in ROD-Block
200 channelNumber=((rodChannelNumber&0xe)<<2) + ((rodChannelNumber&0x1)<<6) + (rodChannelNumber>>4); //channel number of the FEB
201
202 // get information available for all cells
203 // Energy on a 16 bit word and decode ranges
204 uint16_t encodedEnergy; // 16 bits Encoded Energy word
205 int32_t aux;
206 uint16_t range; // 2 bits range
208
209 LARBSDBG("-------->>>> in LArRodBlockPhysicsV5::getNextEnergy : decode energy.....");
210 // decode energy
211 if(m_EnergyIndex & 0x1) encodedEnergy = m_EnergyPointer[m_EnergyIndex-1]; // Big/Little Endien stuff
212 else encodedEnergy = m_EnergyPointer[m_EnergyIndex+1]; // Big/Little Endien stuff
213
214 aux = (int32_t) (encodedEnergy&0x1fff);
215 range = (encodedEnergy & 0xc000) >> 14;
216 if(aux==0 && range>0) aux=0x2000;
217 sign = encodedEnergy & 0x2000;
218 aux <<= 3*range;
219 if(sign) aux = -aux;
220 energy = aux;
221
222 // gain in 2 bits of a 32 bits word
223 if(m_GainPointer) {
224 gain = (uint32_t) ((m_GainPointer[m_EnergyIndex>>4] >> ((m_EnergyIndex&0xf)<<1)) & 0x3);
226 } else gain=0xffffffff;
227
228 // Get Time and Quality if the information is present according to summary block
229 uint32_t hasTQ;
231 hasTQ = (uint32_t) ((m_MaskTimeQualityPointer[m_EnergyIndex>>5] >> (m_EnergyIndex&0x1f)) &0x1);
232 else
233 hasTQ = 0;
235 if (m_TimeQualityPointer && hasTQ) // Data has Time and Quality information
236 {
237 //Time is in 10 ps in ByteStream, hence the factor 10 to convert to ps
238 time = 10*(reinterpret_cast<const int16_t *>(m_TimeQualityPointer))[m_TimeQualityIndex++];
240
241#ifdef LARBSDBGOUTPUT
242 logstr << MYLEVEL <<"This cell has time and Quality information "<<endmsg;
243#endif
244 }
245 else // Data has no Time and Quality information
246 {
247 time=0;
248 quality=-1;
249 }
250
251
252#ifdef LARBSDBGOUTPUT
253 logstr << MYLEVEL <<"Range = "<<range<<endmsg;
254 logstr << MYLEVEL <<"Sign = "<<sign<<endmsg;
255 logstr << MYLEVEL <<" Encoded Energy ="<< MSG::hex << encodedEnergy << MSG::dec << " E=" << energy
256 << " t=" << time
257 << " Q=" << quality
258 << " G=" << gain
259 << " channel Number=" << channelNumber
260 << endmsg;
261#endif
262
263 return 1;
264}
#define LARBSDBG(text)
This class provides decoding/encoding from/to ROD format.
int sign(int a)
const uint16_t * m_TimeQualityPointer
const uint16_t * m_EnergyPointer
time(flags, cells_name, *args, **kw)
setScaleOne setStatusOne setSaturated int16_t

◆ getNextRawData()

int LArRodBlockPhysicsV5::getNextRawData ( int & channelNumber,
std::vector< short > & samples,
uint32_t & gain )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 229 of file LArRodBlockPhysicsV5.cxx.

230{
231#ifdef LARBSDBGOUTPUT
232 MsgStream logstr(Athena::getMessageSvc(), BlockType());
233 //Debug output
234 logstr << MYLEVEL << "Let s go in getNextRawData..." << endmsg;
235 logstr << MYLEVEL << "GetNextRawData for FEB 0x" << MSG::hex << (uint32_t)getHeader32(FEBID) << MSG::dec << endmsg;
236 logstr << MYLEVEL << "m_RawDataPointer=" << m_RawDataPointer << " m_RawDataIndex="<< m_RawDataIndex
237 << " m_channelsPerFEB=" << m_channelsPerFEB << endmsg;
238#endif
239
240 if (m_RawDataIndex>=m_channelsPerFEB) { //Already beyond maximal number of channels
241#ifdef LARBSDBGOUTPUT
242 logstr << MYLEVEL << "Maximum number of channels reached" << endmsg;
243#endif
244 return 0;
245 }
246 if (!m_RawDataPointer) { //Block does not exist
247 // Try to get samples and gain from getNextDigits
248 return getNextDigits(channelNumber,samples,gain);
249 }
250
251 // Get next channel
252 unsigned rodChannelNumber=m_RawDataIndex; // Index of Channel in ROD-Block
253 channelNumber=((rodChannelNumber&0xe)<<2) + ((rodChannelNumber&0x1)<<6) + (rodChannelNumber>>4); //channel number of the FEB
254 uint32_t febgain;
255 const unsigned int nsamples = getHeader16(NSamples) & 0xff;
256 const unsigned int ngains = getHeader16(NGains);
257
258#ifdef LARBSDBGOUTPUT
259 logstr << MYLEVEL << "This FEB has " << nsamples << " samples" << endmsg;
260 logstr << MYLEVEL << "This FEB has " << ngains << " gains" << endmsg;
261#endif
262
263 if(ngains==0 || nsamples==0) return 0;
264 int s_size = nsamples+1;
265 int offset = (10+nsamples)&0xfffc;
266 int index;
267 index = s_size*m_RawDataIndex + offset;
268 uint16_t s[2];
269 if((nsamples+1)&0x7) {
270 s[0] = m_RawDataPointer[index++]>>2;
271 febgain = m_RawDataPointer[index++];
272 samples.push_back(s[0]);
273 for(unsigned int i=0;i<nsamples/2;i++) {
274 s[1] = m_RawDataPointer[index++]>>2;
275 s[0] = m_RawDataPointer[index++]>>2;
276 samples.push_back(s[0]);
277 samples.push_back(s[1]);
278 }
279 } // End of check for 5 samples
280 else {
281 if (!(m_RawDataIndex%2)) {
282 s[0] = m_RawDataPointer[index++]>>2;
283 febgain = m_RawDataPointer[index++];
284 samples.push_back(s[0]);
285 for(unsigned int i=0;i<nsamples/2;i++) {
286 s[1] = m_RawDataPointer[index++]>>2;
287 s[0] = m_RawDataPointer[index++]>>2;
288 samples.push_back(s[0]);
289 samples.push_back(s[1]);
290 }
291 } else {
292 for(unsigned int i=0;i<nsamples/2;i++) {
293 s[1] = m_RawDataPointer[index++]>>2;
294 s[0] = m_RawDataPointer[index++]>>2;
295 samples.push_back(s[0]);
296 samples.push_back(s[1]);
297 }
298 febgain = m_RawDataPointer[index++];
299 s[0] = m_RawDataPointer[index++]>>2;
300 samples.push_back(s[0]);
301 }
302 } // End of >5 check
303 gain=RawToOfflineGain(febgain);
304
305#ifdef LARBSDBGOUTPUT
306 logstr << MYLEVEL << " ===> ROD Channel = " << m_RawDataIndex << endmsg;
307 logstr << MYLEVEL << " ===> FEB Channel = " << channelNumber << endmsg;
308 logstr << MYLEVEL << " ===> Gain = " << gain << endmsg;
309 for(int i=0;i<nsamples;i++)
310 logstr << MYLEVEL << " ===> sample " << i << " = " << samples[i] << endmsg;
311 int n = m_RawDataIndex;
312 int32_t e,t,q;
313 uint32_t g;
315#endif
317 unsigned rearrangeFirstSample=0;
319 rearrangeFirstSample=m_rearrangeFirstSample; //Overwrite by jobOptions
320 else
321 rearrangeFirstSample=getFirstSampleIndex();
322 if (rearrangeFirstSample && rearrangeFirstSample<samples.size()) //FIXME: Very ugly hack! See explanation in LArRodDecoder.h file
323 {//Change e.g. 3 0 1 2 4 to 0 1 2 3 4
324 short movedSample=samples[0];
325 for (unsigned i=1;i<=rearrangeFirstSample;i++)
326 samples[i-1]=samples[i];
327 samples[rearrangeFirstSample]=movedSample;
328 }
329#ifdef LARBSDBGOUTPUT
330 logstr << MYLEVEL << "GetNextRawData for FEB finished 0x" << MSG::hex << (uint32_t)getHeader32(FEBID) << MSG::dec << endmsg;
331#endif
332 return 1;
333}
int getNextDigits(int &channelNumber, std::vector< short > &samples, uint32_t &gain)
virtual int getNextEnergy(int &channelNumber, int32_t &energy, int32_t &time, int32_t &quality, uint32_t &gain)

◆ getNStep()

uint16_t LArRodBlockStructure::getNStep ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockAccumulatedV3, and LArRodBlockCalibrationV3.

Definition at line 474 of file LArRodBlockStructure.h.

475{
476 std::cout << "Error: Function getNStep not implemented in this instance of LArRodBlockStructure!\n";
477 return 0;
478}

◆ getNTrigger()

uint16_t LArRodBlockStructure::getNTrigger ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockAccumulatedV3, LArRodBlockCalibrationV1, LArRodBlockCalibrationV2, and LArRodBlockCalibrationV3.

Definition at line 463 of file LArRodBlockStructure.h.

464{
465 std::cout << "Error: Function getNTrigger not implemented in this instance of LArRodBlockStructure!\n";
466 return 0;
467}

◆ getNumberOfGains()

uint32_t LArRodBlockPhysicsV5::getNumberOfGains ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 487 of file LArRodBlockPhysicsV5.cxx.

488{
489 return getHeader16(NGains);
490}

◆ getNumberOfSamples()

uint32_t LArRodBlockPhysicsV5::getNumberOfSamples ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 482 of file LArRodBlockPhysicsV5.cxx.

483{
484 return getHeader16(NSamples);
485}

◆ getNumberOfWords()

uint32_t LArRodBlockStructure::getNumberOfWords ( ) const
inlineinherited

Definition at line 411 of file LArRodBlockStructure.h.

412{
413 return getHeader32(NWTot);
414}

◆ getPulsed()

bool LArRodBlockStructure::getPulsed ( unsigned channelNumber) const
inlinevirtualinherited

Reimplemented in LArRodBlockCalibrationV0< DSPHEADER >, LArRodBlockCalibrationV1, LArRodBlockCalibrationV2, and LArRodBlockCalibrationV3.

Definition at line 445 of file LArRodBlockStructure.h.

446{
447 std::cout << "Error: Function getPulsed not implemented in this instance of LArRodBlockStructure!\n";
448 return 0;
449}

◆ getRadd()

uint32_t LArRodBlockPhysicsV5::getRadd ( uint32_t adc,
uint32_t sample ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 507 of file LArRodBlockPhysicsV5.cxx.

508{
509 if(!m_RawDataPointer) {
510 if(!m_RaddPointer) return 0;
511 if(sample%2) sample+=2;
512 return m_RaddPointer[sample];
513 }
514 int index;
515 if(sample==0) index=6;
516 else if(sample & 0x1) index=7+sample-1;
517 else index=7+sample+1;
519 if(adc>=8) return x>>8;
520 return x&0xff;
521}

◆ getRawDataSize()

uint16_t LArRodBlockPhysicsV5::getRawDataSize ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 502 of file LArRodBlockPhysicsV5.cxx.

503{
505}

◆ getResults1Size()

uint16_t LArRodBlockPhysicsV5::getResults1Size ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 492 of file LArRodBlockPhysicsV5.cxx.

493{
494 return getHeader16(ResultsDim1);
495}

◆ getResults2Size()

uint16_t LArRodBlockPhysicsV5::getResults2Size ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 497 of file LArRodBlockPhysicsV5.cxx.

498{
499 return getHeader16(ResultsDim2);
500}

◆ getStatus()

uint32_t LArRodBlockPhysicsV5::getStatus ( ) const
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 547 of file LArRodBlockPhysicsV5.cxx.

548{
549 if(getNumberOfWords()<EventStatus/2) return 0;
551 return x;
552}

◆ getStepIndex()

uint16_t LArRodBlockStructure::getStepIndex ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockAccumulatedV3, and LArRodBlockCalibrationV3.

Definition at line 469 of file LArRodBlockStructure.h.

470{
471 std::cout << "Error: Function getStepIndex not implemented in this instance of LArRodBlockStructure!\n";
472 return 0;
473}

◆ getSumE()

int32_t LArRodBlockStructure::getSumE ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockPhysicsV6.

Definition at line 310 of file LArRodBlockStructure.h.

310 {
311 return 0;
312}

◆ getTDCPhase()

uint8_t LArRodBlockStructure::getTDCPhase ( ) const
virtualinherited

Reimplemented in LArRodBlockAccumulatedV3, LArRodBlockCalibrationV1, LArRodBlockCalibrationV2, LArRodBlockCalibrationV3, and LArRodBlockTransparentV0< DSPHEADER >.

Definition at line 147 of file LArRodBlockStructure.cxx.

148{
149 std::cout << "Error: Function getTDCPhase not implemented in this instance of LArRodBlockStructure!\n";
150 return 0;
151}

◆ getVectorHeader16()

uint16_t LArRodBlockStructure::getVectorHeader16 ( const unsigned n) const
inlineprotectedinherited

Definition at line 351 of file LArRodBlockStructure.h.

352{ if (n&0x1) //n is a odd number
353 return (std::as_const(*m_vFragment).at(n>>1) & 0xffff);
354 else //n is a even number
355 return (std::as_const(*m_vFragment).at(n>>1) >> 16);
356}

◆ getVectorHeader32()

uint32_t LArRodBlockStructure::getVectorHeader32 ( const unsigned n) const
inlineprotectedinherited

Definition at line 358 of file LArRodBlockStructure.h.

359{
360 return (*m_vFragment)[n>>1];
361}

◆ getVROBEx()

int32_t LArRodBlockPhysicsV5::getVROBEx ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 295 of file LArRodBlockPhysicsV5.h.

296{
297 const int32_t* p = reinterpret_cast<const int32_t*>(m_virtualROBPointerLocal);
298 if(p) return (p[1]>>9);
299 return 0;
300}
const uint32_t * m_virtualROBPointerLocal

◆ getVROBEy()

int32_t LArRodBlockPhysicsV5::getVROBEy ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 302 of file LArRodBlockPhysicsV5.h.

303{
304 const int32_t* p = reinterpret_cast<const int32_t*>(m_virtualROBPointerLocal);
305 if(p) return (p[2]>>9);
306 return 0;
307}

◆ getVROBEz()

int32_t LArRodBlockPhysicsV5::getVROBEz ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 309 of file LArRodBlockPhysicsV5.h.

310{
311 const int32_t* p = reinterpret_cast<const int32_t*>(m_virtualROBPointerLocal);
312 if(p) return (p[3]>>9);
313 return 0;
314}

◆ getVROBFebId()

uint32_t LArRodBlockPhysicsV5::getVROBFebId ( )
inlinevirtual

◆ getVROBSumE()

int32_t LArRodBlockStructure::getVROBSumE ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockPhysicsV6.

Definition at line 330 of file LArRodBlockStructure.h.

330 {
331return 0;
332}

◆ hasAccumBlock()

virtual uint32_t LArRodBlockStructure::hasAccumBlock ( ) const
inlinevirtualinherited

Reimplemented in LArRodBlockAccumulatedV3.

Definition at line 121 of file LArRodBlockStructure.h.

121{return 0;} ;

◆ hasCalibBlock()

virtual uint32_t LArRodBlockPhysicsV5::hasCalibBlock ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 83 of file LArRodBlockPhysicsV5.h.

83{return 0;} ;

◆ hasControlWords()

virtual uint32_t LArRodBlockPhysicsV5::hasControlWords ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 86 of file LArRodBlockPhysicsV5.h.

86{return getHeader16(RawDataBlkOff);} ;

◆ hasPhysicsBlock()

virtual uint32_t LArRodBlockPhysicsV5::hasPhysicsBlock ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 84 of file LArRodBlockPhysicsV5.h.

84{return getHeader16(ResultsOff1);} ;

◆ hasRawDataBlock()

virtual uint32_t LArRodBlockPhysicsV5::hasRawDataBlock ( ) const
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 85 of file LArRodBlockPhysicsV5.h.

◆ initializeFEB()

void LArRodBlockPhysicsV5::initializeFEB ( const uint32_t id)
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 580 of file LArRodBlockPhysicsV5.cxx.

581{
583 if (m_vFragment->size()<m_iHeadBlockSize) //Got empty or spoiled fragment
584 {
585 m_vFragment->resize(m_iHeadBlockSize,0); //Initialize FEB-Header
586 setHeader32(FEBID,id); //Set Feb ID
587 // At least 10 (head) + 16 (gain/sumblks) + 64 (energies)
588 m_vFragment->reserve(90);
589 }
590
591 m_SumBlkBlockE1.resize(4);
592 for(unsigned int i=0;i<4;i++) m_SumBlkBlockE1[i]=0x0;
593 m_SumBlkBlockE2.resize(4);
594 for(unsigned int i=0;i<4;i++) m_SumBlkBlockE2[i]=0x0;
595 m_GainBlock.resize(8);
596 for(unsigned int i=0;i<8;i++) m_GainBlock[i]=0x0;
597 m_TimeQualityBlock.resize(6);
598 for(unsigned int i=0;i<6;i++) m_TimeQualityBlock[i]=0x0;
599 m_EnergyBlockEncode.resize(128);
600 for(unsigned int i=0;i<128;i++) m_EnergyBlockEncode[i]=0x0;
601 m_DigitsEncode.clear();
602
604
605}

◆ initializeFragment()

void LArRodBlockPhysicsV5::initializeFragment ( std::vector< uint32_t > & fragment)
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 555 of file LArRodBlockPhysicsV5.cxx.

555 {
556 m_pRODblock=&fragment; //remember pointer to fragment
557 if (fragment.size()>m_iHeadBlockSize) { //Got filled fragment
558 unsigned int sizeRead=0;
559 //Store existing data in the FEB-Map
560 while (sizeRead<fragment.size()) {
561 std::vector<uint32_t>::iterator FebIter;
562 FebIter=fragment.begin()+sizeRead; //Store pointer to current Feb-Header
563 m_FebBlock=&(*FebIter); //Set m_FebBlock in order to use getHeader-functions.
564 uint32_t currFEBid=getHeader32(FEBID); //Get this FEB-ID
565 uint16_t currFebSize=getNumberOfWords(); //Size of this FEB-Block
566 if (FebIter+currFebSize>fragment.end()) {
567 fragment.clear(); //Clear existing vector
568 return;
569 }
570 m_mFebBlocks[currFEBid].assign(FebIter,FebIter+currFebSize); //Copy data from ROD-fragment into FEB-Block
571 sizeRead+=currFebSize+m_MiddleHeaderSize; //6 is the middle header size
572 } // end while
573 }
574 fragment.clear(); //Clear existing vector
575 return;
576
577}

◆ LE_getHeader16()

uint16_t LArRodBlockStructure::LE_getHeader16 ( const unsigned n) const
inlineprotectedinherited

Definition at line 388 of file LArRodBlockStructure.h.

389{
390 return (reinterpret_cast<const uint16_t*>(m_FebBlock))[n];
391}

◆ LE_getVectorHeader16()

uint16_t LArRodBlockStructure::LE_getVectorHeader16 ( const unsigned n) const
inlineprotectedinherited

Definition at line 393 of file LArRodBlockStructure.h.

394{
395 const uint32_t* data32 = std::as_const(*m_vFragment).data();
396 const uint16_t* data16 = reinterpret_cast<const uint16_t*> (data32);
397 return data16[n];
398}

◆ LE_setHeader16()

void LArRodBlockStructure::LE_setHeader16 ( const unsigned n,
const uint16_t w )
inlineprotectedinherited

Definition at line 400 of file LArRodBlockStructure.h.

401{
402#ifdef LARBYTESTREAMRODBLOCK_CHCKBOUNDARIES
403 if ((unsigned)n>=m_vFragment->size()*2) {
404 std::cout << "LArRodBlockStructure::LE_setHeader16 Error: WRITE BEYOND ARRAY BONDARY!" << std::endl;
405 std::abort();
406 }
407#endif
408 reinterpret_cast<uint16_t*>(m_vFragment->data())[n] = w;
409}

◆ nextFEB()

bool LArRodBlockStructure::nextFEB ( )
inlineinherited

Definition at line 480 of file LArRodBlockStructure.h.

481{
483
484 const int32_t FebOffset = m_FebBlockSize+m_MiddleHeaderSize;
485 const int32_t LeftSize = m_RodBlockSize-FebOffset+m_RodBlock-m_FebBlock;
486 if ( LeftSize<=0 ) return false; //No next feb
487
488 m_FebBlock += FebOffset; //Jump to the next FEB fragment
489 const int32_t BlockSize = getNumberOfWords();
490 if (BlockSize>LeftSize) {
491 std::cout << "Error while decoding LArByteStream: Found FEB block of size " << BlockSize << " in a ROD block of size " << LeftSize << std::endl;
492 m_error_next_feb = true;
493 return false;
494 }
495
496 m_FebBlockSize=BlockSize;
497
498 setPointers();
499 return true;
500}

◆ offlineCheckSum()

uint32_t LArRodBlockStructure::offlineCheckSum ( ) const
virtualinherited

Reimplemented in LArRodBlockPhysicsV4.

Definition at line 161 of file LArRodBlockStructure.cxx.

162{
163 int end = getNumberOfWords()-3;
164 int start = 1;
165 uint32_t sum = 0;
167 for(int i=start;i<end;i++) {
168 sum += m_FebBlock[i];
169 }
170 return sum & 0x7fffffff;
171}

◆ OfflineToRawGain()

uint32_t LArRodBlockStructure::OfflineToRawGain ( const uint32_t gain) const
inlineinherited

Definition at line 337 of file LArRodBlockStructure.h.

338{return m_OfflineToRawGainMap[gain];} //For efficency, don't check range
static const uint32_t m_OfflineToRawGainMap[3]

◆ onlineCheckSum()

uint32_t LArRodBlockStructure::onlineCheckSum ( ) const
virtualinherited

Reimplemented in LArRodBlockPhysicsV4.

Definition at line 153 of file LArRodBlockStructure.cxx.

154{
155 int offset = getNumberOfWords()-1;
156 if(offset>=m_FebBlockSize) return 0;
157 if(offset<0) return 0;
158 return m_FebBlock[offset];
159}

◆ operator()()

template<class RAWDATA>
bool LArRodBlockPhysicsV5::operator() ( const RAWDATA * ch1,
const RAWDATA * ch2 ) const

Definition at line 875 of file LArRodBlockPhysicsV5.cxx.

877{
878 HWIdentifier id1 = ch1->channelID();
879 HWIdentifier id2 = ch2->channelID();
880
881 HWIdentifier febId1= m_onlineHelper->feb_Id(id1);
882 HWIdentifier febId2= m_onlineHelper->feb_Id(id2);
883
884 if(febId1 == febId2 ){
885 int cId1 = m_onlineHelper->channel(id1);
886 int cId2 = m_onlineHelper->channel(id2);
887 return FebToRodChannel(cId1) < FebToRodChannel(cId2);
888 }
889
890 return febId1 < febId2 ;
891}
virtual int FebToRodChannel(int ch) const

◆ RawToOfflineGain()

uint32_t LArRodBlockStructure::RawToOfflineGain ( const uint32_t gain) const
inlineinherited

Definition at line 334 of file LArRodBlockStructure.h.

335{return m_RawToOfflineGainMap[gain];} //For efficency, don't check range
static const uint32_t m_RawToOfflineGainMap[4]

◆ report_error()

bool LArRodBlockStructure::report_error ( void ) const
inlineinherited

Definition at line 173 of file LArRodBlockStructure.h.

173{ return m_error_next_feb;}

◆ resetPointers()

void LArRodBlockPhysicsV5::resetPointers ( )
privatevirtual

◆ setBit()

void LArRodBlockStructure::setBit ( uint32_t *const p,
const unsigned chan )
inlineprotectedinherited

Definition at line 427 of file LArRodBlockStructure.h.

428{
429 // chan = (0,127)
430 // int a = chan/32;
431 // int r = chan%32;
432 int a = chan>>5;
433 int r = chan&0x1f;
434 // a = (0,3), r = ( 0, 31 )
435 *(p+a) |= (1<<r) ;
436 return;
437}

◆ setDAC()

void LArRodBlockStructure::setDAC ( const uint16_t DACValue)
virtualinherited

Reimplemented in LArRodBlockCalibrationV0< DSPHEADER >.

Definition at line 55 of file LArRodBlockStructure.cxx.

56{
57 std::cout << "Error: Function setDAC not implemented in this instance of LArRodBlockStructure!\n";
58 return;
59}

◆ setDelay()

void LArRodBlockStructure::setDelay ( const uint16_t DelayValue)
virtualinherited

Reimplemented in LArRodBlockCalibrationV0< DSPHEADER >.

Definition at line 61 of file LArRodBlockStructure.cxx.

62{
63 std::cout << "Error: Function setDelay not implemented in this instance of LArRodBlockStructure!\n";
64 return;
65}

◆ setE()

void LArRodBlockPhysicsV5::setE ( unsigned index,
double E )
inlineprivate

Definition at line 322 of file LArRodBlockPhysicsV5.h.

322 {
323
324 union {
325 int32_t i;
326 uint16_t us[2];
327 } conv;
328 conv.i = static_cast<int32_t>(E);
329 // Write as multiplication, not as left shift, since left-shifting
330 // a negative number is undefined in C++.
331 // Compiles to the same code on x86_64.
332 conv.i *= (1<<9);
333 if ( m_TimeQualityBlock.size()>=6 ){
336 }
337}

◆ setEtQ()

void LArRodBlockStructure::setEtQ ( const int channel,
const int32_t energy,
const int32_t time,
const int32_t quality,
const uint32_t gain )
virtualinherited

Definition at line 118 of file LArRodBlockStructure.cxx.

119{
120 std::cout << "Error: Function setEtQ not implemented in this instance of LArRodBlockStructure!\n";
121 return;
122}

◆ setEx()

void LArRodBlockPhysicsV5::setEx ( double Ex)
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 340 of file LArRodBlockPhysicsV5.h.

340 {
341 setE (0, Ex);
342}
void setE(unsigned index, double E)

◆ setEy()

void LArRodBlockPhysicsV5::setEy ( double Ey)
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 345 of file LArRodBlockPhysicsV5.h.

345 {
346 setE (2, Ey);
347}

◆ setEz()

void LArRodBlockPhysicsV5::setEz ( double Ez)
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 350 of file LArRodBlockPhysicsV5.h.

350 {
351 setE (4, Ez);
352}

◆ setFirstSample()

void LArRodBlockStructure::setFirstSample ( const int rearrangeFirstSample)
inlineinherited

Definition at line 166 of file LArRodBlockStructure.h.

167 {m_rearrangeFirstSample=rearrangeFirstSample;}

◆ setFragment()

bool LArRodBlockStructure::setFragment ( const uint32_t * p,
uint32_t n )
inlineinherited

Definition at line 248 of file LArRodBlockStructure.h.

249{
251 m_FebBlock = p;
254 uint32_t BlockSize = 0;
255
256 if (n==0) {
257 std::cout << "Error while decoding LArByteStream: Got Rod block size 0" << std::endl;
258 return false;
259 }
260
261 BlockSize = getNumberOfWords();
262 if (BlockSize>n) {
263 std::cout << "Error while decoding LArByteStream: Found FEB block of size " << BlockSize << " in a ROD block of size " << n << std::endl;
264 return false;
265 }
266 m_FebBlockSize = BlockSize;
267
268 m_error_next_feb = false;
269
270 return setPointers();
271}

◆ setFragmentVirtualROB()

int LArRodBlockStructure::setFragmentVirtualROB ( const uint32_t * p,
uint32_t n )
inlinevirtualinherited

Definition at line 275 of file LArRodBlockStructure.h.

276 {
281 if ( n<2 ) { std::cout << "Error" << std::endl; return 0;}
282 m_virtualROBJump = ((*p)>>16)>>1; // Divide by two (two FEBs-1ROB)
283 m_ROB_to_decode = ( (*p) & 0xFFFF )<<1; // Multiply by two
284 if ( (n - m_virtualROBJump*m_ROB_to_decode-1) ){
285 std::cout << "Error AGAIN" << std::endl;
286 return 0;
287 }
289 return m_ROB_to_decode;
290}

◆ setGain()

virtual int LArRodBlockPhysicsV5::setGain ( const int )
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 98 of file LArRodBlockPhysicsV5.h.

98{ return 1; };

◆ setHeader16()

void LArRodBlockStructure::setHeader16 ( const unsigned n,
const uint16_t w )
inlineprotectedinherited

Definition at line 363 of file LArRodBlockStructure.h.

364{
365#ifdef LARBYTESTREAMRODBLOCK_CHCKBOUNDARIES
366 if ((unsigned)n>=m_vFragment->size()*2) {
367 std::cout << "Error WRITE BEYOND ARRAY BONDARY!" << std::endl;
368 std::abort();
369 }
370#endif
371 if (n&0x1) // n is a odd number
372 m_vFragment->at(n>>1)=((m_vFragment->at(n>>1) & 0xffff0000) | w);
373 else
374 m_vFragment->at(n>>1)=((m_vFragment->at(n>>1) & 0xffff) | (w << 16));
375}

◆ setHeader32()

void LArRodBlockStructure::setHeader32 ( const unsigned n,
const uint32_t w )
inlineprotectedinherited

Definition at line 377 of file LArRodBlockStructure.h.

378{
379#ifdef LARBYTESTREAMRODBLOCK_CHCKBOUNDARIES
380 if ((unsigned)n>=m_vFragment->size()*2) {
381 std::cout << "Error WRITE BEYOND ARRAY BONDARY!" << std::endl;
382 std::abort();
383 }
384#endif
385 m_vFragment->at(n>>1) = w;
386}

◆ setNextEnergy() [1/2]

void LArRodBlockPhysicsV5::setNextEnergy ( const int channel,
const int32_t energy,
const int32_t time,
const int32_t quality,
const uint32_t gain )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 607 of file LArRodBlockPhysicsV5.cxx.

609{
610 int rcNb=FebToRodChannel(channel);
611 //rcNb is supposed to be equal or bigger than m_EIndex.
612 //In the latter case, we fill up the missing channels with zero
613 if (rcNb<m_EnergyIndex) {
614 return;
615 }
616
617 //Fill up missing channels with zeros:
618 while (m_EnergyIndex<rcNb)
619 setNextEnergy((int16_t)0,(int16_t)32767,(int16_t)-32767,(uint32_t)0);
620
621 // transform 32 bits data into 16 bits data
622
623 uint16_t theenergy;
624 uint32_t abse,EncodedE;
625 int16_t thetime,thequality;
626 int32_t sign;
627
628 //Time is in 10 ps in ByteStream, hence the factor 10 to convert from ps
629 thetime = (int16_t) time/10;
630 thequality = (int16_t) quality;
631
632 sign=(energy>=0?1:-1); // get sign of energy
633 abse=(uint32_t)abs(energy);
634
635 EncodedE=abse; // range 0
636
637 if ((abse>8192)&&(abse<65536))
638 {
639 EncodedE=((abse>>3)|0x4000); // range 1 : drop last 3 bits and put range bits (bits 14 and 13 = 01)
640 }
641 else if ((abse>65535)&&(abse<524288))
642 {
643 EncodedE=((abse>>6)|0x8000); // range 2 : drop last 6 bits and put range bits (bits 14 and 13 = 10)
644 }
645 else if ((abse>524288))
646 {
647 EncodedE=((abse>>9)|0xc000); // range 3 : drop last 9 bits and put range bits (bits 14 and 13 = 11)
648 }
649
650 // treat sign now :
651
652 if (sign<0) EncodedE |= 0x2000;
653 theenergy = (uint16_t) EncodedE;
654
655
656 // Add data...
657
658 if (abse> m_EnergyThreshold1)
659 {
660 setNextEnergy(theenergy,thetime,thequality,gain);
661 }
662 else
663 {
664 setNextEnergy(theenergy,(int16_t)32767,(int16_t)-32767,gain);
665 }
666 return;
667}

◆ setNextEnergy() [2/2]

void LArRodBlockPhysicsV5::setNextEnergy ( const uint16_t energy,
const int16_t time,
const int16_t quality,
const uint32_t gain )
private

Definition at line 670 of file LArRodBlockPhysicsV5.cxx.

671{
672 if (m_EnergyIndex>=m_channelsPerFEB) //Use m_EIndex to count total number of channels
673 {
674 return;
675 }
676
677 // Energy
678 int endianindex;
679 if (m_EnergyIndex & 0x1) endianindex = m_EnergyIndex-1;
680 else endianindex = m_EnergyIndex+1;
681 m_EnergyBlockEncode[endianindex] = energy;
682
683 // Find correct position
684
685 // update summary block
686 // Gain is composed of two bits per cell
687 uint16_t gain_idx=m_EnergyIndex>>4;
688 uint16_t gain_bit=(m_EnergyIndex&0xf)*2;
690 m_GainBlock[gain_idx] |= (gain1 << gain_bit);
691
692 // write Time and Chi2 for cells above HighEnergyCellCut threshold
693
694 if (quality!=-32767) // Do write Time and Chi2 information
695 {
696 // count the number of hot cells
698 // count the number of cells offtime
699 if (abs(time)>m_OffTimeCut) m_numberHotCellOffTime++;
700 uint16_t mask_idx=m_EnergyIndex>>5;
701 uint16_t mask_bit=(m_EnergyIndex&0x1f);
702 m_SumBlkBlockE1[mask_idx] |= (0x1 << mask_bit);
703
704 m_TimeQualityBlock.push_back(*((uint16_t*)&time));
705 m_TimeQualityBlock.push_back(*((uint16_t*)&quality));
706 }
707 m_EnergyIndex++; //Use m_EIndex to count the channels put in the Energy block
708
709}
uint32_t OfflineToRawGain(const uint32_t gain) const

◆ setNTrigger()

void LArRodBlockStructure::setNTrigger ( const uint16_t NTrigger)
virtualinherited

Definition at line 73 of file LArRodBlockStructure.cxx.

74{
75 std::cout << "Error: Function setNTrigger not implemented in this instance of LArRodBlockStructure!\n";
76 return;
77}

◆ setNumberOfGains()

void LArRodBlockStructure::setNumberOfGains ( const uint8_t n)
virtualinherited

Reimplemented in LArRodBlockPhysicsV0, LArRodBlockPhysicsV1, LArRodBlockPhysicsV2, LArRodBlockPhysicsV3, and LArRodBlockTransparentV0< DSPHEADER >.

Definition at line 136 of file LArRodBlockStructure.cxx.

137{
138 std::cout << "Error: Function setNumberOfGains not implemented in this instance of LArRodBlockStructure!\n";
139 return;
140}

◆ setNumberOfSamples()

void LArRodBlockStructure::setNumberOfSamples ( const uint8_t n)
virtualinherited

Reimplemented in LArRodBlockPhysicsV0, LArRodBlockPhysicsV1, LArRodBlockPhysicsV2, LArRodBlockPhysicsV3, and LArRodBlockTransparentV0< DSPHEADER >.

Definition at line 130 of file LArRodBlockStructure.cxx.

131{
132 std::cout << "Error: Function setNumberOfSamples not implemented in this instance of LArRodBlockStructure!\n";
133 return;
134}

◆ setPointers()

bool LArRodBlockPhysicsV5::setPointers ( )
privatevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 80 of file LArRodBlockPhysicsV5.cxx.

81{
83 {
84 int off = -8;
85 int ns = getHeader16(NSamples) & 0xff;
86 if (m_requiredNSamples > 0 && m_requiredNSamples != ns) return false;
87 int radd = (ns+1)/2;
88 int dim1 = getHeader16(ResultsDim1);
89 int off1 = getHeader16(ResultsOff1);
90 int off2 = getHeader16(ResultsOff2);
91 int dim2 = getHeader16(ResultsDim2);
92 int off3 = getHeader16(RawDataBlkOff);
93 int dim3 = getHeader16(RawDataBlkDim);
94
95 if (off1 && dim1+off1+off<m_FebBlockSize) {
96 off1 += off;
97 if (dim1>=8)
98 m_GainPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+off1);
99 if (dim1>=12)
100 m_MaskTimeQualityPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+off1+8);
101 if (dim1>=16)
102 m_MaskDigitsPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+off1+12);
103 if (dim1>=16+radd)
104 m_RaddPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off1+16);
105 if (dim1>=80+radd)
106 m_EnergyPointer=reinterpret_cast<const uint16_t*> (m_FebBlock+off1+16+radd);
107 if (dim1>=83+radd)
108 m_SumPointer=reinterpret_cast<const int32_t*>(m_FebBlock+off1+80+radd);
109 if (dim1>83+radd)
110 m_TimeQualityPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off1+83+radd);
111 off1 -= off;
112 }
113 if (off2 && dim2+off2+off<m_FebBlockSize) {
114 m_DigitsPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off2+off);
115 }
116 if (off3 && dim3+off3+off<m_FebBlockSize) {
117 m_RawDataPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off3+off);
118 }
119
120 // Check for offsets problems
121 uint32_t problem = 0;
122 int n1, n2;
123 int n1_tmp, n2_tmp;
124 int off1_tmp, dim1_tmp;
125 int off2_tmp, dim2_tmp;
126 int off3_tmp, dim3_tmp;
127 if(off1==0) {
128 n1 = n2 = 0;
129 n1_tmp = n2_tmp =0;
130 off1_tmp = dim1_tmp = 0;
131 off2_tmp = dim2_tmp = 0;
132 off3_tmp = dim3_tmp = 0;
133 }
134 else {
135 m_RaddPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+26);
136 m_MaskTimeQualityPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+18);
137 m_MaskDigitsPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+22);
138 n1 = getNbSweetCells1();
139 n2 = getNbSweetCells2();
140 n1_tmp = getNbSweetCells1FromMask();
141 n2_tmp = getNbSweetCells2FromMask();
142 off1_tmp = 10-off;
143 dim1_tmp = 83+(ns+1)/2+n1;
144 if ( m_requiredNSamples > 0 )
145 dim1_tmp = 83 +(m_requiredNSamples+1)/2+n1;
146 off2_tmp = off1_tmp+dim1_tmp;
147 dim2_tmp = (n2*ns+1)/2;
148 off3_tmp = off2_tmp+dim2_tmp;
149 dim3_tmp = getNumberOfWords()-3-off3_tmp-off;
150 if(dim2_tmp==0) off2_tmp = 0;
151 if(dim3_tmp==0) off3_tmp = 0;
152 }
153
154 if(off1 != off1_tmp) problem=1;
155 if(dim1 != dim1_tmp) problem=2;
156 if(off2 != off2_tmp) problem=3;
157 if(dim2 != dim2_tmp) problem=4;
158 if(off3 != off3_tmp) problem=5;
159 if(dim3 != dim3_tmp) problem=6;
160 if(n1 != n1_tmp) problem=7;
161 if(n2 != n2_tmp) problem=8;
162 if (m_requiredNSamples > 0 &&
163 getHeader32(NGains) != (uint32_t)0x10000 + m_requiredNSamples) problem=9;
164 if(problem) { // Try to recompute offsets
165 std::cout << "LArByteStreamProblem " << problem << std::endl;
166 std::cout << "NSamples = " << std::dec << ns << std::endl;
167 std::cout << "getHeader32(NGains) = " << std::hex << getHeader32(NGains) << std::endl;
168 std::cout << "NWTot: " << std::hex << getNumberOfWords() << " n1=" << n1 << " (" << n1_tmp << ") n2=" << n2 << " (" << n2_tmp << ")" << std::endl;
169 std::cout << "Found 1: " << off1 << " " << dim1 << std::endl;
170 std::cout << "Found 2: " << off2 << " " << dim2 << std::endl;
171 std::cout << "Found 3: " << off3 << " " << dim3 << std::dec << std::endl;
172
173 if(n1==n1_tmp && n2==n2_tmp) { // Check consistency of cells above threshold
174 off1 = off1_tmp;
175 dim1 = dim1_tmp;
176 off2 = off2_tmp;
177 dim2 = dim2_tmp;
178 off3 = off3_tmp;
179 dim3 = dim3_tmp;
180 std::cout << "Recomputed 1: " << std::hex << off1 << " " << dim1 << std::endl;
181 std::cout << "Recomputed 2: " << off2 << " " << dim2 << std::endl;
182 std::cout << "Recomputed 3: " << off3 << " " << dim3 << std::dec << std::endl;
183
184 if (off1 && dim1+off1+off<m_FebBlockSize) {
185 off1 += off;
186 if (dim1>=8)
187 m_GainPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+off1);
188 if (dim1>=12)
189 m_MaskTimeQualityPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+off1+8);
190 if (dim1>=16)
191 m_MaskDigitsPointer=reinterpret_cast<const uint32_t*>(m_FebBlock+off1+12);
192 if (dim1>=16+radd)
193 m_RaddPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off1+16);
194 if (dim1>=80+radd)
195 m_EnergyPointer=reinterpret_cast<const uint16_t*> (m_FebBlock+off1+16+radd);
196 if (dim1>=83+radd)
197 m_SumPointer=reinterpret_cast<const int32_t*>(m_FebBlock+off1+80+radd);
198 if (dim1>83+radd)
199 m_TimeQualityPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off1+83+radd);
200 }
201 if (off2 && dim2+off2+off<m_FebBlockSize) {
202 m_DigitsPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off2+off);
203 }
204 if (off3 && dim3+off3+off<m_FebBlockSize) {
205 m_RawDataPointer=reinterpret_cast<const uint16_t*>(m_FebBlock+off3+off);
206 }
207 }
208 }
209
210 problem=0;
211 // Recheck offsets
212 if(off1< off2 && off1 + dim1 > off2) problem = 1;
213 if(off1< off3 && off1 + dim1 > off3) problem = 2;
214 if(off2< off1 && off2 + dim2 > off1) problem = 3;
215 if(off2< off3 && off2 + dim2 > off3) problem = 4;
216 if(off3< off1 && off3 + dim3 > off1) problem = 5;
217 if(off3< off2 && off3 + dim3 > off2) problem = 6;
218
219 if(problem) {
221 std::cout << "LArByteStreamProblem " << problem << std::endl;
222 std::cout << "Unrecoverable problem" << std::endl;
223 }
224 }
225
226 return true;
227}
virtual uint16_t getNbSweetCells2() const
uint16_t getNbSweetCells1FromMask() const
uint16_t getNbSweetCells2FromMask() const
virtual uint16_t getNbSweetCells1() const

◆ setPulsed()

void LArRodBlockStructure::setPulsed ( const unsigned channelNumber)
virtualinherited

Reimplemented in LArRodBlockCalibrationV0< DSPHEADER >.

Definition at line 67 of file LArRodBlockStructure.cxx.

68{
69 std::cout << "Error: Function setPulsed not implemented in this instance of LArRodBlockStructure!\n";
70 return;
71}

◆ setRawData()

void LArRodBlockPhysicsV5::setRawData ( const int chIdx,
const std::vector< short > & samples,
const uint32_t  )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 711 of file LArRodBlockPhysicsV5.cxx.

711 {
712
713 // First of all, set the bits
714 int cchIdx = FebToRodChannel(chIdx);
715 uint16_t mask_idx=cchIdx>>5;
716 uint16_t mask_bit=(cchIdx&0x1f);
717 m_SumBlkBlockE2[mask_idx] |= (0x1 << mask_bit);
718 for(std::vector<short>::const_iterator i=samples.begin();i!=samples.end();++i){
719 m_DigitsEncode.push_back((*i)<<2);
720 }
721
722}

◆ setRawDataFixed()

void LArRodBlockStructure::setRawDataFixed ( const int channel,
const std::vector< short > & samples,
const uint32_t gain )
virtualinherited

Reimplemented in LArRodBlockTransparentV0< DSPHEADER >.

Definition at line 124 of file LArRodBlockStructure.cxx.

125{
126 std::cout << "Error: Function setRawDataFixed not implemented in this instance of LArRodBlockStructure!\n";
127 return;
128}

◆ setRequiredNSamples()

void LArRodBlockPhysicsV5::setRequiredNSamples ( unsigned short ns)
inline

Definition at line 112 of file LArRodBlockPhysicsV5.h.

◆ setSumE()

void LArRodBlockStructure::setSumE ( double )
virtualinherited

Reimplemented in LArRodBlockPhysicsV6.

Definition at line 216 of file LArRodBlockStructure.cxx.

217{
218 // Do not set anything
219 return;
220}

◆ setTDCPhase()

void LArRodBlockStructure::setTDCPhase ( const uint8_t n)
virtualinherited

Reimplemented in LArRodBlockTransparentV0< DSPHEADER >.

Definition at line 142 of file LArRodBlockStructure.cxx.

143{
144 std::cout << "Error: Function setTDCPhase not implemented in this instance of LArRodBlockStructure!\n";
145 return;
146}

◆ sortDataVector() [1/5]

void LArRodBlockStructure::sortDataVector ( std::vector< const LArAccumulatedCalibDigit * > & )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 90 of file LArRodBlockStructure.cxx.

238{return;}

◆ sortDataVector() [2/5]

void LArRodBlockStructure::sortDataVector ( std::vector< const LArAccumulatedDigit * > & )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 91 of file LArRodBlockStructure.cxx.

241{return;}

◆ sortDataVector() [3/5]

void LArRodBlockStructure::sortDataVector ( std::vector< const LArCalibDigit * > & )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 89 of file LArRodBlockStructure.cxx.

235{return;}

◆ sortDataVector() [4/5]

void LArRodBlockPhysicsV5::sortDataVector ( std::vector< const LArDigit * > & vDigit)
inlinevirtual

Reimplemented from LArRodBlockStructure.

Definition at line 118 of file LArRodBlockPhysicsV5.h.

119 {std::sort(vDigit.begin(),vDigit.end(),*this); }
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.

◆ sortDataVector() [5/5]

void LArRodBlockStructure::sortDataVector ( std::vector< const LArRawChannel * > & )
virtual

Reimplemented from LArRodBlockStructure.

Definition at line 87 of file LArRodBlockStructure.cxx.

229{return;}

Member Data Documentation

◆ m_channelsPerFEB

int LArRodBlockStructure::m_channelsPerFEB
protectedinherited

Definition at line 219 of file LArRodBlockStructure.h.

◆ m_DigitsChannel

int LArRodBlockPhysicsV5::m_DigitsChannel = 0
private

Definition at line 142 of file LArRodBlockPhysicsV5.h.

◆ m_DigitsEncode

std::vector<uint16_t> LArRodBlockPhysicsV5::m_DigitsEncode
private

Definition at line 137 of file LArRodBlockPhysicsV5.h.

◆ m_DigitsIndex

int LArRodBlockPhysicsV5::m_DigitsIndex = 0
private

Definition at line 141 of file LArRodBlockPhysicsV5.h.

◆ m_DigitsPointer

const uint16_t* LArRodBlockPhysicsV5::m_DigitsPointer = nullptr
private

Definition at line 151 of file LArRodBlockPhysicsV5.h.

◆ m_EnergyBlockEncode

std::vector<uint16_t> LArRodBlockPhysicsV5::m_EnergyBlockEncode
private

Definition at line 136 of file LArRodBlockPhysicsV5.h.

◆ m_EnergyIndex

int LArRodBlockPhysicsV5::m_EnergyIndex = 0
private

Definition at line 139 of file LArRodBlockPhysicsV5.h.

◆ m_EnergyPointer

const uint16_t* LArRodBlockPhysicsV5::m_EnergyPointer = nullptr
private

Definition at line 148 of file LArRodBlockPhysicsV5.h.

◆ m_EnergyThreshold1

uint16_t LArRodBlockPhysicsV5::m_EnergyThreshold1
private

Definition at line 159 of file LArRodBlockPhysicsV5.h.

◆ m_EnergyThreshold2

uint16_t LArRodBlockPhysicsV5::m_EnergyThreshold2
private

Definition at line 160 of file LArRodBlockPhysicsV5.h.

◆ m_error_next_feb

bool LArRodBlockStructure::m_error_next_feb
protectedinherited

Definition at line 237 of file LArRodBlockStructure.h.

◆ m_Ex

int32_t LArRodBlockStructure::m_Ex
protectedinherited

Definition at line 203 of file LArRodBlockStructure.h.

◆ m_Ey

int32_t LArRodBlockStructure::m_Ey
protectedinherited

Definition at line 204 of file LArRodBlockStructure.h.

◆ m_Ez

int32_t LArRodBlockStructure::m_Ez
protectedinherited

Definition at line 205 of file LArRodBlockStructure.h.

◆ m_FebBlock

const uint32_t* LArRodBlockStructure::m_FebBlock
protectedinherited

Definition at line 221 of file LArRodBlockStructure.h.

◆ m_FebBlockSize

int32_t LArRodBlockStructure::m_FebBlockSize
protectedinherited

Definition at line 223 of file LArRodBlockStructure.h.

◆ m_FebInfoBlock

std::vector<uint32_t> LArRodBlockPhysicsV5::m_FebInfoBlock
private

Definition at line 133 of file LArRodBlockPhysicsV5.h.

◆ m_fixedGain

int LArRodBlockPhysicsV5::m_fixedGain
private

Definition at line 155 of file LArRodBlockPhysicsV5.h.

◆ m_GainBlock

std::vector<uint32_t> LArRodBlockPhysicsV5::m_GainBlock
private

Definition at line 132 of file LArRodBlockPhysicsV5.h.

◆ m_GainPointer

const uint32_t* LArRodBlockPhysicsV5::m_GainPointer = nullptr
private

Definition at line 144 of file LArRodBlockPhysicsV5.h.

◆ m_iHeadBlockSize

unsigned short LArRodBlockStructure::m_iHeadBlockSize
protectedinherited

Definition at line 217 of file LArRodBlockStructure.h.

◆ m_logstr

MsgStream LArRodBlockStructure::m_logstr
protectedinherited

Definition at line 244 of file LArRodBlockStructure.h.

◆ m_MaskDigitsPointer

const uint32_t* LArRodBlockPhysicsV5::m_MaskDigitsPointer = nullptr
private

Definition at line 146 of file LArRodBlockPhysicsV5.h.

◆ m_MaskTimeQualityPointer

const uint32_t* LArRodBlockPhysicsV5::m_MaskTimeQualityPointer = nullptr
private

Definition at line 145 of file LArRodBlockPhysicsV5.h.

◆ m_mFebBlocks

FEBMAPTYPE LArRodBlockStructure::m_mFebBlocks
protectedinherited

Definition at line 228 of file LArRodBlockStructure.h.

◆ m_MiddleHeaderSize

int32_t LArRodBlockStructure::m_MiddleHeaderSize
protectedinherited

Definition at line 233 of file LArRodBlockStructure.h.

◆ m_numberHotCell

uint16_t LArRodBlockPhysicsV5::m_numberHotCell = 0U
private

Definition at line 156 of file LArRodBlockPhysicsV5.h.

◆ m_numberHotCellOffTime

uint16_t LArRodBlockPhysicsV5::m_numberHotCellOffTime
private

Definition at line 157 of file LArRodBlockPhysicsV5.h.

◆ m_numberOfEx

uint32_t LArRodBlockStructure::m_numberOfEx
protectedinherited

Definition at line 210 of file LArRodBlockStructure.h.

◆ m_numberOfEy

uint32_t LArRodBlockStructure::m_numberOfEy
protectedinherited

Definition at line 211 of file LArRodBlockStructure.h.

◆ m_numberOfEz

uint32_t LArRodBlockStructure::m_numberOfEz
protectedinherited

Definition at line 212 of file LArRodBlockStructure.h.

◆ m_numberOfSumE

uint32_t LArRodBlockStructure::m_numberOfSumE
protectedinherited

Definition at line 213 of file LArRodBlockStructure.h.

◆ m_OfflineToRawGainMap

const uint32_t LArRodBlockStructure::m_OfflineToRawGainMap ={3,2,1}
staticprotectedinherited

Definition at line 17 of file LArRodBlockStructure.h.

◆ m_OffTimeCut

int16_t LArRodBlockPhysicsV5::m_OffTimeCut
private

Definition at line 161 of file LArRodBlockPhysicsV5.h.

◆ m_onlineHelper

const LArOnlineID* LArRodBlockPhysicsV5::m_onlineHelper
private

Definition at line 164 of file LArRodBlockPhysicsV5.h.

◆ m_pRODblock

std::vector<uint32_t>* LArRodBlockStructure::m_pRODblock
protectedinherited

Definition at line 226 of file LArRodBlockStructure.h.

◆ m_RaddPointer

const uint16_t* LArRodBlockPhysicsV5::m_RaddPointer = nullptr
private

Definition at line 147 of file LArRodBlockPhysicsV5.h.

◆ m_RawDataBlock

std::vector<uint32_t> LArRodBlockPhysicsV5::m_RawDataBlock
private

Definition at line 135 of file LArRodBlockPhysicsV5.h.

◆ m_RawDataIndex

int LArRodBlockPhysicsV5::m_RawDataIndex = 0
private

Definition at line 143 of file LArRodBlockPhysicsV5.h.

◆ m_RawDataPointer

const uint16_t* LArRodBlockPhysicsV5::m_RawDataPointer = nullptr
private

Definition at line 152 of file LArRodBlockPhysicsV5.h.

◆ m_RawToOfflineGainMap

const uint32_t LArRodBlockStructure::m_RawToOfflineGainMap ={0, 2, 1,0}
staticprotectedinherited

Definition at line 16 of file LArRodBlockStructure.h.

◆ m_rearrangeFirstSample

unsigned int LArRodBlockStructure::m_rearrangeFirstSample
protectedinherited

Definition at line 235 of file LArRodBlockStructure.h.

◆ m_requiredNSamples

unsigned short LArRodBlockPhysicsV5::m_requiredNSamples
private

Definition at line 163 of file LArRodBlockPhysicsV5.h.

◆ m_ROB_to_decode

int32_t LArRodBlockStructure::m_ROB_to_decode
protectedinherited

Definition at line 240 of file LArRodBlockStructure.h.

◆ m_RodBlock

const uint32_t* LArRodBlockStructure::m_RodBlock
protectedinherited

Definition at line 222 of file LArRodBlockStructure.h.

◆ m_RodBlockSize

int32_t LArRodBlockStructure::m_RodBlockSize
protectedinherited

Definition at line 224 of file LArRodBlockStructure.h.

◆ m_SumBlkBlockE1

std::vector<uint32_t> LArRodBlockPhysicsV5::m_SumBlkBlockE1
private

Definition at line 130 of file LArRodBlockPhysicsV5.h.

◆ m_SumBlkBlockE2

std::vector<uint32_t> LArRodBlockPhysicsV5::m_SumBlkBlockE2
private

Definition at line 131 of file LArRodBlockPhysicsV5.h.

◆ m_SumE

int32_t LArRodBlockStructure::m_SumE
protectedinherited

Definition at line 206 of file LArRodBlockStructure.h.

◆ m_SumPointer

const int32_t* LArRodBlockPhysicsV5::m_SumPointer = nullptr
private

Definition at line 149 of file LArRodBlockPhysicsV5.h.

◆ m_TimeQualityBlock

std::vector<uint16_t> LArRodBlockPhysicsV5::m_TimeQualityBlock
private

Definition at line 134 of file LArRodBlockPhysicsV5.h.

◆ m_TimeQualityIndex

int LArRodBlockPhysicsV5::m_TimeQualityIndex = 0
private

Definition at line 140 of file LArRodBlockPhysicsV5.h.

◆ m_TimeQualityPointer

const uint16_t* LArRodBlockPhysicsV5::m_TimeQualityPointer = nullptr
private

Definition at line 150 of file LArRodBlockPhysicsV5.h.

◆ m_vFragment

std::vector<uint32_t>* LArRodBlockStructure::m_vFragment
protectedinherited

Definition at line 225 of file LArRodBlockStructure.h.

◆ m_virtualROBJump

uint32_t LArRodBlockStructure::m_virtualROBJump
protectedinherited

Definition at line 239 of file LArRodBlockStructure.h.

◆ m_virtualROBPointer

const uint32_t* LArRodBlockStructure::m_virtualROBPointer
protectedinherited

Definition at line 241 of file LArRodBlockStructure.h.

◆ m_virtualROBPointerLocal

const uint32_t* LArRodBlockStructure::m_virtualROBPointerLocal
protectedinherited

Definition at line 242 of file LArRodBlockStructure.h.


The documentation for this class was generated from the following files: