#include <LArShapeDumper.h>

Inheritance 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.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
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.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

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::ReadHandleKey< LArRawSCContainer >	m_rawRecomputedscKey {this, "RecomputedSCKey", "SC_ET_RECO", "key for recomputed 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)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
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);
 
  m_bcMaskSC.setSC();
}

◆ ~LArShapeDumper()

LArShapeDumper::~LArShapeDumper ( )

Definition at line 83 of file LArShapeDumper.cxx.

84{

85}

Member Function Documentation

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty	(	Gaudi::Property< T, V, H > &	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());
 
  }

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & 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; }

AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ evtStore()

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

AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ execute()

StatusCode LArShapeDumper::execute ( )

overridevirtual

Definition at line 218 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() );
 
    if (hdlSC->empty()) {
      ATH_MSG_WARNING ( "LArRawSCContainer with key=" << m_rawscKey.key() << " is empty!" );
      return StatusCode::SUCCESS;
    }
 
    SG::ReadHandle<LArRawSCContainer> hdlrecoSC(m_rawRecomputedscKey, ctx);
    if(!hdlrecoSC.isValid()) {
       ATH_MSG_WARNING( "Unable to retrieve LArRawSCContainer with key " << m_rawRecomputedscKey << " from EventStore. " );
       return StatusCode::SUCCESS;
    } else
        ATH_MSG_DEBUG( "Got LArRawSCContainer with key " << m_rawRecomputedscKey.key() );
 
    if (hdlrecoSC->empty()) {
      ATH_MSG_WARNING ( "LArRawSCContainer with key=" << m_rawRecomputedscKey.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, std::pair<float,float> > > channelsToKeepSC;
 
    for (const LArRawSC* rawSC : *hdlSC) { // from bytestream
 
      // is SC masked ?
      if (m_bcMaskSC.cellShouldBeMasked(bcCont,rawSC->hardwareID())) continue;
 
      const std::vector<unsigned short>& bcids = rawSC->bcids();
      const std::vector<int>& energies = rawSC->energies();
      const std::vector<bool>& satur = rawSC->satur();
  
      // Look for bcid:
      float scEne = 0;
      float defValue = -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_energyCutSC) 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;
      }    
 
      channelsToKeepSC[hash] = std::make_pair(scEne, std::make_pair(defValue,defValue));
 
      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);
        }      
      }
    }
    ATH_MSG_INFO("SC to keep in this event: "<<channelsToKeepSC.size());
 
    for (const LArRawSC* rawSC : *hdlrecoSC) { //reconstructed
 
      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;
      }    
 
      // do we keep this SC ?
      if (channelsToKeepSC.find(hash) == channelsToKeepSC.end()) continue;
 
      const std::vector<unsigned short>& bcids = rawSC->bcids();
      const std::vector<int>& energies = rawSC->energies();
      const std::vector<int>& tauenergies = rawSC->tauEnergies();
  
      // 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;
  
      scEne = energies[i]; 
 
      if(tauenergies.size() && scEne != 0) scTim = tauenergies[i] / scEne;
 
      channelsToKeepSC[hash].second = 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, 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.second,
                               findChannel->second.second.first,    
                               eventIndex,
                               LArVectorProxy(), 
                     -1, pedestal, 0, 0);
      
     
      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 726 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 88 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_rawRecomputedscKey.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 770 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()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

◆ 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 147 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 719 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, AthHistogramAlgorithm, and PyAthena::Alg.

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.

137{this,"AutocorrKey","LArAutoCorr","SG Key of LArAutoCorr object"};

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

119{this,"ADC2MeVKey","LArADC2MeV","SG Key of ADC2MeV conditions object"};

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

131{this, "BunchCrossingCondDataKey", "BunchCrossingData" ,"SG Key of BunchCrossing CDO"};

◆ m_BCKey

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

private

Definition at line 126 of file LArShapeDumper.h.

126{this, "BadChanKey", "LArBadChannel", "SG bad channels key"};

◆ m_BCKeySC

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

private

Definition at line 127 of file LArShapeDumper.h.

127{this, "BadChanKeySC", "LArBadChannelSC", "SG bad channels key"};

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

123{this,"CablingKey","LArOnOffIdMap","SG Key of LArOnOffIdMapping object"};

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

124{this,"CablingKeySC","LArOnOffIdMapSC","SG Key of LArOnOffIdMapping object for SC"};

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

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

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

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

◆ 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 112 of file LArShapeDumper.h.

112{this, "ChannelsKey", "LArRawChannels", "key for LArRawChannels"};

◆ 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 111 of file LArShapeDumper.h.

111{this, "DigitsKey", "FREE", "key for LArDigitContainer"};

◆ m_digitsKeySC

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

private

Definition at line 114 of file LArShapeDumper.h.

114{this, "DigitsKeySC", "SC_ADC_BAS", "key for LArDigitContainer for SC"};

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

104{this,"LArShapeDumperTool","LArShapeDumperTool"};

◆ m_dumperToolSC

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

private

Definition at line 105 of file LArShapeDumper.h.

105{this,"LArShapeDumperToolSC","LArShapeDumperToolSC"};

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

153{};

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

129{this,"CaloNoiseKey","totalNoise","SG Key of CaloNoise data object"};

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

134{this,"PedestalKey","LArPedestal","SG Key of LArPedestal object"};

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

135{this,"PedestalKeySC","LArPedestalSC","SG Key of LArPedestal object"};

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

108{this,"ProblemsToMask",{}, "Bad-Channel categories to patch"};

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

110{this,"ProblemsToMaskSC",{}, "Bad-Channel categories to mask"};

◆ m_random

TRandom LArShapeDumper::m_random

private

Definition at line 160 of file LArShapeDumper.h.

◆ m_rawRecomputedscKey

SG::ReadHandleKey<LArRawSCContainer> LArShapeDumper::m_rawRecomputedscKey {this, "RecomputedSCKey", "SC_ET_RECO", "key for recomputed LArRawSCContainer"}

private

Definition at line 116 of file LArShapeDumper.h.

116{this, "RecomputedSCKey", "SC_ET_RECO", "key for recomputed LArRawSCContainer"};

◆ m_rawscKey

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

private

Definition at line 115 of file LArShapeDumper.h.

115{this, "RawSCKey", "SC_ET_ID", "key for LArRawSCContainer"};

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

121{this, "TrigDecisionTool", "", "Handle to the TrigDecisionTool"};

◆ 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()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ finalize()

◆ initialize()

◆ inputHandles()

◆ makeEvent()

◆ msg()

◆ 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

◆ m_gains

◆ m_gainSpec

◆ m_maxChannels

◆ m_minADCMax

◆ m_minADCMaxSC

◆ m_nLArError

◆ m_nNoDigits

◆ m_nNoDigitsSC

◆ m_noiseCDOKey

◆ m_noiseSignifCut