MakeLArCellFromRaw Class Reference

Building LArCell objects from LArRawChannel. More...

#include <MakeLArCellFromRaw.h>

Collaboration diagram for MakeLArCellFromRaw:

Classes
struct	CellInfo

Public Member Functions
	MakeLArCellFromRaw ()
	constructor
	~MakeLArCellFromRaw ()
	destructor
void	initialize (const LArRoIMap &roiMap, const LArOnOffIdMapping &onOffMap, const CaloDetDescrManager &man, const std::vector< const CaloCellCorrection * > *pCorr, unsigned int poolMaxSize=190000)
	initialize the internal map
LArCell *	getLArCell (unsigned int feb, unsigned int chan, int e, int t, int q, CaloGain::CaloGain caloGain, unsigned int &ttId) const
	access by Identifier
LArCell *	getLArCell (const HWIdentifier &id, int e, int t, int q, CaloGain::CaloGain g) const
LArCell *	getLArCell (unsigned int feb, unsigned int chan, int e, int t, int q, unsigned int &ttId) const
LArCell *	getLArCell (const HWIdentifier &id, int e, int t, int q) const
void	setThreshold (float t)

Private Types
typedef std::vector< CellInfo >	CELL_VEC
typedef std::map< unsigned int, CELL_VEC >	CELL_MAP

Private Member Functions
double	getCorrection (LArCell *cell, const std::vector< const CaloCellCorrection * > &vCorr, const EventContext &ctx) const

Private Attributes
CELL_MAP	m_cellMap
size_t	m_poolMaxSize
IMessageSvc *	m_msgSvc
const LArOnlineID *	m_onlineID
float	m_ethreshold

Detailed Description

Building LArCell objects from LArRawChannel.

e_threshold:          energy threshold

Definition at line 30 of file MakeLArCellFromRaw.h.

Member Typedef Documentation

CELL_MAP

typedef std::map<unsigned int, CELL_VEC > MakeLArCellFromRaw::CELL_MAP

private

Definition at line 92 of file MakeLArCellFromRaw.h.

CELL_VEC

typedef std::vector<CellInfo> MakeLArCellFromRaw::CELL_VEC

private

Definition at line 89 of file MakeLArCellFromRaw.h.

Constructor & Destructor Documentation

MakeLArCellFromRaw()

MakeLArCellFromRaw::MakeLArCellFromRaw

(

)

constructor

Definition at line 31 of file MakeLArCellFromRaw.cxx.

  : m_poolMaxSize(0),
    m_msgSvc (nullptr),
    m_onlineID (nullptr),
    m_ethreshold(-1.e5) 
{
}

~MakeLArCellFromRaw()

MakeLArCellFromRaw::~MakeLArCellFromRaw

(

)

destructor

Definition at line 40 of file MakeLArCellFromRaw.cxx.

{
 return ; 
}

Member Function Documentation

getCorrection()

double MakeLArCellFromRaw::getCorrection ( LArCell * cell, const std::vector< const CaloCellCorrection * > & vCorr, const EventContext & ctx ) const

private

Definition at line 321 of file MakeLArCellFromRaw.cxx.

{
// LArCell was made with energy = 50. 
    double en= 50. * GeV ; 
//  cell->setEnergy(en); 
 
        // apply corrections. 
        for (const CaloCellCorrection* corr : vCorr)
         {
           corr->MakeCorrection (cell, ctx);
         } 
 
    double c=  cell->energy()/en; 
    return c;        
}

getLArCell() [1/4]

LArCell * MakeLArCellFromRaw::getLArCell	(	const HWIdentifier &	id,
		int	e,
		int	t,
		int	q ) const

Definition at line 244 of file MakeLArCellFromRaw.cxx.

{
  return getLArCell(id, e, t, q, CaloGain::UNKNOWNGAIN);
  
}

getLArCell() [2/4]

LArCell * MakeLArCellFromRaw::getLArCell	(	const HWIdentifier &	id,
		int	e,
		int	t,
		int	q,
		CaloGain::CaloGain	g ) const

Definition at line 252 of file MakeLArCellFromRaw.cxx.

{
// GU to change
 
  HWIdentifier feb   = m_onlineID->feb_Id(id);
  unsigned int chan  = m_onlineID->channel(id); 
  unsigned int tt; // not needed for this method. 
//  unsigned int ifeb = m_onlineID->feb_Hash(feb);
  return getLArCell(feb.get_identifier32().get_compact() ,  chan , e,  t,  q, g, tt) ; 
} 

getLArCell() [3/4]

LArCell * MakeLArCellFromRaw::getLArCell	(	unsigned int	feb,
		unsigned int	chan,
		int	e,
		int	t,
		int	q,
		CaloGain::CaloGain	caloGain,
		unsigned int &	ttId ) const

access by Identifier

Definition at line 272 of file MakeLArCellFromRaw.cxx.

{
 
  CELL_MAP::const_iterator it =m_cellMap.find( feb ); 
 
  if(it == m_cellMap.end()){
    MsgStream log(m_msgSvc, "MakeLArCellFromRaw");
    log << MSG::ERROR <<" MakeLArCellFromRaw ERROR, failed to find existing cells. "<< endmsg; 
    std::cout <<std::hex<<"FEBID = "<< feb <<std::dec<< std::endl ; 
    //DR    assert(0); 
    return nullptr;
  }
 
  const CELL_VEC& cellVec = (*it).second; 
  if(cellVec.size()<= chan){
    MsgStream log(m_msgSvc, "MakeLArCellFromRaw");
    log << MSG::FATAL <<" MakeLArCellFromRaw ERROR, channel number= "<<chan<< endmsg ;
    assert(0);
    return nullptr;
  }
  
  const CellInfo& info = cellVec[chan] ; 
  // do not create cell if it is not connected. 
  if(info.elem == nullptr) return nullptr; 
 
  //remove /1000 (GeV->MeV)
  double de = e * info.eCorr; 
  LArCell* cell =nullptr;
 
  // DR convert time from ps (in LArRawChannel ) to ns
  double time = t/1000.0;
 
  if (m_poolMaxSize > 0) {
    DataPool<LArCell> pool;
    cell = new ( pool.nextElementPtr()) LArCell(info.elem,de,time,q,g); 
  }
  else
    {
    cell = new LArCell(info.elem,de,time,q,g); 
    }
 
  ttId = info.tt; 
 
  return cell; 
 
}

getLArCell() [4/4]

LArCell * MakeLArCellFromRaw::getLArCell	(	unsigned int	feb,
		unsigned int	chan,
		int	e,
		int	t,
		int	q,
		unsigned int &	ttId ) const

Definition at line 264 of file MakeLArCellFromRaw.cxx.

{
  return getLArCell(feb, chan ,e, t, q,  CaloGain::UNKNOWNGAIN, ttId );
  
}

initialize()

void MakeLArCellFromRaw::initialize	(	const LArRoIMap &	roiMap,
		const LArOnOffIdMapping &	onOffMap,
		const CaloDetDescrManager &	man,
		const std::vector< const CaloCellCorrection * > *	pCorr,
		unsigned int	poolMaxSize = 190000 )

initialize the internal map

Definition at line 45 of file MakeLArCellFromRaw.cxx.

{
  const EventContext& ctx = Gaudi::Hive::currentContext();
 
  m_poolMaxSize = poolMaxSize;
 
  MsgStream log(Athena::getMessageSvc(), "MakeLArCellFromRaw");
 
  log << MSG::INFO <<" initializing MakeLArCellFromRaw"<< endmsg ;
 
  if(pCorr) log <<MSG::INFO <<" Number of Corrections "<< pCorr->size()<< endmsg ;
  
 
  SmartIF<StoreGateSvc> detStore{Gaudi::svcLocator( )->service("DetectorStore")};
  if (!detStore)
  {
    log << MSG::ERROR << "MakeLArCellFromRaw ERROR cannot retrieve DetectorStore " << endmsg;
    return;
  }
 
  
  // create DataPool if requested
  if (poolMaxSize!=0) {
    log << MSG::INFO << "MakeLArCellFromRaw Creating DataPool<LArCell> of size " << poolMaxSize << endmsg ;
  }else{
    log << MSG::INFO << "MakeLArCellFromRaw do not use DataPool" << endmsg ;
  }
 
 
  const CaloCell_ID* calo_id = man.getCaloCell_ID();
  const LArEM_ID& em_id = *calo_id->em_idHelper();
  const LArHEC_ID& hec_id = *calo_id->hec_idHelper();
  const LArFCAL_ID& fcal_id = *calo_id->fcal_idHelper();
 
  if (detStore->retrieve(m_onlineID, "LArOnlineID").isFailure()) {
    log <<MSG::ERROR << "cannot find LArOnlineID in MakeLArCellFromRaw " << endmsg;
  }
 
 
  int n_em   = 0 ; 
  int n_hec  = 0 ; 
  int n_fcal = 0 ; 
 
  int n_em_err   = 0 ; 
  int n_hec_err  = 0 ; 
  int n_fcal_err = 0 ; 
 
 
  const std::vector<Identifier>& emIds   = em_id.channel_ids();
  const std::vector<Identifier>& hecIds  = hec_id.channel_ids();
  const std::vector<Identifier>& fcalIds = fcal_id.channel_ids();
 
  double en=50. * GeV;
  double time=0;
  double qu=0;
 
  CellInfo info0; 
  info0.eta=0; 
  info0.phi=0; 
  info0.x=0; 
  info0.y=0; 
  info0.z=0; 
  info0.fcal=false; 
  info0.tt=0; 
  info0.eCorr=1.; 
  info0.elem = nullptr ; 
 
  // EM
  std::vector<Identifier>::const_iterator it = emIds.begin();
  std::vector<Identifier>::const_iterator it2 = emIds.end();
  for(; it!=it2;++it)
  {
    const Identifier& chan_id = *it ;
 
    try{ 
     
      HWIdentifier sigId = onOffMap.createSignalChannelID(chan_id);
      HWIdentifier feb   = m_onlineID->feb_Id(sigId);
      CELL_VEC::size_type chan  = m_onlineID->channel(sigId);   
      CELL_VEC& cellVec= m_cellMap[ feb.get_identifier32().get_compact() ] ; 
      if (cellVec.empty())cellVec.resize(128,info0); 
 
      CellInfo& cell = cellVec[chan] ;
 
      cell.tt  = roiMap.trigTowerID(sigId); 
 
      const CaloDetDescrElement* caloDDE =man.get_element( chan_id);
 
      cell.id  =  chan_id ; 
      cell.eta =  caloDDE->eta();
      cell.phi =  caloDDE->phi(); 
      cell.fcal= false; 
      cell.elem =caloDDE; 
// to compute correction, assume that cell is in mid gain 
      const CaloGain::CaloGain g=CaloGain::LARMEDIUMGAIN;
      
      ++n_em; 
      if(pCorr)  {
    // make a fake LArCell, and get the correction factor. 
       LArCell larCell(caloDDE,en,time,qu,g); 
           cell.eCorr   = getCorrection(&larCell, *pCorr, ctx);
      }
    } catch (LArID_Exception& ex){
                    ++n_em_err; 
    }
  } // end of EM 
 
 
 
  // HEC
  it = hecIds.begin();
  it2 = hecIds.end();
  for(; it!=it2;++it)
  {
    const Identifier& chan_id = *it ;
 
    try{ 
     
      HWIdentifier sigId = onOffMap.createSignalChannelID(chan_id);
// to change using LArOnline methods
      HWIdentifier feb   = m_onlineID->feb_Id(sigId);
      CELL_VEC::size_type chan  = m_onlineID->channel(sigId);  
      CELL_VEC& cellVec= m_cellMap[ feb.get_identifier32().get_compact() ] ; 
      if (cellVec.empty())cellVec.resize(128,info0); 
 
      CellInfo& cell = cellVec[chan] ;
 
      cell.tt  = roiMap.trigTowerID(sigId); 
 
      const CaloDetDescrElement* caloDDE =man.get_element( chan_id);
 
      cell.id  =  chan_id ; 
      cell.eta =  caloDDE->eta();
      cell.phi = caloDDE->phi(); 
      cell.fcal= false; 
      cell.elem =caloDDE; 
      const CaloGain::CaloGain g=CaloGain::LARMEDIUMGAIN ;
 
      if(pCorr)  {
    // make a fake LArCell, and get the correction factor. 
       LArCell larCell(caloDDE,en,time,qu,g); 
           cell.eCorr   = getCorrection(&larCell, *pCorr, ctx);
      }
      ++n_hec; 
    } catch (LArID_Exception& ex){
                    ++n_hec_err; 
    }
  } // end of HEC
 
  // FCAL
  it = fcalIds.begin();
  it2 = fcalIds.end();
  for(; it!=it2;++it)
  {
    const Identifier& chan_id = *it ;
 
    try{ 
     
      HWIdentifier sigId = onOffMap.createSignalChannelID(chan_id);
// GU to change to use LArOnlineID method
      HWIdentifier feb   = m_onlineID->feb_Id(sigId);
      CELL_VEC::size_type chan  = m_onlineID->channel(sigId);  
      CELL_VEC& cellVec= m_cellMap[ feb.get_identifier32().get_compact() ] ; 
      if (cellVec.empty())cellVec.resize(128,info0); 
 
      CellInfo& cell = cellVec[chan] ;
 
      cell.tt  = roiMap.trigTowerID(sigId); 
 
      const CaloDetDescrElement* caloDDE =man.get_element( chan_id);
 
      cell.id  =  chan_id ; 
      cell.eta =  caloDDE->eta();
      cell.phi = caloDDE->phi(); 
      cell.fcal= false; 
      cell.elem =caloDDE; 
      const CaloGain::CaloGain g=CaloGain::LARMEDIUMGAIN ;
 
      if(pCorr)  {
    // make a fake LArCell, and get the correction factor. 
       LArCell larCell(caloDDE,en,time,qu,g); 
           cell.eCorr   = getCorrection(&larCell, *pCorr, ctx);
      }
      ++n_fcal; 
    } catch (LArID_Exception& ex){
                    ++n_fcal_err; 
    }
  } // end of FCAL
 
 
  log <<MSG::INFO <<" number of em, hec, fcal cells="<<n_em<<" "<<n_hec<<" "<<n_fcal<<endmsg;
  log <<MSG::INFO<<" number of exceptions for em,hec,fcal="<<n_em_err<<" "<<n_hec_err<<" "<<n_fcal_err<< endmsg;
  log <<MSG::INFO<<" done with initializing MakeLArCellFromRaw"<< endmsg;
  return; 
}

setThreshold()

void MakeLArCellFromRaw::setThreshold

(

float

t

)

Definition at line 341 of file MakeLArCellFromRaw.cxx.

{
  m_ethreshold = t; 
  return; 
}

Member Data Documentation

m_cellMap

CELL_MAP MakeLArCellFromRaw::m_cellMap

private

Definition at line 94 of file MakeLArCellFromRaw.h.

m_ethreshold

float MakeLArCellFromRaw::m_ethreshold

private

Definition at line 101 of file MakeLArCellFromRaw.h.

m_msgSvc

IMessageSvc* MakeLArCellFromRaw::m_msgSvc

private

Definition at line 97 of file MakeLArCellFromRaw.h.

m_onlineID

const LArOnlineID* MakeLArCellFromRaw::m_onlineID

private

Definition at line 99 of file MakeLArCellFromRaw.h.

m_poolMaxSize

size_t MakeLArCellFromRaw::m_poolMaxSize

private

Definition at line 96 of file MakeLArCellFromRaw.h.

---

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

• MakeLArCellFromRaw.h
• MakeLArCellFromRaw.cxx