LArNoisyROTool Class Reference

#include <LArNoisyROTool.h>

Inheritance diagram for LArNoisyROTool:

Collaboration diagram for LArNoisyROTool:

Classes
class	FEBEvtStat

Public Member Functions
	LArNoisyROTool ()=delete
	LArNoisyROTool (const std::string &type, const std::string &name, const IInterface *parent)
	Constructor with parameters: More...
virtual	~LArNoisyROTool ()=default
	Destructor: More...
virtual StatusCode	initialize ()
virtual std::unique_ptr< LArNoisyROSummary >	process (const EventContext &, const CaloCellContainer *, const std::set< unsigned int > *, const std::vector< HWIdentifier > *, const LArHVNMap *, const CaloDetDescrManager *, const LArHVIdMapping *) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...
virtual StatusCode	sysStart () override
	Handle START transition. More...
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef std::unordered_map< unsigned int, FEBEvtStat >	FEBEvtStatMap
typedef std::unordered_map< unsigned int, FEBEvtStat >::iterator	FEBEvtStatMapIt
typedef std::unordered_map< unsigned int, FEBEvtStat >::const_iterator	FEBEvtStatMapCstIt
typedef std::unordered_map< HWIdentifier, unsigned int >	HVlinesStatMap
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
size_t	partitionNumber (const HWIdentifier) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
std::unordered_map< unsigned int, unsigned int >	m_mapPSFEB
ToolHandle< ILArHVMapTool >	m_hvMapTool
const CaloCell_ID *	m_calo_id {}
const LArOnlineID *	m_onlineID {}
const LArElectrodeID *	m_elecID {}
SG::ReadCondHandleKey< LArOnOffIdMapping >	m_cablingKey {this, "CablingKey", "LArOnOffIdMap", "key to read OnOff mapping"}
unsigned int	m_CellQualityCut {}
bool	m_ignore_masked_cells {}
bool	m_ignore_front_innerwheel_cells {}
unsigned int	m_BadChanPerPA = 0U
unsigned int	m_BadChanPerFEB {}
unsigned int	m_MinBadFEB {}
unsigned int	m_SaturatedCellQualityCut {}
float	m_SaturatedCellEnergyTightCut {}
unsigned int	m_SaturatedCellTightCut {}
unsigned int	m_SaturatedCellTightCutEvents = 0U
bool	m_doHVline {}
float	m_BadChanFracPerHVline {}
unsigned int	m_MinBadHV {}
unsigned int	m_MNBLooseCut {}
unsigned int	m_MNBTightCut {}
std::vector< unsigned int >	m_MNBTight_PsVetoCut
std::array< uint8_t, 4 >	m_partitionMask {}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 42 of file LArNoisyROTool.h.

Member Typedef Documentation

FEBEvtStatMap

typedef std::unordered_map<unsigned int, FEBEvtStat> LArNoisyROTool::FEBEvtStatMap

private

Definition at line 103 of file LArNoisyROTool.h.

FEBEvtStatMapCstIt

typedef std::unordered_map<unsigned int, FEBEvtStat>::const_iterator LArNoisyROTool::FEBEvtStatMapCstIt

private

Definition at line 105 of file LArNoisyROTool.h.

FEBEvtStatMapIt

typedef std::unordered_map<unsigned int, FEBEvtStat>::iterator LArNoisyROTool::FEBEvtStatMapIt

private

Definition at line 104 of file LArNoisyROTool.h.

HVlinesStatMap

typedef std::unordered_map<HWIdentifier, unsigned int> LArNoisyROTool::HVlinesStatMap

private

Definition at line 107 of file LArNoisyROTool.h.

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

LArNoisyROTool() [1/2]

LArNoisyROTool::LArNoisyROTool ( )

delete

LArNoisyROTool() [2/2]

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

Constructor with parameters:

Definition at line 21 of file LArNoisyROTool.cxx.

                                           : 
  ::AthAlgTool  ( type, name, parent   ),m_hvMapTool("LArHVMapTool",this),
  m_calo_id(nullptr), m_onlineID(nullptr), 
  m_partitionMask({{LArNoisyROSummary::EMECAMask,LArNoisyROSummary::EMBAMask,LArNoisyROSummary::EMBCMask,LArNoisyROSummary::EMECCMask}}) //beware: The order matters! 
{
  declareInterface<ILArNoisyROTool >(this);
  declareProperty( "BadChanPerFEB", m_BadChanPerFEB=30 );
  declareProperty( "CellQualityCut", m_CellQualityCut=4000 );
  declareProperty( "IgnoreMaskedCells", m_ignore_masked_cells=false );
  declareProperty( "IgnoreFrontInnerWheelCells", m_ignore_front_innerwheel_cells=true );
  declareProperty( "BadFEBCut", m_MinBadFEB=3 );
 
  declareProperty( "MNBLooseCut",m_MNBLooseCut=5,"Number of cells above CellQualityCut");
  declareProperty( "MNBTightCut",m_MNBTightCut=17,"Number of cells above CellQualityCut");
  declareProperty( "MNBTight_PsVetoCut",m_MNBTight_PsVetoCut={13,3},"Number of cells above CellQualityCut");
  declareProperty( "SaturatedCellQualityCut", m_SaturatedCellQualityCut=65535);
  declareProperty( "SaturatedCellEnergyTightCut", m_SaturatedCellEnergyTightCut=1000.);
  declareProperty( "SaturatedCellTightCut", m_SaturatedCellTightCut=20);
 
  declareProperty( "DoHVflag", m_doHVline=true );
  declareProperty( "BadChanFracPerHVline", m_BadChanFracPerHVline=0.25 );
  declareProperty( "BadHVCut", m_MinBadHV=3 );
}

~LArNoisyROTool()

virtual LArNoisyROTool::~LArNoisyROTool ( )

virtualdefault

Destructor:

Member Function Documentation

declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleKeyArrayType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleKeyType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
  }

declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t, const SG::NotHandleType &    )

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(t);
  }

declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, SG::VarHandleBase &  hndl, const std::string &  doc, const SG::VarHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, SG::VarHandleKey &  hndl, const std::string &  doc, const SG::VarHandleKeyType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, SG::VarHandleKeyArray &  hndArr, const std::string &  doc, const SG::VarHandleKeyArrayType &    )

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc, const SG::NotHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(name, property, doc);
  }

declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc = "none"  )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

  {
    typedef typename SG::HandleClassifier<T>::type htype;
    return declareProperty (name, property, doc, htype());
  }

declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > &  t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                 {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

{ return m_evtStore; }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

initialize()

StatusCode LArNoisyROTool::initialize ( )

virtual

Definition at line 51 of file LArNoisyROTool.cxx.

                                      {
 
  if ( m_CellQualityCut > m_SaturatedCellQualityCut ) {
    ATH_MSG_FATAL( "Configuration problem: LArNoisyROTool assumes that the QFactor cut to declare a channel noisy is softer than the QFactor cut to declare the quality saturated !"  );
    return StatusCode::FAILURE;
  }
 
  if ( m_MNBLooseCut > m_MNBTightCut) {
    ATH_MSG_FATAL( "Configuration problem: LArNoisyROTool assumes that MNBLooseCut is smaller than MNBTightCut"  );
    return StatusCode::FAILURE;
  }
 
  ATH_CHECK(detStore()->retrieve(m_calo_id,"CaloCell_ID"));
  ATH_CHECK(detStore()->retrieve(m_onlineID,"LArOnlineID"));
  ATH_CHECK(detStore()->retrieve(m_elecID,"LArElectrodeID"));
  ATH_CHECK( m_cablingKey.initialize() );
 
  // Fill the map between any EMB FEB and the same FT PS FEB
  // Filled only for EMB so far
  for (std::vector<HWIdentifier>::const_iterator allFeb = m_onlineID->feb_begin(); 
       allFeb != m_onlineID->feb_end(); ++allFeb) {
    
    HWIdentifier febid = HWIdentifier(*allFeb);    
    int FEBIndex = febid.get_identifier32().get_compact();
    int FEBIndex_PS = 0;
    int barrel_ec = m_onlineID->barrel_ec(febid);
    if (barrel_ec == 0){
      int pos_neg   = m_onlineID->pos_neg(febid);
      int ft        = m_onlineID->feedthrough(febid);
      FEBIndex_PS = (m_onlineID->feb_Id(0,pos_neg,ft,1)).get_identifier32().get_compact();      
    }
    m_mapPSFEB[FEBIndex] = FEBIndex_PS;
  }
 
  return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

interfaceID()

static const InterfaceID& LArNoisyROTool::interfaceID ( )

static

msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

partitionNumber()

size_t LArNoisyROTool::partitionNumber ( const HWIdentifier  hwid ) const

inlineprivate

Definition at line 165 of file LArNoisyROTool.h.

                                                                           {
  int pn=m_onlineID->pos_neg(hwid);
  if (m_onlineID->isEMECchannel(hwid)) {
    if (pn) 
      return 0; //positive EMECA side
    else
      return 3; //negative EMECC side
  }
  if (m_onlineID->isEMBchannel(hwid)) {
    if (pn) 
      return 1; //positive EMBA side
    else
      return 2; //negative EMBC side
  }
  return 4;//Anything else
}

process()

std::unique_ptr< LArNoisyROSummary > LArNoisyROTool::process ( const EventContext &  ctx, const CaloCellContainer *  cellContainer, const std::set< unsigned int > *  knownBadFEBs, const std::vector< HWIdentifier > *  knownMNBFEBs, const LArHVNMap *  nCellsperLine, const CaloDetDescrManager *  cddm, const LArHVIdMapping *  hvid  ) const

virtual

Implements ILArNoisyROTool.

Definition at line 89 of file LArNoisyROTool.cxx.

                                                                                                                                                                                                                                                                                                                 {
 
  // sanity check
  bool doHVline=m_doHVline;
  if(doHVline && ( !nCellsperLine || !cddm || !hvid)) {
     ATH_MSG_ERROR("HV line flagging asked, but missing ingrediences, switching off !!");
     doHVline = false;
  }
 
  std::unique_ptr<LArNoisyROSummary> noisyRO(new LArNoisyROSummary);
 
  if(!cellContainer) return noisyRO;
 
  SG::ReadCondHandle<LArOnOffIdMapping> larCablingHdl(m_cablingKey, ctx);
  const LArOnOffIdMapping* cabling=*larCablingHdl;
  
  FEBEvtStatMap FEBStats; //counter per FEB
  HVlinesStatMap HVStats; //counter per HV line
 
  unsigned int NsaturatedTightCutBarrelA = 0;
  unsigned int NsaturatedTightCutBarrelC = 0;
  unsigned int NsaturatedTightCutEMECA = 0;
  unsigned int NsaturatedTightCutEMECC = 0;
 
 
  CaloCellContainer::const_iterator cellItr    = cellContainer->begin();
  CaloCellContainer::const_iterator cellItrEnd = cellContainer->end();
  for ( ; cellItr != cellItrEnd; ++cellItr )
  {
    const CaloCell* cell = (*cellItr);
    if (!cell) continue;
 
    // only cells with a bad enough Quality Factor
    if ( cell->quality() < m_CellQualityCut ) continue;
 
    // cells with zero energy have been masked by previous algorithms
    // they should not matter for physics so don't consider them
    if ( m_ignore_masked_cells && std::abs(cell->e()) < 0.1 ) continue; //Fixme: use provenance
 
    Identifier id = cell->ID();
    if (m_ignore_front_innerwheel_cells && m_calo_id->is_em_endcap_inner(id) && m_calo_id->sampling(id) == 1) continue; // Front inner wheel cells are ignored
 
    // saturated Qfactor ? Tight cuts.
    if ( cell->quality()>=m_SaturatedCellQualityCut && 
     std::abs(cell->e()) > m_SaturatedCellEnergyTightCut )
    {
      bool sideA = cell->eta() > 0.;
      if ( m_calo_id->is_em_barrel(id) )
      {
    if ( sideA ) { NsaturatedTightCutBarrelA++; } 
    else { NsaturatedTightCutBarrelC++; }
      }
      else if ( m_calo_id->is_em_endcap(id) )
      {
    if ( sideA ) { NsaturatedTightCutEMECA++; } 
    else { NsaturatedTightCutEMECC++; }
      }
    }
 
 
    // only LAr EM for bad FEBs
    if ( m_calo_id->is_em(id) ) 
    {
      // get FEB ID and channel number
      HWIdentifier hwid = cabling->createSignalChannelID(id);
      HWIdentifier febid = m_onlineID->feb_Id(hwid);
      unsigned int FEBindex = febid.get_identifier32().get_compact();
      unsigned int channel = m_onlineID->channel(hwid);    
      FEBStats[FEBindex].addBadChannel(channel);
    }
 
    if(doHVline) {
       //HVline, in all calos
       std::vector<HWIdentifier> hvlines;
       m_hvMapTool->GetHVLines(id, cddm,hvlines);
       for(unsigned int i=0; i<hvlines.size(); ++i) {
          if(HVStats.contains(hvlines[i])) HVStats[hvlines[i]] += 1; else HVStats[hvlines[i]]=1;
       }
    }
 
  }
 
  // Store the Saturated flag per partition
  uint8_t SatTightPartitions = 0;
  if ( NsaturatedTightCutBarrelA >= m_SaturatedCellTightCut ) SatTightPartitions |= LArNoisyROSummary::EMBAMask;
  if ( NsaturatedTightCutBarrelC >= m_SaturatedCellTightCut ) SatTightPartitions |= LArNoisyROSummary::EMBCMask;
  if ( NsaturatedTightCutEMECA >= m_SaturatedCellTightCut ) SatTightPartitions |= LArNoisyROSummary::EMECAMask;
  if ( NsaturatedTightCutEMECC >= m_SaturatedCellTightCut ) SatTightPartitions |= LArNoisyROSummary::EMECCMask;
  bool badSaturatedTightCut = (SatTightPartitions != 0);
  if ( badSaturatedTightCut ) noisyRO-> SetSatTightFlaggedPartitions(SatTightPartitions);
 
  // loop on all FEBs and check whether FEB can be declared as bad for the different type of flags:
  // regular noise burst, weighted noise burst, MNB tight and loose
  for ( FEBEvtStatMapCstIt it = FEBStats.begin(); it != FEBStats.end(); ++it ) {
    ATH_MSG_VERBOSE(" candidate FEB " << it->first << " with " << it->second.badChannels() << " bad channels");
    if ( it->second.badChannels() > m_BadChanPerFEB ) {
      noisyRO->add_noisy_feb(HWIdentifier(it->first));
    }
 
    // Loose MNBs
    if ( it->second.badChannels() > m_MNBLooseCut ){
       noisyRO->add_MNBLoose_feb(HWIdentifier(it->first));
       ATH_MSG_DEBUG("Loose bad FEB " << it->first << " " << m_onlineID->channel_name(HWIdentifier(it->first)) << " with " << it->second.badChannels() << " bad channels");
       // Tight_PsVeto MNBs
       if ( it->second.badChannels() > m_MNBTight_PsVetoCut[0] ){
         unsigned int associatedPSFEB = m_mapPSFEB.find(it->first)->second;
         if (associatedPSFEB != 0){ // Check if a PS FEB is associated (TRUE only for EMB FEBs)
           if (FEBStats.count(associatedPSFEB) == 0) noisyRO->add_MNBTight_PsVeto_feb(HWIdentifier(it->first));
           else if (FEBStats[associatedPSFEB].badChannels() < m_MNBTight_PsVetoCut[1]) noisyRO->add_MNBTight_PsVeto_feb(HWIdentifier(it->first));
         }
       }
       // Tight MNBs
       if ( it->second.badChannels() > m_MNBTightCut ){
          noisyRO->add_MNBTight_feb(HWIdentifier(it->first));
       }
    }
 
 
//  // Noisy preamp removed as no used currently
//  // Kept here just in case we may want to revive it
//    const unsigned int* PAcounters = it->second.PAcounters();
//    for ( size_t i = 0; i < 32; i++ ) {
//      if ( PAcounters[i] > m_BadChanPerPA ) {
//  uint64_t PAid = static_cast<uint64_t>(1000000000)*static_cast<uint64_t>(i)+static_cast<uint64_t>(it->first);
//  ATH_MSG_DEBUG(" bad preamp " << i << " in FEB " << it->first << "  ID " << PAid);
//  noisyRO->add_noisy_preamp(HWIdentifier(it->first),4*i);
//   if (m_printSummary) m_badPA_counters[PAid]++;
//      }
//    }
 
  }//end loop over m_FEBats
 
 
  // Count noisy FEB per partition EMEC-EMB - Simple and weighted quantities
  unsigned int NBadFEBEMECA = 0; unsigned int NBadFEBEMECA_W = 0;
  unsigned int NBadFEBEMECC = 0; unsigned int NBadFEBEMECC_W = 0;
  unsigned int NBadFEBEMBA = 0; unsigned int NBadFEBEMBA_W = 0;
  unsigned int NBadFEBEMBC = 0; unsigned int NBadFEBEMBC_W = 0;
 
  const std::vector<HWIdentifier>& badfebs = noisyRO->get_noisy_febs();
  
  //for ( std::vector<HWIdentifier>::const_iterator febit = badfebs.begin();
  //    febit != badfebs.end(); febit++ )
  for (const HWIdentifier& febid : badfebs) 
  {
    // first channel of FEB, as safety since FEBid seem to be the Id of the
    // first channel (no guarantee?)
    HWIdentifier chanID = m_onlineID->channel_Id(febid,0);
 
    int weight = 1;
    // If the FEB is known to be subject to noise burst (list defiend as property)
    // give a weight 2
    const unsigned int int_id =  febid.get_identifier32().get_compact();
    if (knownBadFEBs->find(int_id)!=knownBadFEBs->end()) weight=2;
 
    if ( m_onlineID->isEMBchannel(chanID) ) 
    {
      if ( m_onlineID->pos_neg(chanID) == 1 ){
      NBadFEBEMBA_W = NBadFEBEMBA_W + weight;
      NBadFEBEMBA++;
      }
      else{
    NBadFEBEMBC_W = NBadFEBEMBC_W + weight;
    NBadFEBEMBC++;
      }
    }
    else if ( m_onlineID->isEMECchannel(chanID) ) 
    {
      if ( m_onlineID->pos_neg(chanID) == 1 ){
    NBadFEBEMECA_W = NBadFEBEMECA_W + weight;
    NBadFEBEMECA++;
      }
      else{
    NBadFEBEMECC_W = NBadFEBEMECC_W + weight;
    NBadFEBEMECC++;
      }
    }
  } 
 
  uint8_t BadFEBPartitions = 0;
  if ( NBadFEBEMBA  > m_MinBadFEB )  BadFEBPartitions |= LArNoisyROSummary::EMBAMask;
  if ( NBadFEBEMBC  > m_MinBadFEB )  BadFEBPartitions |= LArNoisyROSummary::EMBCMask;
  if ( NBadFEBEMECA  > m_MinBadFEB )  BadFEBPartitions |= LArNoisyROSummary::EMECAMask;
  if ( NBadFEBEMECC  > m_MinBadFEB )  BadFEBPartitions |= LArNoisyROSummary::EMECCMask;
  if ( BadFEBPartitions != 0  ) noisyRO-> SetBadFEBFlaggedPartitions(BadFEBPartitions);
 
  uint8_t BadFEBPartitions_W = 0;
  if ( NBadFEBEMBA_W  > m_MinBadFEB )  BadFEBPartitions_W |= LArNoisyROSummary::EMBAMask;
  if ( NBadFEBEMBC_W  > m_MinBadFEB )  BadFEBPartitions_W |= LArNoisyROSummary::EMBCMask;
  if ( NBadFEBEMECA_W  > m_MinBadFEB )  BadFEBPartitions_W |= LArNoisyROSummary::EMECAMask;
  if ( NBadFEBEMECC_W  > m_MinBadFEB )  BadFEBPartitions_W |= LArNoisyROSummary::EMECCMask;
  if ( BadFEBPartitions_W != 0 ) noisyRO-> SetBadFEB_WFlaggedPartitions(BadFEBPartitions_W);
 
 
 
  //Check for Mini Noise Bursts:
  uint8_t MNBTightPartition=0;
  uint8_t MNBTight_PsVetoPartition=0;
  uint8_t MNBLoosePartition=0;
  
  std::array<unsigned,5> nTightMNBFEBSperPartition{};
  std::array<unsigned,5> nTight_PsVetoMNBFEBSperPartition{};
  std::array<unsigned,5> nLooseMNBFEBSperPartition{};
  for (HWIdentifier febid: *knownMNBFEBs) { //Loop over known MNB FEBs
    FEBEvtStatMapCstIt statIt=FEBStats.find(febid.get_identifier32().get_compact());
    if (statIt!=FEBStats.end()) {
      if (statIt->second.badChannels()>=m_MNBLooseCut) {
          (nLooseMNBFEBSperPartition[partitionNumber(febid)])++;
        // Tight_PsVeto MNBs
          if ( statIt->second.badChannels() > m_MNBTight_PsVetoCut[0] ){
            auto found = m_mapPSFEB.find(statIt->first);
            if (found != m_mapPSFEB.end()){
              if (unsigned int associatedPSFEB = found->second; associatedPSFEB != 0){
                if (FEBStats.count(associatedPSFEB) == 0) (nTight_PsVetoMNBFEBSperPartition[partitionNumber(febid)])++;
                else if (FEBStats[associatedPSFEB].badChannels() < m_MNBTight_PsVetoCut[1]) (nTight_PsVetoMNBFEBSperPartition[partitionNumber(febid)])++;
              }
            }
          }
        // Tight MNBs
        if (statIt->second.badChannels()>=m_MNBTightCut)
          (nTightMNBFEBSperPartition[partitionNumber(febid)])++;
      }
    }//End FEB in list of bad-Q FEBs
  }//end loop over known MNB Febs
 
 
  for (unsigned iP=0;iP<4;++iP) {
    ATH_MSG_DEBUG( "Partition " << iP << ": Found " << nLooseMNBFEBSperPartition[iP] << " MNB FEBs with more than " <<  m_MNBLooseCut << " bad-Q channels"  );
    ATH_MSG_DEBUG( "Partition " << iP << ": Found " << nTightMNBFEBSperPartition[iP] << " MNB FEBs with more than " <<  m_MNBTightCut << " bad-Q channels"  );
    if (nLooseMNBFEBSperPartition[iP]>0) MNBLoosePartition |= m_partitionMask[iP];
    if (nTightMNBFEBSperPartition[iP]>0) MNBTightPartition |= m_partitionMask[iP];
    if (nTight_PsVetoMNBFEBSperPartition[iP]>0) MNBTight_PsVetoPartition |= m_partitionMask[iP];
  }// end loop over partitions      
 
 
  noisyRO->SetMNBTightFlaggedPartitions(MNBTightPartition);
  noisyRO->SetMNBTight_PsVetoFlaggedPartitions(MNBTight_PsVetoPartition);
  noisyRO->SetMNBLooseFlaggedPartitions(MNBLoosePartition);
 
  if(!hvid || !nCellsperLine ) return noisyRO; // do not have HVcells map
 
  // Count noisy HVlines per partition
  unsigned int NBadHVEMECA = 0;
  unsigned int NBadHVEMECC = 0;
  unsigned int NBadHVEMBA = 0; 
  unsigned int NBadHVEMBC = 0; 
  unsigned int NBadHVHECA = 0; 
  unsigned int NBadHVHECC = 0; 
  unsigned int NBadHVFCALA = 0; 
  unsigned int NBadHVFCALC = 0; 
 
  // loop over HVlines, to check if they are qualified as noisy
  for ( HVlinesStatMap::const_iterator it = HVStats.begin(); it != HVStats.end(); ++it ) {
    ATH_MSG_DEBUG(ctx.eventID().event_number()<<" candidate HVline " << it->first << " with " << it->second << " bad channels, out of "<<nCellsperLine->HVNcell(HWIdentifier(it->first))<<" channels");
    if ( it->second >= m_BadChanFracPerHVline * nCellsperLine->HVNcell(it->first) ) {
      HWIdentifier hwd(it->first); 
      noisyRO->add_noisy_hvline(hwd);
      const std::vector<HWIdentifier> elecVec = hvid->getLArElectrodeIDvec(hwd);
      int side = m_elecID->zside(elecVec[0]);
      int part = m_elecID->detector(elecVec[0]);
      switch(side) {
         case 1: { ATH_MSG_DEBUG("Elec. side: "<<side);
                    switch(part){
                    case 0: case 1:{NBadHVEMBC += 1; break;}
                    case 2: case 3:{NBadHVEMECC += 1; break;}
                    case 4:        {NBadHVHECC += 1; break;}
                    case 5:        {NBadHVFCALC += 1; break;}
                    default:       {ATH_MSG_WARNING("Wrong HV line detector "<<part); break;}
                   };
                   break; 
                 }
         case 0: { switch(part){
                    case 0: case 1:{NBadHVEMBA += 1; break;}
                    case 2: case 3:{NBadHVEMECA += 1; break;}
                    case 4:        {NBadHVHECA += 1; break;}
                    case 5:        {NBadHVFCALA += 1; break;}
                    default:       {ATH_MSG_WARNING("Wrong HV line detector "<<part); break;}
                   };
                   break;               
                 }
      }
    }// HVline is noisy
  }// all hvlines 
 
  uint8_t BadHVPartitions = 0;
  if ( NBadHVEMBA  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::EMBAMask;
  if ( NBadHVEMBC  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::EMBCMask;
  if ( NBadHVEMECA  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::EMECAMask;
  if ( NBadHVEMECC  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::EMECCMask;
  if ( NBadHVHECA  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::HECAMask;
  if ( NBadHVHECC  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::HECCMask;
  if ( NBadHVFCALA  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::FCALAMask;
  if ( NBadHVFCALC  >= m_MinBadHV )  BadHVPartitions |= LArNoisyROSummary::FCALCMask;
 
  if ( BadHVPartitions != 0 ) noisyRO-> SetBadHVlinesPartitions(BadHVPartitions);
 
  return noisyRO;
}

renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T &  h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_BadChanFracPerHVline

float LArNoisyROTool::m_BadChanFracPerHVline {}

private

Definition at line 151 of file LArNoisyROTool.h.

m_BadChanPerFEB

unsigned int LArNoisyROTool::m_BadChanPerFEB {}

private

Definition at line 129 of file LArNoisyROTool.h.

m_BadChanPerPA

unsigned int LArNoisyROTool::m_BadChanPerPA = 0U

private

Definition at line 126 of file LArNoisyROTool.h.

m_cablingKey

SG::ReadCondHandleKey<LArOnOffIdMapping> LArNoisyROTool::m_cablingKey {this, "CablingKey", "LArOnOffIdMap", "key to read OnOff mapping"}

private

Definition at line 114 of file LArNoisyROTool.h.

m_calo_id

const CaloCell_ID* LArNoisyROTool::m_calo_id {}

private

Definition at line 111 of file LArNoisyROTool.h.

m_CellQualityCut

unsigned int LArNoisyROTool::m_CellQualityCut {}

private

Definition at line 117 of file LArNoisyROTool.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_doHVline

bool LArNoisyROTool::m_doHVline {}

private

Definition at line 148 of file LArNoisyROTool.h.

m_elecID

const LArElectrodeID* LArNoisyROTool::m_elecID {}

private

Definition at line 113 of file LArNoisyROTool.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_hvMapTool

ToolHandle<ILArHVMapTool> LArNoisyROTool::m_hvMapTool

private

Definition at line 110 of file LArNoisyROTool.h.

m_ignore_front_innerwheel_cells

bool LArNoisyROTool::m_ignore_front_innerwheel_cells {}

private

Definition at line 123 of file LArNoisyROTool.h.

m_ignore_masked_cells

bool LArNoisyROTool::m_ignore_masked_cells {}

private

Definition at line 120 of file LArNoisyROTool.h.

m_mapPSFEB

std::unordered_map<unsigned int,unsigned int> LArNoisyROTool::m_mapPSFEB

private

Definition at line 106 of file LArNoisyROTool.h.

m_MinBadFEB

unsigned int LArNoisyROTool::m_MinBadFEB {}

private

Definition at line 132 of file LArNoisyROTool.h.

m_MinBadHV

unsigned int LArNoisyROTool::m_MinBadHV {}

private

Definition at line 154 of file LArNoisyROTool.h.

m_MNBLooseCut

unsigned int LArNoisyROTool::m_MNBLooseCut {}

private

Definition at line 155 of file LArNoisyROTool.h.

m_MNBTight_PsVetoCut

std::vector<unsigned int> LArNoisyROTool::m_MNBTight_PsVetoCut

private

Definition at line 157 of file LArNoisyROTool.h.

m_MNBTightCut

unsigned int LArNoisyROTool::m_MNBTightCut {}

private

Definition at line 156 of file LArNoisyROTool.h.

m_onlineID

const LArOnlineID* LArNoisyROTool::m_onlineID {}

private

Definition at line 112 of file LArNoisyROTool.h.

m_partitionMask

std::array<uint8_t,4> LArNoisyROTool::m_partitionMask {}

private

Definition at line 158 of file LArNoisyROTool.h.

m_SaturatedCellEnergyTightCut

float LArNoisyROTool::m_SaturatedCellEnergyTightCut {}

private

Definition at line 139 of file LArNoisyROTool.h.

m_SaturatedCellQualityCut

unsigned int LArNoisyROTool::m_SaturatedCellQualityCut {}

private

Definition at line 136 of file LArNoisyROTool.h.

m_SaturatedCellTightCut

unsigned int LArNoisyROTool::m_SaturatedCellTightCut {}

private

Definition at line 142 of file LArNoisyROTool.h.

m_SaturatedCellTightCutEvents

unsigned int LArNoisyROTool::m_SaturatedCellTightCutEvents = 0U

private

Definition at line 145 of file LArNoisyROTool.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• LArNoisyROTool.h
• LArNoisyROTool.cxx