#include <LArShapeDumper.h>

Inheritance diagram for LArShapeDumper:

Collaboration diagram for LArShapeDumper:

Public Member Functions
	LArShapeDumper (const std::string &name, ISvcLocator *pSvcLocator)

	~LArShapeDumper ()

virtual StatusCode	initialize () override

virtual StatusCode	start () override

virtual StatusCode	execute () override

virtual StatusCode	stop () override

virtual StatusCode	finalize () override

int	makeEvent (LArSamples::EventData *&eventData, int run, int event, int lumiBlock, int bunchXing) const

virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...

virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...

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	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

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 ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
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::string	m_fileName

int	m_prescale

unsigned	m_count

unsigned	m_maxChannels

unsigned	m_nWrongBunchGroup

unsigned	m_nPrescaledAway

unsigned	m_nLArError

unsigned	m_nNoDigits

unsigned	m_nNoDigitsSC

std::string	m_caloType

double	m_energyCut

double	m_noiseSignifCut

double	m_energyCutSC

int	m_minADCMax

int	m_minADCMaxSC

std::string	m_gainSpec

bool	m_dumpDisc

std::vector< std::string >	m_triggerNames

ToolHandle< ILArShapeDumperTool >	m_dumperTool {this,"LArShapeDumperTool","LArShapeDumperTool"}

ToolHandle< ILArShapeDumperTool >	m_dumperToolSC {this,"LArShapeDumperToolSC","LArShapeDumperToolSC"}

LArBadChannelMask	m_bcMask

Gaudi::Property< std::vector< std::string > >	m_problemsToMask {this,"ProblemsToMask",{}, "Bad-Channel categories to patch"}

LArBadChannelMask	m_bcMaskSC

Gaudi::Property< std::vector< std::string > >	m_problemsToMaskSC {this,"ProblemsToMaskSC",{}, "Bad-Channel categories to mask"}

SG::ReadHandleKey< LArDigitContainer >	m_digitsKey {this, "DigitsKey", "FREE", "key for LArDigitContainer"}

SG::ReadHandleKey< LArRawChannelContainer >	m_channelsKey {this, "ChannelsKey", "LArRawChannels", "key for LArRawChannels"}

SG::ReadHandleKey< LArDigitContainer >	m_digitsKeySC {this, "DigitsKeySC", "SC_ADC_BAS", "key for LArDigitContainer for SC"}

SG::ReadHandleKey< LArRawSCContainer >	m_rawscKey {this, "RawSCKey", "SC_ET_ID", "key for LArRawSCContainer"}

SG::ReadCondHandleKey< LArADC2MeV >	m_adc2mevKey {this,"ADC2MeVKey","LArADC2MeV","SG Key of ADC2MeV conditions object"}

PublicToolHandle< Trig::TrigDecisionTool >	m_trigDec {this, "TrigDecisionTool", "", "Handle to the TrigDecisionTool"}

SG::ReadCondHandleKey< LArOnOffIdMapping >	m_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}

SG::ReadCondHandleKey< LArOnOffIdMapping >	m_cablingKeySC {this,"CablingKeySC","LArOnOffIdMapSC","SG Key of LArOnOffIdMapping object for SC"}

SG::ReadCondHandleKey< LArBadChannelCont >	m_BCKey {this, "BadChanKey", "LArBadChannel", "SG bad channels key"}

SG::ReadCondHandleKey< LArBadChannelCont >	m_BCKeySC {this, "BadChanKeySC", "LArBadChannelSC", "SG bad channels key"}

SG::ReadCondHandleKey< CaloNoise >	m_noiseCDOKey {this,"CaloNoiseKey","totalNoise","SG Key of CaloNoise data object"}

SG::ReadCondHandleKey< BunchCrossingCondData >	m_bcDataKey {this, "BunchCrossingCondDataKey", "BunchCrossingData" ,"SG Key of BunchCrossing CDO"}

SG::ReadCondHandleKey< ILArPedestal >	m_pedestalKey {this,"PedestalKey","LArPedestal","SG Key of LArPedestal object"}

SG::ReadCondHandleKey< ILArPedestal >	m_pedestalKeySC {this,"PedestalKeySC","LArPedestalSC","SG Key of LArPedestal object"}

SG::ReadCondHandleKey< ILArAutoCorr >	m_acorrKey {this,"AutocorrKey","LArAutoCorr","SG Key of LArAutoCorr object"}

SG::ReadCondHandleKey< CaloDetDescrManager >	m_caloMgrKey

SG::ReadCondHandleKey< CaloSuperCellDetDescrManager >	m_caloSuperCellMgrKey

const LArOnlineID *	m_onlineHelper

const LArOnline_SuperCellID *	m_onlineHelperSC

bool	m_doStream

bool	m_doTrigger

bool	m_doOFCIter

bool	m_doAllEvents

bool	m_doRoIs

bool	m_doAllLvl1

bool	m_dumpChannelInfos

bool	m_doEM

bool	m_doHEC

bool	m_doFCAL

bool	m_doSC

bool	m_gains [CaloGain::LARNGAIN] {}

bool	m_onlyEmptyBC

LArSamples::DataStore *	m_samples

std::unique_ptr< LArSamples::RunData >	m_runData

std::vector< const Trig::ChainGroup * >	m_triggerGroups

TRandom	m_random

DataObjIDColl	m_extendedExtraObjects

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

Author: Nicol.nosp@m.as.B.nosp@m.erger.nosp@m.@cer.nosp@m.n.ch

Definition at line 66 of file LArShapeDumper.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ LArShapeDumper()

LArShapeDumper::LArShapeDumper	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 42 of file LArShapeDumper.cxx.

                                                                                 : 
   AthAlgorithm(name, pSvcLocator),
   m_count(0),
   m_nWrongBunchGroup(0),
   m_nPrescaledAway(0),
   m_nLArError(0),
   m_nNoDigits(0),
   m_nNoDigitsSC(0),
   m_onlineHelper(nullptr),
   m_onlineHelperSC(nullptr),
   m_doEM(false),
   m_doHEC(false),
   m_doFCAL(false),
   m_doSC(false),
   m_samples(nullptr)
 {
   declareProperty("FileName", m_fileName = "samples.root");
   declareProperty("MaxChannels", m_maxChannels = 200000);
   declareProperty("Prescale", m_prescale = 1);
   declareProperty("CaloType", m_caloType = "EMHECFCALSC");
   declareProperty("EnergyCut", m_energyCut = -1);
   declareProperty("EnergyCutSC", m_energyCutSC = -1);
   declareProperty("NoiseSignifCut", m_noiseSignifCut = 3);
   declareProperty("MinADCMax", m_minADCMax = -1);
   declareProperty("MinADCMaxSC", m_minADCMaxSC = -1);
   declareProperty("Gains", m_gainSpec = "HIGH,MEDIUM,LOW");
   declareProperty("DumpDisconnected", m_dumpDisc = false);
   declareProperty("DoStream", m_doStream = false);
   declareProperty("DoTrigger", m_doTrigger = true);
   declareProperty("DoOFCIter", m_doOFCIter = true);
   declareProperty("DoAllEvents", m_doAllEvents = true);
   declareProperty("DumpChannelInfos", m_dumpChannelInfos = false);
   declareProperty("DoRoIs", m_doRoIs = true);
   declareProperty("TriggerNames", m_triggerNames);
   declareProperty("DoAllLvl1", m_doAllLvl1 = true);
   declareProperty("onlyEmptyBC",m_onlyEmptyBC=false);
 }

◆ ~LArShapeDumper()

LArShapeDumper::~LArShapeDumper ( )

Definition at line 81 of file LArShapeDumper.cxx.

82 {

83 }

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::detStore ( ) const

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ execute()

StatusCode LArShapeDumper::execute ( )

overridevirtual

Definition at line 215 of file LArShapeDumper.cxx.

 {    
   m_count++;
   const EventContext& ctx = Gaudi::Hive::currentContext();
   if ((m_prescale > 1 && m_random.Rndm() > 1.0/m_prescale) || m_prescale <= 0) {
     ATH_MSG_VERBOSE ( "======== prescaling event "<< m_count << " ========" );
     m_nPrescaledAway++;
     return StatusCode::SUCCESS;
   }
  
   ATH_MSG_VERBOSE ( "======== executing event "<< m_count << " ========" );
  
   const xAOD::EventInfo* eventInfo = nullptr;
   ATH_CHECK( evtStore()->retrieve(eventInfo) );
   
   int event     = eventInfo->eventNumber();
   int run       = eventInfo->runNumber();
   int lumiBlock = eventInfo->lumiBlock();
   int bunchId   = eventInfo->bcid();
  
   
   if (eventInfo->errorState(xAOD::EventInfo::LAr)==xAOD::EventInfo::Error) {
     ATH_MSG_DEBUG("Ignoring Event b/c of LAr ERROR");
     m_nLArError++;
     return StatusCode::SUCCESS;
   }
  
  
   SG::ReadCondHandle<BunchCrossingCondData> bccd (m_bcDataKey,ctx);
   const BunchCrossingCondData* bunchCrossing=*bccd;
   if (!bunchCrossing) {
     ATH_MSG_ERROR("Failed to retrieve Bunch Crossing obj");
     return StatusCode::FAILURE;
   }
  
   SG::ReadCondHandle<CaloDetDescrManager> caloMgrHandle{m_caloMgrKey};
   ATH_CHECK(caloMgrHandle.isValid());
   const CaloDetDescrManager* caloMgr = *caloMgrHandle;  
  
   if (m_onlyEmptyBC) {
     if (!bccd->isFilled(bunchId)) {
       ATH_MSG_DEBUG("Ignoring Event with bunch crossing type ");
       m_nWrongBunchGroup++;
       return StatusCode::SUCCESS;
     }
    }
  
  
   EventData* eventData = nullptr;
   int eventIndex = -1;
   if (m_doAllEvents) {
     eventIndex = makeEvent(eventData, run, event, lumiBlock, bunchId);
     if (eventIndex < 0) return StatusCode::FAILURE;
   }
  
   SG::ReadHandle<LArDigitContainer> hdlDigit(m_digitsKey, ctx);
   if(!hdlDigit.isValid()) {
      ATH_MSG_WARNING( "Unable to retrieve LArDigitContainer with key " << m_digitsKey << " from DetectorStore. " );
      return StatusCode::SUCCESS;
   } else
       ATH_MSG_DEBUG( "Got LArDigitContainer with key " << m_digitsKey.key() );
   const LArDigitContainer* larDigitContainer = &(*hdlDigit);
  
   if (larDigitContainer->empty()) {
     ATH_MSG_WARNING ( "LArDigitContainer with key=" << m_digitsKey.key() << " is empty!" );
     m_nNoDigits++;
     return StatusCode::SUCCESS;
   }
  
   SG::ReadHandle<LArRawChannelContainer> hdlRaw(m_channelsKey, ctx);
   if(!hdlRaw.isValid()) {
      ATH_MSG_WARNING( "Unable to retrieve LArRawChannelsContainer with key " << m_channelsKey.key() << " from DetectorStore. " );
      return StatusCode::SUCCESS;
   } else
       ATH_MSG_DEBUG( "Got LArRawChannelsContainer with key " << m_channelsKey.key() );
   const LArRawChannelContainer* rawChannelContainer = &(*hdlRaw);;
  
   if (rawChannelContainer->empty()) {
     ATH_MSG_WARNING ( "LArRawChannelContainer with key=" << m_channelsKey << " is empty!" );
     m_nNoDigits++;
     return StatusCode::SUCCESS;
   }
  
   
   
   SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey};
   const LArOnOffIdMapping* cabling=*cablingHdl;
   if(!cabling) {
      ATH_MSG_ERROR( "Do not have cabling object LArOnOffIdMapping" );
      return StatusCode::FAILURE;
   }
  
   SG::ReadCondHandle<LArADC2MeV> adc2mevHdl(m_adc2mevKey, ctx);
   const LArADC2MeV* adc2MeV=*adc2mevHdl;
   if(!adc2MeV) {
      ATH_MSG_ERROR( "Failed to retreive ADC2MeV cond obj" );
      return StatusCode::FAILURE;
   }
  
   SG::ReadCondHandle<ILArPedestal> pedHdl(m_pedestalKey, ctx);
   const ILArPedestal* pedestals=*pedHdl;
   if (!pedestals) {
     ATH_MSG_ERROR("Failed to retrieve pedestal cond obj");
      return StatusCode::FAILURE;
   }
  
   const ILArAutoCorr* aCorr=nullptr;
   if(m_dumperTool->doShape()) {
      SG::ReadCondHandle<ILArAutoCorr> acorrHdl(m_acorrKey, ctx);
      aCorr=*acorrHdl;
      if (!aCorr) {
         ATH_MSG_ERROR("Failed to retrieve AutoCorr cond obj");
         return StatusCode::FAILURE;
      }
   }
  
   SG::ReadCondHandle<LArBadChannelCont> readHandle{m_BCKey};
   const LArBadChannelCont *bcCont {*readHandle};
   if(!bcCont) {
      ATH_MSG_ERROR( "Do not have Bad chan container " << m_BCKey.key() );
      return StatusCode::FAILURE;
   }
   
   const LArOFIterResultsContainer* ofIterResult = nullptr;
   if (m_doOFCIter) {
     if (evtStore()->contains<LArOFIterResultsContainer> ("LArOFIterResult")) {
        ATH_CHECK( evtStore()->retrieve(ofIterResult, "LArOFIterResult") );
     } else {
        ATH_MSG_WARNING("Do not have LArOFIterResult in this event");
     }
   }
  
   const LArFebErrorSummary* larFebErrorSummary = nullptr;
   ATH_CHECK( evtStore()->retrieve(larFebErrorSummary, "LArFebErrorSummary") );
   const std::map<unsigned int,uint16_t>& febErrorMap = larFebErrorSummary->get_all_febs();
   std::map<unsigned int, const LArRawChannel*> channelsToKeep;
   
   for (LArRawChannelContainer::const_iterator channel = rawChannelContainer->begin(); 
        channel != rawChannelContainer->end(); ++channel) 
   {
     if (m_energyCut > 0 && TMath::Abs(channel->energy()) < m_energyCut) continue;
     if (m_bcMask.cellShouldBeMasked(bcCont,channel->channelID())) continue;
  
     IdentifierHash hash = m_onlineHelper->channel_Hash(channel->channelID());
     
     if (!hash.is_valid()) {
       ATH_MSG_FATAL ( "Found a LArRawChannel whose HWIdentifier (" << channel->channelID()
                       << ") does not correspond to a valid hash -- returning StatusCode::FAILURE." );
       return StatusCode::FAILURE;
     }    
     channelsToKeep[hash] = &*channel;
     if (m_dumpChannelInfos) {
       HistoryContainer* histCont = m_samples->hist_cont(hash);
       CellInfo* info = nullptr;
       if (!histCont) {
         HWIdentifier channelID = channel->hardwareID();
         const Identifier id = cabling->cnvToIdentifier(channelID);
         const CaloDetDescrElement* caloDetElement = caloMgr->get_element(id);
         info = m_dumperTool->makeCellInfo(channelID, id, caloDetElement);
         if (!info) continue;
         m_samples->makeNewHistory(hash, info);
       }      
     }
   }
  
   std::map<HWIdentifier, LArOFIterResultsContainer::const_iterator> ofcResultPosition;
   if (m_doOFCIter && ofIterResult) { 
     for (LArOFIterResultsContainer::const_iterator ofResult = ofIterResult->begin();
          ofResult != ofIterResult->end(); ++ofResult) 
       ofcResultPosition[ofResult->getChannelID()] = ofResult;
   
     ATH_MSG_DEBUG ( "njpbSizes : " << larDigitContainer->size()
                  << " " << (ofIterResult ? ofIterResult->size() : 0) << " " 
                  << rawChannelContainer->size() << " " << channelsToKeep.size() );
   }
   SG::ReadCondHandle<CaloNoise> noiseHdl{m_noiseCDOKey};
   const CaloNoise* noiseCDO=*noiseHdl;
   
   for (LArDigitContainer::const_iterator digit = larDigitContainer->begin();
        digit != larDigitContainer->end(); ++digit) 
   {    
     //Check Energy selection
     IdentifierHash hash = m_onlineHelper->channel_Hash((*digit)->channelID());
     
     std::map<unsigned int, const LArRawChannel*>::const_iterator findChannel = channelsToKeep.find(hash);
     if (findChannel == channelsToKeep.end()) continue;
     const LArRawChannel* rawChannel = findChannel->second;
  
     //Check detector part
     HWIdentifier channelID = (*digit)->hardwareID();
     if ((m_onlineHelper->isEMBchannel(channelID) || m_onlineHelper->isEMECchannel(channelID)) && !m_doEM) continue;
     if (m_onlineHelper->isHECchannel(channelID) && !m_doHEC) continue;
     if (m_onlineHelper->isFCALchannel(channelID) && !m_doFCAL) continue;
  
     //Check gain
     CaloGain::CaloGain gain=(*digit)->gain();
     if (gain >= CaloGain::LARNGAIN || m_gains[gain] == false) continue;
  
     //Check if connected
     const bool connected = cabling->isOnlineConnected(channelID);
     if (!connected && !m_dumpDisc) continue;
    
     // Check ADCMax selection
     float pedestal = pedestals->pedestal(channelID, gain);
     float pedestalRMS = pedestals->pedestalRMS(channelID, gain);
     if (m_minADCMax > 0 || m_noiseSignifCut > 0) {
       const std::vector<short>& samples = (*digit)->samples();
       double maxValue = -1;
       for (short sample : samples)
         if (sample - pedestal > maxValue) maxValue = sample - pedestal;
       if (m_minADCMax > 0 && fabs(maxValue) < m_minADCMax) continue;
       if (m_noiseSignifCut > 0 && fabs(maxValue) < pedestalRMS*m_noiseSignifCut) continue;
     }
    
     const Identifier id = cabling->cnvToIdentifier(channelID);
     const CaloDetDescrElement* caloDetElement = nullptr;
   
     HistoryContainer* histCont = m_samples->hist_cont(hash);
     CellInfo* info = nullptr;
     if (!histCont) {
       if (!caloDetElement) caloDetElement = caloMgr->get_element(id);
       info = m_dumperTool->makeCellInfo(channelID, id, caloDetElement);
       if (!info) continue;
       histCont = m_samples->makeNewHistory(hash, info);
     }
     else 
       info = histCont->cell_info();
  
     float noise = -1;
     unsigned int status = 0xFFFFFFFF;
     if (connected) {
       if (!caloDetElement) caloDetElement = caloMgr->get_element(id);
       noise = noiseCDO->getNoise(id,gain);
       status = bcCont->status(channelID).packedData();
       HWIdentifier febId = m_onlineHelper->feb_Id(m_onlineHelper->feedthrough_Id(channelID), m_onlineHelper->slot(channelID));
       std::map<unsigned int,uint16_t>::const_iterator findError = febErrorMap.find(febId.get_identifier32().get_compact());
       if (findError != febErrorMap.end()) {
     unsigned int febErrWord = findError->second;
     status = status & (febErrWord >> 17);
       }
     }
  
     //std::vector<float> autoCorr;
     ILArAutoCorr::AutoCorrRef_t autoCorr;
     if (m_dumperTool->doShape()) {
       const LArAutoCorrComplete* autoCorrObj = dynamic_cast<const LArAutoCorrComplete*>(aCorr);
       if (!autoCorrObj)
         ATH_MSG_WARNING ( "AutoCorr object is not of type LArAutoCorrComplete!" );
       else
         autoCorr = autoCorrObj->autoCorr(channelID, gain);
     }
  
     if (!info->shape((*digit)->gain())) // this happens if doAllShapes is off
       info->setShape((*digit)->gain(), m_dumperTool->retrieveShape(channelID, gain));
     
     if (!eventData) {
       eventIndex = makeEvent(eventData, run, event, lumiBlock, bunchId); // this happens if doAllEvents is off
       if (eventIndex < 0) return StatusCode::FAILURE;
     }
     
     DataContainer* data = 
         new DataContainer((*digit)->gain(), (*digit)->samples(),
                              rawChannel->energy(),
                              rawChannel->time()/double(1000),
                              rawChannel->quality(),    
                              eventIndex,
                              autoCorr, 
                  noise, pedestal, pedestalRMS, status);
     
    //  std::map<HWIdentifier, LArOFIterResultsContainer::const_iterator>::const_iterator findResult = ofcResultPosition.find(channelID);    
 //     if (findResult != ofcResultPosition.end()) {
 //       LArOFIterResultsContainer::const_iterator ofResult = findResult->second;
 //       if (ofResult->getValid() && ofResult->getConverged())
 //         data->setADCMax(ofResult->getAmplitude());
 //     }
 //     else
 //       msg() << MSG::INFO << "OFResult for channel 0x" << MSG::hex << channelID << MSG::dec 
 //           << " not found. (size was " << ofcResultPosition.size() << ")" << endmsg;
     
    
     const auto ramp=adc2MeV->ADC2MEV(channelID,gain); //dudu
     data->setADCMax(rawChannel->energy()/ramp[1]); //pow(ADCPeak,i); //dudu
         
  
     histCont->add(data);
   }
   
   if(m_doSC) {
     std::map<unsigned int, const LArRawSC*> scToKeep;
     
     SG::ReadHandle<LArRawSCContainer> hdlSC(m_rawscKey, ctx);
     if(!hdlSC.isValid()) {
        ATH_MSG_WARNING( "Unable to retrieve LArRawSCContainer with key " << m_rawscKey << " from EventStore. " );
        return StatusCode::SUCCESS;
     } else
         ATH_MSG_DEBUG( "Got LArRawSCContainer with key " << m_rawscKey.key() );
     const LArRawSCContainer*  etcontainer = &(*hdlSC);
  
     if (etcontainer->empty()) {
       ATH_MSG_WARNING ( "LArRawSCContainer with key=" << m_rawscKey.key() << " is empty!" );
       return StatusCode::SUCCESS;
     }
  
     SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKeySC, ctx};
     const LArOnOffIdMapping* cablingSC=*cablingHdl;
     if(!cablingSC) {
        ATH_MSG_ERROR( "Do not have cabling object LArOnOffIdMapping" );
        return StatusCode::FAILURE;
     }
  
     SG::ReadCondHandle<CaloSuperCellDetDescrManager> caloMgrHandle{m_caloSuperCellMgrKey};
     ATH_CHECK(caloMgrHandle.isValid());
     const CaloSuperCellDetDescrManager* caloMgrSC = *caloMgrHandle;  
  
     SG::ReadCondHandle<ILArPedestal> pedHdl(m_pedestalKeySC, ctx);
     const ILArPedestal* pedestalsSC=*pedHdl;
     if (!pedestals) {
       ATH_MSG_ERROR("Failed to retrieve pedestal cond obj for SC");
        return StatusCode::FAILURE;
     }
  
     std::map<unsigned int, std::pair<float,float> > channelsToKeepSC;
  
     for (const LArRawSC* rawSC : *hdlSC) {
  
       const std::vector<unsigned short>& bcids = rawSC->bcids();
       const std::vector<int>& energies = rawSC->energies();
       const std::vector<int>& tauenergies = rawSC->tauEnergies();
       const std::vector<bool>& satur = rawSC->satur();
   
       // Look for bcid:
       float scEne = 0;
       float scTim = -99999999.;
   
       const size_t nBCIDs = bcids.size();
       size_t i = 0;
       for (i = 0; i < nBCIDs && bcids[i] != bunchId; i++)
         ;
       if(i==nBCIDs) continue;
       if (satur[i]) continue;
   
       scEne = energies[i]; 
       if (m_energyCutSC > 0 && TMath::Abs(scEne) < m_energyCut) continue;
       if (m_bcMaskSC.cellShouldBeMasked(bcCont,rawSC->hardwareID())) continue;
  
       IdentifierHash hash = m_onlineHelperSC->channel_Hash(rawSC->hardwareID());
       
       if (!hash.is_valid()) {
         ATH_MSG_FATAL ( "Found a LArRawSC whose HWIdentifier (" << rawSC->hardwareID()
                         << ") does not correspond to a valid hash -- returning StatusCode::FAILURE." );
         return StatusCode::FAILURE;
       }    
  
       if(tauenergies.size() && scEne != 0) scTim = tauenergies[i] / scEne;
       channelsToKeepSC[hash] = std::make_pair(scEne, scTim);
  
       if (m_dumpChannelInfos) {
         HistoryContainer* histCont = m_samples->hist_cont_sc(hash);
         CellInfo* info = nullptr;
         if (!histCont) {
           HWIdentifier channelID = rawSC->hardwareID();
           const Identifier id = cablingSC->cnvToIdentifier(channelID);
           const CaloDetDescrElement* caloDetElement = caloMgrSC->get_element(id);
           info = m_dumperToolSC->makeCellInfo(channelID, id, caloDetElement);
           if (!info) continue;
           m_samples->makeNewHistorySC(hash, info);
         }      
       }
     }
  
     
     SG::ReadHandle<LArDigitContainer> hdlSCDigit(m_digitsKeySC, ctx);
     if(!hdlSCDigit.isValid()) {
        ATH_MSG_WARNING( "Unable to retrieve LArDigitContainer with key " << m_digitsKeySC << " from DetectorStore. " );
        return StatusCode::SUCCESS;
     } else
         ATH_MSG_DEBUG( "Got LArDigitContainer with key " << m_digitsKeySC.key() );
     const LArDigitContainer* larSCDigitContainer = &(*hdlSCDigit);
  
     if (larSCDigitContainer->empty()) {
       ATH_MSG_WARNING ( "LArDigitContainer with key=" << m_digitsKeySC.key() << " is empty!" );
       m_nNoDigitsSC++;
       return StatusCode::SUCCESS;
     }
  
     for (LArDigitContainer::const_iterator digit = larSCDigitContainer->begin();
          digit != larSCDigitContainer->end(); ++digit) 
     {    
       //Check Energy selection
       IdentifierHash hash = m_onlineHelperSC->channel_Hash((*digit)->channelID());
       
       std::map<unsigned int, std::pair<float,float> >::const_iterator findChannel = channelsToKeepSC.find(hash);
       if (findChannel == channelsToKeepSC.end()) continue;
  
       // Check ADCMax selection
       float pedestal = pedestalsSC->pedestal((*digit)->channelID(), 0);
       if (m_minADCMaxSC > 0 ) {
         const std::vector<short>& samples = (*digit)->samples();
         double maxValue = -1;
         for (short sample : samples)
           if (sample - pedestal > maxValue) maxValue = sample - pedestal;
         if (fabs(maxValue) < m_minADCMax) continue;
       }
      
       const Identifier id = cablingSC->cnvToIdentifier((*digit)->channelID());
       const CaloDetDescrElement* caloDetElement = nullptr;
     
       HistoryContainer* histCont = m_samples->hist_cont_sc(hash);
       CellInfo* info = nullptr;
       if (!histCont) {
         if (!caloDetElement) caloDetElement = caloMgrSC->get_element(id);
         info = m_dumperToolSC->makeCellInfo((*digit)->channelID(), id, caloDetElement);
         if (!info) continue;
         histCont = m_samples->makeNewHistorySC(hash, info);
       }
       else 
         info = histCont->cell_info();
  
       if (!eventData) {
         eventIndex = makeEvent(eventData, run, event, lumiBlock, bunchId); // this happens if doAllEvents is off
         if (eventIndex < 0) return StatusCode::FAILURE;
       }
       
       DataContainer* data = 
           new DataContainer((*digit)->gain(), (*digit)->samples(),
                                findChannel->second.first,
                                findChannel->second.second/double(1000),
                                -1,    
                                eventIndex,
                                LArVectorProxy(), 
                      -1, pedestal);
       
      
       histCont->add(data);
     }
     
   } //m_doSC
   
   return StatusCode::SUCCESS;
 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

 {
   // If we didn't find any symlinks to add, just return the collection
   // from the base class.  Otherwise, return the extended collection.
   if (!m_extendedExtraObjects.empty()) {
     return m_extendedExtraObjects;
   }
   return Algorithm::extraOutputDeps();
 }

◆ finalize()

StatusCode LArShapeDumper::finalize ( )

overridevirtual

Definition at line 664 of file LArShapeDumper.cxx.

 {
   ATH_MSG_DEBUG ("in finalize() ");
  
   if (m_prescale>1) ATH_MSG_INFO("Prescale dropped " <<  m_nPrescaledAway << "/" << m_count << " events"); 
   if (m_onlyEmptyBC) ATH_MSG_INFO("Dropped " << m_nWrongBunchGroup << "/" << m_count << " events b/c of wrong bunch group"); 
   ATH_MSG_INFO("Dropped " << m_nLArError << "/" <<  m_count << " Events b/c of LAr Veto (Noise burst or corruption)");
  
  
   int n = 0;
   for (unsigned int i = 0; i < m_samples->nChannels(); i++)
     if (m_samples->historyContainer(i)) {
       if (m_samples->historyContainer(i)->cellInfo() == nullptr)
     ATH_MSG_INFO ( "Cell with no cellInfo at index " << i << " !!" );
       //else if (m_samples->historyContainer(i)->cellInfo()->shape() == 0)
     //msg() << MSG::INFO << "Cell with no ShapeInfo at index " << i << " !!" << endmsg;
       //msg() << MSG::INFO << "Non-zero cell at index " << i << " " << m_samples->shape(i)->size() << endmsg;
       n++;
     }
  
   //for (unsigned int i = 0; i < m_samples->nEvents(); i++) {
   //   msg() << MSG::INFO << "Event " << i << " = " 
   //            << m_samples->eventData(i)->run() << " " << m_samples->eventData(i)->event()
   //            << "trigger = " << m_samples->eventData(i)->triggers() << ", nRoIs = " << m_samples->eventData(i)->nRoIs() << endmsg;
   // }
   ATH_MSG_INFO ( "Non-zero cells = " << n << ", footprint = " << m_samples->footprint() );
   ATH_MSG_INFO ( "Writing..." );
  
   if (!m_doStream) {
     m_samples->writeTrees(m_fileName.c_str());
 /*    TFile* f = TFile::Open(m_fileName.c_str(), "RECREATE");
     msg() << MSG::INFO << "Writing (2)..." << endmsg;
     f->WriteObjectAny(m_samples, "Container", "LArSamples");
     msg() << MSG::INFO << "Closing..." << endmsg;
     f->Close();
     msg() << MSG::INFO << "Deleting..." << endmsg;
     delete m_samples;*/
     msg() << MSG::INFO << "Done!" << endmsg;
   }
  
   return StatusCode::SUCCESS;
 }

◆ initialize()

StatusCode LArShapeDumper::initialize ( )

overridevirtual

Get bad-channel mask (only if jO IgnoreBadChannels is true)

Definition at line 86 of file LArShapeDumper.cxx.

 {
   ATH_MSG_DEBUG ("in initialize()");
  
   m_samples = new DataStore();
  
   std::transform(m_caloType.begin(), m_caloType.end(), m_caloType.begin(), toupper);
   m_doEM   = (m_caloType.find("EM")   != std::string::npos);
   m_doHEC  = (m_caloType.find("HEC")  != std::string::npos);
   m_doFCAL = (m_caloType.find("FCAL") != std::string::npos);
   m_doSC   = (m_caloType.find("SC") != std::string::npos);
  
   ATH_CHECK( m_cablingKey.initialize() );
   ATH_CHECK( m_BCKey.initialize() );
   ATH_CHECK( m_BCKeySC.initialize(m_doSC) );
   ATH_CHECK( m_noiseCDOKey.initialize() );
   ATH_CHECK( m_adc2mevKey.initialize() );
   ATH_CHECK( m_pedestalKey.initialize() );
   ATH_CHECK( m_bcDataKey.initialize() );
  
   ATH_CHECK( m_digitsKey.initialize() );
   ATH_CHECK( m_channelsKey.initialize() );
  
   ATH_CHECK( m_digitsKeySC.initialize(m_doSC) );
   ATH_CHECK( m_rawscKey.initialize(m_doSC) );
   ATH_CHECK( m_cablingKeySC.initialize(m_doSC) );
   ATH_CHECK( m_pedestalKeySC.initialize(m_doSC) );
  
   ATH_CHECK( detStore()->retrieve(m_onlineHelper, "LArOnlineID") );
   ATH_CHECK(m_caloMgrKey.initialize());
   ATH_CHECK(m_caloSuperCellMgrKey.initialize(m_doSC));
   if(m_doSC) ATH_CHECK( detStore()->retrieve(m_onlineHelperSC, "LArOnline_SuperCellID") );
  
   ATH_CHECK(m_bcMask.buildBitMask(m_problemsToMask,msg()));
   if(m_doSC) ATH_CHECK(m_bcMaskSC.buildBitMask(m_problemsToMaskSC,msg()));
  
  
   if (m_doTrigger) {
     CHECK( m_trigDec.retrieve() );
   }
  
   ATH_CHECK( m_dumperTool.retrieve() );
  
   ATH_CHECK(m_acorrKey.initialize(m_dumperTool->doShape()) );
   
   std::transform(m_gainSpec.begin(), m_gainSpec.end(), m_gainSpec.begin(), toupper);
   m_gains[CaloGain::LARHIGHGAIN]   = (m_gainSpec.find("HIGH")   != std::string::npos);
   m_gains[CaloGain::LARMEDIUMGAIN] = (m_gainSpec.find("MEDIUM") != std::string::npos);
   m_gains[CaloGain::LARLOWGAIN]    = (m_gainSpec.find("LOW")    != std::string::npos);
   
   //if (m_onlyEmptyBC)
   //ATH_CHECK(m_bcidTool.retrieve());
  
   return StatusCode::SUCCESS; 
 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::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.

◆ makeEvent()

int LArShapeDumper::makeEvent	(	LArSamples::EventData *&	eventData,
		int	run,
		int	event,
		int	lumiBlock,
		int	bunchXing
	)		const

Definition at line 708 of file LArShapeDumper.cxx.

 {
   std::vector<unsigned int> triggerWords;
   if (m_doTrigger) {
     const ROIB::RoIBResult* l1Result = nullptr;
     if (evtStore()->retrieve(l1Result).isFailure() || !l1Result) {
       ATH_MSG_FATAL ( "Could not retrieve RoIBResult!" );
       return -1;
     }
     const std::vector<ROIB::CTPRoI> tav = l1Result->cTPResult().TAV();
     for (const ROIB::CTPRoI& word : tav)
       triggerWords.push_back(word.roIWord());
  
     for (const std::pair<const TString, unsigned int>& p : m_runData->triggerConfig()) {
       while (triggerWords.size() <= p.second/32) triggerWords.push_back(0);
       if (m_trigDec->isPassed(p.first.Data())) {
     triggerWords[p.second/32] |= (0x1 << (p.second % 32));
       //msg() << MSG::INFO << "Trigger line " << p.first.Data() << " passed" << endmsg;
       }
     }
     //msg() << MSG::INFO << "Trigger words : ";
     //for (unsigned int i = 0; i < triggerWords.size(); i++) msg() << MSG::INFO << triggerWords[i] << " ";
     //msg() << MSG::INFO << endmsg;
   }
   
   eventData = new EventData(event, 0, lumiBlock, bunchXing);
   if (m_runData->run() == 0) m_runData->setRun(run);
   eventData->setRunData(m_runData.get());
   eventData->setTriggerData(triggerWords);
   if (m_doRoIs) {
     //msg() << MSG::INFO << "Filling RoI list" << endmsg;
     for (const Trig::ChainGroup* group : m_triggerGroups) {
       std::vector<Trig::Feature<TrigRoiDescriptor> > roIs = group->features().get<TrigRoiDescriptor>();
       for (const Trig::Feature<TrigRoiDescriptor>& roI : roIs) {
     //msg() << MSG::INFO << "Found an roi for chain ";
         //for (unsigned int i = 0; i < group->getListOfTriggers().size(); i++) cout << group->getListOfTriggers()[i] << " ";
         //cout << "@ " << roI.cptr()->eta() << ", " << roI.cptr()->phi() << ", TE = " 
     //   << roI.te()->getId() << " " << Trig::getTEName(*roI.te()) << " with label " << roI.label() << endmsg;
     eventData->addRoI(roI.cptr()->eta(), roI.cptr()->phi(), group->getListOfTriggers()[0].c_str(), roI.label().c_str());
     //msg() << MSG::INFO << "nRoIs so far = " << eventData->nRoIs() << endmsg;
       }
     }
   }
   return m_samples->addEvent(eventData);
 }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< Algorithm >::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< Algorithm > >::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.

◆ 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< Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ start()

StatusCode LArShapeDumper::start ( )

overridevirtual

Definition at line 144 of file LArShapeDumper.cxx.

 {
   m_runData = std::make_unique<RunData>(0);
  
   if (m_doTrigger) {
     std::vector<std::regex> regexs;
     regexs.reserve(m_triggerNames.size());
 for (const std::string& name : m_triggerNames) {
       regexs.push_back(std::regex(name));
     }
  
     std::vector<std::string> chains = m_trigDec->getListOfTriggers();
     std::vector<std::string> myChains;
     std::cmatch match;
  
     for (const std::string& chain : chains) {
       ATH_MSG_INFO ( "Configured chain : " << chain );
       for (const std::regex& regex : regexs) {
         if (std::regex_match(chain.c_str(), match, regex)) myChains.push_back(chain);
       }
     }
     for (const std::string& group : m_trigDec->getListOfGroups())
       ATH_MSG_INFO ( "Configured group : " << group );
     const Trig::ChainGroup* calibStreamGroup = m_trigDec->getChainGroup("Calibration"); 
     if (calibStreamGroup) {
       std::vector<std::string> chains = calibStreamGroup->getListOfTriggers();
       ATH_MSG_INFO ( "Chains for Calibration group:" );
       for (const std::string& chain : chains)
         ATH_MSG_INFO ( "Calib chain : " << chain );
     }
  
     unsigned int idx = 0;
   
     if (m_doAllLvl1) {
       const Trig::ChainGroup* group = m_trigDec->getChainGroup("L1_.*");
       for (const std::string& l1Item : group->getListOfTriggers()) {
         const TrigConf::TriggerItem* confItem = m_trigDec->ExperimentalAndExpertMethods().getItemConfigurationDetails(l1Item);
         if (!confItem) {
           ATH_MSG_WARNING ( "LVL1 item " << l1Item << ", obtained from TrigConfig, cannot be retrieved!" );
           continue;
         }
         int pos = confItem->ctpId();
         if (pos < 0 || pos >= 256) {
           ATH_MSG_WARNING ( "LVL1 item " << l1Item << "has out-of-range ctpId " << pos );
           continue;
         }
         m_runData->addBit(l1Item.c_str(), pos);
         ATH_MSG_INFO ( "Adding LVL1 trigger bit for " << l1Item << " at position " << pos );
       }
       idx = 256;
     }
  
     for (const std::string& name : myChains) {
       if (m_trigDec->getListOfTriggers(name).empty()) {
         ATH_MSG_WARNING ( "Requested trigger name " << name << " is not configured in this run" );
         continue;
       }
       const Trig::ChainGroup* group = m_trigDec->getChainGroup(name);
       if (!group) {
         ATH_MSG_WARNING ( "Could not retrieve chain group for trigger " << name );
         continue;
       }
       m_runData->addBit(name.c_str(), idx++);
       m_triggerGroups.push_back(group);
       ATH_MSG_INFO ( "Adding trigger bit for " << name << " at position " << idx-1 );
     }
   }
   return StatusCode::SUCCESS;
 }

◆ stop()

StatusCode LArShapeDumper::stop ( )

overridevirtual

Definition at line 657 of file LArShapeDumper.cxx.

 {
   m_samples->addRun(m_runData.release());
   return StatusCode::SUCCESS;
 }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

                                        {
   StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
  
   if (sc.isFailure()) {
     return sc;
   }
   ServiceHandle<ICondSvc> cs("CondSvc",name());
   for (auto h : outputHandles()) {
     if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
       // do this inside the loop so we don't create the CondSvc until needed
       if ( cs.retrieve().isFailure() ) {
         ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
         return StatusCode::SUCCESS;
       }
       if (cs->regHandle(this,*h).isFailure()) {
         sc = StatusCode::FAILURE;
         ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                       << " with CondSvc");
       }
     }
   }
   return sc;
 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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_acorrKey

SG::ReadCondHandleKey<ILArAutoCorr> LArShapeDumper::m_acorrKey {this,"AutocorrKey","LArAutoCorr","SG Key of LArAutoCorr object"}

private

Definition at line 137 of file LArShapeDumper.h.

◆ m_adc2mevKey

SG::ReadCondHandleKey<LArADC2MeV> LArShapeDumper::m_adc2mevKey {this,"ADC2MeVKey","LArADC2MeV","SG Key of ADC2MeV conditions object"}

private

Definition at line 119 of file LArShapeDumper.h.

◆ m_bcDataKey

SG::ReadCondHandleKey<BunchCrossingCondData> LArShapeDumper::m_bcDataKey {this, "BunchCrossingCondDataKey", "BunchCrossingData" ,"SG Key of BunchCrossing CDO"}

private

Definition at line 131 of file LArShapeDumper.h.

◆ m_BCKey

SG::ReadCondHandleKey<LArBadChannelCont> LArShapeDumper::m_BCKey {this, "BadChanKey", "LArBadChannel", "SG bad channels key"}

private

Definition at line 126 of file LArShapeDumper.h.

◆ m_BCKeySC

SG::ReadCondHandleKey<LArBadChannelCont> LArShapeDumper::m_BCKeySC {this, "BadChanKeySC", "LArBadChannelSC", "SG bad channels key"}

private

Definition at line 127 of file LArShapeDumper.h.

◆ m_bcMask

LArBadChannelMask LArShapeDumper::m_bcMask

private

Definition at line 107 of file LArShapeDumper.h.

◆ m_bcMaskSC

LArBadChannelMask LArShapeDumper::m_bcMaskSC

private

Definition at line 109 of file LArShapeDumper.h.

◆ m_cablingKey

SG::ReadCondHandleKey<LArOnOffIdMapping> LArShapeDumper::m_cablingKey {this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"}

private

Definition at line 123 of file LArShapeDumper.h.

◆ m_cablingKeySC

SG::ReadCondHandleKey<LArOnOffIdMapping> LArShapeDumper::m_cablingKeySC {this,"CablingKeySC","LArOnOffIdMapSC","SG Key of LArOnOffIdMapping object for SC"}

private

Definition at line 124 of file LArShapeDumper.h.

◆ m_caloMgrKey

SG::ReadCondHandleKey<CaloDetDescrManager> LArShapeDumper::m_caloMgrKey

private

Initial value:

{ this
      , "CaloDetDescrManager"
      , "CaloDetDescrManager"
      , "SG Key for CaloDetDescrManager in the Condition Store" }

Definition at line 138 of file LArShapeDumper.h.

◆ m_caloSuperCellMgrKey

SG::ReadCondHandleKey<CaloSuperCellDetDescrManager> LArShapeDumper::m_caloSuperCellMgrKey

private

Initial value:

{this, 
     "CaloSuperCellDetDescrManager", "CaloSuperCellDetDescrManager", 
     "SG key of the resulting CaloSuperCellDetDescrManager" }

Definition at line 142 of file LArShapeDumper.h.

◆ m_caloType

std::string LArShapeDumper::m_caloType

private

Definition at line 96 of file LArShapeDumper.h.

◆ m_channelsKey

SG::ReadHandleKey<LArRawChannelContainer> LArShapeDumper::m_channelsKey {this, "ChannelsKey", "LArRawChannels", "key for LArRawChannels"}

private

Definition at line 113 of file LArShapeDumper.h.

◆ m_count

unsigned LArShapeDumper::m_count

private

Definition at line 87 of file LArShapeDumper.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_digitsKey

SG::ReadHandleKey<LArDigitContainer> LArShapeDumper::m_digitsKey {this, "DigitsKey", "FREE", "key for LArDigitContainer"}

private

Definition at line 112 of file LArShapeDumper.h.

◆ m_digitsKeySC

SG::ReadHandleKey<LArDigitContainer> LArShapeDumper::m_digitsKeySC {this, "DigitsKeySC", "SC_ADC_BAS", "key for LArDigitContainer for SC"}

private

Definition at line 115 of file LArShapeDumper.h.

◆ m_doAllEvents

bool LArShapeDumper::m_doAllEvents

private

Definition at line 151 of file LArShapeDumper.h.

◆ m_doAllLvl1

bool LArShapeDumper::m_doAllLvl1

private

Definition at line 151 of file LArShapeDumper.h.

◆ m_doEM

bool LArShapeDumper::m_doEM

private

Definition at line 152 of file LArShapeDumper.h.

◆ m_doFCAL

bool LArShapeDumper::m_doFCAL

private

Definition at line 152 of file LArShapeDumper.h.

◆ m_doHEC

bool LArShapeDumper::m_doHEC

private

Definition at line 152 of file LArShapeDumper.h.

◆ m_doOFCIter

bool LArShapeDumper::m_doOFCIter

private

Definition at line 150 of file LArShapeDumper.h.

◆ m_doRoIs

bool LArShapeDumper::m_doRoIs

private

Definition at line 151 of file LArShapeDumper.h.

◆ m_doSC

bool LArShapeDumper::m_doSC

private

Definition at line 152 of file LArShapeDumper.h.

◆ m_doStream

bool LArShapeDumper::m_doStream

private

Definition at line 150 of file LArShapeDumper.h.

◆ m_doTrigger

bool LArShapeDumper::m_doTrigger

private

Definition at line 150 of file LArShapeDumper.h.

◆ m_dumpChannelInfos

bool LArShapeDumper::m_dumpChannelInfos

private

Definition at line 151 of file LArShapeDumper.h.

◆ m_dumpDisc

bool LArShapeDumper::m_dumpDisc

private

Definition at line 101 of file LArShapeDumper.h.

◆ m_dumperTool

ToolHandle<ILArShapeDumperTool> LArShapeDumper::m_dumperTool {this,"LArShapeDumperTool","LArShapeDumperTool"}

private

Definition at line 104 of file LArShapeDumper.h.

◆ m_dumperToolSC

ToolHandle<ILArShapeDumperTool> LArShapeDumper::m_dumperToolSC {this,"LArShapeDumperToolSC","LArShapeDumperToolSC"}

private

Definition at line 105 of file LArShapeDumper.h.

◆ m_energyCut

double LArShapeDumper::m_energyCut

private

Definition at line 97 of file LArShapeDumper.h.

◆ m_energyCutSC

double LArShapeDumper::m_energyCutSC

private

Definition at line 98 of file LArShapeDumper.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_fileName

std::string LArShapeDumper::m_fileName

private

Definition at line 84 of file LArShapeDumper.h.

◆ m_gains

bool LArShapeDumper::m_gains[CaloGain::LARNGAIN] {}

private

Definition at line 153 of file LArShapeDumper.h.

◆ m_gainSpec

std::string LArShapeDumper::m_gainSpec

private

Definition at line 100 of file LArShapeDumper.h.

◆ m_maxChannels

unsigned LArShapeDumper::m_maxChannels

private

Definition at line 88 of file LArShapeDumper.h.

◆ m_minADCMax

int LArShapeDumper::m_minADCMax

private

Definition at line 99 of file LArShapeDumper.h.

◆ m_minADCMaxSC

int LArShapeDumper::m_minADCMaxSC

private

Definition at line 99 of file LArShapeDumper.h.

◆ m_nLArError

unsigned LArShapeDumper::m_nLArError

private

Definition at line 91 of file LArShapeDumper.h.

◆ m_nNoDigits

unsigned LArShapeDumper::m_nNoDigits

private

Definition at line 92 of file LArShapeDumper.h.

◆ m_nNoDigitsSC

unsigned LArShapeDumper::m_nNoDigitsSC

private

Definition at line 93 of file LArShapeDumper.h.

◆ m_noiseCDOKey

SG::ReadCondHandleKey<CaloNoise> LArShapeDumper::m_noiseCDOKey {this,"CaloNoiseKey","totalNoise","SG Key of CaloNoise data object"}

private

Definition at line 129 of file LArShapeDumper.h.

◆ m_noiseSignifCut

double LArShapeDumper::m_noiseSignifCut

private

Definition at line 97 of file LArShapeDumper.h.

◆ m_nPrescaledAway

unsigned LArShapeDumper::m_nPrescaledAway

private

Definition at line 90 of file LArShapeDumper.h.

◆ m_nWrongBunchGroup

unsigned LArShapeDumper::m_nWrongBunchGroup

private

Definition at line 89 of file LArShapeDumper.h.

◆ m_onlineHelper

const LArOnlineID* LArShapeDumper::m_onlineHelper

private

Definition at line 145 of file LArShapeDumper.h.

◆ m_onlineHelperSC

const LArOnline_SuperCellID* LArShapeDumper::m_onlineHelperSC

private

Definition at line 146 of file LArShapeDumper.h.

◆ m_onlyEmptyBC

bool LArShapeDumper::m_onlyEmptyBC

private

Definition at line 155 of file LArShapeDumper.h.

◆ m_pedestalKey

SG::ReadCondHandleKey<ILArPedestal> LArShapeDumper::m_pedestalKey {this,"PedestalKey","LArPedestal","SG Key of LArPedestal object"}

private

Definition at line 134 of file LArShapeDumper.h.

◆ m_pedestalKeySC

SG::ReadCondHandleKey<ILArPedestal> LArShapeDumper::m_pedestalKeySC {this,"PedestalKeySC","LArPedestalSC","SG Key of LArPedestal object"}

private

Definition at line 135 of file LArShapeDumper.h.

◆ m_prescale

int LArShapeDumper::m_prescale

private

Definition at line 86 of file LArShapeDumper.h.

◆ m_problemsToMask

Gaudi::Property<std::vector<std::string> > LArShapeDumper::m_problemsToMask {this,"ProblemsToMask",{}, "Bad-Channel categories to patch"}

private

Definition at line 108 of file LArShapeDumper.h.

◆ m_problemsToMaskSC

Gaudi::Property<std::vector<std::string> > LArShapeDumper::m_problemsToMaskSC {this,"ProblemsToMaskSC",{}, "Bad-Channel categories to mask"}

private

Definition at line 110 of file LArShapeDumper.h.

◆ m_random

TRandom LArShapeDumper::m_random

private

Definition at line 160 of file LArShapeDumper.h.

◆ m_rawscKey

SG::ReadHandleKey<LArRawSCContainer> LArShapeDumper::m_rawscKey {this, "RawSCKey", "SC_ET_ID", "key for LArRawSCContainer"}

private

Definition at line 116 of file LArShapeDumper.h.

◆ m_runData

std::unique_ptr<LArSamples::RunData> LArShapeDumper::m_runData

private

Definition at line 158 of file LArShapeDumper.h.

◆ m_samples

LArSamples::DataStore* LArShapeDumper::m_samples

private

Definition at line 157 of file LArShapeDumper.h.

◆ m_trigDec

PublicToolHandle< Trig::TrigDecisionTool > LArShapeDumper::m_trigDec {this, "TrigDecisionTool", "", "Handle to the TrigDecisionTool"}

private

Definition at line 121 of file LArShapeDumper.h.

◆ m_triggerGroups

std::vector<const Trig::ChainGroup*> LArShapeDumper::m_triggerGroups

private

Definition at line 159 of file LArShapeDumper.h.

◆ m_triggerNames

std::vector<std::string> LArShapeDumper::m_triggerNames

private

Definition at line 102 of file LArShapeDumper.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ LArShapeDumper()

◆ ~LArShapeDumper()

Member Function Documentation

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ makeEvent()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ start()

◆ stop()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_acorrKey

◆ m_adc2mevKey

◆ m_bcDataKey

◆ m_BCKey

◆ m_BCKeySC

◆ m_bcMask

◆ m_bcMaskSC

◆ m_cablingKey

◆ m_cablingKeySC

◆ m_caloMgrKey

◆ m_caloSuperCellMgrKey

◆ m_caloType

◆ m_channelsKey

◆ m_count

◆ m_detStore

◆ m_digitsKey

◆ m_digitsKeySC

◆ m_doAllEvents

◆ m_doAllLvl1

◆ m_doEM

◆ m_doFCAL

◆ m_doHEC

◆ m_doOFCIter

◆ m_doRoIs

◆ m_doSC

◆ m_doStream

◆ m_doTrigger

◆ m_dumpChannelInfos

◆ m_dumpDisc

◆ m_dumperTool

◆ m_dumperToolSC

◆ m_energyCut

◆ m_energyCutSC

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_fileName