TileCellBuilderFromHit Class Reference

This class creates Cells from RawChannels and stores them in a container. More...

#include <TileCellBuilderFromHit.h>

Inheritance diagram for TileCellBuilderFromHit:

Collaboration diagram for TileCellBuilderFromHit:

Public Member Functions
	TileCellBuilderFromHit (const std::string &type, const std::string &name, const IInterface *parent)
	Contructor. More...
virtual	~TileCellBuilderFromHit ()
	Destructor. More...
virtual StatusCode	initialize () override
	Initializer. More...
virtual StatusCode	finalize () override
virtual StatusCode	process (CaloCellContainer *theCellContainer, const EventContext &ctx) const override
	method to process all raw channels and store them in container More...

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Private Types
typedef TileDrawerEvtStatus	TileDrawerEvtStatusArray[5][64]
	status of every drawer More...

Private Member Functions
template<class ITERATOR , class COLLECTION >
void	build (const CaloNoise *caloNoise, TileDrawerEvtStatusArray &drawerEvtStatus, const ITERATOR &begin, const ITERATOR &end, COLLECTION *coll, TileCellContainer *MBTSCells, TileCellContainer *E4prCells, const TileSamplingFraction *samplingFraction) const
bool	maskBadChannel (TileDrawerEvtStatusArray &drawerEvtStatus, TileCell *pCell) const
	method to check if channels are good or bad. More...
bool	maskBadChannels (TileDrawerEvtStatusArray &drawerEvtStatus, TileCell *pCell, bool single_PMT_C10, bool Ecell) const
void	correctCell (TileCell *pCell, int correction, int pmt, int gain, float ener, float time, unsigned char iqual, unsigned char qbit) const
	Compute calibrated energy, time, etc. More...
unsigned char	iquality (float qual) const
unsigned char	qbits (TileDrawerEvtStatusArray &drawerEvtStatus, int ros, int drawer, bool count_over, bool good_time, bool good_ener, bool overflow, bool underflow, bool good_overflowfit) const
	method to compute the cell quality bits More...
int	mbts_index (int side, int phi, int eta) const
int	e4pr_index (int phi) const

Private Attributes
SG::ReadCondHandleKey< CaloNoise >	m_caloNoiseKey
SG::ReadHandleKey< TileHitContainer >	m_hitContainerKey
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfoKey
SG::WriteHandleKey< TileCellContainer >	m_MBTSContainerKey
SG::WriteHandleKey< TileCellContainer >	m_E4prContainerKey
SG::ReadCondHandleKey< TileSamplingFraction >	m_samplingFractionKey
	Name of TileSamplingFraction in condition store. More...
std::string	m_infoName
float	m_eneForTimeCut
	keep time for channels with energy above cut More...
float	m_eneForTimeCutMBTS
	similar cut for MBTS in pC More...
float	m_zeroEnergy
	energy to store in every PMT if both PMT are bad More...
int	m_qualityCut
	cut on channel quality (set energy to m_zeroEnergy for them) More...
float	m_maxTime
	maximum time for the PMTs in the cels More...
float	m_minTime
	minimum time for the PMTs in the cels More...
bool	m_maskBadChannels
	if true=> bad channels are masked More...
float	m_noiseSigma
	cell electronic noise if CaloNoise is switched off More...
const TileID *	m_tileID {nullptr}
	Pointer to TileID. More...
const TileTBID *	m_tileTBID {nullptr}
	Pointer to TileTBID. More...
const TileHWID *	m_tileHWID {nullptr}
	Pointer to TileHWID. More...
const TileCablingService *	m_cabling {nullptr}
	Pointer to TileCabling. More...
ServiceHandle< IAthRNGSvc >	m_rndmSvc { this, "RndmSvc", "AthRNGSvc", "Random Number Service used in TileCellBuildetFromHit" }
	< Random number service to use More...
ToolHandle< ITileBadChanTool >	m_tileBadChanTool
SG::ReadCondHandleKey< TileEMScale >	m_emScaleKey
	Name of TileEMScale in condition store. More...
ServiceHandle< TileCablingSvc >	m_cablingSvc
	Name of Tile cabling service. More...
const TileDetDescrManager *	m_tileMgr {nullptr}
	Pointer to TileDetDescrManager. More...
const MbtsDetDescrManager *	m_mbtsMgr {nullptr}
	Pointer to MbtsDetDescrManager. More...
TileFragHash::TYPE	m_RChType
	Type of TileRawChannels (Fit, OF2, etc.) More...
std::vector< CaloAffectedRegionInfo >	m_affectedRegionInfo_global
	TileDrawerRunStatus m_drawerRunStatus[5][64]; //!< overall status of drawer in whole run. More...
std::vector< CaloAffectedRegionInfo >	m_affectedRegionInfo_current_run
	method to process raw channels from a given vector and store them in collection More...
bool	m_RUN2
bool	m_RUN2plus
int	m_E1_TOWER

Static Private Attributes
static const int	NSIDE = 2
static const int	NPHI = 8
static const int	NETA = 2
static const int	NCELLMBTS = NSIDE * NPHI * NETA
static const int	E4SIDE = -1
static const int	E4ETA = 2
static const int	E4NPHI = 4
static const int	NCELLE4PR = E4NPHI

Detailed Description

This class creates Cells from RawChannels and stores them in a container.

Definition at line 80 of file TileCellBuilderFromHit.h.

Member Typedef Documentation

TileDrawerEvtStatusArray

typedef TileDrawerEvtStatus TileCellBuilderFromHit::TileDrawerEvtStatusArray[5][64]

private

status of every drawer

Definition at line 102 of file TileCellBuilderFromHit.h.

Constructor & Destructor Documentation

TileCellBuilderFromHit()

TileCellBuilderFromHit::TileCellBuilderFromHit	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Contructor.

Definition at line 71 of file TileCellBuilderFromHit.cxx.

  : base_class(type, name, parent)
  , m_eneForTimeCut(35. * MeV) // keep time only for cells above 70 MeV (more than 35 MeV in at least one PMT to be precise)
  , m_eneForTimeCutMBTS(0.03675) // the same cut for MBTS, but in pC, corresponds to 3 ADC counts or 35 MeV
  , m_qualityCut(254) // cut on overflow in quality (if quality is 255 - assume that channel is bad)
  , m_maxTime (25.)
  , m_minTime(-25.)
  , m_maskBadChannels(false)
  , m_noiseSigma(20.*MeV)
  , m_RChType(TileFragHash::Default)
  , m_RUN2(false)
  , m_RUN2plus(false)
  , m_E1_TOWER(10)
{
  declareInterface<TileCellBuilderFromHit>( this );
 
  //memset(m_drawerRunStatus, 0, sizeof(m_drawerRunStatus));
  //memset(m_eventErrorCounter, 0, sizeof(m_eventErrorCounter));
 
  // never set energy to zero, but set it to some small number
  // this will help TopoCluster to assign proper weight to the cell if needed
  m_zeroEnergy = 0.5 * MeV; // half a MeV in both PMTs i.e. one MeV in a cell
 
  // Noise Sigma
  declareProperty("NoiseSigma",m_noiseSigma);
 
  // Maximum and minimum time for a cell to be included:
  declareProperty("MaxTime", m_maxTime);
  declareProperty("MinTime", m_minTime);
 
  // put zero energy in bad channels and recover from single-channel failure using second PMT is a cell
  declareProperty("maskBadChannels", m_maskBadChannels);
 
  // PMT energy will be set to this value if channel is bad
  declareProperty("BadChannelZeroEnergy", m_zeroEnergy);
  // PMT with energy above cut will preserve time info in ESD
  declareProperty("EneForTimeCut", m_eneForTimeCut);
  declareProperty("EneForTimeCutMBTS", m_eneForTimeCutMBTS);
  // PMT with quality greater than this cut will be masked
  declareProperty("QualityCut", m_qualityCut);
}

~TileCellBuilderFromHit()

TileCellBuilderFromHit::~TileCellBuilderFromHit ( )

virtual

Destructor.

Definition at line 117 of file TileCellBuilderFromHit.cxx.

                                               { 
}

Member Function Documentation

build()

template<class ITERATOR , class COLLECTION >

void TileCellBuilderFromHit::build ( const CaloNoise *  caloNoise, TileDrawerEvtStatusArray &  drawerEvtStatus, const ITERATOR &  begin, const ITERATOR &  end, COLLECTION *  coll, TileCellContainer *  MBTSCells, TileCellContainer *  E4prCells, const TileSamplingFraction *  samplingFraction  ) const

private

< vector to of pointers to MBTS cells

< vector to of pointers to E4' cells

Definition at line 733 of file TileCellBuilderFromHit.cxx.

{
  static const std::string rngname = name() + "-" + ClassName<COLLECTION>::name();
  const EventContext& ctx = Gaudi::Hive::currentContext();
  ATHRNG::RNGWrapper* wrapper = m_rndmSvc->getEngine(this);
  wrapper->setSeed (rngname, ctx);
  CLHEP::HepRandomEngine* engine = wrapper->getEngine (ctx);
 
  SG::ReadCondHandle<TileEMScale> emScale(m_emScaleKey, ctx);
 
  /* zero all counters and sums */
  int nTwo = 0;
  int nCell = 0;
  int nFake = 0;
  int nMBTS = 0;
  int nE4pr = 0;
  int nChan = 0;
  float eCh = 0.0;
  float eCellTot = 0.0;
  float eMBTSTot = 0.0;
  float eE4prTot = 0.0;
  bool EBdrawerPresent[128];
  memset(EBdrawerPresent, 0, sizeof(EBdrawerPresent));
#ifdef USE_TILECELLS_DATAPOOL
  DataPool<TileCell> tileCellsP(5217);
#endif
  //**
  //* Iterate over hits, creating new TileCells (or incrementing
  //* existing ones). Add each new TileCell to the output collection
  //**
 
  bool overflow = false;
  bool underflow = false;
  bool overfit = false;
  float ener_min = (!m_caloNoiseKey.empty()) ? 1.0E-30F : 0.0;
 
  std::vector<TileCell*> allCells (m_tileID->cell_hash_max(), nullptr);
  std::vector<TileCell*> MBTSVec;   
  if (MBTSCells) {
    MBTSVec.resize (NCELLMBTS);
  }
  std::vector<TileCell*> E4prVec;   
  if (E4prCells) {
    E4prVec.resize (NCELLE4PR);
  }
 
  for (ITERATOR hitItr = begin; hitItr != end; ++hitItr) {
 
    const TileHit* pHit = (*hitItr);
 
    // Get hit id (logical pmt id)
    Identifier cell_id = pHit->identify();
 
    double ehit = 0.;
    double thit = 0.;
 
    int hitsize = pHit->size();
    for (int ind = 0; ind < hitsize; ++ind) { // loop over sub-hits 
      float tim = pHit->time(ind);
      if (m_minTime < tim && tim < m_maxTime) {
        float ene = pHit->energy(ind);
        ehit += ene;
        thit += tim*ene;
      }
    }
 
    float ener = ener_min; 
    float time = 0.;
    float qual = 0.;
    int gain = TileID::HIGHGAIN;
    bool good_time = false;
    bool non_zero_time = (ehit!=0.0);
 
    HWIdentifier channel_id = m_cabling->s2h_channel_id(cell_id);
    int ros = m_tileHWID->ros(channel_id);
    int drawer = m_tileHWID->drawer(channel_id);
    int drawerIdx = TileCalibUtils::getDrawerIdx(ros, drawer);
    int channel = m_tileHWID->channel(channel_id);
 
    if (non_zero_time) {
      time = thit / ehit;
      // Convert hit energy deposit to cell energy
      ener = ehit * samplingFraction->getSamplingFraction(drawerIdx, channel);
      if (ener > 10000. * MeV) { // very simple check - if above 10 GeV assume low gain
        gain = TileID::LOWGAIN;
      }
      good_time = (fabs(time) < 25.);
    }
    
    bool MBTS = false, E4pr = false, E1_CELL = false;
    int section,side,tower,sample,pmt,index;
 
    // Get logical ID of cell
    if (m_tileTBID->is_tiletb(cell_id)) {
 
      pmt = 0;
      side = m_tileTBID->type(cell_id);
      int phi = m_tileTBID->module(cell_id);
      int eta = m_tileTBID->channel(cell_id);
      if (eta<2) {
        MBTS = true;
        index = mbts_index(std::max(0,side),phi,eta);
      } else {
        E4pr = true;
        index = e4pr_index(phi);
      }
 
    } else {
 
      MBTS = false;
      E4pr = false;
 
      pmt = m_tileID->pmt(cell_id); // 0 or 1 - first or second PMT of the cell
 
      cell_id = m_tileID->cell_id(cell_id); // remove pmt field from cell_id
      index = m_tileID->cell_hash(cell_id);
 
      section = m_tileID->section(cell_id);
      side = m_tileID->side(cell_id);
      drawer = m_tileID->module(cell_id);
      tower = m_tileID->tower(cell_id);
      sample = m_tileID->sample(cell_id);
 
      // instead of using cabling (which might be wrong for upgrade geometry
      // try to figure out ros/drawer/channel 
      ros = ( TileID::POSITIVE == side ) ? 1 : 2;
      if ( TileID::BARREL != section ) ros += 2;
      else if (tower == 0 && sample == TileID::SAMP_D && pmt == 1) ros = 2; // Cell D0, second PMT in LBC
      channel = -1; // don't know channel
      E1_CELL = (section == TileID::GAPDET && sample == TileID::SAMP_E && tower == m_E1_TOWER);
    }
 
    ++nChan;
    eCh += ener;
 
    if (E4pr) { // E4' cells
 
      if (E4prCells) { // do something with them only if contaier existst
        ++nE4pr;
 
        eE4prTot += ener;
        unsigned char iqual = iquality(qual);
        // for E4' cell qbit use only non_zero_time flag and check that energy is above standatd energy threshold in MeV
        unsigned char qbit = qbits(drawerEvtStatus,
                                   ros, drawer, true, non_zero_time, (fabs(ener) > m_eneForTimeCut)
                                   , overflow, underflow, overfit);
        CaloGain::CaloGain cgain = (gain == TileID::HIGHGAIN)
                                   ? CaloGain::TILEONEHIGH
                                   : CaloGain::TILEONELOW;
 
        TileCell* pCell = NEWTILECELL();
        // no CaloDDE
        // Cell ID is set explicitly
        pCell->set(NULL, cell_id);
        pCell->setEnergy_nonvirt(ener, 0, cgain, 3);
        pCell->setTime_nonvirt(time);
        pCell->setQual1(iqual);
        pCell->setQual2(0);
        pCell->setQbit1(qbit);
        pCell->setQbit2(0);
 
        if (msgLvl(MSG::VERBOSE)) {
          msg(MSG::VERBOSE) << " E4'  cell_id=" << m_tileTBID->to_string(cell_id)
                            << " adc_id=" << "5/0/" << ros << "/" << drawer << "/" << channel << "/" << gain 
                            << " ene= " << ener
                            << " ehit= " << ehit
                            << " time= " << time
                            << " iqual= " << (int) iqual
                            << " qbit = 0x" << MSG::hex << (int) qbit << MSG::dec << endmsg;
        }
 
        if (E4prVec[index]) {
          msg(MSG::WARNING) << " double E4' cell_id=" << m_tileTBID->to_string(cell_id)
                            << "  ignoring previous value" << endmsg;
        }
        E4prVec[index] = pCell;
      }
 
    } else if (MBTS) { // MBTS cells
 
      if (MBTSCells) { // do something with them only if contaier existst
        ++nMBTS;
 
        // convert energy to pCb
        ener /= emScale->calibrateChannel(drawerIdx, channel, gain, 1.,
                                          TileRawChannelUnit::PicoCoulombs,
                                          TileRawChannelUnit::MegaElectronVolts);
        gain = (ener>12.) ? TileID::LOWGAIN : TileID::HIGHGAIN; // assume low gain above 12 pC
        eMBTSTot += ener;
        unsigned char iqual = iquality(qual);
        // for MBTS qbit use AND of good_time and non_zero_time and check that energy is above MBTS energy threshold in pC
        unsigned char qbit = qbits(drawerEvtStatus,
                                   ros, drawer, false, (good_time && non_zero_time)
           , (fabs(ener) > m_eneForTimeCutMBTS), overflow, underflow, overfit);
        CaloGain::CaloGain cgain = (gain == TileID::HIGHGAIN)
                                   ? CaloGain::TILEONEHIGH
                                   : CaloGain::TILEONELOW;
 
        TileCell* pCell = NEWTILECELL();
        // MBTS CaloDDE
        // Cell ID is set explicitly
        pCell->set((m_mbtsMgr) ? m_mbtsMgr->get_element(cell_id) : NULL, cell_id);
        pCell->setEnergy_nonvirt(ener, 0, cgain, 3);
        pCell->setTime_nonvirt(time);
        pCell->setQual1(iqual);
        pCell->setQual2(0);
        pCell->setQbit1(qbit);
        pCell->setQbit2(0);
 
        if (msgLvl(MSG::VERBOSE)) {
          msg(MSG::VERBOSE) << " MBTS cell_id=" << m_tileTBID->to_string(cell_id)
                            << " adc_id=" << "5/0/" << ros << "/" << drawer << "/" << channel << "/" << gain 
                            << " ene= " << ener
                            << " ehit= " << ehit
                            << " time= " << time
                            << " iqual= " << (int) iqual
                            << " qbit = 0x" << MSG::hex << (int) qbit << MSG::dec << endmsg;
        }
 
        if (MBTSVec[index]) {
          msg(MSG::WARNING) << " double MBTS cell_id=" << m_tileTBID->to_string(cell_id)
                            << "  ignoring previous value" << endmsg;
        }
        MBTSVec[index] = pCell;
      }
    } else {
 
      eCellTot += ener;
 
      unsigned char iqual = iquality(qual);
      // for normal cell qbit use only non_zero_time flag and check that energy is above standard energy threshold in MeV
     unsigned char qbit = qbits(drawerEvtStatus, ros, drawer, true, non_zero_time,
                                (fabs(ener) > m_eneForTimeCut), overflow, underflow, overfit);
 
      if (E1_CELL && m_RUN2plus) {
 
        int drawer2 = m_cabling->E1_merged_with_run2plus(ros,drawer);
        if (drawer2 != 0) { // Raw channel split into two E1 cells for Run 2.
          Identifier cell_id2 = m_tileID->cell_id(TileID::GAPDET, side, drawer2, m_E1_TOWER, TileID::SAMP_E);
          int index2 = m_tileID->cell_hash(cell_id2);
          TileCell* pCell2 = NEWTILECELL();
          ++nCell;
          allCells[index2] = pCell2;
          const CaloDetDescrElement* dde2 = m_tileMgr->get_cell_element(index2);
          pCell2->set(dde2, cell_id2);
          pCell2->setEnergy_nonvirt(0.0, 0.0, CaloGain::INVALIDGAIN, CaloGain::INVALIDGAIN);
          int pmt2(0);
          ener *= 0.5F;
          correctCell(pCell2, 1, pmt2, gain, ener, time, iqual, qbit);
 
          ATH_MSG_DEBUG("E1 cell Id => " << m_tileID->to_string(cell_id)
                       << " split into " << m_tileID->to_string(cell_id2));
        }
      }
 
      TileCell* pCell = allCells[index];
      if (pCell) {
        ++nTwo;
        correctCell(pCell, 2, pmt, gain, ener, time, iqual, qbit); // correct & merge 2 PMTs in one cell
      } else {
        ++nCell;
        allCells[index] = pCell = NEWTILECELL();
        const CaloDetDescrElement* dde = m_tileMgr->get_cell_element(index);
        pCell->set(dde, cell_id);
        pCell->setEnergy_nonvirt(0, 0, CaloGain::INVALIDGAIN, CaloGain::INVALIDGAIN);
        correctCell(pCell, 1, pmt, gain, ener, time, iqual, qbit); // correct & save e,t,q in new cell
      }
 
      if (msgLvl(MSG::VERBOSE)) {
        float calib1 = (ehit != 0.0) ? ener / ehit : 0.0;
        msg(MSG::VERBOSE) << " cell_id=" << m_tileID->to_string(cell_id, -2)
                          << " adc_id=" << "5/0/" << ros << "/" << drawer << "/" << channel << "/" << gain 
                          << " calib=" << calib1
                          << " nCell=" << nCell
                          << " energy=" << ener << " (" << pCell->energy() << ", " << pCell->eneDiff() << ")" << endmsg;
 
        msg(MSG::VERBOSE) << " ehit= " << ehit
                          << " time= " << time
                          << " iqual= " << (int) iqual
                          << " qbit = 0x" << MSG::hex << (int) qbit << MSG::dec << endmsg;
      }
    }
  }
 
 
 
  unsigned char iqual = 0;
  unsigned char qbit = m_RChType;
 
  //**
  // Now store all TileCells
  //**
  for (unsigned int index = 0; index < allCells.size(); ++index) {
 
    TileCell * pCell = allCells[index];
    const CaloDetDescrElement* dde = m_tileMgr->get_cell_element(index);
    if(!dde)
    {
      ATH_MSG_VERBOSE ("skipping dummy tower with id=" << m_tileID->to_string(m_tileID->cell_id((IdentifierHash)index), -2) << " hash " << index);
      continue;
    }
    Identifier cell_id = dde->identify();
 
    int section = m_tileID->section(cell_id);
    int side = m_tileID->side(cell_id);
    int module = m_tileID->module(cell_id);
    int sample = m_tileID->sample(cell_id);
 
    bool single_PMT_C10 = ( section == TileID::GAPDET && sample == TileID::SAMP_B && !m_cabling->C10_connected(module) );
    bool missing_D4 = ( section == TileID::GAPDET && sample == TileID::SAMP_D && (module == ((side>0)?14:17)) );
    bool Ecell = (sample == TileID::SAMP_E);
    bool single_PMT = Ecell || single_PMT_C10;
    
    if (!pCell) {
 
      if ( ! missing_D4 ) {
 
        ++nCell;
        if (!single_PMT) ++nTwo;
 
        allCells[index] = pCell = NEWTILECELL();
        pCell->set(dde, cell_id);
 
        pCell->setEnergy_nonvirt(ener_min, 0.0, TileID::HIGHGAIN, (Ecell)?3:TileID::HIGHGAIN); // reset energy completely
        pCell->setTime_nonvirt(0.0); // reset time completely
        pCell->setQuality(iqual, qbit, 0); // reset quality flag for first pmt
        pCell->setQuality(iqual, qbit, 1); // reset quality flag for second pmt
 
        ATH_MSG_VERBOSE ( "adding cell with no signal with id=" << m_tileID->to_string(pCell->ID(), -2));
      }
    }
    
    if (pCell) {      // cell exists
 
      if (m_maskBadChannels)
        if (maskBadChannels (drawerEvtStatus, pCell,single_PMT_C10,Ecell))
          ATH_MSG_VERBOSE ( "cell with id=" << m_tileID->to_string(pCell->ID(), -2)
                           << " bad channels masked, new energy=" << pCell->energy() );
 
      bool bad1=pCell->badch1();
      bool bad2=pCell->badch2();
      
      if (! (bad1 && bad2) ) {
 
        float noiseSigma = m_noiseSigma;
        if (caloNoise) {
          noiseSigma = caloNoise->getEffectiveSigma (pCell->ID(), pCell->gain(), pCell->energy());
        }
 
        if (bad1 || bad2 || single_PMT) {
            
          float ene = RandGaussQ::shoot(engine, 0.0, noiseSigma);
 
          ATH_MSG_VERBOSE ( "Cell noise for cell id=" << m_tileID->to_string(pCell->ID(), -2)
                            << " sigma " << noiseSigma << " noise " << ene );
          if (Ecell)
            pCell->addEnergy(ene,0,pCell->gain1());
          else
            pCell->addEnergy(ene);
          
        } else {
 
          int gain1=pCell->gain1();
          int gain2=pCell->gain2();
            
          if (gain1!=gain2) {
 
            float ene = RandGaussQ::shoot(engine, 0.0, noiseSigma);
            ATH_MSG_VERBOSE ( "Cell noise for cell id=" << m_tileID->to_string(pCell->ID(), -2)
                              << " sigma " << noiseSigma << " noise " << ene );
 
            if (gain1==TileID::LOWGAIN)
              pCell->addEnergy(ene,0,gain1);
            else
              pCell->addEnergy(ene,1,gain2);
 
          } else {
 
            noiseSigma *= M_SQRT1_2;
            float ene1 = RandGaussQ::shoot(engine, 0.0, noiseSigma);
            float ene2 = RandGaussQ::shoot(engine, 0.0, noiseSigma);
 
            ATH_MSG_VERBOSE ( "Cell noise for cell id=" << m_tileID->to_string(pCell->ID(), -2)
                              << " sigma " << noiseSigma*M_SQRT2 << " noise " << ene1+ene2 << " noise1 " << ene1 << " noise2 " << ene2 );
 
            pCell->addEnergy(ene1,0,gain1);
            pCell->addEnergy(ene2,0,gain2);
          }
        }
      }
      coll->push_back(pCell); // store cell in container
      allCells[index] = 0; // clear pointer for next event
    }
  }
 
 
  if (MBTSCells) {
 
    for (int side = 0; side < NSIDE; ++side) {
      for (int phi = 0; phi < NPHI; ++phi) {
        for (int eta = 0; eta < NETA; ++eta) {
 
          int index=mbts_index(side,phi,eta);
          TileCell * pCell = MBTSVec[index];
 
          bool merged_MBTS = ( eta == 1 && (phi&1) == 1 && m_RUN2); // in RUN2 every second outer MBTS does not exist
 
          if (!pCell && !merged_MBTS) {
 
            ++nMBTS;
            MBTSVec[index] = pCell = NEWTILECELL();
 
            Identifier cell_id = m_tileTBID->channel_id((side > 0) ? 1 : -1, phi, eta);
            pCell->set((m_mbtsMgr) ? m_mbtsMgr->get_element(cell_id) : NULL, cell_id);
            pCell->setEnergy_nonvirt(0.0, 0.0, CaloGain::TILEONEHIGH, CaloGain::INVALIDGAIN); // reset energy completely
            pCell->setTime_nonvirt(0.0); // reset time completely
            pCell->setQuality(iqual, qbit, 0); // reset quality flag for first pmt
            pCell->setQuality(iqual, qbit, 1); // reset quality flag for second pmt
          }
          
          if (pCell) {
            if (m_maskBadChannels && maskBadChannel (drawerEvtStatus, pCell))
                ATH_MSG_VERBOSE ( "MBTS cell with id=" << m_tileTBID->to_string(pCell->ID())
                                  << " bad channel masked, new energy=" << pCell->energy() );
 
            MBTSCells->push_back(pCell); // store cell in container
            MBTSVec[index] = 0; // clear pointer for next event
          }
        }
      }
    }
  }
 
  if (E4prCells) {
 
    for (int phi = 0; phi < E4NPHI; ++phi) {
 
      int index = e4pr_index(phi);
      TileCell * pCell = E4prVec[index];
 
      if (!pCell) {
 
        ++nE4pr;
        E4prVec[index] = pCell = NEWTILECELL();
 
        pCell->set(NULL, m_tileTBID->channel_id(E4SIDE, phi, E4ETA));
        pCell->setEnergy_nonvirt(0.0, 0.0, CaloGain::TILEONEHIGH, CaloGain::INVALIDGAIN); // reset energy completely
        pCell->setTime_nonvirt(0.0); // reset time completely
        pCell->setQuality(iqual, qbit, 0); // reset quality flag for first pmt
        pCell->setQuality(iqual, qbit, 1); // reset quality flag for second pmt
      }
 
      if (pCell) {
        if (m_maskBadChannels && maskBadChannel (drawerEvtStatus, pCell))
          ATH_MSG_VERBOSE ( "E4pr cell with id=" << m_tileTBID->to_string(pCell->ID())
                             << " bad channel masked, new energy=" << pCell->energy() );
 
        E4prCells->push_back(pCell); // store cell in container
        E4prVec[index] = 0; // clear pointer for next event
      }
    }
  }
  
  if (msgLvl(MSG::DEBUG)) {
    msg(MSG::DEBUG) << " nChan=" << nChan
                    << " RawChSum=" << eCh
                    << " nCell=" << nCell
                    << " n2=" << nTwo
                    << " nFake=" << nFake
                    << " eneTot=" << eCellTot;
 
    if (MBTSCells)
      msg(MSG::DEBUG) << " nMBTS=" << nMBTS
                      << " eMBTS=" << eMBTSTot;
    if (E4prCells)
      msg(MSG::DEBUG) << " nE4pr=" << nE4pr
                      << " eE4pr=" << eE4prTot;
 
    msg(MSG::DEBUG) << endmsg;
  }
}

correctCell()

void TileCellBuilderFromHit::correctCell ( TileCell *  pCell, int  correction, int  pmt, int  gain, float  ener, float  time, unsigned char  iqual, unsigned char  qbit  ) const

private

Compute calibrated energy, time, etc.

for TileCell and adjust it.

Definition at line 409 of file TileCellBuilderFromHit.cxx.

                                                                                                                  {
//************************************************************************
 
// Merge two pmts in one cell if needed
// and apply corrections
 
  switch (correction) {
    case 1: // first pmt for this cell
      pCell->addEnergy(ener, pmt, gain);
      pCell->setTime(time); // overwrite time completely
      pCell->setQuality(iqual, qbit, pmt);
      pCell->setQuality(0, 0, 1 - pmt);
      break;
    case 2: // second pmt for this cell
      pCell->addEnergy(ener, pmt, gain);
      pCell->setTime(time, pmt); // calculate average time and timeDiff
      pCell->setQuality(iqual, qbit, pmt);
      break;
  }
}

e4pr_index()

int TileCellBuilderFromHit::e4pr_index ( int  phi ) const

inlineprivate

Definition at line 229 of file TileCellBuilderFromHit.h.

                                         {
      return  phi;
    }

finalize()

StatusCode TileCellBuilderFromHit::finalize ( )

overridevirtual

Definition at line 190 of file TileCellBuilderFromHit.cxx.

                                            {
 
  ATH_MSG_INFO( "Finalizing" );
 
  return StatusCode::SUCCESS;
}

initialize()

StatusCode TileCellBuilderFromHit::initialize ( )

overridevirtual

Initializer.

Definition at line 123 of file TileCellBuilderFromHit.cxx.

                                              {
  
  // retrieve MBTS and Tile detector manager, TileID helper and TileIfno from det store
  // retrieve MBTS and Tile detector manager, TileID helper and TileIfno from det store
  if (m_MBTSContainerKey.empty()) {
    m_mbtsMgr = nullptr;
  } else {
 
    ATH_CHECK( m_MBTSContainerKey.initialize() );
    ATH_MSG_INFO( "Storing MBTS cells in " << m_MBTSContainerKey.key() );
    
    if (detStore()->retrieve(m_mbtsMgr).isFailure()) {
      ATH_MSG_WARNING( "Unable to retrieve MbtsDetDescrManager from DetectorStore" );
      m_mbtsMgr = nullptr;
    }
  }
 
  ATH_CHECK( m_eventInfoKey.initialize() );
  ATH_CHECK( m_hitContainerKey.initialize() );
  ATH_CHECK( m_samplingFractionKey.initialize() );
  ATH_CHECK( m_emScaleKey.initialize() );
 
  ATH_CHECK( detStore()->retrieve(m_tileMgr) );
  ATH_CHECK( detStore()->retrieve(m_tileID) );
  ATH_CHECK( detStore()->retrieve(m_tileTBID) );
  ATH_CHECK( detStore()->retrieve(m_tileHWID) );
 
  //=== get TileBadChanTool
  ATH_CHECK( m_tileBadChanTool.retrieve() );
 
  //---- retrieve the noise ----------------
  if (!m_caloNoiseKey.empty()) {
    ATH_CHECK( m_caloNoiseKey.initialize() );
    ATH_MSG_INFO( "Reading electronic noise from DB" );
  } else {
    ATH_MSG_INFO( "Noise Sigma " << m_noiseSigma << " MeV is selected!" );
  }
 
  ATH_CHECK( m_rndmSvc.retrieve());
  (void) m_rndmSvc->getEngine(this); // get this created early.
 
  ATH_MSG_INFO( "max time thr  " << m_maxTime << " ns" );
  ATH_MSG_INFO( "min time thr  " << m_minTime << " ns" );
  
  m_E1_TOWER = (m_tileID->cell_hash_max() < 10000) ? 10 : 40;
 
  ATH_MSG_INFO( "taking hits from '" << m_hitContainerKey.key() << "'" );
 
  ATH_CHECK( m_cablingSvc.retrieve() );
  m_cabling = m_cablingSvc->cablingService();
 
  m_RUN2 = (m_cabling->getCablingType() == TileCablingService::RUN2Cabling);
  m_RUN2plus = m_cabling->isRun2PlusCabling();
 
  if (m_RUN2 && !m_E4prContainerKey.empty()) {
    ATH_CHECK( m_E4prContainerKey.initialize() );
    ATH_MSG_INFO( "Storing E4'  cells in " << m_E4prContainerKey.key() );
  } else {
    m_E4prContainerKey = ""; // no E4' container for RUN1
  }
 
 
  ATH_MSG_INFO( "TileCellBuilderFromHit initialization completed" );
 
  return StatusCode::SUCCESS;
}

interfaceID()

const InterfaceID & TileCellBuilderFromHit::interfaceID ( )

static

Definition at line 66 of file TileCellBuilderFromHit.cxx.

                                                        {
  return IID_ITileCellBuilderFromHit;    
}

iquality()

unsigned char TileCellBuilderFromHit::iquality ( float  qual ) const

inlineprivate

< method to compute the cell quality

Definition at line 206 of file TileCellBuilderFromHit.h.

                                              {
         return std::min(255, abs((int) qual));
    } // keep quality within 8 bits make it "unsigned char"

maskBadChannel()

bool TileCellBuilderFromHit::maskBadChannel ( TileDrawerEvtStatusArray &  drawerEvtStatus, TileCell *  pCell  ) const

private

method to check if channels are good or bad.

Puts zero if both channels are bad or recovers from single-channel failure. It returns true if cell was changed, false otherwise

Definition at line 464 of file TileCellBuilderFromHit.cxx.

{
  Identifier cell_id = pCell->ID();
 
  HWIdentifier channel_id = m_cabling->s2h_channel_id(cell_id);
  int ros = m_tileHWID->ros(channel_id);
  int drawer = m_tileHWID->drawer(channel_id);
  int channel = m_tileHWID->channel(channel_id);
  int drawerIdx = TileCalibUtils::getDrawerIdx(ros, drawer);
  int gain = pCell->gain1();
 
  TileBchStatus chStatus = m_tileBadChanTool->getAdcStatus(drawerIdx, channel, gain);
 
  // check quality first
  bool bad = ((int) pCell->qual1() > m_qualityCut);
  if (bad) {
    ++drawerEvtStatus[ros][drawer].nBadQuality;
 
  } else {
    // check bad status in DB
    bad = chStatus.isBad();
 
  }
 
  if (bad) {
    // only one channel in this cell and it is bad
    ++drawerEvtStatus[ros][drawer].nMaskedChannels;
 
    //pCell->setEnergy(m_zeroEnergy,0.0,TileID::LOWGAIN,CaloGain::INVALIDGAIN); // reset energy completely, indicate problem putting low gain
    //pCell->setTime(0.0); // reset time completely
    //pCell->setQuality(255,TileCell::MASK_BADCH,0); // reset quality flag for first pmt
 
    if (gain == CaloGain::INVALIDGAIN) {
      pCell->setEnergy(0.0, 0.0, TileID::LOWGAIN, CaloGain::INVALIDGAIN); // reset energy completely, indicate problem putting low gain
    } else {
      pCell->setEnergy(0.0, 0.0); // reset energy completely without changing the gain
    }
    pCell->setTime(-100.0); // reset time to big negative number to distinguish this bad cell from good cells
    pCell->setQuality(255, (TileCell::MASK_BADCH | (pCell->qbit1() & TileCell::MASK_ALGO)), 0); // reset quality flag for first pmt
    pCell->setQuality(0, TileCell::MASK_BADCH, 1); // reset quality flag for second pmt
 
    return true;
 
  } else if (chStatus.isBadTiming()) {
    pCell->setTime(0.0); // channel with bad timing - no cell time
    uint8_t qbit1 = pCell->qbit1() & (~(TileCell::MASK_TIME)); // clear time bit for first pmt
    pCell->setQuality(pCell->qual1(), qbit1, 0); // update qbits for first pmt
  }
 
  return false;
}

maskBadChannels()

bool TileCellBuilderFromHit::maskBadChannels ( TileDrawerEvtStatusArray &  drawerEvtStatus, TileCell *  pCell, bool  single_PMT_C10, bool  Ecell  ) const

private

Definition at line 519 of file TileCellBuilderFromHit.cxx.

                                                                                                 {
 
  const CaloDetDescrElement* caloDDE = pCell->caloDDE();
 
  IdentifierHash hash1 = caloDDE->onl1();
  IdentifierHash hash2 = caloDDE->onl2();
 
  int gain1 = pCell->gain1();
 
  HWIdentifier ch_id1 = m_tileHWID->channel_id(hash1);
 
  int ros1 = m_tileHWID->ros(ch_id1);
  int drawer1 = m_tileHWID->drawer(ch_id1);
  int chan1 = m_tileHWID->channel(ch_id1);
  int drawerIdx1 = TileCalibUtils::getDrawerIdx(ros1, drawer1);
  const TileBchStatus& chStatus1 = m_tileBadChanTool->getAdcStatus(drawerIdx1, chan1, (gain1 < 0) ? 1 : gain1);
  
  // check quality first
  bool bad1 = ((int) pCell->qual1() > m_qualityCut);
  if (bad1) {
    ++drawerEvtStatus[ros1][drawer1].nBadQuality;
 
  } else {
    // check bad status in DB
    bad1 = (gain1 < 0) || chStatus1.isBad();
 
  }
 
  if ( Ecell || hash2 == TileHWID::NOT_VALID_HASH) {
    // gap/crack scintillators with one PMT per cell
 
    if (bad1) {
      // only one channel in this cell and it is bad
      ++drawerEvtStatus[ros1][drawer1].nMaskedChannels;
 
      if (gain1 == CaloGain::INVALIDGAIN) {
        pCell->setEnergy(m_zeroEnergy, 0.0, TileID::LOWGAIN, CaloGain::INVALIDGAIN); // reset energy completely, indicate problem putting low gain
      } else {
        pCell->setEnergy(m_zeroEnergy, 0.0); // reset energy completely without changing gain
      }
      pCell->setTime(0.0); // reset time completely
      pCell->setQuality(255, (TileCell::MASK_BADCH | (pCell->qbit1() & TileCell::MASK_ALGO)), 0); // reset quality flag for first pmt
      pCell->setQuality(0, TileCell::MASK_BADCH, 1); // reset quality flag for second pmt
 
      return true;
 
    } else if (chStatus1.isBadTiming()) {
      pCell->setTime(0.0); // channel with bad timing - no cell time
      uint8_t qbit1 = pCell->qbit1() & (~(TileCell::MASK_TIME)); // clear time bit for first pmt
      pCell->setQuality(pCell->qual1(), qbit1, 0); // update qbits for first pmt
    }
    
  } else { //cell has both PMTs
 
    int gain2 = pCell->gain2();
    
    HWIdentifier ch_id2 = m_tileHWID->channel_id(hash2);
 
    int ros2 = m_tileHWID->ros(ch_id2);
    int drawer2 = m_tileHWID->drawer(ch_id2);
    int chan2 = m_tileHWID->channel(ch_id2);
    int drawerIdx2 = TileCalibUtils::getDrawerIdx(ros2, drawer2);
    const TileBchStatus& chStatus2 = m_tileBadChanTool->getAdcStatus(drawerIdx2, chan2, (gain2 < 0) ? 1 : gain2);
 
    // check quality first
    bool bad2 = ((int) pCell->qual2() > m_qualityCut);
    if (bad2) {
      ++drawerEvtStatus[ros2][drawer2].nBadQuality;
 
    } else {
      // check bad status in DB
      bad2 = (gain2 < 0) || chStatus2.isBad();
 
    }
 
    if (single_PMT_C10) {
      // for special C10 disconnected channel might be masked in DB
      // and energy of good channel is taken twice with correct weight
      // but if this channel is not masked in DB - set its bad status
      // equal to bad status of real channel, so that cell is masked correctly
      // if real channel connected to a cell is bad
#ifdef ALLOW_DEBUG_COUT
      static int cnt=0;
      if (++cnt < 17) {
        std::cout << "special C10 in " << ((ros2==TileHWID::EXTBAR_POS) ? "EBA" : "EBC")
        << drawer2+1 << " status " << chan1 << "/" << chan2 << " "
        << (chStatus1.isBad()?"bad":"good") << "/"
        << (chStatus2.isBad()?"bad":"good") << "/"
        << ((m_RUN2plus)?" RUN2+ cabling": "RUN1 cabling")
        << std::endl;
      }
#endif
      if (chan1 == 4) {
        if (m_RUN2plus || !chStatus1.isBad()) {
#ifdef ALLOW_DEBUG_COUT
          if (cnt < 17) {
            std::cout << "Ene of chan1 was " << pCell->ene1() << " changing to half of " << pCell->ene2()
            << " and setting bad1=true" << std::endl;
          }
#endif
          pCell->setEnergy(pCell->ene2()/2., pCell->ene2()/2., gain2, gain2);
          //bad1 = bad2;
          bad1 = true;
          --drawerEvtStatus[ros1][drawer1].nMaskedChannels; // since it's fake masking, decrease counter by 1 in advance
        }
      } else {
        if (m_RUN2plus || !chStatus2.isBad()) {
#ifdef ALLOW_DEBUG_COUT
          if (cnt < 17) {
            std::cout << "Ene of chan2 was " << pCell->ene2() << " changing to half of " << pCell->ene1()
            << " and setting bad2=true" << std::endl;
          }
#endif
          pCell->setEnergy(pCell->ene1()/2., pCell->ene1()/2., gain1, gain1);
          //bad2 = bad1;
          bad2 = true;
          --drawerEvtStatus[ros2][drawer2].nMaskedChannels; // since it's fake masking, decrease counter by 1 in advance
        }
      }
    }
    if (bad1 && bad2) {
      // both channels are bad
      ++drawerEvtStatus[ros1][drawer1].nMaskedChannels;
      ++drawerEvtStatus[ros2][drawer2].nMaskedChannels;
 
      if (gain1 == CaloGain::INVALIDGAIN || gain2 == CaloGain::INVALIDGAIN) {
        if (gain1 == CaloGain::INVALIDGAIN) gain1 = 0; // this is TileID::LOWGAIN; - commented out to make Coverity happy
        if (gain2 == CaloGain::INVALIDGAIN) gain2 = 0; // this is TileID::LOWGAIN; - commented out to make Coverity happy
        pCell->setEnergy(m_zeroEnergy, m_zeroEnergy, gain1, gain2); // reset energy completely, indicate problem putting low gain
      } else {
        pCell->setEnergy(m_zeroEnergy, m_zeroEnergy); // reset energy completely without changing gain
      }
      pCell->setTime(0.0);  // reset time completely
      pCell->setQuality(255, (TileCell::MASK_BADCH | (pCell->qbit1() & TileCell::MASK_ALGO)), 0); // reset quality flag for first pmt
      pCell->setQuality(255, (TileCell::MASK_BADCH | (pCell->qbit2() & TileCell::MASK_ALGO)), 1); // reset quality flag for second pmt
 
      return true;
 
    } else if (bad1 && !bad2) {
      // first channel is bad
      ++drawerEvtStatus[ros1][drawer1].nMaskedChannels;
 
      float ene2 = pCell->ene2();
      pCell->setEnergy(ene2, ene2, gain2, gain2); // use energy/gain from second pmt for both pmts
 
      if (chStatus2.isBadTiming()) {
        pCell->setTime(0.0); // time in second pmt is bad - no cell time
        uint8_t qual2 = pCell->qual2();
        uint8_t qbit2 = pCell->qbit2() & (~(TileCell::MASK_TIME)); // clear time bit for second pmt
        uint8_t qbit1 = qbit2 | TileCell::MASK_BADCH; // set bad channel bit for first pmt
        uint8_t qual1 = qual2 + (gain1 - gain2); // if gains are different, qua11 and qual2 will be different
        if (qual1 > m_qualityCut && gain1 > gain2) qual1 = qual2 - (gain1 - gain2); // new feature in release 17.2
        pCell->setQuality(qual1, qbit1, 0); // change quality and qbits for first pmt
        pCell->setQuality(qual2, qbit2, 1); // update qbits for second pmt
      } else {
        pCell->setTime(pCell->time2());  // use time from second pmt as cell time
        pCell->setQuality(pCell->qual2(), (pCell->qbit2() | TileCell::MASK_BADCH), 0); // change quality flag for first pmt
      }
      
      return true;
 
    } else if (!bad1 && bad2) {
      // second channel is bad
      ++drawerEvtStatus[ros2][drawer2].nMaskedChannels;
 
      float ene1 = pCell->ene1();
      pCell->setEnergy(ene1, ene1, gain1, gain1);  // use energy/gain from first pmt for both pmts
 
      if (chStatus1.isBadTiming()) {
        pCell->setTime(0.0); // time in first pmt is bad - no cell time
        uint8_t qual1 = pCell->qual1();
        uint8_t qbit1 = pCell->qbit1() & (~(TileCell::MASK_TIME)); // clear time bit for first pmt
        uint8_t qbit2 = qbit1 | TileCell::MASK_BADCH; // set bad channel bit for second pmt
        uint8_t qual2 = qual1 + (gain2 - gain1); // if gains are different, qua11 and qual2 will be different
        if (qual2 > m_qualityCut && gain2 > gain1) qual2 = qual1 - (gain2 - gain1); // new feature in release 17.2
        pCell->setQuality(qual1, qbit1, 0); // update qbits for first pmt
        pCell->setQuality(qual2, qbit2, 1); // change quality and qbits for second pmt
      } else {
        pCell->setTime(pCell->time1());  // use time from first pmt as cell time
        pCell->setQuality(pCell->qual1(), (pCell->qbit1() | TileCell::MASK_BADCH), 1); // change quality flag for second pmt
      }
      
      return true;
 
    } else {
 
      if (chStatus1.isBadTiming()) {
 
        if (chStatus2.isBadTiming()) {
          pCell->setTime(0.0); // time in both pmts is bad - no cell time
          uint8_t qbit2 = pCell->qbit2() & (~(TileCell::MASK_TIME)); // clear time bit for second pmt
          pCell->setQuality(pCell->qual2(), qbit2, 1); // update qbits for second pmt
        } else {
          pCell->setTime(pCell->time2());  // use time from second pmt as cell time
        }
        uint8_t qbit1 = pCell->qbit1() & (~(TileCell::MASK_TIME)); // clear time bit for first pmt
        pCell->setQuality(pCell->qual1(), qbit1, 0); // update qbits for first pmt
 
      } else if (chStatus2.isBadTiming()) {
 
        pCell->setTime(pCell->time1());  // use time from first pmt as cell time
        uint8_t qbit2 = pCell->qbit2() & (~(TileCell::MASK_TIME)); // clear time bit for second pmt
        pCell->setQuality(pCell->qual2(), qbit2, 1); // update qbits for second pmt
      }
    }
    
  }
  
  return single_PMT_C10;
}

mbts_index()

int TileCellBuilderFromHit::mbts_index ( int  side, int  phi, int  eta  ) const

inlineprivate

Definition at line 222 of file TileCellBuilderFromHit.h.

                                                            {
      return (side * NPHI + phi) * NETA + eta;
    }

process()

StatusCode TileCellBuilderFromHit::process ( CaloCellContainer *  theCellContainer, const EventContext &  ctx  ) const

overridevirtual

method to process all raw channels and store them in container

Definition at line 197 of file TileCellBuilderFromHit.cxx.

{
  //**
  //* Get TileHits
  //**
 
  TileDrawerEvtStatusArray drawerEvtStatus;
 
  SG::ReadCondHandle<TileSamplingFraction> samplingFraction(m_samplingFractionKey, ctx);
  ATH_CHECK( samplingFraction.isValid() );
 
  SG::ReadHandle<TileHitContainer> hitContainer(m_hitContainerKey, ctx);
 
  if (!hitContainer.isValid()) {
 
    ATH_MSG_WARNING( " Could not find container " << m_hitContainerKey.key() );
    ATH_MSG_WARNING( " do not fill CaloCellContainer " );
 
  } else {
    
    ATH_MSG_DEBUG( "Container " << m_hitContainerKey.key() << " with TileHits found ");
 
    std::unique_ptr<TileCellContainer> MBTSCells;
    if (!m_MBTSContainerKey.empty()) {
      MBTSCells = std::make_unique<TileCellContainer>(SG::VIEW_ELEMENTS);
    }
 
    std::unique_ptr<TileCellContainer> E4prCells;
    if (!m_E4prContainerKey.empty()) {
      E4prCells = std::make_unique<TileCellContainer>(SG::VIEW_ELEMENTS);
    }
    
    SelectAllObject<TileHitContainer> selAll(hitContainer.cptr());
    SelectAllObject<TileHitContainer>::const_iterator begin = selAll.begin();
    SelectAllObject<TileHitContainer>::const_iterator end = selAll.end();
 
    const CaloNoise* caloNoise = nullptr;
    if (!m_caloNoiseKey.empty()) {
      SG::ReadCondHandle<CaloNoise> noiseH (m_caloNoiseKey, ctx);
      caloNoise = noiseH.cptr();
    }
 
    if (begin != end) {
      ATH_MSG_DEBUG( " Calling build() method for hits from " << m_hitContainerKey.key() );
      build (caloNoise, drawerEvtStatus, begin, end, theCellContainer,
             MBTSCells.get(), E4prCells.get(), *samplingFraction);
    }
    
    if (!m_MBTSContainerKey.empty()) {
      SG::WriteHandle<TileCellContainer> MBTSContainer(m_MBTSContainerKey, ctx);
      ATH_CHECK( MBTSContainer.record(std::move(MBTSCells)) );
    }
    
    if (!m_E4prContainerKey.empty()) {
      SG::WriteHandle<TileCellContainer> E4prContainer(m_E4prContainerKey, ctx);
      ATH_CHECK( E4prContainer.record(std::move(E4prCells)) );
    }
 
    CaloCell_ID::SUBCALO caloNum = CaloCell_ID::TILE;
    //specify that a given calorimeter has been filled
    if (theCellContainer->hasCalo(caloNum)) {
      // log << MSG::WARNING << "CaloCellContainer has already been filled with TileCells (caloNum = " 
      //  <<    caloNum << ")" << endreq ;    
    }
    theCellContainer->setHasCalo(caloNum);
  }
  
  //enum EventFlagErrorState { NotSet, Warning, Error };
  xAOD::EventInfo::EventFlagErrorState error = xAOD::EventInfo::NotSet;
  // flag will contain status of a given event
  // every 4 bits - status of partitions LBA,LBC,EBA,EBC
  // bits 0-3   - there is a signal above threshold in partitions
  // bits 4-7   - there are channels with underflow (sample=0) in partition (since rel 17.2.6.4)
  // bits 8-11  - there are channels with overflow (sample=m_tileInfo->ADCmax()) in partition (since rel 17.2.6.4)
  // bits 12-15 - there are at least 16 drawers with bad quality in partition
  // bits 16-19 - maximal length of consecutive bad area (since rel 17.2.6.5)
  // bits 20-23 - there are at least 16 drawers which are completely masked in partition
  // bits 24-27 - there are at least 16 drawers which do not send data in partition
  // bits 28-31 - reserved for global good/warning/bad status
  // bits 20-27 are also used for module number which gives warning status (since release 17.2.6.5)
  //            in case of warning we are sure that bits which indicates error are not filled
  unsigned int flag = 0;
  
  int drConsecMaxMax = 0;
  int drConsecNum = 0;
 
  for (int p = 1; p < 5; ++p) {
    TileDrawerEvtStatus * evt = drawerEvtStatus[p];
    //TileDrawerRunStatus * run = m_drawerRunStatus[p];
    int drAbsent = 0;
    int drMasked = 0;
    int drConsec = 0;
    int drConsecMax = 0;
    int hasBadQ = 0;
    int hasOver = 0;
    int hasUnder = 0;
    int hasSig = 0;
    for (int d = 0; d < 64; ++d) {
      if (evt[d].nChannels == 0) {
        ++drConsec;
        ++drAbsent;
        //++(run[d].drawerAbsent);
      } else if (evt[d].nMaskedChannels >= evt[d].nChannels) {
        ++drConsec;
        ++drMasked;
        //++(run[d].drawerMasked);
      } else {
        if (drConsec > drConsecMax) {
          drConsecMax = drConsec;
          if (drConsecMax > drConsecMaxMax) {
            drConsecMaxMax = drConsecMax;
            drConsecNum = ((p - 1) << 6) | (d - drConsec);
          }
        }
        drConsec = 0;
        //if (evt[d].nMaskedChannels > 0) {
        // ++(run[d].channelsMasked);
        //}
        if (evt[d].nBadQuality) ++hasBadQ;
        if (evt[d].nOverflow) ++hasOver;
        if (evt[d].nUnderflow) ++hasUnder;
        if (evt[d].nSomeSignal) ++hasSig;
      }
    }
    if (drConsec != 0 && drConsecMax < 64) { // 64th drawer is bad - check transition from 64th to 1st drawer
      for (int d = 0; d < drConsecMax; ++d) {
        if (evt[d].nChannels == 0 || evt[d].nMaskedChannels >= evt[d].nChannels) {
          ++drConsec;
        } else {
          break;
        }
      }
      if (drConsec > drConsecMax) {
        drConsecMax = drConsec;
      }
    }
    unsigned int fl = 0;
    if (drAbsent > 15) {
      fl |= 0x01000000;
      error = xAOD::EventInfo::Error;
    }
    if (drMasked > 15) {
      fl |= 0x00100000;
      error = xAOD::EventInfo::Error;
    }
    //if (drConsecMax > 1)fl |= 0x00010000; // want to use these bits for length of consecutive area
    if (hasBadQ > 15) fl |= 0x00001000;
    if (hasOver) fl |= 0x00000100;
    if (hasUnder) fl |= 0x00000010;
    if (hasSig) fl |= 0x00000001;
 
#ifdef ALLOW_DEBUG_COUT
    std::cout<<"partition "<<p<<" drAbsent "<<drAbsent<<" drMasked "<<drMasked<<" drConsec "<<drConsecMax
    <<" hasBadQ "<<hasBadQ<<" hasOver "<<hasOver<<" hasUnder "<<hasUnder<<" hasSig "<<hasSig<<std::endl;
#endif
    flag |= fl << (p - 1);
  }
 
  // number of consecutively masked modules (if it's > 15 we have error already set)
  flag |= (std::min(15, drConsecMaxMax) << 16);
 
  if (drConsecMaxMax > 1 && error < xAOD::EventInfo::Warning) {
    // setting warning flag
    error = xAOD::EventInfo::Warning;
    // putting starting module number of consecutive bad area
    // instead of bits which indicates 16 masked or 16 absent modules in partition
    flag |= (drConsecNum << 20);
#ifdef ALLOW_DEBUG_COUT
    std::cout<<"warning in partition " << (drConsecNum>>6)+1 << " for modules "
    <<(drConsecNum)%64 <<" - " <<(drConsecNum+drConsecMaxMax-1)%64 <<std::endl;
#endif
  }
  
#ifdef ALLOW_DEBUG_COUT
  std::cout<<"partition flag 0x0"<<std::hex<<flag<<std::dec<<" error "<<error<<std::endl;
#endif
 
  //++m_eventErrorCounter[error]; // error index is 0 or 1 or 2 here
  //++m_eventErrorCounter[3]; // count separately total number of events
  
  // retrieve EventInfo
  SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInfoKey, ctx);
 
  if (eventInfo.isValid()) {
 
    if (flag != 0) {
      ATH_MSG_DEBUG( " set eventInfo for Tile for this event to 0x" << MSG::hex << flag << MSG::dec );
      if (!eventInfo->updateEventFlags(xAOD::EventInfo::Tile, flag)) {
        ATH_MSG_WARNING( " cannot set eventInfo for Tile " );
      }
    }
    
    if (error != xAOD::EventInfo::NotSet) {
      ATH_MSG_DEBUG( " set error bits for Tile for this event to " << error );
      if (!eventInfo->updateErrorState(xAOD::EventInfo::Tile, error)) {
        ATH_MSG_WARNING( " cannot set error state for Tile " );
      }
    }
 
  }
  else {
    ATH_MSG_WARNING( " cannot retrieve EventInfo, will not set Tile information " );
  }
  
  // Execution completed.
  ATH_MSG_DEBUG( "TileCellBuilderFromHit execution completed." );
 
  return StatusCode::SUCCESS;
}

qbits()

unsigned char TileCellBuilderFromHit::qbits ( TileDrawerEvtStatusArray &  drawerEvtStatus, int  ros, int  drawer, bool  count_over, bool  good_time, bool  good_ener, bool  overflow, bool  underflow, bool  good_overflowfit  ) const

private

method to compute the cell quality bits

Definition at line 431 of file TileCellBuilderFromHit.cxx.

                                                                                         {
 
  ++drawerEvtStatus[ros][drawer].nChannels;
  // new feature in rel 17.2.7 - count underflows and overflows
  if (count_over) {
    if (overflow) ++drawerEvtStatus[ros][drawer].nOverflow;
    if (underflow) ++drawerEvtStatus[ros][drawer].nUnderflow;
  }
#ifdef ALLOW_DEBUG_COUT
  if (overflow)  std::cout << "channel with overflow " << ((count_over)?"":"MBTS") << std::endl;
  if (underflow) std::cout << "channel with underflow " << ((count_over)?"":"MBTS") << std::endl;
  if (overfit)   std::cout << "channel with corrected overflow " << ((count_over)?"":"MBTS") << std::endl;
#endif
 
  unsigned char qbit = (overfit) ? (TileFragHash::FitFilter & TileCell::MASK_ALGO)
                                 : (m_RChType & TileCell::MASK_ALGO);
  if (good_time) qbit |= TileCell::MASK_TIME;
  if (overflow || underflow) qbit |= TileCell::MASK_OVER;
      
  if (good_ener) {
    qbit |= TileCell::MASK_AMPL;
    if (count_over) {
      ++drawerEvtStatus[ros][drawer].nSomeSignal;
    }
  }
 
  return qbit;
}

Member Data Documentation

E4ETA

const int TileCellBuilderFromHit::E4ETA = 2

staticprivate

Definition at line 226 of file TileCellBuilderFromHit.h.

E4NPHI

const int TileCellBuilderFromHit::E4NPHI = 4

staticprivate

Definition at line 227 of file TileCellBuilderFromHit.h.

E4SIDE

const int TileCellBuilderFromHit::E4SIDE = -1

staticprivate

Definition at line 225 of file TileCellBuilderFromHit.h.

m_affectedRegionInfo_current_run

std::vector<CaloAffectedRegionInfo> TileCellBuilderFromHit::m_affectedRegionInfo_current_run

private

method to process raw channels from a given vector and store them in collection

Definition at line 181 of file TileCellBuilderFromHit.h.

m_affectedRegionInfo_global

std::vector<CaloAffectedRegionInfo> TileCellBuilderFromHit::m_affectedRegionInfo_global

private

TileDrawerRunStatus m_drawerRunStatus[5][64]; //!< overall status of drawer in whole run.

Definition at line 180 of file TileCellBuilderFromHit.h.

m_cabling

const TileCablingService* TileCellBuilderFromHit::m_cabling {nullptr}

private

Pointer to TileCabling.

Definition at line 146 of file TileCellBuilderFromHit.h.

m_cablingSvc

ServiceHandle<TileCablingSvc> TileCellBuilderFromHit::m_cablingSvc

private

Initial value:

{ this,
        "TileCablingSvc", "TileCablingSvc", "The Tile cabling service"}

Name of Tile cabling service.

Definition at line 163 of file TileCellBuilderFromHit.h.

m_caloNoiseKey

SG::ReadCondHandleKey<CaloNoise> TileCellBuilderFromHit::m_caloNoiseKey

private

Initial value:

{this, "CaloNoise",
                                                    "electronicNoise",
                                                    "CaloNoise object to read"}

Definition at line 105 of file TileCellBuilderFromHit.h.

m_E1_TOWER

int TileCellBuilderFromHit::m_E1_TOWER

private

Definition at line 216 of file TileCellBuilderFromHit.h.

m_E4prContainerKey

SG::WriteHandleKey<TileCellContainer> TileCellBuilderFromHit::m_E4prContainerKey

private

Initial value:

{this, "E4prContainer", 
                                                             "E4prContainer",
                                                             "Output Tile E4 prime container key"}

Definition at line 121 of file TileCellBuilderFromHit.h.

m_emScaleKey

SG::ReadCondHandleKey<TileEMScale> TileCellBuilderFromHit::m_emScaleKey

private

Initial value:

{this,
         "TileEMScale", "TileEMScale", "Input Tile EMS calibration constants"}

Name of TileEMScale in condition store.

Definition at line 157 of file TileCellBuilderFromHit.h.

m_eneForTimeCut

float TileCellBuilderFromHit::m_eneForTimeCut

private

keep time for channels with energy above cut

Definition at line 133 of file TileCellBuilderFromHit.h.

m_eneForTimeCutMBTS

float TileCellBuilderFromHit::m_eneForTimeCutMBTS

private

similar cut for MBTS in pC

Definition at line 134 of file TileCellBuilderFromHit.h.

m_eventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> TileCellBuilderFromHit::m_eventInfoKey

private

Initial value:

{this, "EventInfo", 
                                                      "EventInfo", 
                                                      "EventInfo key"}

Definition at line 112 of file TileCellBuilderFromHit.h.

m_hitContainerKey

SG::ReadHandleKey<TileHitContainer> TileCellBuilderFromHit::m_hitContainerKey

private

Initial value:

{this, "TileHitContainer", 
                                                          "TileHitCnt", 
                                                          "Input Tile hit container key"}

Definition at line 108 of file TileCellBuilderFromHit.h.

m_infoName

std::string TileCellBuilderFromHit::m_infoName

private

Definition at line 131 of file TileCellBuilderFromHit.h.

m_maskBadChannels

bool TileCellBuilderFromHit::m_maskBadChannels

private

if true=> bad channels are masked

Definition at line 140 of file TileCellBuilderFromHit.h.

m_maxTime

float TileCellBuilderFromHit::m_maxTime

private

maximum time for the PMTs in the cels

Definition at line 138 of file TileCellBuilderFromHit.h.

m_MBTSContainerKey

SG::WriteHandleKey<TileCellContainer> TileCellBuilderFromHit::m_MBTSContainerKey

private

Initial value:

{this, "MBTSContainer", 
                                                             "MBTSContainer", 
                                                             "Output Tile MBTS container key"}

Definition at line 117 of file TileCellBuilderFromHit.h.

m_mbtsMgr

const MbtsDetDescrManager* TileCellBuilderFromHit::m_mbtsMgr {nullptr}

private

Pointer to MbtsDetDescrManager.

Definition at line 168 of file TileCellBuilderFromHit.h.

m_minTime

float TileCellBuilderFromHit::m_minTime

private

minimum time for the PMTs in the cels

Definition at line 139 of file TileCellBuilderFromHit.h.

m_noiseSigma

float TileCellBuilderFromHit::m_noiseSigma

private

cell electronic noise if CaloNoise is switched off

Definition at line 141 of file TileCellBuilderFromHit.h.

m_qualityCut

int TileCellBuilderFromHit::m_qualityCut

private

cut on channel quality (set energy to m_zeroEnergy for them)

Definition at line 136 of file TileCellBuilderFromHit.h.

m_RChType

TileFragHash::TYPE TileCellBuilderFromHit::m_RChType

private

Type of TileRawChannels (Fit, OF2, etc.)

Definition at line 170 of file TileCellBuilderFromHit.h.

m_rndmSvc

ServiceHandle<IAthRNGSvc> TileCellBuilderFromHit::m_rndmSvc { this, "RndmSvc", "AthRNGSvc", "Random Number Service used in TileCellBuildetFromHit" }

private

< Random number service to use

Definition at line 148 of file TileCellBuilderFromHit.h.

m_RUN2

bool TileCellBuilderFromHit::m_RUN2

private

Definition at line 214 of file TileCellBuilderFromHit.h.

m_RUN2plus

bool TileCellBuilderFromHit::m_RUN2plus

private

Definition at line 215 of file TileCellBuilderFromHit.h.

m_samplingFractionKey

SG::ReadCondHandleKey<TileSamplingFraction> TileCellBuilderFromHit::m_samplingFractionKey

private

Initial value:

{this,
        "TileSamplingFraction", "TileSamplingFraction", "Input Tile sampling fraction"}

Name of TileSamplingFraction in condition store.

Definition at line 127 of file TileCellBuilderFromHit.h.

m_tileBadChanTool

ToolHandle<ITileBadChanTool> TileCellBuilderFromHit::m_tileBadChanTool

private

Initial value:

{this,
        "TileBadChanTool", "TileBadChanTool", "Tile bad channel tool"}

Definition at line 151 of file TileCellBuilderFromHit.h.

m_tileHWID

const TileHWID* TileCellBuilderFromHit::m_tileHWID {nullptr}

private

Pointer to TileHWID.

Definition at line 145 of file TileCellBuilderFromHit.h.

m_tileID

const TileID* TileCellBuilderFromHit::m_tileID {nullptr}

private

Pointer to TileID.

Definition at line 143 of file TileCellBuilderFromHit.h.

m_tileMgr

const TileDetDescrManager* TileCellBuilderFromHit::m_tileMgr {nullptr}

private

Pointer to TileDetDescrManager.

Definition at line 167 of file TileCellBuilderFromHit.h.

m_tileTBID

const TileTBID* TileCellBuilderFromHit::m_tileTBID {nullptr}

private

Pointer to TileTBID.

Definition at line 144 of file TileCellBuilderFromHit.h.

m_zeroEnergy

float TileCellBuilderFromHit::m_zeroEnergy

private

energy to store in every PMT if both PMT are bad

Definition at line 135 of file TileCellBuilderFromHit.h.

NCELLE4PR

const int TileCellBuilderFromHit::NCELLE4PR = E4NPHI

staticprivate

Definition at line 228 of file TileCellBuilderFromHit.h.

NCELLMBTS

const int TileCellBuilderFromHit::NCELLMBTS = NSIDE * NPHI * NETA

staticprivate

Definition at line 221 of file TileCellBuilderFromHit.h.

NETA

const int TileCellBuilderFromHit::NETA = 2

staticprivate

Definition at line 220 of file TileCellBuilderFromHit.h.

NPHI

const int TileCellBuilderFromHit::NPHI = 8

staticprivate

Definition at line 219 of file TileCellBuilderFromHit.h.

NSIDE

const int TileCellBuilderFromHit::NSIDE = 2

staticprivate

Definition at line 218 of file TileCellBuilderFromHit.h.

---

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

• TileCellBuilderFromHit.h
• TileCellBuilderFromHit.cxx