LArRodEncoder Class Reference

This class provides conversion from LArRawChannel and LArDigit to ROD format. More...

#include <LArRodEncoder.h>

Collaboration diagram for LArRodEncoder:

Classes
struct	FebData_t

Public Member Functions
	LArRodEncoder (const LArOnlineID &onlineHelper, const CaloDetDescrManager &calodd, const LArOnOffIdMapping &onOffIdMapping, LArRodBlockStructure *BlStruct)
	~LArRodEncoder ()
void	add (const LArRawChannel *rc)
void	add (const LArDigit *digit)
void	add (const LArDigit *dg, const int gain)
void	add (const LArCalibDigit *dg, const int gain)
void	clear ()
void	fillROD (std::vector< uint32_t > &v, MsgStream &logstr, const CaloNoise &noise, double nsigma)

Private Attributes
std::map< uint32_t, FebData_t >	m_mFEB
const LArOnlineID &	m_onlineHelper
const CaloDetDescrManager &	m_CaloDetDescrManager
const LArOnOffIdMapping &	m_onOffIdMapping
LArRodBlockStructure *	m_BlStruct

Detailed Description

This class provides conversion from LArRawChannel and LArDigit to ROD format.

Author

H. Ma

Version

0-0-1 , Oct 7, 2002

Modified, Jan 02, 2003 Split from LArROD_Decoder.

Modified, Jul 08 2003 Walter Lampl Read ByteStream for LArDigits and LArRawChannel

Modified, Jul 17 2003 Remi Lafaye Split with LArRODBlockStructure

• LArROD_Encoder = Athena world only
• LArRODBlockStructure = ROD format only It might be less CPU efficient Althought it will be easier to update

Modified, Nov 10 2005 D.O. Damazio New LArRodBlockPhysics (Bertrand Laforge format).

Definition at line 54 of file LArRodEncoder.h.

Constructor & Destructor Documentation

LArRodEncoder()

LArRodEncoder::LArRodEncoder	(	const LArOnlineID &	onlineHelper,
		const CaloDetDescrManager &	calodd,
		const LArOnOffIdMapping &	onOffIdMapping,
		LArRodBlockStructure *	BlStruct
	)

Definition at line 19 of file LArRodEncoder.cxx.

  : m_onlineHelper (onlineHelper),
    m_CaloDetDescrManager (calodd),
    m_onOffIdMapping (onOffIdMapping),
    m_BlStruct (BlStruct)
{
}

~LArRodEncoder()

LArRodEncoder::~LArRodEncoder

(

)

Definition at line 31 of file LArRodEncoder.cxx.

{
  clear();
}

Member Function Documentation

add() [1/4]

void LArRodEncoder::add	(	const LArCalibDigit *	dg,
		const int	gain
	)

Definition at line 61 of file LArRodEncoder.cxx.

{
 uint32_t FEB_ID = (m_onlineHelper.feb_Id(dg->hardwareID()).get_identifier32().get_compact());
 m_mFEB[FEB_ID].vLArCalibDigit[gain].push_back(dg);
}

add() [2/4]

void LArRodEncoder::add	(	const LArDigit *	dg,
		const int	gain
	)

Definition at line 55 of file LArRodEncoder.cxx.

{
 uint32_t FEB_ID = (m_onlineHelper.feb_Id(dg->hardwareID()).get_identifier32().get_compact());
 m_mFEB[FEB_ID].vLArDigitFixed[gain].push_back(dg);
}

add() [3/4]

void LArRodEncoder::add

(

const LArDigit *

digit

)

Definition at line 49 of file LArRodEncoder.cxx.

{
 uint32_t FEB_ID = (m_onlineHelper.feb_Id(dg->hardwareID()).get_identifier32().get_compact());
 m_mFEB[FEB_ID].vLArDigit.push_back(dg);
}

add() [4/4]

void LArRodEncoder::add

(

const LArRawChannel *

rc

)

Definition at line 37 of file LArRodEncoder.cxx.

{
 uint32_t FEB_ID = (m_onlineHelper.feb_Id(rc->channelID()).get_identifier32().get_compact());
 if ( m_mFEB[FEB_ID].vLArRC.empty() ){
    m_mFEB[FEB_ID].vLArRC.resize(128,0);
    //for(int i=0;i<128;i++) m_mFEB[FEB_ID].vLArRC.push_back(0);
 }
 uint32_t chan = m_BlStruct->FebToRodChannel(m_onlineHelper.channel( rc->channelID() ) );
 m_mFEB[FEB_ID].vLArRC[chan]=rc;
}

clear()

void LArRodEncoder::clear

(

)

Definition at line 71 of file LArRodEncoder.cxx.

{
  m_mFEB.clear();
}

fillROD()

void LArRodEncoder::fillROD	(	std::vector< uint32_t > &	v,
		MsgStream &	logstr,
		const CaloNoise &	noise,
		double	nsigma
	)

Definition at line 78 of file LArRodEncoder.cxx.

{
 if (!m_BlStruct)
  {logstr << MSG::ERROR << "No LArRodBlockStructure defined! Can't encode fragment!" << endmsg; 
   return;
  }
 //int febcounter=0; 
 //std::cout << "Going to init Fragment..." << std::endl;
 m_BlStruct->initializeFragment(v); //Makes new Fragment or splits v into existing Feb-Blocks
 std::map<uint32_t,FebData_t>::iterator it=m_mFEB.begin();
 std::map<uint32_t,FebData_t>::iterator it_end=m_mFEB.end();
 // The sort alorithm for the vectors inside the data object of the map fails
 // when I use const_iterator at this point. Don't ask me why
 //std::cout << "Have " << m_mFEB.size() << " FEBs in this ROD" << std::endl;
 for(;it!=it_end;++it) {
   //std::cout << "Setting values for feb # " << febcounter << " ID=" << it->first << std::endl;
    m_BlStruct->initializeFEB(it->first);                                     // Set FEB id
    //std::cout << "Setting Energy" << std::endl;
    // ************* Energy Block *************
    if (m_BlStruct->canSetEnergy() && it->second.vLArRC.size()>0) {
      //int j=0;
//      m_BlStruct->sortDataVector(it->second.vLArRC);
      // compute Ex, Ey for this FEB : initialisation phase
      
      {
    double Ex=0;
    double Ey=0;
    double Ez=0;
    double SumE=0;
    Identifier myofflineID;
 
        for (const LArRawChannel *theChannel : it->second.vLArRC) {
      if ( theChannel != nullptr ){
      int cId =  m_onlineHelper.channel(theChannel->hardwareID());
 
      int e=theChannel->energy();
          uint32_t quality = theChannel->quality();
      m_BlStruct->setNextEnergy(cId,e, theChannel->time(), 
                    quality,theChannel->gain());
      
      // you convert from hardwareID to offline channle ID hash (???)
      myofflineID = m_onOffIdMapping.cnvToIdentifier(theChannel->hardwareID()) ;
      //std::cout << "Got Offile id 0x" << std::hex << myofflineID.get_compact() << std::dec << std::endl;
          const CaloDetDescrElement* caloDDE = m_CaloDetDescrManager.get_element(myofflineID);
      // This is probably NOT what one wants. You want the cell gain!
          double cellnoise = noise.getNoise(myofflineID,theChannel->gain());
          if( e > (nsigma*cellnoise) && (quality != 65535 ) ){
        double aux=caloDDE->sinTh();
        double aux_z=tanh(caloDDE->eta());
            Ex += e*aux*caloDDE->cosPhi(); 
        Ey += e*aux*caloDDE->sinPhi();
        Ez += e*aux_z;
        SumE += e;
          }
     } // End of if to check for the LArRawChannel existence
        } // end of chit loop
    // Set after the loop finishes
    m_BlStruct->setEx(Ex);
    m_BlStruct->setEy(Ey);
    m_BlStruct->setEz(Ez);
    m_BlStruct->setSumE(SumE);
      }
    }//End canSetEnergy
    // try to include samples together with data
    if (m_BlStruct->canIncludeRawData() && it->second.vLArDigit.size()>0) {
      // Order channels according to ROD numbering
      m_BlStruct->sortDataVector(it->second.vLArDigit);
      for (const LArDigit* digit : it->second.vLArDigit) {
        int cId =  m_onlineHelper.channel(digit->hardwareID());
        m_BlStruct->setRawData(cId, digit->samples(), digit->gain());
      } // end of for digits
    } // End of can Include Raw Data check
    //std::cout << "Setting Raw Data" << std::endl;
    // ************** Raw Data Block ***********
    if (m_BlStruct->canSetRawData() && it->second.vLArDigit.size()>0) {
      //sort(it->second.vLArDigit.begin(), it->second.vLArDigit.end(), *m_BlStruct);  // Order channels according to ROD numbering
      m_BlStruct->sortDataVector(it->second.vLArDigit);
      std::vector<const LArDigit*>::const_iterator digit_it = it->second.vLArDigit.begin();
      std::vector<const LArDigit*>::const_iterator digit_it_end=it->second.vLArDigit.end();
      if(digit_it!=digit_it_end) { //Container not empty
    m_BlStruct->setNumberOfSamples((*digit_it)->samples().size());
    for (;digit_it!=digit_it_end;++digit_it)  {
      int cId =  m_onlineHelper.channel((*digit_it)->hardwareID());
      //cId = m_BlStruct->FebToRodChannel(cId);
      m_BlStruct->setRawData(cId, (*digit_it)->samples(), (*digit_it)->gain());
    } 
      }// end container not empty
    }// end if can set raw data    
 
    //std::cout << "Setting fixed gain data" << std::endl;
    // ************** Fixed Gain Data Block ***********
    if (m_BlStruct->canSetRawDataFixed()) {
      for (int i=0;i<3;i++) { //Loop over three gains
    //sort(it->second.vLArDigitFixed[i].begin(), it->second.vLArDigitFixed[i].end(), *m_BlStruct);
    m_BlStruct->sortDataVector(it->second.vLArDigitFixed[i]);
    std::vector<const LArDigit*>::const_iterator digit_it = it->second.vLArDigitFixed[i].begin();
    std::vector<const LArDigit*>::const_iterator digit_it_end=it->second.vLArDigitFixed[i].end();
    if(digit_it!=digit_it_end) {//Container not empty
      m_BlStruct->setNumberOfSamples((*digit_it)->samples().size());
      for (;digit_it!=digit_it_end;++digit_it)  {
        int cId =  m_onlineHelper.channel((*digit_it)->hardwareID()); 
        //cId = m_BlStruct->FebToRodChannel(cId);
        m_BlStruct->setRawDataFixed(cId, (*digit_it)->samples(), (*digit_it)->gain());
      }
    } //end Container not empty
      } // end loop over three gains
    } // end can set raw data fixed
  
    // ************** Calibration Digit Block ***********
    if (m_BlStruct->canSetCalibration()) { 
        for (int i=0;i<3;i++) { //Loop over three gains
      //sort(it->second.vLArCalibDigit[i].begin(), it->second.vLArCalibDigit[i].end(), *m_BlStruct);
      m_BlStruct->sortDataVector(it->second.vLArCalibDigit[i]);
      std::vector<const LArCalibDigit*>::const_iterator digit_it = it->second.vLArCalibDigit[i].begin();
      std::vector<const LArCalibDigit*>::const_iterator digit_it_end=it->second.vLArCalibDigit[i].end();
      if(digit_it!=digit_it_end) {//Container not empty
        m_BlStruct->setNumberOfSamples((*digit_it)->samples().size());
        m_BlStruct->setDelay((*digit_it)->delay());
        m_BlStruct->setDAC((*digit_it)->DAC());
        for (;digit_it!=digit_it_end;++digit_it) { 
          int cId =  m_onlineHelper.channel((*digit_it)->hardwareID()); 
          //cId = m_BlStruct->FebToRodChannel(cId);
          m_BlStruct->setRawDataFixed(cId, (*digit_it)->samples(), (*digit_it)->gain());
          if ((*digit_it)->isPulsed())
        m_BlStruct->setPulsed(cId);
        }
      } //End Container not empty
    }// end loop over three gains
    } //end can set calibration
    
    /*
    // ************** Averaged Calibration Digit Block ***********
    if (m_BlStruct->canSetAverages()) { 
      for (int i=0;i<3;i++) { //Loop over three gains
    sort(it->second.vLArCalibDigit[i].begin(), it->second.vLArCalibDigit[i].end(), *m_BlStruct);
    std::vector<const LArAverageDigit*>::const_iterator digit_it = it->second.vLArAverageDigit.begin();
    std::vector<const LArAverageDigit*>::const_iterator digit_it_end=it->second.vLArAverageDigit.end();
    if(digit_it!=digit_it_end) {//Container not empty
      m_BlStruct->setNumberOfSamples((*digit_it)->samples().size());
      m_BlStruct->setDAC((*digit_it)->DAC());
      m_BlStruct->setDelay((*digit_it)->delay());
      for (;digit_it!=digit_it_end;digit_it++) { 
        int cId =  m_onOffIdMapping.channel((*digit_it)->channelID()); 
        cId = m_BlStruct->FebToRodChannel(cId);
        m_BlStruct->setRawData(cId, (*digit_it)->samples(), (*digit_it)->gain());
        m_BlStruct->setIsPulsed(cId);
      }
    } //End Container not empty
      }// end loop over three gains
    }// can set averages
    */ 
    //std::cout << "Going to finalize FEB" << std::endl;
    m_BlStruct->finalizeFEB();
 }// end iterating over FEB-map
 //std::cout << "Going to finalize ROD" << std::endl;
 m_BlStruct->concatinateFEBs();
 //std::cout << "final size of fragment " << v.size() << std::endl;
 //m_BlStruct->dumpFragment(); // For testing purpose
 return; 
}

Member Data Documentation

m_BlStruct

LArRodBlockStructure* LArRodEncoder::m_BlStruct

private

Definition at line 118 of file LArRodEncoder.h.

m_CaloDetDescrManager

const CaloDetDescrManager& LArRodEncoder::m_CaloDetDescrManager

private

Definition at line 116 of file LArRodEncoder.h.

m_mFEB

std::map<uint32_t,FebData_t> LArRodEncoder::m_mFEB

private

Definition at line 113 of file LArRodEncoder.h.

m_onlineHelper

const LArOnlineID& LArRodEncoder::m_onlineHelper

private

Definition at line 115 of file LArRodEncoder.h.

m_onOffIdMapping

const LArOnOffIdMapping& LArRodEncoder::m_onOffIdMapping

private

Definition at line 117 of file LArRodEncoder.h.

---

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

• LArRodEncoder.h
• LArRodEncoder.cxx