TBECLArRawChannelBuilder Class Reference

#include <TBECLArRawChannelBuilder.h>

Inheritance diagram for TBECLArRawChannelBuilder:

Collaboration diagram for TBECLArRawChannelBuilder:

Public Member Functions
	TBECLArRawChannelBuilder (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) override
	Execute method.
virtual StatusCode	finalize () override
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
	Get filter decision:
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
	Set filter decision:
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
virtual bool	isReEntrant () const override final
	Legacy algorithms are not thread-safe.
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
const LArOnlineID *	m_onlineHelper
const CaloCell_ID *	m_calo_id
SG::ReadCondHandleKey< LArOnOffIdMapping >	m_cablingKey {this,"keyCabling", "LArOnOffIdMap", "Input key for Id mapping"}
SG::ReadCondHandleKey< ILArHVScaleCorr >	m_offlineHVScaleCorrKey {this, "keyOfflineHVCorr", "LArHVScaleCorrRecomputed","Key for LArHVScaleCorr"}
std::string	m_DataLocation
std::string	m_ChannelContainerName
bool	m_useTDC
bool	m_useRamp
bool	m_useShape
bool	m_ConvertADCToHighGain
float	m_Ecut
int	m_initialTimeSampleShift
int	m_NOFCPhases
int	m_NOFCTimeBins
bool	m_useIntercept [4]
bool	m_useOFCPhase
bool	m_phaseInv
float	m_ramp_max [3]
double	m_SamplingPeriode
float	m_OFCTimeBin
bool	m_binHalfOffset
bool	m_allowTimeJump
int	m_pedestalFallbackMode
unsigned int	m_iPedestal
unsigned int	m_shapeMode
unsigned	m_skipSaturCells
short	m_AdcMax
int	m_noEnergy
int	m_noTime
int	m_noShape
int	m_noShapeDer
int	m_saturation
int	m_lastNoEnergy
int	m_lastNoTime
int	m_lastNoShape
int	m_lastNoShapeDer
float	m_aveNoEnergy
float	m_aveNoTime
float	m_aveNoShape
float	m_aveNoShapeDer
float	m_aveSaturCells
int	m_nEvents
float	m_aveChannels
bool	m_hvcorr
double	m_SamplingPeriodeUpperLimit
double	m_SamplingPeriodeLowerLimit
const LArEM_ID *	m_emId
float	m_adc2mev [30]
SG::ReadCondHandleKey< ILArOFC >	m_ofcKey {this, "OFCKey", "LArOFC", "SG Key of OFC conditions object" }
	Property: OFC coefficients (conditions input).
SG::ReadCondHandleKey< LArADC2MeV >	m_adc2mevKey { this, "ADC2MeVKey", "LArADC2MeV", "SG Key of the LArADC2MeV CDO" }
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
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

Definition at line 25 of file TBECLArRawChannelBuilder.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

TBECLArRawChannelBuilder()

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

Definition at line 26 of file TBECLArRawChannelBuilder.cxx.

                                                                                                  :
  AthAlgorithm(name, pSvcLocator),
  m_onlineHelper(0),
  m_calo_id(0),
  m_DataLocation("FREE"),
  m_ChannelContainerName("LArRawChannels"),
  m_useTDC(false),
  m_useRamp(true),
  m_useShape(true),
  m_ConvertADCToHighGain(false),
  m_Ecut(256*MeV),
  m_initialTimeSampleShift(0),
  m_ramp_max(),
  m_noEnergy(0),
  m_noTime(0),
  m_noShape(0),
  m_noShapeDer(0),
  m_saturation(0),
  m_lastNoEnergy(0),
  m_lastNoTime(0),
  m_lastNoShape(0),
  m_lastNoShapeDer(0),
  m_aveNoEnergy(0),
  m_aveNoTime(0),
  m_aveNoShape(0),
  m_aveNoShapeDer(0),
  m_aveSaturCells(0),
  m_nEvents(0),
  m_aveChannels(0),
  m_SamplingPeriodeUpperLimit(0),
  m_SamplingPeriodeLowerLimit(0),
  m_emId(0),
  m_adc2mev()
  //m_NOFCPhases(24),
  //m_NOFCTimeBins(24)
 {
   //m_useIntercept={false,false,false,false};
 declareProperty("LArRawChannelContainerName",m_ChannelContainerName);
 declareProperty("DataLocation",              m_DataLocation );
 declareProperty("UseTDC",                    m_useTDC);
 declareProperty("UseRamp",m_useRamp);
 declareProperty("UseShape",m_useShape);
 declareProperty("ConvertADCToHighGain",m_ConvertADCToHighGain);
 declareProperty("Ecut",                      m_Ecut);
 declareProperty("UseHighGainRampIntercept",  m_useIntercept[CaloGain::LARHIGHGAIN]=false);
 declareProperty("UseMedGainRampIntercept",   m_useIntercept[CaloGain::LARMEDIUMGAIN]=false);
 declareProperty("UseLowGainRampIntercept",   m_useIntercept[CaloGain::LARLOWGAIN]=false);
 declareProperty("InitialTimeSampleShift",    m_initialTimeSampleShift);
 declareProperty("NOFCTimeBins",              m_NOFCTimeBins=25); //Number of OFC time bins in a sampling periode
 m_NOFCPhases = m_NOFCTimeBins; 
 declareProperty("NOFCPhases",                m_NOFCPhases);      //Total number of available OFC sets
 declareProperty("UseOFCPhase",               m_useOFCPhase=false);
 declareProperty("PhaseInversion",            m_phaseInv=false);
 declareProperty("SamplingPeriod",            m_SamplingPeriode=1/(40.08*megahertz));
 declareProperty("OFCTimeBin",                m_OFCTimeBin=m_SamplingPeriode/m_NOFCTimeBins);
 declareProperty("BinHalfOffset",             m_binHalfOffset=false);
 declareProperty("AllowTimeSampleJump",       m_allowTimeJump=true);
 declareProperty("PedestalFallbackMode",      m_pedestalFallbackMode=0); // 0=only DB, 1=Only if missing,
 declareProperty("PedestalSample",            m_iPedestal=0);            // 2=Low, 3=Low+Me dium, 4=All LAr
 declareProperty("ShapeMode",                 m_shapeMode=0); 
 declareProperty("SkipSaturCellsMode",        m_skipSaturCells=0);
 declareProperty("ADCMax",                    m_AdcMax=4095);
 declareProperty("HVcorr",                    m_hvcorr=false);
}

Member Function Documentation

declareGaudiProperty()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode TBECLArRawChannelBuilder::execute ( const EventContext & ctx )

overridevirtual

Execute method.

Implements AthAlgorithm.

Definition at line 180 of file TBECLArRawChannelBuilder.cxx.

{
 
  //Counters for errors & warnings per event
  int noEnergy   = 0; // Number of completly failed channels in a given event
  int BadTiming  = 0; // Number of channels with bad timing in a given event
  int noTime     = 0; // Number of channels without time info in a given event
  int noShape    = 0; // Number of channels without Shape (= with no quality factor) in a given event 
  int noShapeDer = 0; // Number of channels without ShapeDerivative in a given event 
  int highE      = 0; // Number of channels with 'high' (above threshold) energy in a given event 
  int saturation = 0; // Number of saturating channels in a given event   
  
  const ILArHVScaleCorr *oflHVCorr=nullptr;
  if(m_hvcorr) {
     SG::ReadCondHandle<ILArHVScaleCorr> oflHVCorrHdl(m_offlineHVScaleCorrKey, ctx);
     oflHVCorr = *oflHVCorrHdl;
     if(!oflHVCorr) {
        ATH_MSG_ERROR( "Could not get the HVScaleCorr from key " << m_offlineHVScaleCorrKey.key() );
        return StatusCode::FAILURE;
     }
  }
  SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey, ctx};
  const LArOnOffIdMapping* cabling{*cablingHdl};
  if(!cabling) {
        ATH_MSG_ERROR( "Could not get the cabling mapping from key " << m_cablingKey.key() );
        return StatusCode::FAILURE;
     }
 
  //Pointer to input data container
  const LArDigitContainer* digitContainer=NULL;//Pointer to LArDigitContainer
  //const TBPhase* theTBPhase; //Pointer to Testbeam TDC-Phase object (if needed)
  float PhaseTime=0;                //Testbeam TDC phase (if needed)
  float globalTimeOffset=0;
  //Pointer to conditions data objects 
  const ILArFEBTimeOffset* larFebTimeOffset=NULL;
  const ILArShape* larShape=NULL;
  //Retrieve Digit Container
 
  ATH_CHECK( evtStore()->retrieve(digitContainer,m_DataLocation) );
  
  //Retrieve calibration data
  const ILArPedestal* larPedestal=nullptr;
  ATH_CHECK( detStore()->retrieve(larPedestal) );
  
  if (m_useShape) {
    ATH_MSG_DEBUG ( "Retrieving LArShape object" );
    StatusCode sc=detStore()->retrieve(larShape);
    if (sc.isFailure()) {
      ATH_MSG_WARNING ( "Can't retrieve LArShape from Conditions Store" << std::endl
                        << "Quality factor will not be caluclated." );
      larShape=NULL;
    }
  }
 
  ATH_MSG_DEBUG ( "Retrieving LArOFC object" );
  SG::ReadCondHandle<ILArOFC> larOFC (m_ofcKey, ctx);
 
  const LArADC2MeV* adc2mev = nullptr;
  if (m_useRamp) {
    SG::ReadCondHandle<LArADC2MeV> adc2mevH (m_adc2mevKey, ctx);
    adc2mev = *adc2mevH;
  }
 
  //retrieve TDC
  if (m_useTDC) { //All this timing business is only necessary if the readout and the beam are not in phase (Testbeam)
    const TBPhase* theTBPhase = nullptr;
    const ILArGlobalTimeOffset* larGlobalTimeOffset = nullptr;
    ATH_CHECK( evtStore()->retrieve(theTBPhase,"TBPhase") );
    //Get Phase in nanoseconds
    PhaseTime = theTBPhase->getPhase();
    // ###
    if (m_phaseInv) PhaseTime = m_SamplingPeriode - PhaseTime ;
    ATH_MSG_DEBUG ( " *** Phase = " << PhaseTime );
    // ###
    
    //Get Global Time Offset
    StatusCode sc=detStore()->retrieve(larGlobalTimeOffset);
    if (sc.isSuccess()) globalTimeOffset = larGlobalTimeOffset->TimeOffset();
 
    //Get FEB time offset
    sc=detStore()->retrieve(larFebTimeOffset);
    if (sc.isFailure()) larFebTimeOffset=NULL;
  }
 
 
  LArRawChannelContainer* larRawChannelContainer=new LArRawChannelContainer();
  larRawChannelContainer->reserve(digitContainer->size());
  StatusCode sc = evtStore()->record(larRawChannelContainer,m_ChannelContainerName);
  if(sc.isFailure()) {
    ATH_MSG_ERROR ( "Can't record LArRawChannelContainer in StoreGate" );
  }
 
  // Average number of LArDigits per event
  m_nEvents++;
  m_aveChannels += digitContainer->size();
 
  bool debugPrint=false;
  if (msgLvl(MSG::DEBUG) ) debugPrint=true;
 
  // Now all data is available, start loop over Digit Container
  int ntot_raw=0;
 
  for (const LArDigit* digit : *(digitContainer)) {
 
    //Data that goes into RawChannel:
    float energy=0;
    float time=0;
    float quality=0;
 
    int OFCTimeBin=0;
    int timeSampleShift=m_initialTimeSampleShift;
 
    //Get data from LArDigit
    const std::vector<short>& samples=digit->samples();
    const unsigned nSamples=samples.size(); 
    const HWIdentifier chid=digit->channelID();
    const CaloGain::CaloGain gain=digit->gain();
    
    // to be used in case of DEBUG output
    int layer  = -99999 ;
    int eta    = -99999 ; 
    int phi    = -99999 ; 
    int region = -99999 ;    
    if (msgLvl(MSG::DEBUG) ) {
      Identifier id ;
      try {
        id = cabling->cnvToIdentifier(chid);
      } catch ( LArID_Exception & except ) {
        ATH_MSG_DEBUG ( "A Cabling exception was caught for channel 0x!" 
                        << MSG::hex << chid.get_compact() << MSG::dec );
        continue ;
      }
      layer  = m_emId->sampling(id);
      eta    = m_emId->eta(id); 
      phi    = m_emId->phi(id);
      region = m_emId->region(id);    
      ATH_MSG_VERBOSE ( "Channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
                        << " [ Layer = " << layer << " - Eta = " << eta 
                        << " - Phi = " << phi << " - Region = " << region << " ] " );
    }
        
    // check for saturation, in case skip channel
    int nSatur=-1 ;
    for (unsigned iSample=0;iSample<samples.size();iSample++) {
      if (samples[iSample]>=m_AdcMax) {
        nSatur++;
        break ;
      }
    }
    if ( nSatur>-1 ) {
      msg() << MSG::DEBUG << "Saturation on channel 0x" << MSG::hex << chid.get_compact() << MSG::dec ;         
      saturation++;
    }
    if ( m_skipSaturCells && nSatur>-1 ) {
      msg() << ". Skipping channel." << endmsg; 
      continue; // Ignore this cell, saturation on at least one sample
    } else if ( nSatur>-1 ) {
      msg() << "." << endmsg;
    }   
    
    //Get conditions data for this channel:
 
    // Pedestal
    float pedestal=larPedestal->pedestal(chid,gain);
 
    float pedestalAverage;
    if (pedestal < (1.0+LArElecCalib::ERRORCODE)) {
      if( m_pedestalFallbackMode >= 1 ) {
        ATH_MSG_DEBUG ( "No pedestal found for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
                        << " Gain " << gain <<".  Using time sample " << m_iPedestal );
        pedestalAverage=samples[m_iPedestal];
      } else {              
    ATH_MSG_DEBUG ( noEnergy << ". No pedestal found for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
            << " [ Layer = " << layer << " - Eta = " << eta << " - Phi = " << phi << " - Region = " << region << " ]"
                        << " Gain = " << gain << ". Skipping channel." );
    noEnergy++;
        continue;
      }
    } else {
      if( ( m_pedestalFallbackMode>=2 && gain==CaloGain::LARLOWGAIN )   ||
          ( m_pedestalFallbackMode>=3 && gain==CaloGain::LARMEDIUMGAIN ) ||
          ( m_pedestalFallbackMode>=4 && gain==CaloGain::LARHIGHGAIN )       ) {
        ATH_MSG_DEBUG ( "Forcing pedestal fallback for  channel 0x" << MSG::hex << chid.get_compact()
                        << MSG::dec   << " Gain=" << gain << ". Using time sample " << m_iPedestal );
        pedestalAverage=samples[m_iPedestal];
      } else {
        pedestalAverage=pedestal;
      }      
    }
 
    // Optimal Filtering Coefficients
    ILArOFC::OFCRef_t ofc_a;
    ILArOFC::OFCRef_t ofc_b;
    {// get OFC from Conditions Store
      float febTimeOffset=0;
      const HWIdentifier febid=m_onlineHelper->feb_Id(chid);
      if (larFebTimeOffset)
    febTimeOffset=larFebTimeOffset->TimeOffset(febid);
      double timeShift=PhaseTime+globalTimeOffset+febTimeOffset;
      if (debugPrint)
        msg() << MSG::VERBOSE << "Channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
       << " phase=" << PhaseTime  << " Feb=" << febTimeOffset 
       << " Global=" << globalTimeOffset;
 
      if (m_useOFCPhase) {
    const double ofcTimeOffset=larOFC->timeOffset(chid,gain);
    timeShift+=ofcTimeOffset;
    if (debugPrint) msg() << MSG::VERBOSE << " OFC=" << ofcTimeOffset;
      }
 
      if (debugPrint) msg() << MSG::VERBOSE << " Total=" << timeShift << endmsg;
      
      if (m_allowTimeJump && timeShift >= m_NOFCPhases*m_OFCTimeBin ) {
    if (debugPrint) ATH_MSG_VERBOSE ( "Time Sample jump: -1" );
    timeSampleShift -= 1;
    //timeShift       -= m_NOFCTimeBins*m_OFCTimeBin ;
    timeShift       -= m_SamplingPeriode ;
      }
      else if (m_allowTimeJump && timeShift < 0 ) {
    if (debugPrint) ATH_MSG_VERBOSE ( "Time Sample jump: +1" );
    timeSampleShift += 1;
    //timeShift       += m_NOFCTimeBins*m_OFCTimeBin ;
        timeShift       += m_SamplingPeriode ;
      }
 
      if (m_allowTimeJump && ( timeShift > m_NOFCPhases*m_OFCTimeBin || timeShift < 0 ) ) {
    BadTiming++;
    noEnergy++;
    ATH_MSG_ERROR ( noEnergy << ". Time offset out of range for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
                        << " Found " << timeShift << ",  expected ( 0 - " << m_NOFCPhases*m_OFCTimeBin << ") ns. Skipping channel." );
    continue;
      }
      
      if (m_allowTimeJump && timeSampleShift < 0) {
    BadTiming++;
    noEnergy++;
    ATH_MSG_ERROR ( noEnergy << ". Negative time sample (" << timeSampleShift << ") shift for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
                        << " Found. Skipping channel." );
    continue;
      }
 
      OFCTimeBin = (int) ( timeShift / m_OFCTimeBin );
      
      if ( !m_phaseInv ) // if not done with PhaseTime at the beginning, invert time bin for OFC bin selection
         OFCTimeBin = ( m_NOFCTimeBins - 1 ) - OFCTimeBin;
      // do not use the following: 24<PhaseTime<25 you always get OFCTimeBin = -1!
      //else 
      //   OFCTimeBin -= 1 ; 
      
      if (debugPrint) ATH_MSG_VERBOSE ( "OFC bin width = " << m_OFCTimeBin << " - OFCBin = " << OFCTimeBin << " - timeShift = " << timeShift );
      
      if ( OFCTimeBin < 0 ) {
        ATH_MSG_ERROR ( "Channel " << MSG::hex << chid.get_compact() << MSG::dec << " asks for OFC bin = " << OFCTimeBin << ". Set to 0." );
        OFCTimeBin=0;
      } else if ( OFCTimeBin >= m_NOFCPhases ) {
        ATH_MSG_ERROR ( "Channel " << MSG::hex << chid.get_compact() << MSG::dec << " asks for OFC bin = " << OFCTimeBin << ". Set to (NOFCPhases-1) =" << m_NOFCTimeBins-1 );
        OFCTimeBin = m_NOFCPhases-1;
      }
 
      ofc_a=larOFC->OFC_a(chid,gain,OFCTimeBin);
      //ofc_b=&(larOFC->OFC_b(chid,gain,OFCTimeBin)); retrieve only when needed
    }
    
    //Check if we have OFC for this channel and time bin
    if (ofc_a.size()==0) {
      noEnergy++;
      ATH_MSG_DEBUG ( noEnergy << ". No OFC's found for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
          << " [ Layer = " << layer << " - Eta = " << eta << " - Phi = " << phi << " - Region = " << region << " ]"
                      << " Time bin = " << OFCTimeBin << ", Gain = " << gain << ". Skipping channel." );
      continue;
    } 
    if (ofc_a.size()+timeSampleShift>nSamples) {
      BadTiming++;
      noEnergy++;
      if (timeSampleShift==0)
    ATH_MSG_DEBUG ( "Found LArDigit with " << nSamples << " samples, but OFCs for " 
                        << ofc_a.size() << " samples. Skipping Channel ");
      else //have time sample shift
    ATH_MSG_DEBUG ( "After time sample shift of " << timeSampleShift << ", " << nSamples-timeSampleShift
                        << " samples left, but have OFCs for " << ofc_a.size() << " samples. Skipping Channel ");
      continue;
    } 
 
    //Now apply Optimal Filtering to get ADC peak
    float ADCPeak=0;
    for (unsigned i=0;i<(ofc_a.size());i++) 
      ADCPeak+=(samples[i+timeSampleShift]-pedestalAverage)*ofc_a.at(i);
    
    if (debugPrint) ATH_MSG_VERBOSE ( "ADC Height calculated " << ADCPeak << " TimeBin=" << OFCTimeBin  );
     
    if (m_useRamp) {
      //ADC2MeV (a.k.a. Ramp)   
      LArVectorProxy ramp = adc2mev->ADC2MEV(chid,gain);
      //Check ramp coefficents
      if (ramp.size()==0) {
    noEnergy++;
    ATH_MSG_DEBUG ( noEnergy << ". No ADC2MeV data found for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec
        << " [ Layer = " << layer << " - Eta = " << eta << " - Phi = " << phi << " - Region = " << region << " ]"
                        << " Gain = " << gain << ". Skipping channel." );
    continue;
      } 
      
      // temporary fix for bad ramps... should be done in the DB
      if(ramp[1]>m_ramp_max[gain] || ramp[1]<0) {
    noEnergy++;
    ATH_MSG_DEBUG ( "Bad ramp for channel " << chid << " (ramp[1] = " << ramp[1] << "): skip this channel" );
    continue;
      } 
      
      float ADCPeakPower=ADCPeak;
      
      if (m_useIntercept[gain])
    energy=ramp[0];
      //otherwise ignore intercept, E=0;
      for (unsigned i=1;i<ramp.size();i++)
    {energy+=ramp[i]*ADCPeakPower; //pow(ADCPeak,i);
    ADCPeakPower*=ADCPeak;
    }
    } else {
      energy = ADCPeak;
      if (m_ConvertADCToHighGain && gain == CaloGain::LARMEDIUMGAIN) 
    energy *= 9.5;
      Identifier id = cabling->cnvToIdentifier(chid);
      int is = m_calo_id->calo_sample(id);
      energy *= m_adc2mev[is];                // Ramp for h.g. scale
    }
    
// HV correction
 
    if (m_hvcorr) {
// HV tool
       float hvCorr = oflHVCorr-> HVScaleCorr(chid);
       energy = energy*hvCorr;
    }
  
    //Check if energy is above threshold for time & quality calculation
    if (energy>m_Ecut) {
      highE++;
      ofc_b=larOFC->OFC_b(chid,gain,OFCTimeBin);
      if (ofc_b.size() != ofc_a.size()) {//don't have proper number of coefficients
    if (ofc_b.size()==0)
      ATH_MSG_DEBUG ( "No time-OFC's found for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
                          << " Gain "<< gain << " found. Time not calculated." );
    else
      ATH_MSG_DEBUG ( "OFC for time size " << ofc_b.size() 
              << " not equal to OFC for energy size " << ofc_a.size() 
                          << "   Time not calculated " );
    noTime++;
      }else{
    for (unsigned i=0;i<(ofc_b.size());i++) 
      time+=(samples[i+timeSampleShift]-pedestalAverage)*ofc_b.at(i);
    time/=ADCPeak;
    // !! Time is now in ns with respect to calibration pulse shape
    // Used to calculate quality factor
      }
      if (debugPrint) ATH_MSG_VERBOSE ( "Time calculated " << time << " TimeBin=" << OFCTimeBin  );
 
      //Calculate Quality factor
      if (larShape) { //Have shape object
    //Get Shape & Shape Derivative for this channel
 
    //const std::vector<float>& shape=larShape->Shape(chid,gain,OFCTimeBin);
    //const std::vector<float>& shapeDer=larShape->ShapeDer(chid,gain,OFCTimeBin);
        // ###
        ILArShape::ShapeRef_t shape=larShape->Shape(chid,gain,OFCTimeBin,m_shapeMode);
        ILArShape::ShapeRef_t shapeDer=larShape->ShapeDer(chid,gain,OFCTimeBin,m_shapeMode);
        // ###
    
    //Check Shape
        if (shape.size() < ofc_a.size()) {
      if (shape.size()==0) 
        ATH_MSG_DEBUG ( "No Shape found for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
                            << " Gain "<< gain << ". Quality factor not calculated." );
          else
            ATH_MSG_DEBUG ( "Shape size " << shape.size() 
                << "smaller than OFC size " << ofc_a.size() 
                << "for channel 0x" << MSG::hex << chid.get_compact() 
                            << MSG::dec << ". Quality factor not calculated." );
      quality=0;  //Can't calculate chi^2, assume good hit.
      noShape++;
    }
    else {//Shape ok
      if (time!=0 && shapeDer.size()!=shape.size()) { 
        //Send warning
        ATH_MSG_DEBUG ( "Shape-Derivative has different size than Shape for channel 0x" << MSG::hex << chid.get_compact() << MSG::dec 
                            << ". Derivative not taken into accout for quality factor." );
        noShapeDer++;
      }//end-if 
      if (time==0 || shapeDer.size()!=shape.size() ) { //Calculate Q without time info
        for (unsigned i=0;i<(ofc_a.size());i++)
          quality+=((samples[i+timeSampleShift]-pedestalAverage)-shape[i]*ADCPeak)*
               ((samples[i+timeSampleShift]-pedestalAverage)-shape[i]*ADCPeak);
      }
      else { //All input data ok, calculate Q with time info
        for (unsigned i=0;i<(ofc_a.size());i++) 
          quality+=((samples[i+timeSampleShift]-pedestalAverage)-((shape[i]-shapeDer[i]*time)*ADCPeak))*
               ((samples[i+timeSampleShift]-pedestalAverage)-((shape[i]-shapeDer[i]*time)*ADCPeak));
      }
    } // end else (Shape ok)
      } //end if larShape
      else { //No Shape found at all 
    quality=0; //Can't calculate chi^2, assume good hit.
    noShape++;
      }
    //   if (m_useTDC) //Correct time according to EMTB definition (do we really want this?)
    //  time= -time+24.5-tbin;
    }// end-if energy>Ecut
    else 
      quality=-1; //in case E<Ecut
    //time*=1000.0; 
    time=time/picosecond; //Convert time to ps
    //Make LArRawChannel Object with new data
 
    uint16_t iquality=0; 
    uint16_t iprovenance=0xA5; 
    if (quality>=0) {  
       iquality = ((int)(quality) ) & 0xFFFF; 
       iprovenance=iprovenance | 0x2000; 
    }
 
 
    LArRawChannel larRawChannel(chid,(int)energy,(int)time,iquality,iprovenance, gain);
    larRawChannelContainer->push_back(larRawChannel); //Add to container 
    ntot_raw++;
    if (debugPrint)
      ATH_MSG_VERBOSE ( "Got LArRawChannel #" << ntot_raw << ", chid=0x" << MSG::hex << chid.get_compact() << MSG::dec  
                        << " e=" << energy << " t=" << time << " Q=" << quality );
  } // End loop over LArDigits
 
  ATH_MSG_DEBUG (  ntot_raw << " channels successfully processed, (" << highE << " with high energy)" );
 
  // deal with bad timing
  if(BadTiming>=128){
    ATH_MSG_ERROR ( "Too many channels (" <<BadTiming<<  " !) have a bad timing !!" );
    ATH_MSG_ERROR ( "OFC time constants should be revisited !!!" );
    ATH_MSG_ERROR ( "Event is skipped" );
    larRawChannelContainer->clear();
    //return StatusCode::SUCCESS;
  }
  
  // in case of at least one saturating cell, skip all event (if selected) 
  if ( saturation && m_skipSaturCells == 2 ) {
    ATH_MSG_ERROR ( saturation << " saturating channels were found. Event is skipped." );
    larRawChannelContainer->clear();
  }
  
  //Put this LArRawChannel container in the transient store
  //sc = evtStore()->record(m_larRawChannelContainer, m_ChannelContainerName);
  //if(sc.isFailure()) {
  // log << MSG::ERROR << "Can't record LArRawChannelContainer in StoreGate" << endmsg;
  //}
  //else
  //  std::cout << "Successfully recorded LArRawChannelContainer to StoreGate" << std::endl;
  
  /*
   
   Error & Warning summary *per event*
   
   Strategy: 'No Energy' is an ERROR, no time or no quality is a WARNING
   
   Missing calibration constants are most likly missing for an entire run, threfore:
   In DEBUG: Print summary for each event if something is missing
   otherwise: Print summary only for new problems (different number of missing channels) 
  
   Saturatin cells summary is shown in any case, WARNING if not skipped, ERROR if skipped
   
  */
  
  if (noEnergy)   m_noEnergy++;
  if (noTime)     m_noTime++;
  if (noShape)    m_noShape++;
  if (noShapeDer) m_noShapeDer++;
  if (saturation) m_saturation++;
  
  m_aveNoEnergy   += noEnergy;  
  m_aveNoTime     += noTime;
  m_aveNoShape    += noShape;
  m_aveNoShapeDer += noShapeDer;
  m_aveSaturCells += saturation;
  
  if ( (   noEnergy!=m_lastNoEnergy 
        || noTime!=m_lastNoTime 
        || noShape>m_lastNoShape 
        || noShapeDer>m_lastNoShapeDer 
        || saturation>0 ) 
       || ( msgSvc()->outputLevel(name()) <= MSG::DEBUG && ( noEnergy || noTime || noShape || noShapeDer || saturation ) )
      ) {
    
    m_lastNoEnergy = noEnergy;
    m_lastNoTime   = noTime;
    if (noShape>m_lastNoShape) m_lastNoShape=noShape;
    if (noShapeDer>m_lastNoShapeDer) m_lastNoShapeDer=noShapeDer;
    //m_lastSaturCells = saturation ;
 
    MSG::Level msglvl;
    if (noEnergy) 
      msglvl=MSG::ERROR;
    else
      msglvl=MSG::WARNING;
    msg() << msglvl << " *** Error & Warning summary for this event *** " << std::endl;
    
    if ( noEnergy ) {
      msg() << msglvl << "   " << noEnergy << " out of " 
      << digitContainer->size() << " channel(s) skipped due to a lack of basic calibration constants." 
      << std::endl;       
    }
    if ( noTime ) {
      msg() << msglvl << "   " << noTime << " out of " 
      << highE << " high-enegy channel(s) have no time-info due to a lack of Optimal Filtering Coefficients." 
      << std::endl;
    }
    if ( noShape ) {
      msg() << msglvl << "   " << noShape << " out of " 
      << highE << " high-enegy channel(s) have no quality factor due to a lack of shape." 
      << std::endl;
    }
    if ( noShapeDer ) { 
      msg() << msglvl << "   " << noShapeDer << " out of " 
      << highE << " high-enegy channel(s) lack the derivative of the shape. Not taken into accout for Quality factor." 
      << std::endl;
    }
    if ( saturation ) {
      if ( m_skipSaturCells == 2 )
        msg() << MSG::ERROR << "   " << saturation << " out of " 
          << digitContainer->size() << " channel(s) showed saturations. The complete event was skipped." << std::endl;
      else if ( m_skipSaturCells == 1 )
        msg() << MSG::ERROR << "   " << saturation << " out of " 
          << digitContainer->size() << " channel(s) showed saturations and were skipped." << std::endl;
      else
        msg() << MSG::WARNING << "   " << saturation << " out of " 
          << digitContainer->size() << " channel(s) showed saturations." << std::endl;
    }
    msg() << endmsg;
  }
    
  // lock raw channel container
  ATH_CHECK( evtStore()->setConst(larRawChannelContainer) );
 
  return StatusCode::SUCCESS;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::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 & AthCommonAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 89 of file AthCommonAlgorithm.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 BaseAlg::extraOutputDeps();
}

filterPassed()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

finalize()

StatusCode TBECLArRawChannelBuilder::finalize ( )

overridevirtual

Definition at line 721 of file TBECLArRawChannelBuilder.cxx.

{ 
 
  m_aveNoEnergy   /= m_noEnergy;
  m_aveNoTime     /= m_noTime;
  m_aveNoShape    /= m_noShape;
  m_aveNoShapeDer /= m_noShapeDer;
  m_aveSaturCells /= m_saturation;
 
  m_aveChannels   /= m_nEvents;
 
  // Error and Warning Summary for this job:
  
  ATH_MSG_DEBUG ( "  TBECLArRawChannelBuilder::finalize " 
                  << m_noEnergy << " " << m_noTime << " " << m_noShape << " " << m_noShapeDer << " " << m_saturation );
  
  if ( m_noEnergy || m_noTime || m_noShape || m_noShapeDer || m_saturation ) {
    MSG::Level msglvl;
    if ( m_noEnergy || m_skipSaturCells ) 
      msglvl=MSG::ERROR;
    else
      msglvl=MSG::WARNING;
    msg() << msglvl << " *** Error & Warning Summary for all events *** " << std::endl ;
    
    if (m_noEnergy)
      msg() << msglvl << "   " << m_noEnergy << " events had on average " << (int)round(m_aveNoEnergy) 
          << " channels out of " << (int)round(m_aveChannels) << " without basic calibration constants." 
      << std::endl;
    
    if (m_noTime) 
      msg() << msglvl << "   " << m_noTime  << " events had on average " << (int)round(m_aveNoTime) 
          << " channels out of " << (int)round(m_aveChannels) << " without OFCs for timing." 
      << std::endl ;
   
    if (m_noShape)
      msg() << msglvl << "   " << m_noShape << " events had on average " << (int)round(m_aveNoShape) 
          << " channels out of " << (int)round(m_aveChannels) << " without shape information." 
      << std::endl;
    
    if (m_noShapeDer)
      msg() << msglvl << "   " << m_noShapeDer << " events had on average " << (int)round(m_aveNoShapeDer) 
          << " channels out of " << (int)round(m_aveChannels) << " without shape derivative." 
      << std::endl;
    
    if ( m_saturation )
      msg() << msglvl << "   " << m_saturation << " events had on average " << (int)round(m_aveSaturCells) 
          << " out of " << (int)round(m_aveChannels) << " saturating channels."  
      << std::endl ;
    
    msg() << endmsg;
  } 
  else
    ATH_MSG_INFO ( "TBECLArRawChannelBuilder finished without errors or warnings." );
 
  //if (m_larRawChannelContainer) {
    //m_larRawChannelContainer->release();
    //m_larRawChannelContainer = 0;
  //}
 
  return StatusCode::SUCCESS;
}

initialize()

StatusCode TBECLArRawChannelBuilder::initialize ( )

overridevirtual

Definition at line 92 of file TBECLArRawChannelBuilder.cxx.

                                               {
 
  ATH_CHECK( detStore()->retrieve(m_onlineHelper, "LArOnlineID") );
 
  ATH_CHECK( m_ofcKey.initialize() );
  ATH_CHECK( m_adc2mevKey.initialize (m_useRamp) );
 
  if (!m_useRamp)
  {
    // pointer to CaloCell ID helper:
    ATH_CHECK( detStore()->retrieve (m_calo_id, "CaloCell_ID") );
    
    for (int i=0; i<30; i++) {
      m_adc2mev[i] = 0;
      if (i == 6)  m_adc2mev[i] = 0.041*637;    // EMEC2
      if (i == 7)  m_adc2mev[i] = 0.030*637;    // EMEC3
      if (i == 8)  m_adc2mev[i] = 0.00360*3270; // HEC0
      if (i == 9)  m_adc2mev[i] = 0.00380*3270; // HEC1
      if (i == 10) m_adc2mev[i] = 0.00186*6540; // HEC2
      if (i == 21 || i == 22) m_adc2mev[i] = 0.1087*767;  // FCal1,2
      if (i == 23) m_adc2mev[i] = 0.1087*1508;  // FCal3
    }
  }
  
  // ***
 
  m_emId=m_calo_id->em_idHelper();
  
  // translate offline ID into online ID
  ATH_CHECK( m_cablingKey.initialize() );
  
  // ***
 
  ATH_CHECK( m_offlineHVScaleCorrKey.initialize(m_hvcorr) );
 
 
  //m_larRawOrdering.setMap(&(*m_roiMap)); 
 
  //Set counters for errors and warnings to zero
  m_noEnergy   = 0; // Number of events with at least completly failed channel
  m_noTime     = 0; // Number of events with at least one channel without time info
  m_noShape    = 0; // Number of events with at least one channel without Shape (=with not quality factor);
  m_noShapeDer = 0; // Number of events with at least one channel without ShapeDerivative (=not taken into accout for quality factor);
  m_saturation = 0; // Number of events with at least one saturating channel
  
  m_lastNoEnergy   = -1; // Number of completly failed channels in previous event
  m_lastNoTime     = -1; // Number of channels without time info in previous event
  m_lastNoShape    = -1; // Number of channels without Shape (=with not quality factor) in previous event
  m_lastNoShapeDer = -1; // Number of channels without ShapeDerivative in previous event
  
  //m_lastSaturCells = -1; // Number of saturating channels without in previous event (not used)
  
  m_aveNoEnergy    = 0.; // Average number of completly failed channels per event
  m_aveNoTime      = 0.; // Average number of channels without time info per event
  m_aveNoShape     = 0.; // Average number of channels without Shape (=with not quality factor) per event
  m_aveNoShapeDer  = 0.; // Average number of channels without ShapeDerivative per event
  m_aveSaturCells  = 0.; // Average number of saturating channels without per event
  
  m_nEvents        = 0 ; // Total number of processed events ;
  m_aveChannels    = 0 ; // Average number of readout channels per event
 
  if ( m_skipSaturCells > 2 ) m_skipSaturCells = 0 ;
 
  m_ramp_max[CaloGain::LARHIGHGAIN]   = 500.;
  m_ramp_max[CaloGain::LARMEDIUMGAIN] = 5000.;
  m_ramp_max[CaloGain::LARLOWGAIN]    = 50000.;
 
  // Validity range for a set of OFC's. If the time shift is larger than this number,
  // we make a ADC sample jump (e.g. from [0,5] to [1,6]. The second half of the uppermost 
  // bin should already be rounded to the 0th bin of the following ADC sample.  
  if ( m_binHalfOffset ) {
    m_SamplingPeriodeUpperLimit = m_SamplingPeriode-m_SamplingPeriode/(2*m_NOFCTimeBins);
    m_SamplingPeriodeLowerLimit = -m_SamplingPeriode/(2*m_NOFCTimeBins);
  } else {
    m_SamplingPeriodeUpperLimit = m_SamplingPeriode;
    m_SamplingPeriodeLowerLimit = 0;
  }
 
  ATH_MSG_DEBUG ( "Number of OFC time bins per sampling periode=" << m_NOFCTimeBins );
  ATH_MSG_DEBUG ( "Sampling Periode=" << m_SamplingPeriode << "ns" );
  ATH_MSG_DEBUG ( "Sampling Periode Limits: (" << m_SamplingPeriodeLowerLimit
                  << "," << m_SamplingPeriodeUpperLimit << ") ns" );
 
  return StatusCode::SUCCESS;
}

inputHandles()

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

isClonable()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

                                           {
    return true;
  }

isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const

inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 47 of file AthAlgorithm.h.

{ return false; }

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

virtual void AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::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< Gaudi::Algorithm > >.

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

Definition at line 60 of file AthCommonAlgorithm.cxx.

                                                      {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::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< Gaudi::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< Gaudi::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_adc2mev

float TBECLArRawChannelBuilder::m_adc2mev[30]

private

Definition at line 99 of file TBECLArRawChannelBuilder.h.

m_adc2mevKey

SG::ReadCondHandleKey<LArADC2MeV> TBECLArRawChannelBuilder::m_adc2mevKey { this, "ADC2MeVKey", "LArADC2MeV", "SG Key of the LArADC2MeV CDO" }

private

Definition at line 105 of file TBECLArRawChannelBuilder.h.

{ this, "ADC2MeVKey", "LArADC2MeV", "SG Key of the LArADC2MeV CDO" };

m_AdcMax

short TBECLArRawChannelBuilder::m_AdcMax

private

Definition at line 63 of file TBECLArRawChannelBuilder.h.

m_allowTimeJump

bool TBECLArRawChannelBuilder::m_allowTimeJump

private

Definition at line 57 of file TBECLArRawChannelBuilder.h.

m_aveChannels

float TBECLArRawChannelBuilder::m_aveChannels

private

Definition at line 86 of file TBECLArRawChannelBuilder.h.

m_aveNoEnergy

float TBECLArRawChannelBuilder::m_aveNoEnergy

private

Definition at line 79 of file TBECLArRawChannelBuilder.h.

m_aveNoShape

float TBECLArRawChannelBuilder::m_aveNoShape

private

Definition at line 81 of file TBECLArRawChannelBuilder.h.

m_aveNoShapeDer

float TBECLArRawChannelBuilder::m_aveNoShapeDer

private

Definition at line 82 of file TBECLArRawChannelBuilder.h.

m_aveNoTime

float TBECLArRawChannelBuilder::m_aveNoTime

private

Definition at line 80 of file TBECLArRawChannelBuilder.h.

m_aveSaturCells

float TBECLArRawChannelBuilder::m_aveSaturCells

private

Definition at line 83 of file TBECLArRawChannelBuilder.h.

m_binHalfOffset

bool TBECLArRawChannelBuilder::m_binHalfOffset

private

Definition at line 56 of file TBECLArRawChannelBuilder.h.

m_cablingKey

SG::ReadCondHandleKey<LArOnOffIdMapping> TBECLArRawChannelBuilder::m_cablingKey {this,"keyCabling", "LArOnOffIdMap", "Input key for Id mapping"}

private

Definition at line 38 of file TBECLArRawChannelBuilder.h.

{this,"keyCabling", "LArOnOffIdMap", "Input key for Id mapping"} ;

m_calo_id

const CaloCell_ID* TBECLArRawChannelBuilder::m_calo_id

private

Definition at line 36 of file TBECLArRawChannelBuilder.h.

m_ChannelContainerName

std::string TBECLArRawChannelBuilder::m_ChannelContainerName

private

Definition at line 42 of file TBECLArRawChannelBuilder.h.

m_ConvertADCToHighGain

bool TBECLArRawChannelBuilder::m_ConvertADCToHighGain

private

Definition at line 44 of file TBECLArRawChannelBuilder.h.

m_DataLocation

std::string TBECLArRawChannelBuilder::m_DataLocation

private

Definition at line 42 of file TBECLArRawChannelBuilder.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

m_Ecut

float TBECLArRawChannelBuilder::m_Ecut

private

Definition at line 45 of file TBECLArRawChannelBuilder.h.

m_emId

const LArEM_ID* TBECLArRawChannelBuilder::m_emId

private

Definition at line 96 of file TBECLArRawChannelBuilder.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

m_hvcorr

bool TBECLArRawChannelBuilder::m_hvcorr

private

Definition at line 87 of file TBECLArRawChannelBuilder.h.

m_initialTimeSampleShift

int TBECLArRawChannelBuilder::m_initialTimeSampleShift

private

Definition at line 46 of file TBECLArRawChannelBuilder.h.

m_iPedestal

unsigned int TBECLArRawChannelBuilder::m_iPedestal

private

Definition at line 59 of file TBECLArRawChannelBuilder.h.

m_lastNoEnergy

int TBECLArRawChannelBuilder::m_lastNoEnergy

private

Definition at line 73 of file TBECLArRawChannelBuilder.h.

m_lastNoShape

int TBECLArRawChannelBuilder::m_lastNoShape

private

Definition at line 75 of file TBECLArRawChannelBuilder.h.

m_lastNoShapeDer

int TBECLArRawChannelBuilder::m_lastNoShapeDer

private

Definition at line 76 of file TBECLArRawChannelBuilder.h.

m_lastNoTime

int TBECLArRawChannelBuilder::m_lastNoTime

private

Definition at line 74 of file TBECLArRawChannelBuilder.h.

m_nEvents

int TBECLArRawChannelBuilder::m_nEvents

private

Definition at line 85 of file TBECLArRawChannelBuilder.h.

m_noEnergy

int TBECLArRawChannelBuilder::m_noEnergy

private

Definition at line 66 of file TBECLArRawChannelBuilder.h.

m_NOFCPhases

int TBECLArRawChannelBuilder::m_NOFCPhases

private

Definition at line 47 of file TBECLArRawChannelBuilder.h.

m_NOFCTimeBins

int TBECLArRawChannelBuilder::m_NOFCTimeBins

private

Definition at line 48 of file TBECLArRawChannelBuilder.h.

m_noShape

int TBECLArRawChannelBuilder::m_noShape

private

Definition at line 68 of file TBECLArRawChannelBuilder.h.

m_noShapeDer

int TBECLArRawChannelBuilder::m_noShapeDer

private

Definition at line 69 of file TBECLArRawChannelBuilder.h.

m_noTime

int TBECLArRawChannelBuilder::m_noTime

private

Definition at line 67 of file TBECLArRawChannelBuilder.h.

m_ofcKey

SG::ReadCondHandleKey<ILArOFC> TBECLArRawChannelBuilder::m_ofcKey {this, "OFCKey", "LArOFC", "SG Key of OFC conditions object" }

private

Property: OFC coefficients (conditions input).

Definition at line 102 of file TBECLArRawChannelBuilder.h.

{this, "OFCKey", "LArOFC", "SG Key of OFC conditions object" };

m_OFCTimeBin

float TBECLArRawChannelBuilder::m_OFCTimeBin

private

Definition at line 55 of file TBECLArRawChannelBuilder.h.

m_offlineHVScaleCorrKey

SG::ReadCondHandleKey<ILArHVScaleCorr> TBECLArRawChannelBuilder::m_offlineHVScaleCorrKey {this, "keyOfflineHVCorr", "LArHVScaleCorrRecomputed","Key for LArHVScaleCorr"}

private

Definition at line 39 of file TBECLArRawChannelBuilder.h.

{this, "keyOfflineHVCorr", "LArHVScaleCorrRecomputed","Key for LArHVScaleCorr"};

m_onlineHelper

const LArOnlineID* TBECLArRawChannelBuilder::m_onlineHelper

private

Definition at line 35 of file TBECLArRawChannelBuilder.h.

m_pedestalFallbackMode

int TBECLArRawChannelBuilder::m_pedestalFallbackMode

private

Definition at line 58 of file TBECLArRawChannelBuilder.h.

m_phaseInv

bool TBECLArRawChannelBuilder::m_phaseInv

private

Definition at line 51 of file TBECLArRawChannelBuilder.h.

m_ramp_max

float TBECLArRawChannelBuilder::m_ramp_max[3]

private

Definition at line 52 of file TBECLArRawChannelBuilder.h.

m_SamplingPeriode

double TBECLArRawChannelBuilder::m_SamplingPeriode

private

Definition at line 54 of file TBECLArRawChannelBuilder.h.

m_SamplingPeriodeLowerLimit

double TBECLArRawChannelBuilder::m_SamplingPeriodeLowerLimit

private

Definition at line 90 of file TBECLArRawChannelBuilder.h.

m_SamplingPeriodeUpperLimit

double TBECLArRawChannelBuilder::m_SamplingPeriodeUpperLimit

private

Definition at line 90 of file TBECLArRawChannelBuilder.h.

m_saturation

int TBECLArRawChannelBuilder::m_saturation

private

Definition at line 70 of file TBECLArRawChannelBuilder.h.

m_shapeMode

unsigned int TBECLArRawChannelBuilder::m_shapeMode

private

Definition at line 60 of file TBECLArRawChannelBuilder.h.

m_skipSaturCells

unsigned TBECLArRawChannelBuilder::m_skipSaturCells

private

Definition at line 62 of file TBECLArRawChannelBuilder.h.

m_useIntercept

bool TBECLArRawChannelBuilder::m_useIntercept[4]

private

Definition at line 49 of file TBECLArRawChannelBuilder.h.

m_useOFCPhase

bool TBECLArRawChannelBuilder::m_useOFCPhase

private

Definition at line 50 of file TBECLArRawChannelBuilder.h.

m_useRamp

bool TBECLArRawChannelBuilder::m_useRamp

private

Definition at line 43 of file TBECLArRawChannelBuilder.h.

m_useShape

bool TBECLArRawChannelBuilder::m_useShape

private

Definition at line 43 of file TBECLArRawChannelBuilder.h.

m_useTDC

bool TBECLArRawChannelBuilder::m_useTDC

private

Definition at line 43 of file TBECLArRawChannelBuilder.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• TBECLArRawChannelBuilder.h
• TBECLArRawChannelBuilder.cxx