ATLAS Offline Software
Loading...
Searching...
No Matches
TileRawChannelBuilder Class Reference

#include <TileRawChannelBuilder.h>

Inherits AthAlgTool.

Inherited by TileRawChannelBuilderFitFilter, TileRawChannelBuilderFitFilterCool, TileRawChannelBuilderFlatFilter, TileRawChannelBuilderMF, TileRawChannelBuilderManyAmps, TileRawChannelBuilderOpt2Filter, TileRawChannelBuilderQIEFilter, TileRawChannelBuilderTest, and TileRawChannelBuilderWienerFilter.

Public Member Functions

 TileRawChannelBuilder (const std::string &type, const std::string &name, const IInterface *parent)
 Constructor.
virtual ~TileRawChannelBuilder ()
 Destructor.
virtual StatusCode initialize ()
 Initializer.
virtual StatusCode finalize ()
virtual StatusCode createContainer (const EventContext &ctx)
 Create container in SG with name given by parameter (m_rawChannelContainerKey).
virtual TileRawChannelrawChannel (const TileDigits *digits, const EventContext &ctx)
 Builder virtual method to be implemented by subclasses.
virtual StatusCode commitContainer (const EventContext &ctx)
 Commit RawChannelContiner in SG and make const.
void initLog (const EventContext &ctx)
void endLog ()
StatusCode build (const TileDigitsCollection *collection, const EventContext &ctx)
void resetDrawer ()
void resetOverflows (void)
Overflows_tgetOverflowedChannels (void)
std::string getTileRawChannelContainerID (void)
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 sysInitialize () override
 Perform system initialization for an algorithm.
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

Static Public Member Functions

static const InterfaceID & interfaceID ()
 AlgTool InterfaceID.
static double correctAmp (double phase, bool of2=true)
 Amplitude correction factor according to the time when using weights for tau=0 without iterations.
static double correctTime (double phase, bool of2=true)
 Time correction factor.
static int CorruptedData (int ros, int drawer, int channel, int gain, const std::vector< float > &digits, float &dmin, float &dmax, float ADCmaxMinusEps, float ADCmaskValueMinusEps)
static const char * BadPatternName (float ped)

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.

Protected Attributes

SG::ReadHandleKey< TileDQstatusm_DQstatusKey
SG::ReadHandleKey< TileRawChannelContainerm_DSPContainerKey {this, "DSPContainer", "", "DSP Container key"}
SG::WriteHandleKey< TileRawChannelContainerm_rawChannelContainerKey
std::unique_ptr< TileMutableRawChannelContainerm_rawChannelCnt
TileFragHash::TYPE m_rChType
TileRawChannelUnit::UNIT m_rChUnit
unsigned int m_bsflags
int m_firstSample
bool m_calibrateEnergy
bool m_correctTime
bool m_useDSP
float m_ampMinThresh
 correct amplitude if it's above amplitude threshold (in ADC counts)
float m_timeMinThresh
 correct amplitude is time is above time min threshold
float m_timeMaxThresh
 correct amplitude is time is below time max threshold
int m_runType
const TileIDm_tileID
const TileHWIDm_tileHWID
ToolHandleArray< ITileRawChannelToolm_noiseFilterTools
ToolHandle< TileCondToolEmscalem_tileToolEmscale
ToolHandle< TileCondToolTimingm_tileToolTiming
ToolHandle< TileCondIdTransformsm_tileIdTransforms
ServiceHandle< TileCablingSvcm_cablingSvc
 Name of Tile cabling service.
const TileCablingServicem_cabling
 TileCabling instance.
Gaudi::Property< std::vector< int > > m_demoFragIDs
int m_trigType
bool m_idophys
bool m_idolas
bool m_idoped
bool m_idocis
int m_cischan
double m_capdaq
unsigned int m_evtCounter
unsigned int m_chCounter
int m_nChL
int m_nChH
double m_RChSumL
double m_RChSumH
Overflows_t m_overflows
int m_dataPoollSize
int m_error [MAX_CHANNELS]
int m_lastDrawer = -1
bool m_badDrawer = false
bool m_notUpgradeCabling
std::string m_infoName
const TileInfom_tileInfo
int m_i_ADCmax
float m_f_ADCmax
int m_i_ADCmaxPlus1
float m_f_ADCmaxPlus1
float m_ADCmaxMinusEps
float m_ADCmaskValueMinusEps
 indicates channels which were masked in background dataset

Static Protected Attributes

static const int MAX_CHANNELS = 48
static const int MAX_DMUS = 16

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

void fill_drawer_errors (const EventContext &ctx, const TileDigitsCollection *collection)
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

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

Friends

class TileHid2RESrcID

Detailed Description

Definition at line 59 of file TileRawChannelBuilder.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TileRawChannelBuilder()

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

Constructor.

Definition at line 48 of file TileRawChannelBuilder.cxx.

50 : AthAlgTool(type, name, parent)
53 , m_bsflags(0)
54 , m_tileID(nullptr)
55 , m_tileHWID(nullptr)
56 , m_trigType(0)
57 , m_idophys(false)
58 , m_idolas(false)
59 , m_idoped(false)
60 , m_idocis(false)
61 , m_cischan(-1)
62 , m_capdaq(100)
63 , m_evtCounter(0)
64 , m_chCounter(0)
65 , m_nChL(0)
66 , m_nChH(0)
67 , m_RChSumL(0.0)
68 , m_RChSumH(0.0)
70 , m_tileInfo(nullptr)
71{
73 memset(m_error, 0, sizeof(m_error));
74
75 declareProperty("calibrateEnergy", m_calibrateEnergy = false);
76 declareProperty("correctTime", m_correctTime = false);
77 declareProperty("AmpMinForAmpCorrection", m_ampMinThresh = 15.0);
78 declareProperty("TimeMinForAmpCorrection", m_timeMinThresh = -12.5);
79 declareProperty("TimeMaxForAmpCorrection", m_timeMaxThresh = 12.5);
80 declareProperty("RunType", m_runType = 0);
81 declareProperty("DataPoolSize", m_dataPoollSize = -1);
82 declareProperty("UseDSPCorrection", m_useDSP = true);
83 declareProperty("TileInfoName", m_infoName = "TileInfo");
84 declareProperty("FirstSample",m_firstSample = 0);
85
86}
AthAlgTool()
Default constructor:
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
TileRawChannelUnit::UNIT m_rChUnit
float m_timeMinThresh
correct amplitude is time is above time min threshold
float m_ampMinThresh
correct amplitude if it's above amplitude threshold (in ADC counts)
float m_timeMaxThresh
correct amplitude is time is below time max threshold

◆ ~TileRawChannelBuilder()

TileRawChannelBuilder::~TileRawChannelBuilder ( )
virtual

Destructor.

Definition at line 91 of file TileRawChannelBuilder.cxx.

91 {
92}

Member Function Documentation

◆ BadPatternName()

const char * TileRawChannelBuilder::BadPatternName ( float ped)
static

Definition at line 457 of file TileRawChannelBuilder.cxx.

457 {
458 static const char * const errname[26] = {
459 "-10 - good signal",
460 "-9 - underflow",
461 "-8 - overflow",
462 "-7 - underflow and overflow",
463 "-6 - constant signal",
464 "-5 - disconnected channel",
465 "-4 - half a drawer masked",
466 "-3 - bad DQ status",
467 "-2 - underflow in all samples",
468 "-1 - overflow in all samples",
469 "0 - unknown error",
470 "1 - jump from zero to saturation",
471 "2 - samples with zeros",
472 "3 - at least two saturated. others - close to pedestal",
473 "4 - two distinct levels with at least 2 samples each",
474 "5 - pedestal with jump up in one sample",
475 "6 - pedestal with jump down in one sample",
476 "7 - signal with jump up in one sample",
477 "8 - signal with jump down in one sample",
478 "9 - base line above threshold in low gain",
479 "10 - jump down in first sample in low gain",
480 "11 - jump down in last sample in low gain",
481 "12 - jump up in one sample above const",
482 "13 - jump down in one sample below const",
483 "14 - unrecoverable timing jump",
484 "15 - unknown error"
485 };
486
487 return errname[std::min(25, std::max(0, int((ped + 500) * 1e-4)))];
488}

◆ build()

StatusCode TileRawChannelBuilder::build ( const TileDigitsCollection * collection,
const EventContext & ctx )

Definition at line 491 of file TileRawChannelBuilder.cxx.

492{
493
494 int frag = coll->identify();
495
496 // make sure that error array is up-to-date
497 if (frag != m_lastDrawer && m_notUpgradeCabling) {
498 fill_drawer_errors(ctx, coll);
499 }
500
501 // Iterate over all digits in this collection
502 TileDigitsCollection::const_iterator digitItr = coll->begin();
503 TileDigitsCollection::const_iterator lastDigit = coll->end();
504
505 for (; digitItr != lastDigit; ++digitItr) {
506
507 TileRawChannel* rch = rawChannel((*digitItr), ctx);
508
510
511 int err = m_error[m_tileHWID->channel(rch->adc_HWID())];
512
513 if (err) {
514 if (err == -8 || err == -7) m_overflows.push_back(std::make_pair(rch, (*digitItr)));
515 float ped = rch->pedestal() + 100000 + 10000 * err;
516 rch->setPedestal(ped);
517 if (msgLvl(MSG::VERBOSE) && !m_badDrawer) {
518 if (err < -5) {
519 msg(MSG::VERBOSE) << "BadCh " << m_tileHWID->to_string(rch->adc_HWID())
520 << " warning = " << BadPatternName(ped) << endmsg;
521 } else {
522 msg(MSG::VERBOSE) << "BadCh " << m_tileHWID->to_string(rch->adc_HWID())
523 << " error = " << BadPatternName(ped) << endmsg;
524 }
525 }
526 }
527
528 }
529
530 ATH_CHECK( m_rawChannelCnt->push_back (rch) );
531 }
532
533 IdentifierHash hash = m_rawChannelCnt->hashFunc().hash(coll->identify());
534 TileRawChannelCollection* rawChannelCollection = m_rawChannelCnt->indexFindPtr(hash);
535 rawChannelCollection->setLvl1Id(coll->getLvl1Id());
536 rawChannelCollection->setLvl1Type(coll->getLvl1Type());
537 rawChannelCollection->setDetEvType(coll->getDetEvType());
538 rawChannelCollection->setRODBCID(coll->getRODBCID());
539
540 return StatusCode::SUCCESS;
541}
#define endmsg
#define ATH_CHECK
Evaluate an expression and check for errors.
bool msgLvl(const MSG::Level lvl) const
MsgStream & msg() const
DataModel_detail::const_iterator< DataVector > const_iterator
Standard const_iterator.
Definition DataVector.h:838
static const char * BadPatternName(float ped)
std::unique_ptr< TileMutableRawChannelContainer > m_rawChannelCnt
void fill_drawer_errors(const EventContext &ctx, const TileDigitsCollection *collection)
virtual TileRawChannel * rawChannel(const TileDigits *digits, const EventContext &ctx)
Builder virtual method to be implemented by subclasses.
float pedestal(void) const
void setPedestal(float ped)
void setLvl1Type(uint32_t lvl1Type)
Setter for level1 type from ROD header.
void setLvl1Id(uint32_t lvl1Id)
Setter for level1 id from ROD header.
void setRODBCID(uint32_t rodBCID)
Setter for BCID from ROD header.
void setDetEvType(uint32_t detEvType)
Setter for detector event type from ROD header.
HWIdentifier adc_HWID(void) const
Definition TileRawData.h:53

◆ commitContainer()

StatusCode TileRawChannelBuilder::commitContainer ( const EventContext & ctx)
virtual

Commit RawChannelContiner in SG and make const.

Definition at line 543 of file TileRawChannelBuilder.cxx.

544{
545
546 const TileDQstatus* DQstatus = SG::makeHandle (m_DQstatusKey, ctx).get();
547
548 ToolHandleArray<ITileRawChannelTool>::iterator itrTool = m_noiseFilterTools.begin();
549 ToolHandleArray<ITileRawChannelTool>::iterator endTool = m_noiseFilterTools.end();
550
551 if ( m_useDSP && !m_DSPContainerKey.key().empty() &&
552 (DQstatus->incompleteDigits() || m_chCounter<12288) && itrTool!=endTool )
553 {
554 const TileRawChannelContainer * dspCnt = SG::makeHandle (m_DSPContainerKey, ctx).get();
555 ATH_MSG_DEBUG( "Incomplete container - use noise filter corrections from DSP container" );
556
557 uint32_t bsFlags = dspCnt->get_bsflags();
558 std::vector<IdentifierHash> hashes = m_rawChannelCnt->GetAllCurrentHashes();
559 std::vector<IdentifierHash> dspHashes = dspCnt->GetAllCurrentHashes();
560 if (bsFlags == 0) {
561 ATH_MSG_WARNING("Problem in applying noise corrections: DSP container ("
562 << m_DSPContainerKey.key() << ") seems to be emtpy!");
563 } else if (hashes != dspHashes) {
564 ATH_MSG_ERROR( " Error in applying noise corrections; "
565 "hash vectors do not match.");
566 } else {
567 // Go through all TileRawChannelCollections
568 for (IdentifierHash hash : hashes) {
569 TileRawChannelCollection* coll = m_rawChannelCnt->indexFindPtr (hash);
570 const TileRawChannelCollection* dcoll = dspCnt->indexFindPtr (hash);
571
572 if (coll->identify() != dcoll->identify()) {
573
574 ATH_MSG_ERROR( " Error in applying noise corrections " << MSG::hex
575 << " collection IDs 0x" << coll->identify() << " and 0x" << dcoll->identify()
576 << " do not match " << MSG::dec );
577 break;
578 }
579
580 // iterate over all channels in a collection
583
584 for (TileRawChannel* rch : *coll) {
585 HWIdentifier adc_id = rch->adc_HWID();
586 while (dspItr != dspLast && adc_id != (*dspItr)->adc_HWID()) {
587 ++dspItr;
588 }
589 if (dspItr != dspLast) {
590 float corr = (*dspItr)->pedestal();
591 ATH_MSG_VERBOSE( "Ch "<<m_tileHWID->to_string(adc_id)
592 <<" amp " << rch->amplitude() << " ped " << rch->pedestal()
593 << " corr " << corr );
594 if (corr<10000.) {
595 rch->setAmplitude (rch->amplitude() - corr); // just baseline shift
596 rch->setPedestal (rch->pedestal() + corr); // just baseline shift
597 } else {
598 float ped = rch->pedestal();
599 if (corr > ped) {
600 rch->setPedestal (fmod(ped,10000.) + int(corr)/10000 * 10000); // changing error status
601 ATH_MSG_VERBOSE( "New error status in ped "<<rch->pedestal());
602 }
603 }
604 } else {
605 ATH_MSG_WARNING(" Problem in applying noise corrections "
606 << " can not find channel in DSP container with HWID "
607 << m_tileHWID->to_string(adc_id) );
608 dspItr = dcoll->begin();
609 }
610 }
611 }
612 }
613
614 } else {
615
616 for (ToolHandle<ITileRawChannelTool>& noiseFilterTool : m_noiseFilterTools) {
617 if (noiseFilterTool->process(*m_rawChannelCnt.get(), ctx).isFailure()) {
618 ATH_MSG_ERROR( " Error status returned from noise filter " );
619 } else {
620 ATH_MSG_DEBUG( "Noise filter applied to the container" );
621 }
622 }
623
624 }
625
626 ATH_MSG_DEBUG( " nCh=" << m_chCounter
627 << " nChH/L=" << m_nChH << "/" << m_nChL
628 << " RChSumH/L=" << m_RChSumH << "/" << m_RChSumL );
629
630 SG::WriteHandle<TileRawChannelContainer> rawChannelsContainer(m_rawChannelContainerKey, ctx);
631 ATH_CHECK( rawChannelsContainer.record(std::move(m_rawChannelCnt)) );
632
633 endLog();
634
635 return StatusCode::SUCCESS;
636}
#define ATH_MSG_ERROR(x)
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
Athena::TPCnvVers::Old Athena::TPCnvVers::Old Athena::TPCnvVers::Current Athena::TPCnvVers::Old Athena::TPCnvVers::Current TileRawChannelContainer
const_iterator end() const noexcept
Return a const_iterator pointing past the end of the collection.
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
virtual std::vector< IdentifierHash > GetAllCurrentHashes() const override final
Returns a collection of all hashes availiable in this IDC.
virtual const T * indexFindPtr(IdentifierHash hashId) const override final
return pointer on the found entry or null if out of range using hashed index - fast version,...
bool incompleteDigits() const
A few extra items (from TileBeamInfoProvider).
SG::ReadHandleKey< TileDQstatus > m_DQstatusKey
SG::ReadHandleKey< TileRawChannelContainer > m_DSPContainerKey
ToolHandleArray< ITileRawChannelTool > m_noiseFilterTools
SG::WriteHandleKey< TileRawChannelContainer > m_rawChannelContainerKey
uint32_t get_bsflags() const
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
setEventNumber uint32_t

◆ correctAmp()

double TileRawChannelBuilder::correctAmp ( double phase,
bool of2 = true )
static

Amplitude correction factor according to the time when using weights for tau=0 without iterations.

Definition at line 645 of file TileRawChannelBuilder.cxx.

645 {
646
647 double corr = 1.0;
648 if (of2) {
649 // estimation from Belen for rel 14.0.0
650 /*double a,b,c;
651 if(fabs(phase)<5.){
652 a=0.137; b=0.0877; c=0.0865;
653 }else{
654 a=0.565; b=0.116; c=0.0751;
655 }
656 corr=(1+(a+b*phase+c*phase*phase)/100.);
657 */
658
659 // estimation from Vakhtang for rel 14.4.0
660 /*double k = (phase < 0.0 ? 0.0009400 : 0.0010160);
661 corr = (1.0 + k * phase * phase);
662 */
663
664 // Parabolic correction from Tigran
665 double a1,a2,b,c;
666 a1 = phase < 0.0 ? 0.000940774 : 0.00102111;
667 a2 = phase < 0.0 ? 0.000759051 : 0.000689625;
668 b = phase < 0.0 ? -2.0 * 7.0 * (a1 - a2) : 2.0 * 12.5 * (a1 - a2);
669 c = phase < 0.0 ? 1.0 - 7.0 * 7.0 * (a1-a2) : 1.0 - 12.5 * 12.5 * (a1-a2);
670 if (phase < 12.5 && phase > -7.0) corr = a1 * phase * phase + 1.0;
671 else corr = phase * ( a2 * phase + b) + c;
672
673
674 } else {
675 /*double a,b,c;
676 if(phase<0){
677 a=1.0002942; b=0.0003528; c=0.0005241;
678 }else{
679 a=1.0001841; b=-0.0004182; c=0.0006167;
680 }
681 corr = a + phase * ( b + c * phase);
682 */
683
684 /*double k = (phase < 0.0 ? 0.0005241 : 0.0006167);
685 corr = (1.0 + k * phase * phase);
686 */
687
688 // 4th degree polynomial correction from Tigran
689 double k1 = (phase < 0.0 ? -0.0000326707:0.000380336);
690 double k2 = (phase < 0.0 ? -0.000560962:-0.000670487);
691 double k3 = (phase < 0.0 ? -0.00000807869:0.00000501773);
692 double k4 = (phase < 0.0 ? -0.000000145008:0.0000000584647);
693
694 corr = 1.0 / (1.0 + (k1 + (k2 + (k3 + k4 *phase)*phase)*phase)*phase);
695
696
697 }
698
699 return corr;
700}

◆ correctTime()

double TileRawChannelBuilder::correctTime ( double phase,
bool of2 = true )
static

Time correction factor.

Definition at line 704 of file TileRawChannelBuilder.cxx.

704 {
705
706 double correction = 0.0;
707
708 if (of2) {
709 if(phase < 0) {
710 correction = (-0.00695743 + (0.0020673 - (0.0002976 + 0.00000361305 * phase) * phase) * phase) * phase;
711 } else {
712 correction = (0.0130013 + (0.00128769 + (-0.000550218 + 0.00000755344 * phase) * phase) * phase) * phase;
713 }
714 }
715 // OF1 does not need correction
716
717 return correction;
718}
correction(mu, runmode, campaign, run=None)
Definition zlumi_mc_cf.py:4

◆ CorruptedData()

int TileRawChannelBuilder::CorruptedData ( int ros,
int drawer,
int channel,
int gain,
const std::vector< float > & digits,
float & dmin,
float & dmax,
float ADCmaxMinusEps,
float ADCmaskValueMinusEps )
static

Definition at line 722 of file TileRawChannelBuilder.cxx.

723 {
724 bool eb = (ros > 2);
725 bool ebsp = ((ros == 3 && drawer == 14) || (ros == 4 && drawer == 17));
726 bool empty = ((eb && ((channel > 23 && channel < 30) || channel > 41)) || (ebsp && channel < 3));
727 bool not_gap = !(empty || (eb && (channel == 0 || channel == 1 || channel == 12 || channel == 13))
728 || (ebsp && (channel == 18 || channel == 19)));
729
730 const float epsilon = 4.1; // allow +/- 2 counts fluctuations around const value
731 const float delta[4] = { 29.9, 29.9, 49.9, 99.9 }; // jump levels between constLG, constHG, non-constLG, non-constHG
732 const float level1 = 99.9; // jump from this level to m_i_ADCmax is bad
733 const float level2 = 149.9; // base line at this level in low gain is bad
734 const float narrowLevel[2] = { 29.9, 49.9 }; // minimal amplitude for narrow pulses
735 const float delt = std::min(std::min(std::min(delta[0], delta[1]), std::min(delta[2], delta[3])),
736 std::min(narrowLevel[0], narrowLevel[1]));
737 const float secondMaxLevel = 0.3;
738
739 int error = 0;
740
741 unsigned int nSamp = digits.size();
742 if (nSamp) {
743 dmin = dmax = digits[0];
744 unsigned int pmin = 0;
745 unsigned int pmax = 0;
746 unsigned int nzero = (dmin < 0.01) ? 1 : 0;
747 unsigned int nover = (dmax > ADCmaxMinusEps) ? 1 : 0;
748
749 for (unsigned int i = 1; i < nSamp; ++i) {
750 float dig = digits[i];
751 if (dig > dmax) {
752 dmax = dig;
753 pmax = i;
754 } else if (dig < dmin) {
755 dmin = dig;
756 pmin = i;
757 }
758 if (dig < 0.01) ++nzero;
759 else if (dig > ADCmaxMinusEps) ++nover;
760 }
761
762 float dmaxmin = dmax - dmin;
763 //std::cout << " ros " << ros << " drawer " << drawer << " channel " << channel << " not_gap " << not_gap << " nzero " << nzero << " nover " << nover << std::endl;
764
765 if (dmin > ADCmaxMinusEps) { // overflow in all samples
766 error = (dmin > ADCmaskValueMinusEps) ? -3 : -1; // dmin=m_tileInfo->ADCmaskValue() - masking in overlay job (set in TileDigitsMaker)
767
768 } else if (dmax < 0.01) { // underflow in all samples
769 error = (empty) ? -5 : -2; // set different type of errors for exsiting and non-existing channels
770
771 } else if (dmaxmin < 0.01) { // constant value in all samples
772 error = -6;
773
774 } else if (nzero && nover) { // jump from zero to saturation
775 error = 1;
776
777 } else if ((nzero && (not_gap || empty)) || nzero > 1) { // one sample at zero in normal channel
778 error = 2; // or 2 samples at zero in gap/crack/MBTS
779
780 } else if (gain == 0 && dmin > level2) { // baseline above threshold in low gain is bad
781 error = 9;
782
783 } else if (dmaxmin > delt) { // check that max-min is above minimal allowed jump
784
785 float abovemin = dmax;
786 float belowmax = dmin;
787 unsigned int nmin = 0;
788 unsigned int nmax = 0;
789 for (unsigned int i = 0; i < nSamp; ++i) {
790 float smp = digits[i];
791 if (smp - dmin < epsilon) {
792 ++nmin;
793 }
794 if (dmax - smp < epsilon) {
795 ++nmax;
796 }
797 if (smp < abovemin && smp > dmin) {
798 abovemin = smp;
799 }
800 if (smp > belowmax && smp < dmax) {
801 belowmax = smp;
802 }
803 }
804 // more than two different values - shift index by 2, i.e. use thresholds for non-const levels
805 int gainInd = (abovemin != dmax || belowmax != dmin) ? gain + 2 : gain;
806 bool big_jump = (dmaxmin > delta[gainInd]);
807 bool max_in_middle = (pmax > 0 && pmax < nSamp - 1);
808 bool min_in_middle = (pmin > 0 && pmin < nSamp - 1);
809
810 if (nover > 1 && belowmax < level1) { // at least two saturated. others - close to pedestal
811 error = 3;
812 } else if (nmax + nmin == nSamp && big_jump) {
813 if (nmax > 1 && nmin > 1) { // at least 2 samples at two distinct levels
814 error = 4;
815 } else if (nmax == 1) {
816 if (max_in_middle) { // jump up in one sample, but not at the edge
817 error = 5;
818 }
819 } else if (nmin == 1) { // jump down in one sample
820 error = 6;
821 }
822 }
823 if (error == 0 && dmaxmin > narrowLevel[gain]) {
824 float secondMax = dmaxmin * secondMaxLevel;
825 float dminPlus = dmin + secondMax;
826 float dmaxMinus = dmax - secondMax;
827 if (not_gap) { // jumps above two (or one) neighbour samples
828 if (max_in_middle && std::max(digits[pmax - 1], digits[pmax + 1]) < dminPlus) {
829 error = 7; // jump up in one sample in the middle, which is much higher than two neighbours
830 } else if (min_in_middle && std::min(digits[pmin - 1], digits[pmin + 1]) > dmaxMinus) {
831 error = 8; // jump down in one sample, which is much lower than two neighbours
832 } else if (big_jump && gain == 0) { // check first and last sample only in low gain
833 if (pmin == 0 && digits[1] > dmax - secondMax) {
834 error = 10; // jump down in first sample. which is much lower than next one
835 } else if (pmin == nSamp - 1 && digits[pmin - 1] > dmax - secondMax) {
836 error = 11; // jump down in last sample. which is much lower than previous one
837 }
838 }
839 }
840 if (!error && big_jump) { // jumps above all samples
841 if ((max_in_middle || gain == 0) && nmax == 1 && belowmax < dminPlus) {
842 error = 12; // jump up in one sample in the middle, which is much higher than all others
843 } else if ((min_in_middle || gain == 0) && nmin == 1 && abovemin > dmaxMinus) {
844 error = 13; // jump down in one sample, which is much lower than all others (
845 }
846 }
847 }
848 }
849
850 } else {
851 dmin = dmax = 0.0;
852 }
853
854 return error;
855}
const int nmax(200)
static const Attributes_t empty

◆ createContainer()

StatusCode TileRawChannelBuilder::createContainer ( const EventContext & ctx)
virtual

Create container in SG with name given by parameter (m_rawChannelContainerKey).

Definition at line 232 of file TileRawChannelBuilder.cxx.

232 {
233 initLog(ctx);
234
235 // create TRC container
236 m_rawChannelCnt = std::make_unique<TileMutableRawChannelContainer>(true, m_rChType, m_rChUnit, SG::VIEW_ELEMENTS);
237 ATH_CHECK( m_rawChannelCnt->status() );
238 m_rawChannelCnt->set_bsflags(m_bsflags);
239
240 ATH_MSG_DEBUG( "Created TileRawChannelContainer '" << m_rawChannelContainerKey.key() << "'" );
241
242 return StatusCode::SUCCESS;
243}
void initLog(const EventContext &ctx)
@ VIEW_ELEMENTS
this data object is a view, it does not own its elmts

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::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.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }

◆ declareProperty()

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

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

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

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

Definition at line 95 of file AthCommonDataStore.h.

◆ endLog()

void TileRawChannelBuilder::endLog ( )

Definition at line 638 of file TileRawChannelBuilder.cxx.

638 {
639 m_chCounter = 0;
640 m_nChL = m_nChH = 0;
641 m_RChSumL = m_RChSumH = 0.0;
642
643}

◆ evtStore()

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

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

Definition at line 85 of file AthCommonDataStore.h.

◆ extraDeps_update_handler()

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

Add StoreName to extra input/output deps as needed.

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

◆ fill_drawer_errors()

void TileRawChannelBuilder::fill_drawer_errors ( const EventContext & ctx,
const TileDigitsCollection * collection )
private

Definition at line 299 of file TileRawChannelBuilder.cxx.

301{
302 const TileDQstatus* DQstatus = SG::makeHandle (m_DQstatusKey, ctx).get();
303
304 int frag = coll->identify();
305 int ros = (frag >> 8);
306 int drawer = (frag & 0xff);
307
308 m_lastDrawer = frag;
309
310 memset(m_error, 0, sizeof(m_error));
311 int dmuerr[MAX_DMUS] = {0};
312 int nch = 0;
313 bool bigain = DQstatus->isBiGain();
314 if (!bigain) { // in bigain runs we don't have DQ status fragment
315 for (int ch = 0; ch < MAX_CHANNELS; ch += 3) {
316 if (!DQstatus->isAdcDQgood(ros, drawer, ch, 0)) {
317 m_error[ch + 2] = m_error[ch + 1] = m_error[ch] = -3;
318 dmuerr[ch / 3] = 3;
319 nch += 3;
320 }
321 }
322 }
323 if (nch == MAX_CHANNELS) { // all bad - nothing to do
324 m_badDrawer = true;
325 ATH_MSG_VERBOSE( "Drawer 0x" << MSG::hex << frag << MSG::dec
326 << " is bad - skipping bad patterns check " );
327 return;
328 } else {
329 m_badDrawer = false;
330 ATH_MSG_VERBOSE( "Drawer 0x" << MSG::hex << frag << MSG::dec
331 << " looking for bad patterns in digits" );
332 }
333
334 float mindig, maxdig;
335 int nchbad[2] = { 0, 0 };
336
337 // Iterate over all digits in this collection
338 TileDigitsCollection::const_iterator digitItr = coll->begin();
339 TileDigitsCollection::const_iterator lastDigit = coll->end();
340
341 for (; digitItr != lastDigit; ++digitItr) {
342 const TileDigits * pDigits = (*digitItr);
343 HWIdentifier adcId = pDigits->adc_HWID();
344 int channel = m_tileHWID->channel(adcId);
345 int gain = m_tileHWID->adc(adcId);
346
347 if (m_error[channel]) {
348 ATH_MSG_VERBOSE( "BadCh " << ros
349 << "/" << drawer
350 << "/" << channel
351 << "/" << gain << " BAD DQ STATUS ");
352
353 } else {
354
355 int err = CorruptedData(ros, drawer, channel, gain, pDigits->samples(), mindig, maxdig, m_ADCmaxMinusEps, m_ADCmaskValueMinusEps);
356
357 if (err) {
358
360 if (err > -5) {
361 ++dmuerr[channel / 3];
362 ++nchbad[channel / 24];
363 }
364
365 if (msgLvl(MSG::VERBOSE)) {
366
367 msg(MSG::VERBOSE) << "BadCh " << ros
368 << "/" << drawer
369 << "/" << channel
370 << "/" << gain;
371 if (err < -5) msg(MSG::VERBOSE) << " Warning " << err;
372 else msg(MSG::VERBOSE) << " Error " << err;
373 if (mindig > m_ADCmaskValueMinusEps) msg(MSG::VERBOSE) << " BADDQ";
374 if (maxdig > m_ADCmaxMinusEps) msg(MSG::VERBOSE) << " Overflow";
375 if (mindig < 0.1) msg(MSG::VERBOSE) << " Underflow";
376 if (err < 0) msg(MSG::VERBOSE) << " Const";
377
378 msg(MSG::VERBOSE) << " samp=";
379 std::vector<float> digits = pDigits->samples();
380 for (unsigned int i = 0; i < digits.size(); ++i) {
381 msg(MSG::VERBOSE) << " " << digits[i];
382 }
383 msg(MSG::VERBOSE) << endmsg;
384 }
385
386 } else {
387 if (mindig < 0.01) err += 1;
388 if (maxdig > m_ADCmaxMinusEps) err += 2;
389 if (err) m_error[channel] = err - 10;
390 }
391 }
392 }
393
394 // check if we want to mask half a drawer
395 // in this case set error = -4 for channels which were good before
396
397 int ndmubad[2] = { 0, 0 };
398 int dmu = 0;
399 for (; dmu < MAX_DMUS / 2; ++dmu) { // first half
400 if (dmuerr[dmu] > 1)
401 ++ndmubad[0]; // count DMUs with at least two bad channels
402 }
403 for (; dmu < MAX_DMUS; ++dmu) { // second half
404 if (dmuerr[dmu] > 1)
405 ++ndmubad[1]; // count DMUs with at least two bad channels
406 }
407
408 int ndmulimit[2] = { 3, 3 }; // max number of bad DMUs when half-drawer is not yet masked
409 // if 4 DMUs will be bad - mask whole half-drawer
410 if (frag > 0x2ff) { // if extended barrel
411 if (frag == 0x30e || frag == 0x411)
412 ndmulimit[0] = 4; // in EB special one DMU is always bad (missing)
413 ndmulimit[1] = 5; // in second half of EB 4 DMUs ara always bad (missing)
414 // only if 7 DMUs are bad, mask whole half-drawer
415 }
416
417 bool printall = true;
418 for (int p = 0; p < 2; ++p) {
419 if (ndmubad[p] > ndmulimit[p] && nchbad[p] > 0) {
420 if (msgLvl(MSG::VERBOSE)) {
421 msg(MSG::VERBOSE) << "Drawer 0x" << MSG::hex << frag << MSG::dec
422 << " masking whole " << ((p) ? "second" : "first")
423 << " half" << endmsg;
424 if (printall) {
425 msg(MSG::VERBOSE) << "nDMuErr ";
426 for (int d = 0; d < MAX_DMUS; ++d) {
427 msg(MSG::VERBOSE) << " " << dmuerr[d];
428 }
429 msg(MSG::VERBOSE) << " total " << ndmubad[p] << " errors" << endmsg;
430
431 msg(MSG::VERBOSE) << "ChErr ";
432 int ch = 0;
433 while (ch < MAX_CHANNELS-2) {
434 msg(MSG::VERBOSE) << " " << m_error[ch++];
435 msg(MSG::VERBOSE) << " " << m_error[ch++];
436 msg(MSG::VERBOSE) << " " << m_error[ch++];
437 msg(MSG::VERBOSE) << " ";
438 }
439
440 msg(MSG::VERBOSE) << " total " << nchbad[p]
441 << " bad patterns" << endmsg;
442 printall = false;
443 }
444 }
445 int ch = (p) ? MAX_CHANNELS / 2 : 0;
446 int chmax = (p) ? MAX_CHANNELS : MAX_CHANNELS / 2;
447 for (; ch < chmax; ++ch) {
448 if (m_error[ch] == 0 || m_error[ch] < -5) { // channel was good before
449 m_error[ch] = -4;
450 }
451 }
452 }
453 }
454
455}
bool isAdcDQgood(int partition, int drawer, int ch, int gain) const
returns status of single ADC returns False if there are any errors
bool isBiGain() const
returns gain mode of run
const std::vector< float > & samples() const
Definition TileDigits.h:58
float m_ADCmaskValueMinusEps
indicates channels which were masked in background dataset
static int CorruptedData(int ros, int drawer, int channel, int gain, const std::vector< float > &digits, float &dmin, float &dmax, float ADCmaxMinusEps, float ADCmaskValueMinusEps)

◆ finalize()

StatusCode TileRawChannelBuilder::finalize ( )
virtual

◆ getOverflowedChannels()

Overflows_t & TileRawChannelBuilder::getOverflowedChannels ( void )

Definition at line 37 of file TileRawChannelBuilder.cxx.

37 {
38 return m_overflows;
39}

◆ getTileRawChannelContainerID()

std::string TileRawChannelBuilder::getTileRawChannelContainerID ( void )

Definition at line 41 of file TileRawChannelBuilder.cxx.

41 {
42 return m_rawChannelContainerKey.key();
43}

◆ initialize()

StatusCode TileRawChannelBuilder::initialize ( )
virtual

Initializer.

Reimplemented in TileRawChannelBuilderFitFilter, TileRawChannelBuilderFitFilterCool, TileRawChannelBuilderFlatFilter, TileRawChannelBuilderManyAmps, TileRawChannelBuilderMF, TileRawChannelBuilderOpt2Filter, TileRawChannelBuilderQIEFilter, and TileRawChannelBuilderWienerFilter.

Definition at line 97 of file TileRawChannelBuilder.cxx.

97 {
98
99 ATH_MSG_INFO( "TileRawChannelBuilder::initialize()" );
100
102 m_idophys = ((m_trigType == 0) || (m_trigType == 1));
103 m_idolas = ((m_trigType == 2) || (m_trigType == 3));
104 m_idoped = ((m_trigType == 4) || (m_trigType == 5));
105 m_idocis = ((m_trigType == 8) || (m_trigType == 9));
106 m_cischan = -1;
107 m_capdaq = 0;
108 m_chCounter = 0;
109 m_evtCounter = 0;
110 m_rawChannelCnt = nullptr;
111 m_nChL = m_nChH = 0;
112 m_RChSumL = m_RChSumH = 0.0;
113 m_evtCounter = -1;
114 // retrieve TileID helpers and TileIfno from det store
115 ATH_CHECK( detStore()->retrieve(m_tileID, "TileID") );
116 ATH_CHECK( detStore()->retrieve(m_tileHWID, "TileHWID") );
117
118 ATH_CHECK( detStore()->retrieve(m_tileInfo, m_infoName) );
119 m_i_ADCmax = m_tileInfo->ADCmax();
124 m_ADCmaskValueMinusEps = m_tileInfo->ADCmaskValue() - 0.01; // indicates channels which were masked in background dataset
125
126 // access tools and store them
127 ATH_CHECK( m_noiseFilterTools.retrieve() );
128 ATH_MSG_DEBUG( "Successfully retrieve NoiseFilterTools: " << m_noiseFilterTools );
129
131 ATH_MSG_DEBUG( "Obsolete calibrateEnergy flag is set to True in jobOptions - disabling it" );
132 m_calibrateEnergy = false;
133 }
134
135 // check if we want to keep ADC counts or convert them to pCb
137 : TileRawChannelUnit::ADCcounts;
138
139 // if unit is not pCb, but already MeV one can use method TileRawChannelContainer::set_unit()
140 // later to declare that
141
142 // 8 upper bits of bsflags:
143 // UUPPSTTT
144 // 31,30 - units
145 // 29,28 - pulse type = 3 for offline reco
146 // 27 - 7(=0) or 9(=1) samples
147 // 24,25,26 - TileFragHash::TYPE - OF algorithm type
148 int nsamp = (m_tileInfo->NdigitSamples() > 7) ? 1 : 0;
149 m_bsflags = (m_rChUnit << 30) | (3 << 28) | (nsamp << 27) | (m_rChType << 24);
150
151 // bits 12-15 - various options
152 if (m_correctTime) m_bsflags |= 0x1000;
153
154 if (msgLvl(MSG::DEBUG)) {
155 msg(MSG::DEBUG) << "TileRawChannelBuilder created, storing rc in '"
156 << m_rawChannelContainerKey.key() << "'" << endmsg;
157 msg(MSG::DEBUG) << " calibrate energy = " << m_calibrateEnergy << endmsg;
158 msg(MSG::DEBUG) << " correct time = " << m_correctTime << endmsg;
159 msg(MSG::DEBUG) << " run type = " << m_runType << endmsg;
160 }
161
162 if (m_dataPoollSize < 0) m_dataPoollSize = m_tileHWID->channel_hash_max();
163
164 ATH_CHECK( m_cablingSvc.retrieve());
165
166 m_cabling = m_cablingSvc->cablingService();
167 if (!m_cabling) {
168 ATH_MSG_ERROR( "Unable to retrieve TileCablingService" );
169 return StatusCode::FAILURE;
170 }
171
173
174 int runPeriod = m_cabling->runPeriod();
175 std::ostringstream os;
176 if (runPeriod==3) {
177 if ( m_demoFragIDs.empty() ) {
178 std::vector<int> v = { 0x10d }; // LBA14 is demonstrator in RUN3
180 }
181 os << " in RUN3";
182 }
183
184 if ( !m_demoFragIDs.empty() ) {
186 os << " (frag IDs):";
187 for (int fragID : m_demoFragIDs) {
188 if (fragID>0)
189 os << " 0x" << std::hex << fragID << std::dec;
190 else
191 os << " " << fragID;
192 }
193 ATH_MSG_INFO("Enable special treatment for demonstrator modules" << os.str());
194 }
195
196 if (m_calibrateEnergy) {
197 ATH_CHECK( m_tileToolEmscale.retrieve() );
198 } else {
199 m_tileToolEmscale.disable();
200 }
201
202 if (m_correctTime) {
203 ATH_CHECK( m_tileToolTiming.retrieve() );
204 } else {
205 m_tileToolTiming.disable();
206 }
207
209 ATH_CHECK( m_tileIdTransforms.retrieve() );
210 } else {
211 m_tileIdTransforms.disable();
212 }
213
216
217 if (m_useDSP && !m_DSPContainerKey.key().empty()) {
218 ATH_CHECK( m_DSPContainerKey.initialize() );
219 }
220 else {
222 }
223
224 return StatusCode::SUCCESS;
225}
const ServiceHandle< StoreGateSvc > & detStore() const
ToolHandle< TileCondToolTiming > m_tileToolTiming
ServiceHandle< TileCablingSvc > m_cablingSvc
Name of Tile cabling service.
ToolHandle< TileCondIdTransforms > m_tileIdTransforms
const TileCablingService * m_cabling
TileCabling instance.
ToolHandle< TileCondToolEmscale > m_tileToolEmscale
Gaudi::Property< std::vector< int > > m_demoFragIDs
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.

◆ initLog()

void TileRawChannelBuilder::initLog ( const EventContext & ctx)

Definition at line 245 of file TileRawChannelBuilder.cxx.

245 {
246
247 const TileDQstatus* DQstatus = SG::makeHandle (m_DQstatusKey, ctx).get();
248
249 // update only if there is new event
250 if (m_evtCounter != ctx.evt()) {
251
252 m_evtCounter = ctx.evt();
253 if (m_runType != 0) m_trigType = m_runType;
254 else m_trigType = DQstatus->trigType();
255
256 if (0 == m_trigType) {
257 m_idophys = (DQstatus->calibMode() == 0);
258 m_idolas = false;
259 m_idoped = false;
260 m_idocis = (DQstatus->calibMode() == 1);
261 } else {
262 m_idophys = (m_trigType <= 1);
263 m_idolas = ((m_trigType == 2) || (m_trigType == 3));
264 m_idoped = ((m_trigType == 4) || (m_trigType == 5));
265 m_idocis = ((m_trigType == 8) || (m_trigType == 9));
266 }
267
268 const unsigned int *cispar = DQstatus->cispar();
269 if (0 == cispar[7]) { // if capdaq not set, it can't be CIS event
270 if (m_idocis) { // cis flag was set incorrectly, change to ped
271 m_idoped = true;
272 m_idocis = false;
273 }
274 m_capdaq = 0.0;
275 } else {
276 m_capdaq = (cispar[7] < 10) ? 5.2 : 100.0;
277 }
278 m_cischan = cispar[8] - 1; // channel where CIS is fired (-1 = all channels)
279
280 ATH_MSG_DEBUG( "Trig type is " << m_trigType
281 << "; dophys is " << ((m_idophys) ? "true" : "false")
282 << "; dolas is " << ((m_idolas) ? "true" : "false")
283 << "; doped is " << ((m_idoped) ? "true" : "false")
284 << "; docis is " << ((m_idocis) ? "true" : "false")
285 << "; capacitor is " << m_capdaq
286 << "; cis chan is " << m_cischan );
287 }
288}
int trigType() const
Trigger type.
uint32_t calibMode() const
Calibration mode.
const uint32_t * cispar() const
CIS parameters.

◆ inputHandles()

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

Return this algorithm's input handles.

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

◆ interfaceID()

const InterfaceID & TileRawChannelBuilder::interfaceID ( )
static

AlgTool InterfaceID.

Definition at line 23 of file TileRawChannelBuilder.cxx.

23 {
25}
static const InterfaceID IID_ITileRawChannelBuilder("TileRawChannelBuilder", 1, 0)

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ outputHandles()

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

Return this algorithm's output handles.

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

◆ rawChannel()

TileRawChannel * TileRawChannelBuilder::rawChannel ( const TileDigits * digits,
const EventContext & ctx )
virtual

Builder virtual method to be implemented by subclasses.

Parameters
digitsPointer to TileDigitsContainer

Reimplemented in TileRawChannelBuilderFitFilter, TileRawChannelBuilderFitFilterCool, TileRawChannelBuilderFlatFilter, TileRawChannelBuilderManyAmps, TileRawChannelBuilderMF, TileRawChannelBuilderOpt2Filter, TileRawChannelBuilderQIEFilter, TileRawChannelBuilderTest, and TileRawChannelBuilderWienerFilter.

Definition at line 290 of file TileRawChannelBuilder.cxx.

290 {
291 ++m_chCounter;
292 ATH_MSG_WARNING( "Default constructor for rawChannel!" );
293 DataPool<TileRawChannel> tileRchPool (100);
294 TileRawChannel *rawCh = tileRchPool.nextElementPtr();
295 rawCh->assign (digits->adc_HWID(), 0.0, 0.0, 0.0);
296 return rawCh;
297}
void assign(const HWIdentifier &id, float amplitude, float time, float quality, float ped=0.0)

◆ renounce()

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

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

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

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ resetDrawer()

void TileRawChannelBuilder::resetDrawer ( )

Definition at line 28 of file TileRawChannelBuilder.cxx.

28 {
29 m_lastDrawer = -1;
30 m_badDrawer = false;
31}

◆ resetOverflows()

void TileRawChannelBuilder::resetOverflows ( void )

Definition at line 33 of file TileRawChannelBuilder.cxx.

33 {
34 m_overflows.clear();
35}

◆ sysInitialize()

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

Perform system initialization for an algorithm.

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

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

◆ sysStart()

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

Handle START transition.

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

◆ updateVHKA()

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

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }
std::vector< SG::VarHandleKeyArray * > m_vhka

◆ TileHid2RESrcID

friend class TileHid2RESrcID
friend

Definition at line 116 of file TileRawChannelBuilder.h.

Member Data Documentation

◆ m_ADCmaskValueMinusEps

float TileRawChannelBuilder::m_ADCmaskValueMinusEps
protected

indicates channels which were masked in background dataset

Definition at line 222 of file TileRawChannelBuilder.h.

◆ m_ADCmaxMinusEps

float TileRawChannelBuilder::m_ADCmaxMinusEps
protected

Definition at line 221 of file TileRawChannelBuilder.h.

◆ m_ampMinThresh

float TileRawChannelBuilder::m_ampMinThresh
protected

correct amplitude if it's above amplitude threshold (in ADC counts)

Definition at line 152 of file TileRawChannelBuilder.h.

◆ m_badDrawer

bool TileRawChannelBuilder::m_badDrawer = false
protected

Definition at line 210 of file TileRawChannelBuilder.h.

◆ m_bsflags

unsigned int TileRawChannelBuilder::m_bsflags
protected

Definition at line 138 of file TileRawChannelBuilder.h.

◆ m_cabling

const TileCablingService* TileRawChannelBuilder::m_cabling
protected

TileCabling instance.

Definition at line 182 of file TileRawChannelBuilder.h.

◆ m_cablingSvc

ServiceHandle<TileCablingSvc> TileRawChannelBuilder::m_cablingSvc
protected
Initial value:
{ this,
"TileCablingSvc", "TileCablingSvc", "The Tile cabling service"}

Name of Tile cabling service.

Definition at line 179 of file TileRawChannelBuilder.h.

179 { this,
180 "TileCablingSvc", "TileCablingSvc", "The Tile cabling service"};

◆ m_calibrateEnergy

bool TileRawChannelBuilder::m_calibrateEnergy
protected

Definition at line 142 of file TileRawChannelBuilder.h.

◆ m_capdaq

double TileRawChannelBuilder::m_capdaq
protected

Definition at line 193 of file TileRawChannelBuilder.h.

◆ m_chCounter

unsigned int TileRawChannelBuilder::m_chCounter
protected

Definition at line 196 of file TileRawChannelBuilder.h.

◆ m_cischan

int TileRawChannelBuilder::m_cischan
protected

Definition at line 192 of file TileRawChannelBuilder.h.

◆ m_correctTime

bool TileRawChannelBuilder::m_correctTime
protected

Definition at line 145 of file TileRawChannelBuilder.h.

◆ m_dataPoollSize

int TileRawChannelBuilder::m_dataPoollSize
protected

Definition at line 204 of file TileRawChannelBuilder.h.

◆ m_demoFragIDs

Gaudi::Property<std::vector<int> > TileRawChannelBuilder::m_demoFragIDs
protected
Initial value:
{this,
"DemoFragIDs", {}, "List of Tile frag IDs with new electronics (demonstrator)"}

Definition at line 184 of file TileRawChannelBuilder.h.

184 {this,
185 "DemoFragIDs", {}, "List of Tile frag IDs with new electronics (demonstrator)"};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

◆ m_DQstatusKey

SG::ReadHandleKey<TileDQstatus> TileRawChannelBuilder::m_DQstatusKey
protected
Initial value:
{this, "TileDQstatus",
"TileDQstatus",
"TileDQstatus key"}

Definition at line 120 of file TileRawChannelBuilder.h.

120 {this, "TileDQstatus",
121 "TileDQstatus",
122 "TileDQstatus key"};

◆ m_DSPContainerKey

SG::ReadHandleKey<TileRawChannelContainer> TileRawChannelBuilder::m_DSPContainerKey {this, "DSPContainer", "", "DSP Container key"}
protected

Definition at line 125 of file TileRawChannelBuilder.h.

126{this, "DSPContainer", "", "DSP Container key"};

◆ m_error

int TileRawChannelBuilder::m_error[MAX_CHANNELS]
protected

Definition at line 208 of file TileRawChannelBuilder.h.

◆ m_evtCounter

unsigned int TileRawChannelBuilder::m_evtCounter
protected

Definition at line 195 of file TileRawChannelBuilder.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

◆ m_f_ADCmax

float TileRawChannelBuilder::m_f_ADCmax
protected

Definition at line 218 of file TileRawChannelBuilder.h.

◆ m_f_ADCmaxPlus1

float TileRawChannelBuilder::m_f_ADCmaxPlus1
protected

Definition at line 220 of file TileRawChannelBuilder.h.

◆ m_firstSample

int TileRawChannelBuilder::m_firstSample
protected

Definition at line 140 of file TileRawChannelBuilder.h.

◆ m_i_ADCmax

int TileRawChannelBuilder::m_i_ADCmax
protected

Definition at line 217 of file TileRawChannelBuilder.h.

◆ m_i_ADCmaxPlus1

int TileRawChannelBuilder::m_i_ADCmaxPlus1
protected

Definition at line 219 of file TileRawChannelBuilder.h.

◆ m_idocis

bool TileRawChannelBuilder::m_idocis
protected

Definition at line 191 of file TileRawChannelBuilder.h.

◆ m_idolas

bool TileRawChannelBuilder::m_idolas
protected

Definition at line 189 of file TileRawChannelBuilder.h.

◆ m_idoped

bool TileRawChannelBuilder::m_idoped
protected

Definition at line 190 of file TileRawChannelBuilder.h.

◆ m_idophys

bool TileRawChannelBuilder::m_idophys
protected

Definition at line 188 of file TileRawChannelBuilder.h.

◆ m_infoName

std::string TileRawChannelBuilder::m_infoName
protected

Definition at line 215 of file TileRawChannelBuilder.h.

◆ m_lastDrawer

int TileRawChannelBuilder::m_lastDrawer = -1
protected

Definition at line 209 of file TileRawChannelBuilder.h.

◆ m_nChH

int TileRawChannelBuilder::m_nChH
protected

Definition at line 199 of file TileRawChannelBuilder.h.

◆ m_nChL

int TileRawChannelBuilder::m_nChL
protected

Definition at line 198 of file TileRawChannelBuilder.h.

◆ m_noiseFilterTools

ToolHandleArray<ITileRawChannelTool> TileRawChannelBuilder::m_noiseFilterTools
protected
Initial value:
{this,
"NoiseFilterTools", {}, "Tile noise filter tools"}

Definition at line 163 of file TileRawChannelBuilder.h.

163 {this,
164 "NoiseFilterTools", {}, "Tile noise filter tools"};

◆ m_notUpgradeCabling

bool TileRawChannelBuilder::m_notUpgradeCabling
protected

Definition at line 212 of file TileRawChannelBuilder.h.

◆ m_overflows

Overflows_t TileRawChannelBuilder::m_overflows
protected

Definition at line 203 of file TileRawChannelBuilder.h.

◆ m_rawChannelCnt

std::unique_ptr<TileMutableRawChannelContainer> TileRawChannelBuilder::m_rawChannelCnt
protected

Definition at line 133 of file TileRawChannelBuilder.h.

◆ m_rawChannelContainerKey

SG::WriteHandleKey<TileRawChannelContainer> TileRawChannelBuilder::m_rawChannelContainerKey
protected
Initial value:
{this,"TileRawChannelContainer","TileRawChannelFiltered",
"Output Tile raw channels container key"}

Definition at line 129 of file TileRawChannelBuilder.h.

129 {this,"TileRawChannelContainer","TileRawChannelFiltered",
130 "Output Tile raw channels container key"};

◆ m_RChSumH

double TileRawChannelBuilder::m_RChSumH
protected

Definition at line 201 of file TileRawChannelBuilder.h.

◆ m_RChSumL

double TileRawChannelBuilder::m_RChSumL
protected

Definition at line 200 of file TileRawChannelBuilder.h.

◆ m_rChType

TileFragHash::TYPE TileRawChannelBuilder::m_rChType
protected

Definition at line 136 of file TileRawChannelBuilder.h.

◆ m_rChUnit

TileRawChannelUnit::UNIT TileRawChannelBuilder::m_rChUnit
protected

Definition at line 137 of file TileRawChannelBuilder.h.

◆ m_runType

int TileRawChannelBuilder::m_runType
protected

Definition at line 157 of file TileRawChannelBuilder.h.

◆ m_tileHWID

const TileHWID* TileRawChannelBuilder::m_tileHWID
protected

Definition at line 161 of file TileRawChannelBuilder.h.

◆ m_tileID

const TileID* TileRawChannelBuilder::m_tileID
protected

Definition at line 160 of file TileRawChannelBuilder.h.

◆ m_tileIdTransforms

ToolHandle<TileCondIdTransforms> TileRawChannelBuilder::m_tileIdTransforms
protected
Initial value:
{this,
"TileCondIdTransforms", "TileCondIdTransforms",
"Tile tool to tranlate hardware identifier to the drawerIdx, channel, and adc"}

Definition at line 172 of file TileRawChannelBuilder.h.

172 {this,
173 "TileCondIdTransforms", "TileCondIdTransforms",
174 "Tile tool to tranlate hardware identifier to the drawerIdx, channel, and adc"};

◆ m_tileInfo

const TileInfo* TileRawChannelBuilder::m_tileInfo
protected

Definition at line 216 of file TileRawChannelBuilder.h.

◆ m_tileToolEmscale

ToolHandle<TileCondToolEmscale> TileRawChannelBuilder::m_tileToolEmscale
protected
Initial value:
{this,
"TileCondToolEmscale", "TileCondToolEmscale", "Tile EM scale calibration tool"}

Definition at line 166 of file TileRawChannelBuilder.h.

166 {this,
167 "TileCondToolEmscale", "TileCondToolEmscale", "Tile EM scale calibration tool"};

◆ m_tileToolTiming

ToolHandle<TileCondToolTiming> TileRawChannelBuilder::m_tileToolTiming
protected
Initial value:
{this,
"TileCondToolTiming", "TileCondToolTiming", "Tile timing tool"}

Definition at line 169 of file TileRawChannelBuilder.h.

169 {this,
170 "TileCondToolTiming", "TileCondToolTiming", "Tile timing tool"};

◆ m_timeMaxThresh

float TileRawChannelBuilder::m_timeMaxThresh
protected

correct amplitude is time is below time max threshold

Definition at line 154 of file TileRawChannelBuilder.h.

◆ m_timeMinThresh

float TileRawChannelBuilder::m_timeMinThresh
protected

correct amplitude is time is above time min threshold

Definition at line 153 of file TileRawChannelBuilder.h.

◆ m_trigType

int TileRawChannelBuilder::m_trigType
protected

Definition at line 187 of file TileRawChannelBuilder.h.

◆ m_useDSP

bool TileRawChannelBuilder::m_useDSP
protected

Definition at line 149 of file TileRawChannelBuilder.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ MAX_CHANNELS

const int TileRawChannelBuilder::MAX_CHANNELS = 48
staticprotected

Definition at line 206 of file TileRawChannelBuilder.h.

◆ MAX_DMUS

const int TileRawChannelBuilder::MAX_DMUS = 16
staticprotected

Definition at line 207 of file TileRawChannelBuilder.h.


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