LArRawChannelSimpleBuilder Class Reference

CTB simplified LArRawChannelBuilder using hardcoded calib. constants. More...

#include <LArRawChannelSimpleBuilder.h>

Inheritance diagram for LArRawChannelSimpleBuilder:

Public Member Functions
	LArRawChannelSimpleBuilder (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode	initialize () override
virtual StatusCode	execute (const EventContext &ctx) const override
virtual StatusCode	finalize () override
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
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
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
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
SG::ReadCondHandleKey< LArOnOffIdMapping >	m_cablingKey {this, "OnOffMap", "LArOnOffIdMap", "SG key for mapping object"}
SG::ReadCondHandleKey< LArADC2MeV >	m_adc2mevKey { this, "ADC2MeVKey", "LArADC2MeV", "SG Key of the LArADC2MeV CDO" }
const LArEM_ID *	m_emId
const LArFCAL_ID *	m_fcalId
const LArHEC_ID *	m_hecId
const LArOnlineID *	m_onlineHelper
ToolHandle< LArParabolaPeakRecoTool >	m_peakParabolaTool
SG::ReadHandleKey< LArDigitContainer >	m_DataLocation { this, "DataLocation", "FREE", "" }
SG::WriteHandleKey< LArRawChannelContainer >	m_ChannelContainerName { this, "LArRawChannelContainerName", "LArRawChannels", "" }
int	m_imaxSamp
std::string	m_mode
std::string	m_FCALmodeTime
float	m_cubicAdcCut
bool	m_useRampDB
bool	m_usePedestalDB
int	m_averageSamplesEM
int	m_averageSamplesHEC
int	m_averageSamplesFCAL
float	m_averageScaleEM
float	m_averageScaleHEC
float	m_averageScaleFCAL
float	m_ADCtoMeVFCAL [3] {}
float	m_ADCtoMeVHEC [2] {}
float	m_ADCtoMeVEMECInner [2] {}
float	m_ADCtoMeVEMECOuter [4] {}
float	m_ADCtoMeVEMB [4] {}
unsigned int	m_iPedestal
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

CTB simplified LArRawChannelBuilder using hardcoded calib. constants.

Author

Walter Lampl

Rob McPherson

Sandrine Laplace

Definition at line 35 of file LArRawChannelSimpleBuilder.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

LArRawChannelSimpleBuilder()

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

Definition at line 23 of file LArRawChannelSimpleBuilder.cxx.

                                                                                                      :
  AthReentrantAlgorithm(name, pSvcLocator),
  m_emId(nullptr),
  m_fcalId(nullptr),
  m_hecId(nullptr),
  m_onlineHelper(nullptr),
  m_peakParabolaTool("LArParabolaPeakTool"),
  m_iPedestal(0)// jobO ?
{
   //m_useIntercept={false,false,false,false};
 declareProperty("maxSamp",m_imaxSamp=8);
 declareProperty("RecoMode",m_mode="CUBIC");
 declareProperty("CubicRecoTimeModeFCAL",m_FCALmodeTime="LINEAR");
 declareProperty("CubicAdcCut",m_cubicAdcCut=15.0);
 declareProperty("PedestalSample",m_iPedestal=0);
 declareProperty("UsePedestalDB",m_usePedestalDB=false);
 declareProperty("UseRampDB",m_useRampDB=false);
 declareProperty("AverageSamplesEM",m_averageSamplesEM=5); 
 declareProperty("AverageSamplesHEC",m_averageSamplesHEC=5); 
 declareProperty("AverageSamplesFCAL",m_averageSamplesFCAL=3); 
 declareProperty("AverageScaleEM",m_averageScaleEM=2.6); 
 declareProperty("AverageScaleHEC",m_averageScaleHEC=2.6); 
 declareProperty("AverageScaleFCAL",m_averageScaleFCAL=1.8); 
 
 // m_peakParabolaTool = NULL ; // FIXME RS use empty ToolHandle - python
 
  m_ADCtoMeVFCAL[0]      =   87.0 * MeV;  // FCAL1 High gain
  m_ADCtoMeVFCAL[1]      =  117.0 * MeV;  // FCAL2 High gain
  m_ADCtoMeVFCAL[2]      =  193.0 * MeV;  // FCAL3 High gain
  m_ADCtoMeVHEC[0]       =  136.0 / 9.8 * MeV;  // HEC 1 Medium gain from Monika.  Need / 9.8  ?? 
  m_ADCtoMeVHEC[1]       =  136.0 / 9.8 * MeV;  // HEC 2 Medium gain from Monika.  Need / 9.8  ?? 
  // m_ADCtoMeVEMECInner[0] =   25.22 * MeV; // EMEC  High gain from Pascal, approximate
  // m_ADCtoMeVEMECInner[1] =   19.4  * MeV;  // EMEC  High gain from Pascal, approximate
  m_ADCtoMeVEMECInner[0] =   20.0 * MeV;  // EMEC IW s=1  High gain : fixed 18/8/2004 RMcP
  m_ADCtoMeVEMECInner[1] =   20.0 * MeV;  // EMEC IW s=2 High gain : fixed 18/8/2004 RMcP
  m_ADCtoMeVEMECOuter[0] =   16.0 * MeV;  // EMEC OW pre, s=0 
  m_ADCtoMeVEMECOuter[1] =   16.0 * MeV;  // EMEC OW s=1 High gain from Monika, approximate
  m_ADCtoMeVEMECOuter[2] =   16.0 * MeV;  // EMEC OW s=2 High gain from Monika, approximate
  m_ADCtoMeVEMECOuter[3] =   16.0 * MeV;  // EMEC OW s=3 High gain from Monika, approximate
  m_ADCtoMeVEMB[0]       =   7.0  * MeV;  // EMB pre, s=0 High gain from Isabelle, approximate
  m_ADCtoMeVEMB[1]       =   2.5  * MeV;  // EMB s=1 High gain from Isabelle, approximate
  m_ADCtoMeVEMB[2]       =   18.0 * MeV;  // EMB s=2 High gain from Isabelle, approximate
  m_ADCtoMeVEMB[3]       =   9.0  * MeV;  // EMB s=3 High gain from Isabelle, approximate
 
 
}

Member Function Documentation

cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

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 LArRawChannelSimpleBuilder::execute ( const EventContext & ctx ) const

overridevirtual

Definition at line 103 of file LArRawChannelSimpleBuilder.cxx.

{
  ATH_MSG_DEBUG( "In execute"  );
  
  //Retrieve Digit Container
  ATH_MSG_DEBUG( "About to retrieve LArDigitContainer with key " << m_DataLocation  );
  SG::ReadHandle<LArDigitContainer> digitContainer (m_DataLocation, ctx);
  ATH_MSG_DEBUG( "1) LArDigitContainer container size = " <<  digitContainer->size()  );
 
  ATH_MSG_DEBUG( "2) LArDigitContainer container size = " <<  digitContainer->size()  );
  if( digitContainer->empty() ) {
    ATH_MSG_INFO( "Empty LArDigitContainer container."  );
    return StatusCode::SUCCESS;
  }
  
  auto larRawChannelContainer = std::make_unique<LArRawChannelContainer>();
 
  //Pointer to conditions data objects 
  const ILArPedestal* larPedestal=nullptr;
  if (m_usePedestalDB) {
    if (detStore()->retrieve(larPedestal).isFailure()) {
      larPedestal=nullptr;
      ATH_MSG_DEBUG( "No pedestal found in database. Use default values."  );
    }
  }
 
  SG::ReadCondHandle<LArOnOffIdMapping> cablingHdl{m_cablingKey, ctx};
  const LArOnOffIdMapping* cabling{*cablingHdl};
  if(!cabling) {
     ATH_MSG_ERROR("Do not have mapping object " << m_cablingKey.key() );
     return StatusCode::FAILURE;
  }
 
  const LArADC2MeV* adc2mev = nullptr;
  if (m_useRampDB) {
    SG::ReadCondHandle<LArADC2MeV> adc2mevH (m_adc2mevKey, ctx);
    adc2mev = *adc2mevH;
  }
 
  //loop twice over the digits. In the first pass the best window is
  //found for the averaged signal. In the second pass the full raw
  //channel is created and the averaging in the window found in the
  //first path is used in case the cubic fit fails or the max ADC is
  //below threshold
 
  int nMinEM(0),nMinHEC(0),nMinFCAL(0);
  std::vector<float> fSumEM,fSumHEC,fSumFCAL;
  std::vector<float> *pSum = nullptr;
  
  fSumEM.resize(0);
  fSumHEC.resize(0);
  fSumFCAL.resize(0);
 
  for(int iloop=0;iloop<2;iloop++) {
   
    //Now all data is available, start loop over Digit Container
    ATH_MSG_DEBUG( "Loop over Digit Container "  );
 
    for (const LArDigit* digit : *digitContainer) {
    
      //Get data from LArDigit
      const std::vector<short>& samples  = digit->samples();
      unsigned int              nSamples = samples.size(); 
      const HWIdentifier        chid     = digit->channelID();
      const CaloGain::CaloGain  gain     = digit->gain();
      
      
      float thePedestal=-1;    
      if (larPedestal) {
    float DBpedestal =larPedestal->pedestal(chid,gain);
    if (DBpedestal <= (1.+LArElecCalib::ERRORCODE))
      thePedestal=DBpedestal;
      }
      if (thePedestal<0) {
    thePedestal = (float)samples[m_iPedestal];
    ATH_MSG_DEBUG( "No pedestal found for this cell. Use default value " << thePedestal  );
      }
      //>>>> PL June 20, 2004: subtract pedestal first, assume first sample ->JO?
      std::vector<float> mySamples;
      mySamples.resize(samples.size());
      // FIXME!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! only for 5 leading noise samples!
      
      //    log << MSG::INFO
      //  << "pedestal " << thePedestal << endmsg;
      for ( unsigned int i=0; i<samples.size(); i++ )
    {
      mySamples[i] = ((float)samples[i]) - thePedestal;
    }
      //<<<<
      
      float GainFactor;
      if( gain == CaloGain::LARLOWGAIN ) {
    GainFactor = 9.8*9.8;
      } else if (gain == CaloGain::LARMEDIUMGAIN ) {
    GainFactor = 9.8;
      } else if (gain == CaloGain::LARHIGHGAIN ) {
    GainFactor = 1.0;
      } else {
    GainFactor = 1.0;
    ATH_MSG_ERROR( "Channel " << chid << "unknown gain: " << gain  );
      }
      
      // Get hardware identifier for later use.
      HWIdentifier FebID = m_onlineHelper->feb_Id(chid);
      unsigned int channel  = m_onlineHelper->channel(chid);
      ATH_MSG_DEBUG( std::hex << " FebID / chid / channel = " << FebID << " / " << chid << " / " << channel );
      
      // And note if this is an FCAL channel with fast pulses (cubic fails)
      // EM and HEC are also tested to adjust the signal range for averaging
      bool isFCAL = false;
      bool isEM   = false;
      bool isHEC  = false;
      int nMin = 0;
      unsigned int nAverage = 1;
      float myScale = 1;
      try {
    const Identifier id = cabling->cnvToIdentifier(chid);    
    if (m_fcalId->is_lar_fcal(id)) {
      isFCAL = true;
      nMin = nMinFCAL;
      pSum = &fSumFCAL;
      nAverage = m_averageSamplesFCAL;
      myScale = m_averageScaleFCAL;
    }
    else if (m_emId->is_lar_em(id)) {
      isEM = true;
      nMin = nMinEM;
      pSum = &fSumEM;
      nAverage = m_averageSamplesEM;
      myScale = m_averageScaleEM;
    }
    else if (m_hecId->is_lar_hec(id)) {
      isHEC = true;
      nMin = nMinHEC;
      pSum = &fSumHEC;
      nAverage = m_averageSamplesHEC;
      myScale = m_averageScaleHEC;
    }
      }
      //catch( LArOnlID_Exception & except) 
      catch (LArID_Exception & execpt) {
    ATH_MSG_DEBUG (" is disconnected.");
        // The question now being:  do we want to skip this channel???
        // Probably yes, so do so.  RMcP 9 June 2006
        continue;
      }
      // In the unrealistic case, that pSum is still not defined
      if(!pSum) continue;
 
      // Find peak time sample and scaled average for selected sample range 
      // (i.e. poor man's digital filter)
 
      // sanity checks
      if ( nAverage > nSamples ) {
    ATH_MSG_WARNING( " Number of samples to average (" 
                         << nAverage << ") is larger than total number of samples (" 
                         << nSamples << ")! adjusting nAverage ... "  );
    nAverage = nSamples;
      }
 
 
      if ( iloop == 0 ) {
 
    if ( pSum->empty()) 
//    pSum->resize(nSamples-nAverage+1,0);
      pSum->resize(nSamples,0);
 
    for( unsigned i=0;i<nSamples-nAverage+1;i++ ) {
      for( unsigned j=0;j<nAverage;j++ ) {
        (*pSum)[i] += mySamples[i+j];
      }
    } 
      }
      else {
    float        maxADCPeak = 0.;
    unsigned int iPeakSamp  = 0;
    float        averagedADC = 0;
    for( unsigned i=0;i<nSamples;i++ ) {
      if ( maxADCPeak < mySamples[i] ) 
        {
          maxADCPeak = mySamples[i];
          iPeakSamp  = i;
        }
      if ( (int)i >= nMin && i < nMin+nAverage ) 
        averagedADC += mySamples[i];
    }
    averagedADC /= myScale;
 
    bool CubicFailed = false;
    float ADCPeak=0.;
    float time=0;
      
      
    // maximum amplitude/variable time slice method
    if( m_mode == "MAX" ) 
      {
        ADCPeak = maxADCPeak;
        time    = ((float)iPeakSamp) * 25.0 * nanosecond;
      }
      
      
    // other choices are different for FCAL and others
    else {
    
      // first, deal with FCAL: only one reconstruction is available
      if(isFCAL &&   m_mode == "CUBIC" &&
         ! ( CubicFailed = maxADCPeak <= m_cubicAdcCut ) ) {
      
        ATH_MSG_DEBUG( " Special reconstruction for FCAL."  );
          
        unsigned int it0;
      
        const float invT[3][3] 
          = { {    1,       0,       0 },
          { -1.5,       2,    -0.5 },
          {  0.5,      -1,     0.5 }  };
      
        // peak slice very early
        if ( iPeakSamp <= 1 ) {
          it0 = 1;
        } else if ( iPeakSamp >= nSamples - 1 ) {  // peak is late
          it0 = nSamples - 3;
        } else { // peak in safe region
          it0 = iPeakSamp - 1;
        }
      
        // Quadratic interpolation using
        // 3 samples to be used start at 0 <= t0 <= nsamples-3
        float A[3] = {0, 0, 0};
        float dtmax = 0.0;
        // S = TA -> A = inv(T)S
        for (int ia = 0; ia < 3; ia++) 
          for (int it = 0; it < 3; it++) 
        A[ia] += invT[ia][it] * mySamples[it0+it];
      
        // fit parameters
        if ( not ( CubicFailed = ( A[2] == 0 ) ) )   {
          dtmax = -1.0 * A[1] / 2.0 / A[2];
          if ( ! ( CubicFailed = ( dtmax < 0 || dtmax > 2 ) ) ) {
        //time = (float(it0) + dtmax) * 25.0 * nanosecond; // nsec
        time=dtmax*25.0*ns;
        for(int ia = 0; ia < 3; ia++)
          ADCPeak += A[ia] * pow(dtmax, ia);
          }
        }
      
        // Now use jobOptions to pick time and height of FCAL pulses
        if( m_FCALmodeTime == "LINEAR" ) {
          float weightSum = 0.;
          float timeSum   = 0.;
          for( int it=0; it<3; it++) {
        timeSum   += float(mySamples[it0+it]) * float(it);
        weightSum += float(mySamples[it0+it]);
          }
          time    = (float(it0) + timeSum / weightSum) * 25.0 * nanosecond;
          // ADCPeak = mySamples[iPeakSamp];
          ADCPeak = mySamples[m_imaxSamp];
          CubicFailed = false;
        }
      
        // then, deal with non-FCAL
      } else if(not isFCAL ||  m_mode != "CUBIC" ) {
      
        // bias-corrected parabola extrapolation for selected channels
        if ( m_mode == "PARABOLA" &&
         ! ( CubicFailed = maxADCPeak <= m_cubicAdcCut ) ) {
          if( m_peakParabolaTool) {
        int layer = 0;
        try {
          const Identifier id = cabling->cnvToIdentifier(chid);
          ATH_MSG_DEBUG( std::hex << "     id = " << id );
          if (m_emId->is_em_barrel(id)) {
            layer= m_emId->sampling(id);
          }
        }
        catch (LArID_Exception & execpt){
                  ATH_MSG_DEBUG( std::hex
                                 << " Cannot get offline identifier from online ID = " << chid );
        }
          
        std::vector<float> peak =m_peakParabolaTool->peak(samples,layer,thePedestal); 
        if(peak.size() >1){
          ADCPeak = peak[0]-thePedestal;
          if (peak.size()==2)time = peak[1];
          else time = peak[2];
        }else{
          ATH_MSG_DEBUG( "No pic is computed from Parabola. Use scaled average of selected samples"  );
          ADCPeak = averagedADC;
          time    = -99.;
        }
          }else{
        ATH_MSG_FATAL( "No parabola tool available ! Choose another mode"  );
          }
        }
      
        // cubic extrapolation for selected channels
        else if ( m_mode == "CUBIC" &&
              ! ( CubicFailed = maxADCPeak <= m_cubicAdcCut ) ) {
        
          unsigned int it0;
          const float invT[4][4] 
        = { {        1,       0,       0,        0},
            { -1.83333,       3,    -1.5, 0.333333},
            {        1,    -2.5,       2,     -0.5},
            {-0.166666,     0.5,    -0.5, 0.166666}  };
        
          // peak slice very early
          if ( iPeakSamp <= 1 ) {
        it0 = m_iPedestal + 1;
          } else if ( iPeakSamp >= nSamples - 2 ) {  // peak is late
        it0 = nSamples - 4;
          } else { // peak in safe region
        it0 = ( mySamples[iPeakSamp-2] > mySamples[iPeakSamp+2] )
          ? iPeakSamp - 2
          : iPeakSamp - 1;
          }
        
          // 4 samples to be used start at 0 <= t0 <= nsamples-4
          float A[4] = {0, 0, 0, 0};
          float dtmax = 0.0;
          float disc;
          // S = TA -> A = inv(T)S
          for (int ia = 0; ia < 4; ia++) 
        for (int it = 0; it < 4; it++) 
          A[ia] += invT[ia][it] * mySamples[it0+it];
        
          // fit parameters
          disc = A[2]*A[2] - 3*A[1]*A[3];
          if ( ! ( CubicFailed = ( disc < 0 || A[3] == 0 ) ) )   {
        dtmax = (-A[2]-std::sqrt(disc))/(A[3]*3);
        if ( ! ( CubicFailed = ( dtmax < 0 || dtmax > 3 ) ) ) {
          time = (float(it0) + dtmax) * 25.0 * nanosecond; // nsec
          for(int ia = 0; ia < 4; ia++)
            ADCPeak += A[ia] * pow(dtmax, ia);
        }
          }
        }
      }
    }
      
    // fixed time slice or insufficient signal for cubic fit : use scaled average of selected samples  
    if(m_mode == "FIXED" || CubicFailed ) {
      ADCPeak = averagedADC;
      time    = -99. ;
    }
      
    ATH_MSG_DEBUG( "Flag: "
                       << CubicFailed
                       << " Detector: "
                       << (isEM?"EM":(isHEC?"HEC":(isFCAL?"FCAL":"none")))
                       << " Mode: "
                       << m_mode
                       << " Signal: "
                       <<  ADCPeak
                       << " Peak: "
                       << maxADCPeak
                       << " PeakSample: "
                       << iPeakSamp
                       );
      
    float energy=-9999;
    if (m_useRampDB) {
      float ADCPeakPower=ADCPeak;
      //ADC2MeV (a.k.a. Ramp)   
      LArVectorProxy ramp = adc2mev->ADC2MEV(chid,gain);
      //Check ramp coefficents
      if (ramp.size()>1 && ramp[1]<500 && ramp[1]>0) {
        energy=0;
        for (unsigned i=1;i<ramp.size();i++)
          {energy+=ramp[i]*ADCPeakPower; //pow(ADCPeak,i);
          //std::cout << "Step "<< i <<":" << ramp[i] << " * " << pow(ADCPeak,i) << "Sum=" << energy << std::endl;
          ADCPeakPower*=ADCPeak;
          }
      }
    }
    if (energy==-9999) {
      ATH_MSG_DEBUG( "No Ramp found for this cell. Use default values"  );
      //Apply default values
      // Now must get subdetector ID and feed in here ...
      float ADCtoMeV      = 10.0;
    
      //     HWIdentifier FebID = m_onlineHelper->feb_Id(chid);
      //     unsigned int channel  = m_onlineHelper->channel(chid);
      //     MSG::hex(log) << MSG::DEBUG
      //                   << " FebID / chid / channel = " << FebID << " / " << chid << " / " << channel << endmsg;
    
      try {
        const Identifier id = cabling->cnvToIdentifier(chid);
      
        ATH_MSG_DEBUG( std::hex << "     id = " << id );
      
        if (m_emId->is_em_barrel(id)) {
          const int layer= m_emId->sampling(id);
          const int eta=m_emId->eta(id);
          ADCtoMeV = m_ADCtoMeVEMB[layer];
          if (layer==2 && eta<32)
        ADCtoMeV=ADCtoMeV*(12./18.); //Correct for lead thickness
          ATH_MSG_DEBUG( " in EMB s="<< layer <<", using ADCtoMeV = " << ADCtoMeV  );
        } else if (m_emId->is_em_endcap_inner(id)) {
          // m_emId->sampling(id);
          ADCtoMeV = m_ADCtoMeVEMECInner[m_emId->sampling(id)-1];
          ATH_MSG_DEBUG( " in EMEC inner s="<<m_emId->sampling(id)<<", using ADCtoMeV = " << ADCtoMeV  );
        } else if (m_emId->is_em_endcap_outer(id)) {
          // m_emId->sampling(id);
          // ADCtoMeV = m_ADCtoMeVEMECOuter[m_emId->sampling(id)-1];
              ADCtoMeV = m_ADCtoMeVEMECOuter[m_emId->sampling(id)];
          ATH_MSG_DEBUG( " in EMEC outer s="<<m_emId->sampling(id)<<", using ADCtoMeV = " << ADCtoMeV  );
        } else if (m_fcalId->is_lar_fcal(id)) {
          // m_fcalId->module(chid);
          ADCtoMeV = m_ADCtoMeVFCAL[m_fcalId->module(id)-1];
          ATH_MSG_DEBUG( " in FCAL m=" << m_fcalId->module(id)<<", using ADCtoMeV = " << ADCtoMeV  );
        } else if (m_hecId->is_lar_hec(id)) {
          // m_layer[m_cellIndex]=m_hecId->sampling(id);
          ADCtoMeV = m_ADCtoMeVHEC[0];
          ATH_MSG_DEBUG( " in HEC s="<<m_hecId->sampling(id)<<", using ADCtoMeV = " << ADCtoMeV  );
        }
      }
      //catch( LArOnlID_Exception & except) 
      catch (LArID_Exception & execpt)
        {
              ATH_MSG_DEBUG( " is disconnected. Set ADCtoMeV to 1" );
              ADCtoMeV=1;
        }
      energy  = ADCPeak * ADCtoMeV * GainFactor;
    }
        uint16_t iquality=0;
        uint16_t iprovenance=0;
      
    time*=1000; //Convert to picoseconds
    //Make LArRawChannel Object with new data
    LArRawChannel larRawChannel(chid,(int)energy,(int)time,iquality,iprovenance, gain);   
    larRawChannelContainer->add(larRawChannel); //Add to container
      }
    }// End loop over LArDigits
    if ( iloop == 0 ) {
      // select the best windows for the event
      unsigned int i;
      float tmpSum = 0;
 
      if ( !fSumEM.empty() ) {
        for(i=0;i<fSumEM.size();i++) {
      if (i == 0 || fSumEM[i] > tmpSum ) {
        nMinEM = i;
        tmpSum = fSumEM[i];
      }
        }
    ATH_MSG_DEBUG( "Found best EM window starting at sample <" << nMinEM << ">"  );
      }
 
      for(i=0;i<fSumHEC.size();i++) {
    if (i == 0 || fSumHEC[i] > tmpSum ) {
      nMinHEC = i;
      tmpSum = fSumHEC[i];
    }
      }
      if ( !fSumHEC.empty() ) {
    ATH_MSG_DEBUG( "Found best HEC window starting at sample <" << nMinHEC << ">"  );
      }
 
      for(i=0;i<fSumFCAL.size();i++) {
    if (i == 0 || fSumFCAL[i] > tmpSum ) {
      nMinFCAL = i;
      tmpSum = fSumFCAL[i];
    }
      }
      if ( !fSumFCAL.empty() ) {
    ATH_MSG_DEBUG( "Found best FCAL window starting at sample <" << nMinFCAL << ">"  );
      }
 
    }
  }// End of double loop over Digit containers
 
  ATH_MSG_DEBUG( "Finished loop over Digit Container "  );
 
  // Put this LArRawChannel container in the transient store
  ATH_CHECK( SG::makeHandle (m_ChannelContainerName, ctx).record (std::move (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 & AthCommonReentrantAlgorithm< 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 94 of file AthCommonReentrantAlgorithm.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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

finalize()

StatusCode LArRawChannelSimpleBuilder::finalize ( )

overridevirtual

Definition at line 583 of file LArRawChannelSimpleBuilder.cxx.

{//Error and Warning Summary for this job:
  
  ATH_MSG_INFO( "LArRawChannelSimpleBuilder finished."  );
 
  return StatusCode::SUCCESS;
}

initialize()

StatusCode LArRawChannelSimpleBuilder::initialize ( )

overridevirtual

Definition at line 70 of file LArRawChannelSimpleBuilder.cxx.

{
  ATH_MSG_DEBUG( "In Initialize."  );
 
  ATH_CHECK( m_adc2mevKey.initialize (m_useRampDB) );
  
  if ( m_mode == "PARABOLA"){
    if (m_peakParabolaTool.retrieve().isFailure())
    {
      ATH_MSG_ERROR( "Can't get LArParabolaPeakRecoTool"  );
      return StatusCode::SUCCESS;
    }
    ATH_MSG_DEBUG( "LArParabolaPeakRecoTool retrieved with success!"  );
  }
 
  const CaloCell_ID* idHelper = nullptr;
  ATH_CHECK( detStore()->retrieve (idHelper, "CaloCell_ID") );
  m_emId=idHelper->em_idHelper();
  m_fcalId=idHelper->fcal_idHelper();
  m_hecId=idHelper->hec_idHelper();
 
  ATH_CHECK(m_cablingKey.initialize());
 
  ATH_CHECK( detStore()->retrieve(m_onlineHelper, "LArOnlineID") );
 
 
  ATH_CHECK( m_DataLocation.initialize() );
  ATH_CHECK( m_ChannelContainerName.initialize() );
  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()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, 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, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

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

sysExecute()

StatusCode AthCommonReentrantAlgorithm< 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.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonReentrantAlgorithm< 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 HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.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_adc2mevKey

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

private

Definition at line 49 of file LArRawChannelSimpleBuilder.h.

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

m_ADCtoMeVEMB

float LArRawChannelSimpleBuilder::m_ADCtoMeVEMB[4] {}

private

Definition at line 75 of file LArRawChannelSimpleBuilder.h.

{};

m_ADCtoMeVEMECInner

float LArRawChannelSimpleBuilder::m_ADCtoMeVEMECInner[2] {}

private

Definition at line 73 of file LArRawChannelSimpleBuilder.h.

{};

m_ADCtoMeVEMECOuter

float LArRawChannelSimpleBuilder::m_ADCtoMeVEMECOuter[4] {}

private

Definition at line 74 of file LArRawChannelSimpleBuilder.h.

{};

m_ADCtoMeVFCAL

float LArRawChannelSimpleBuilder::m_ADCtoMeVFCAL[3] {}

private

Definition at line 71 of file LArRawChannelSimpleBuilder.h.

{};

m_ADCtoMeVHEC

float LArRawChannelSimpleBuilder::m_ADCtoMeVHEC[2] {}

private

Definition at line 72 of file LArRawChannelSimpleBuilder.h.

{};

m_averageSamplesEM

int LArRawChannelSimpleBuilder::m_averageSamplesEM

private

Definition at line 68 of file LArRawChannelSimpleBuilder.h.

m_averageSamplesFCAL

int LArRawChannelSimpleBuilder::m_averageSamplesFCAL

private

Definition at line 68 of file LArRawChannelSimpleBuilder.h.

m_averageSamplesHEC

int LArRawChannelSimpleBuilder::m_averageSamplesHEC

private

Definition at line 68 of file LArRawChannelSimpleBuilder.h.

m_averageScaleEM

float LArRawChannelSimpleBuilder::m_averageScaleEM

private

Definition at line 69 of file LArRawChannelSimpleBuilder.h.

m_averageScaleFCAL

float LArRawChannelSimpleBuilder::m_averageScaleFCAL

private

Definition at line 69 of file LArRawChannelSimpleBuilder.h.

m_averageScaleHEC

float LArRawChannelSimpleBuilder::m_averageScaleHEC

private

Definition at line 69 of file LArRawChannelSimpleBuilder.h.

m_cablingKey

SG::ReadCondHandleKey<LArOnOffIdMapping> LArRawChannelSimpleBuilder::m_cablingKey {this, "OnOffMap", "LArOnOffIdMap", "SG key for mapping object"}

private

Definition at line 48 of file LArRawChannelSimpleBuilder.h.

{this, "OnOffMap", "LArOnOffIdMap", "SG key for mapping object"};

m_ChannelContainerName

SG::WriteHandleKey<LArRawChannelContainer> LArRawChannelSimpleBuilder::m_ChannelContainerName { this, "LArRawChannelContainerName", "LArRawChannels", "" }

private

Definition at line 61 of file LArRawChannelSimpleBuilder.h.

{ this, "LArRawChannelContainerName", "LArRawChannels", "" };

m_cubicAdcCut

float LArRawChannelSimpleBuilder::m_cubicAdcCut

private

Definition at line 66 of file LArRawChannelSimpleBuilder.h.

m_DataLocation

SG::ReadHandleKey<LArDigitContainer> LArRawChannelSimpleBuilder::m_DataLocation { this, "DataLocation", "FREE", "" }

private

Definition at line 59 of file LArRawChannelSimpleBuilder.h.

{ this, "DataLocation", "FREE", "" };

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_emId

const LArEM_ID* LArRawChannelSimpleBuilder::m_emId

private

Definition at line 51 of file LArRawChannelSimpleBuilder.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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_fcalId

const LArFCAL_ID* LArRawChannelSimpleBuilder::m_fcalId

private

Definition at line 52 of file LArRawChannelSimpleBuilder.h.

m_FCALmodeTime

std::string LArRawChannelSimpleBuilder::m_FCALmodeTime

private

Definition at line 65 of file LArRawChannelSimpleBuilder.h.

m_hecId

const LArHEC_ID* LArRawChannelSimpleBuilder::m_hecId

private

Definition at line 53 of file LArRawChannelSimpleBuilder.h.

m_imaxSamp

int LArRawChannelSimpleBuilder::m_imaxSamp

private

Definition at line 63 of file LArRawChannelSimpleBuilder.h.

m_iPedestal

unsigned int LArRawChannelSimpleBuilder::m_iPedestal

private

Definition at line 77 of file LArRawChannelSimpleBuilder.h.

m_mode

std::string LArRawChannelSimpleBuilder::m_mode

private

Definition at line 64 of file LArRawChannelSimpleBuilder.h.

m_onlineHelper

const LArOnlineID* LArRawChannelSimpleBuilder::m_onlineHelper

private

Definition at line 54 of file LArRawChannelSimpleBuilder.h.

m_peakParabolaTool

ToolHandle<LArParabolaPeakRecoTool> LArRawChannelSimpleBuilder::m_peakParabolaTool

private

Definition at line 56 of file LArRawChannelSimpleBuilder.h.

m_usePedestalDB

bool LArRawChannelSimpleBuilder::m_usePedestalDB

private

Definition at line 67 of file LArRawChannelSimpleBuilder.h.

m_useRampDB

bool LArRawChannelSimpleBuilder::m_useRampDB

private

Definition at line 67 of file LArRawChannelSimpleBuilder.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:

• LArRawChannelSimpleBuilder.h
• LArRawChannelSimpleBuilder.cxx