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

#include <TileRawChannelBuilderManyAmps.h>

Inheritance diagram for TileRawChannelBuilderManyAmps:
Collaboration diagram for TileRawChannelBuilderManyAmps:

Public Member Functions

 TileRawChannelBuilderManyAmps (const std::string &type, const std::string &name, const IInterface *parent)
 Standard constructor.
 ~TileRawChannelBuilderManyAmps ()
 Destructor.
virtual StatusCode initialize () override
 Initialize.
virtual StatusCode finalize () override
 Finalize.
virtual TileRawChannelrawChannel (const TileDigits *digits, const EventContext &ctx) override
 Builder virtual method to be implemented by subclasses.
virtual StatusCode createContainer (const EventContext &ctx)
 Create container in SG with name given by parameter (m_rawChannelContainerKey).
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

int m_digitFilterMode
int m_digitFilterLevel
int m_digitFilterTest
TileFilterManagerm_tileFilterManagerHi
TileFilterManagerm_tileFilterManagerLo
ToolHandle< TileCondToolNoiseSamplem_tileToolNoiseSample
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default).
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 31 of file TileRawChannelBuilderManyAmps.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

◆ TileRawChannelBuilderManyAmps()

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

Standard constructor.

Definition at line 36 of file TileRawChannelBuilderManyAmps.cxx.

38 : TileRawChannelBuilder(type, name, parent)
41{
42 //declare interfaces
43 declareInterface< TileRawChannelBuilder >( this );
44 declareInterface< TileRawChannelBuilderManyAmps >(this);
45
46 m_rawChannelContainerKey = "TileRawChannelManyAmps";
47
48 //declare properties
49 declareProperty("FilterMode", m_digitFilterMode = 2); // can be 2 or 3
50 declareProperty("FilterLevel", m_digitFilterLevel = 5); // number of parameters for fit (3-9 for mode 2)
51 declareProperty("FilterTester", m_digitFilterTest = 0); // non-zero means call Tester (during initialization phase).
52}
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
SG::WriteHandleKey< TileRawChannelContainer > m_rawChannelContainerKey
TileRawChannelBuilder(const std::string &type, const std::string &name, const IInterface *parent)
Constructor.

◆ ~TileRawChannelBuilderManyAmps()

TileRawChannelBuilderManyAmps::~TileRawChannelBuilderManyAmps ( )

Destructor.

Definition at line 58 of file TileRawChannelBuilderManyAmps.cxx.

58 {
59}

Member Function Documentation

◆ BadPatternName()

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

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

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

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

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

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

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

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}
TileRawChannelUnit::UNIT m_rChUnit
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 ( )
inherited

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

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 TileRawChannelBuilderManyAmps::finalize ( )
overridevirtual

Finalize.

Reimplemented from TileRawChannelBuilder.

Definition at line 137 of file TileRawChannelBuilderManyAmps.cxx.

137 {
138
139 ATH_MSG_DEBUG( "Finalizing" );
142 return StatusCode::SUCCESS;
143}

◆ getOverflowedChannels()

Overflows_t & TileRawChannelBuilder::getOverflowedChannels ( void )
inherited

Definition at line 37 of file TileRawChannelBuilder.cxx.

37 {
38 return m_overflows;
39}

◆ getTileRawChannelContainerID()

std::string TileRawChannelBuilder::getTileRawChannelContainerID ( void )
inherited

Definition at line 41 of file TileRawChannelBuilder.cxx.

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

◆ initialize()

StatusCode TileRawChannelBuilderManyAmps::initialize ( )
overridevirtual

Initialize.

Reimplemented from TileRawChannelBuilder.

Definition at line 65 of file TileRawChannelBuilderManyAmps.cxx.

65 {
66
67 ATH_MSG_INFO( "TileRawChannelBuilderManyAmps::initialize()" );
68
70
71 // init in superclass
73
74
75 // **************************************************************************************
76 ATH_MSG_INFO( " TileRawChannelBuilderManyAmps init:"
77 << " FilterMode=" << m_digitFilterMode
78 << ", FilterLevel=" << m_digitFilterLevel
79 << ", FilterTest=" << m_digitFilterTest );
80
81 if (!(m_digitFilterMode == 2) && !(m_digitFilterMode == 3)) { // filter mode different from 2 or 3
82
83 ATH_MSG_WARNING( "Wrong Filter mode defined: FilterMode=" << m_digitFilterMode );
84 ATH_MSG_WARNING( "Switching to FilterMode=2" );
86 }
87
88 // if( (m_digitFilterMode == 2) || (m_digitFilterMode == 3) ) {
89 /* Get information needed to initialize TileFilterManager. */
90 int Nsamp = m_tileInfo->NdigitSamples();
91 // Set maximum number of parameters to fit, using FilterLevel
92 int nParamMax = m_digitFilterLevel;
93 if (nParamMax > Nsamp + 1) nParamMax = Nsamp;
94 if (nParamMax < 3) nParamMax = 3;
95 int InTsamp = m_tileInfo->ItrigSample();
96 int jBsamp = 0;
97 int jEsamp = Nsamp - 1;
98 int jBcross = 0;
99 int jEcross = Nsamp - 1;
100 int Nshape = m_tileInfo->NdigitSamples();
101 int InTshape = m_tileInfo->ItrigSample(); // need new method to give this !!!
102 std::vector<double> ShapeXHi = m_tileInfo->digitsShapeHi();
103 std::vector<double> ShapeXLo = m_tileInfo->digitsShapeLo();
104
105 bool lVerbose = msgLvl(MSG::VERBOSE);
106 m_tileFilterManagerHi = new TileFilterManager(m_digitFilterMode, m_digitFilterLevel, nParamMax,
107 Nsamp, InTsamp, jBsamp, jEsamp, jBcross, jEcross, Nshape, InTshape, ShapeXHi, lVerbose);
108
109 m_tileFilterManagerLo = new TileFilterManager(m_digitFilterMode, m_digitFilterLevel, nParamMax,
110 Nsamp, InTsamp, jBsamp, jEsamp, jBcross, jEcross, Nshape, InTshape, ShapeXLo, lVerbose);
111
112 // Run TileFilterTester if m_digitFilterTest in non-zero).
113
114 if (m_digitFilterTest > 0) {
115
116 TileFilterTester* tFilterTestHi = new TileFilterTester(m_tileFilterManagerHi
118
119 tFilterTestHi->genEvents(10);
120 delete tFilterTestHi;
121
122 TileFilterTester* tFilterTestLo = new TileFilterTester(m_tileFilterManagerLo
124
125 tFilterTestLo->genEvents(10);
126 delete tFilterTestLo;
127 } // end FilterTest
128
129 ATH_MSG_DEBUG( "TileRawChannelBuilderManyAmps::initialize() completed successfully" );
130
131 return StatusCode::SUCCESS;
132}
#define ATH_MSG_INFO(x)
#define CHECK(...)
Evaluate an expression and check for errors.
void genEvents(int nEvent)
virtual StatusCode initialize()
Initializer.

◆ initLog()

void TileRawChannelBuilder::initLog ( const EventContext & ctx)
inherited

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

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 * TileRawChannelBuilderManyAmps::rawChannel ( const TileDigits * digits,
const EventContext & ctx )
overridevirtual

Builder virtual method to be implemented by subclasses.

Parameters
digitsPointer to TileDigitsContainer

Reimplemented from TileRawChannelBuilder.

Definition at line 145 of file TileRawChannelBuilderManyAmps.cxx.

145 {
146
147 ++m_chCounter;
148
149 const HWIdentifier adcId = tiledigits->adc_HWID();
150 unsigned int drawerIdx(0), channel(0), gain(0);
151 m_tileIdTransforms->getIndices(adcId, drawerIdx, channel, gain);
152
153 bool lVerbose(false);
154 if (msgLvl(MSG::VERBOSE)) {
155 msg(MSG::VERBOSE) << "Running ManyAmps Fit for TileRawChannel with HWID "
156 << m_tileHWID->to_string(adcId) << endmsg;
157 lVerbose = true;
158 }
159
160
161 /* get pedestal value assumed for reconstruction */
162 //double digSigma = m_tileInfo->DigitsResolution(gain);
163 // new way to get channel-dependent sigma:
164 // but we lost difference between sigma used in digitization and
165 // sigma assumed in reconstruction - it's the same sigma now
166 double digSigma = m_tileToolNoiseSample->getHfn(drawerIdx, channel, gain, TileRawChannelUnit::ADCcounts, ctx);
167
168 /* Get vector of time-slice amplitudes. */
169 std::vector<float> digits = tiledigits->samples();
170 double amp_ch = 0; // Fitted amplitude of RC (to be returned from Filtering code).
171 double err_ch; // Error in amp_ch (from Filtering code)
172 double ped_ch; // Fitted pedestal of raw channel (from Filtering).
173 double chisq_ch = 0.; // Chisq resulting from Filtering.
174 double t_ch = 0.; // Fitted time to be supplied by Filtering code
175 // ---------------------------------------------------------------------------------------
176
177 //int icode = 0;
178
179 // Instantiate tResult, which will collect results from Filtering code.
180 TileFilterResult tResult(digits, digSigma);
181 // Call Fitter to extract the in-time pulse height and related info.
182 if (TileID::HIGHGAIN == gain)
183 /*icode =*/m_tileFilterManagerHi->fitDigits(tResult, lVerbose);
184
185 if (TileID::LOWGAIN == gain)
186 /*icode =*/m_tileFilterManagerLo->fitDigits(tResult, lVerbose);
187
188 tResult.getInTime(amp_ch, err_ch, ped_ch, chisq_ch, t_ch);
189
190 ATH_MSG_VERBOSE( " TileRawChannelBuilderManyAmps: return from FilterManager/Fitter."
191 << ", chisq_ch=" << chisq_ch );
192
193 // convert to pCb (if needed)
194 if (m_calibrateEnergy) {
195 amp_ch = m_tileToolEmscale->doCalibCis(drawerIdx, channel, gain, amp_ch);
196 }
197 // we know that time is zero here, put negative chi^2 to indicate that
198 chisq_ch = -fabs(chisq_ch);
199
200 // TileRawChannel *rawCh = new TileRawChannel(adcId,amp_ch,t_ch,chisq_ch);
201 DataPool<TileRawChannel> tileRchPool(m_dataPoollSize);
202 TileRawChannel *rawCh = tileRchPool.nextElementPtr();
203 rawCh->assign (adcId, amp_ch, t_ch, chisq_ch, 0);
204 ATH_MSG_VERBOSE( "Creating RawChannel"
205 << " a=" << amp_ch
206 << " t=" << t_ch
207 << " q=" << chisq_ch );
208
209 if (TileID::HIGHGAIN == gain) {
210 ++m_nChH;
211 m_RChSumH += amp_ch;
212 } else {
213 ++m_nChL;
214 m_RChSumL += amp_ch;
215 }
216
217 return rawCh;
218}
ToolHandle< TileCondToolNoiseSample > m_tileToolNoiseSample
ToolHandle< TileCondIdTransforms > m_tileIdTransforms
ToolHandle< TileCondToolEmscale > m_tileToolEmscale
void assign(const HWIdentifier &id, float amplitude, float time, float quality, float ped=0.0)
MsgStream & msg
Definition testRead.cxx:32

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

Definition at line 28 of file TileRawChannelBuilder.cxx.

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

◆ resetOverflows()

void TileRawChannelBuilder::resetOverflows ( void )
inherited

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

Member Data Documentation

◆ m_ADCmaskValueMinusEps

float TileRawChannelBuilder::m_ADCmaskValueMinusEps
protectedinherited

indicates channels which were masked in background dataset

Definition at line 222 of file TileRawChannelBuilder.h.

◆ m_ADCmaxMinusEps

float TileRawChannelBuilder::m_ADCmaxMinusEps
protectedinherited

Definition at line 221 of file TileRawChannelBuilder.h.

◆ m_ampMinThresh

float TileRawChannelBuilder::m_ampMinThresh
protectedinherited

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
protectedinherited

Definition at line 210 of file TileRawChannelBuilder.h.

◆ m_bsflags

unsigned int TileRawChannelBuilder::m_bsflags
protectedinherited

Definition at line 138 of file TileRawChannelBuilder.h.

◆ m_cabling

const TileCablingService* TileRawChannelBuilder::m_cabling
protectedinherited

TileCabling instance.

Definition at line 182 of file TileRawChannelBuilder.h.

◆ m_cablingSvc

ServiceHandle<TileCablingSvc> TileRawChannelBuilder::m_cablingSvc
protectedinherited
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
protectedinherited

Definition at line 142 of file TileRawChannelBuilder.h.

◆ m_capdaq

double TileRawChannelBuilder::m_capdaq
protectedinherited

Definition at line 193 of file TileRawChannelBuilder.h.

◆ m_chCounter

unsigned int TileRawChannelBuilder::m_chCounter
protectedinherited

Definition at line 196 of file TileRawChannelBuilder.h.

◆ m_cischan

int TileRawChannelBuilder::m_cischan
protectedinherited

Definition at line 192 of file TileRawChannelBuilder.h.

◆ m_correctTime

bool TileRawChannelBuilder::m_correctTime
protectedinherited

Definition at line 145 of file TileRawChannelBuilder.h.

◆ m_dataPoollSize

int TileRawChannelBuilder::m_dataPoollSize
protectedinherited

Definition at line 204 of file TileRawChannelBuilder.h.

◆ m_demoFragIDs

Gaudi::Property<std::vector<int> > TileRawChannelBuilder::m_demoFragIDs
protectedinherited
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_digitFilterLevel

int TileRawChannelBuilderManyAmps::m_digitFilterLevel
private

Definition at line 50 of file TileRawChannelBuilderManyAmps.h.

◆ m_digitFilterMode

int TileRawChannelBuilderManyAmps::m_digitFilterMode
private

Definition at line 49 of file TileRawChannelBuilderManyAmps.h.

◆ m_digitFilterTest

int TileRawChannelBuilderManyAmps::m_digitFilterTest
private

Definition at line 51 of file TileRawChannelBuilderManyAmps.h.

◆ m_DQstatusKey

SG::ReadHandleKey<TileDQstatus> TileRawChannelBuilder::m_DQstatusKey
protectedinherited
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"}
protectedinherited

Definition at line 125 of file TileRawChannelBuilder.h.

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

◆ m_error

int TileRawChannelBuilder::m_error[MAX_CHANNELS]
protectedinherited

Definition at line 208 of file TileRawChannelBuilder.h.

◆ m_evtCounter

unsigned int TileRawChannelBuilder::m_evtCounter
protectedinherited

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
protectedinherited

Definition at line 218 of file TileRawChannelBuilder.h.

◆ m_f_ADCmaxPlus1

float TileRawChannelBuilder::m_f_ADCmaxPlus1
protectedinherited

Definition at line 220 of file TileRawChannelBuilder.h.

◆ m_firstSample

int TileRawChannelBuilder::m_firstSample
protectedinherited

Definition at line 140 of file TileRawChannelBuilder.h.

◆ m_i_ADCmax

int TileRawChannelBuilder::m_i_ADCmax
protectedinherited

Definition at line 217 of file TileRawChannelBuilder.h.

◆ m_i_ADCmaxPlus1

int TileRawChannelBuilder::m_i_ADCmaxPlus1
protectedinherited

Definition at line 219 of file TileRawChannelBuilder.h.

◆ m_idocis

bool TileRawChannelBuilder::m_idocis
protectedinherited

Definition at line 191 of file TileRawChannelBuilder.h.

◆ m_idolas

bool TileRawChannelBuilder::m_idolas
protectedinherited

Definition at line 189 of file TileRawChannelBuilder.h.

◆ m_idoped

bool TileRawChannelBuilder::m_idoped
protectedinherited

Definition at line 190 of file TileRawChannelBuilder.h.

◆ m_idophys

bool TileRawChannelBuilder::m_idophys
protectedinherited

Definition at line 188 of file TileRawChannelBuilder.h.

◆ m_infoName

std::string TileRawChannelBuilder::m_infoName
protectedinherited

Definition at line 215 of file TileRawChannelBuilder.h.

◆ m_lastDrawer

int TileRawChannelBuilder::m_lastDrawer = -1
protectedinherited

Definition at line 209 of file TileRawChannelBuilder.h.

◆ m_nChH

int TileRawChannelBuilder::m_nChH
protectedinherited

Definition at line 199 of file TileRawChannelBuilder.h.

◆ m_nChL

int TileRawChannelBuilder::m_nChL
protectedinherited

Definition at line 198 of file TileRawChannelBuilder.h.

◆ m_noiseFilterTools

ToolHandleArray<ITileRawChannelTool> TileRawChannelBuilder::m_noiseFilterTools
protectedinherited
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
protectedinherited

Definition at line 212 of file TileRawChannelBuilder.h.

◆ m_overflows

Overflows_t TileRawChannelBuilder::m_overflows
protectedinherited

Definition at line 203 of file TileRawChannelBuilder.h.

◆ m_rawChannelCnt

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

Definition at line 133 of file TileRawChannelBuilder.h.

◆ m_rawChannelContainerKey

SG::WriteHandleKey<TileRawChannelContainer> TileRawChannelBuilder::m_rawChannelContainerKey
protectedinherited
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
protectedinherited

Definition at line 201 of file TileRawChannelBuilder.h.

◆ m_RChSumL

double TileRawChannelBuilder::m_RChSumL
protectedinherited

Definition at line 200 of file TileRawChannelBuilder.h.

◆ m_rChType

TileFragHash::TYPE TileRawChannelBuilder::m_rChType
protectedinherited

Definition at line 136 of file TileRawChannelBuilder.h.

◆ m_rChUnit

TileRawChannelUnit::UNIT TileRawChannelBuilder::m_rChUnit
protectedinherited

Definition at line 137 of file TileRawChannelBuilder.h.

◆ m_runType

int TileRawChannelBuilder::m_runType
protectedinherited

Definition at line 157 of file TileRawChannelBuilder.h.

◆ m_tileFilterManagerHi

TileFilterManager* TileRawChannelBuilderManyAmps::m_tileFilterManagerHi
private

Definition at line 55 of file TileRawChannelBuilderManyAmps.h.

◆ m_tileFilterManagerLo

TileFilterManager* TileRawChannelBuilderManyAmps::m_tileFilterManagerLo
private

Definition at line 56 of file TileRawChannelBuilderManyAmps.h.

◆ m_tileHWID

const TileHWID* TileRawChannelBuilder::m_tileHWID
protectedinherited

Definition at line 161 of file TileRawChannelBuilder.h.

◆ m_tileID

const TileID* TileRawChannelBuilder::m_tileID
protectedinherited

Definition at line 160 of file TileRawChannelBuilder.h.

◆ m_tileIdTransforms

ToolHandle<TileCondIdTransforms> TileRawChannelBuilder::m_tileIdTransforms
protectedinherited
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
protectedinherited

Definition at line 216 of file TileRawChannelBuilder.h.

◆ m_tileToolEmscale

ToolHandle<TileCondToolEmscale> TileRawChannelBuilder::m_tileToolEmscale
protectedinherited
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_tileToolNoiseSample

ToolHandle<TileCondToolNoiseSample> TileRawChannelBuilderManyAmps::m_tileToolNoiseSample
private
Initial value:
{this,
"TileCondToolNoiseSample", "TileCondToolNoiseSample", "Tile sample noise tool"}

Definition at line 58 of file TileRawChannelBuilderManyAmps.h.

58 {this,
59 "TileCondToolNoiseSample", "TileCondToolNoiseSample", "Tile sample noise tool"};

◆ m_tileToolTiming

ToolHandle<TileCondToolTiming> TileRawChannelBuilder::m_tileToolTiming
protectedinherited
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
protectedinherited

correct amplitude is time is below time max threshold

Definition at line 154 of file TileRawChannelBuilder.h.

◆ m_timeMinThresh

float TileRawChannelBuilder::m_timeMinThresh
protectedinherited

correct amplitude is time is above time min threshold

Definition at line 153 of file TileRawChannelBuilder.h.

◆ m_trigType

int TileRawChannelBuilder::m_trigType
protectedinherited

Definition at line 187 of file TileRawChannelBuilder.h.

◆ m_useDSP

bool TileRawChannelBuilder::m_useDSP
protectedinherited

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
staticprotectedinherited

Definition at line 206 of file TileRawChannelBuilder.h.

◆ MAX_DMUS

const int TileRawChannelBuilder::MAX_DMUS = 16
staticprotectedinherited

Definition at line 207 of file TileRawChannelBuilder.h.


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