InDet::TRT_StrawStatus Class Reference

This algorithm finds dead or hot TRT straws / chips / boards based on
occupancy and hits on track. More...

#include <TRT_StrawStatus.h>

Inheritance diagram for InDet::TRT_StrawStatus:

Public Member Functions
	TRT_StrawStatus (const std::string &name, ISvcLocator *pSvcLocator)
	Standard Athena-Algorithm Constructor.
	~TRT_StrawStatus ()
	Default Destructor.
StatusCode	initialize ()
	standard Athena-Algorithm method
StatusCode	execute ()
	standard Athena-Algorithm method
StatusCode	finalize ()
	standard Athena-Algorithm method
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef std::array< std::array< std::array< std::array< int, 6 >, 5482 >, 32 >, 2 >	ACCHITS_t
	accumulate hits, last index: 0 - all hits, 1 - hits on track, 2 - all HT (TR) hits, 3 - HT (TR) hits on track
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	clear ()
StatusCode	reportResults ()
void	printDetailedInformation ()
void	myStrawIndex (Identifier id, int *index)
	function that returns straw index (in range 0-5481; 0-1641 for barrel, the rest for endcap) same convention as for TRT_monitoring (copied from there)
int	barrelStrawNumber (int strawNumber, int strawlayerNumber, int LayerNumber)
int	endcapStrawNumber (int strawNumber, int strawLayerNumber, int LayerNumber)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ServiceHandle< ITRT_HWMappingSvc >	m_mapSvc {this,"HWMapSvc","TRT_HWMappingSvc","" }
ServiceHandle< ITRT_DCS_ConditionsSvc >	m_DCSSvc {this,"InDetTRT_DCS_ConditionsSvc","TRT_DCS_ConditionsSvc","" }
ServiceHandle< ITRT_StrawNeighbourSvc >	m_TRTStrawNeighbourSvc {this,"TRT_StrawNeighbourSvc","TRT_StrawNeighbourSvc","retrieve barrel and end-cap straw number later on, as well as DTMROC" }
ToolHandle< ITRT_StrawStatusSummaryTool >	m_TRTStrawStatusSummaryTool {this, "TRT_StrawStatusSummaryTool", "ITRT_StrawStatusSummaryTool", ""}
ToolHandle< Trk::ITrackHoleSearchTool >	m_trt_hole_finder {this, "trt_hole_finder", "Trk::ITrackHoleSearchTool", ""}
PublicToolHandle< Trk::IUpdator >	m_updator {this, "KalmanUpdator", "Trk::KalmanUpdator/TrkKalmanUpdator",""}
Gaudi::Property< double >	m_locR_cut {this, "locR_cut", 1.4, ""}
Gaudi::Property< int >	m_skipBusyEvents {this, "skipBusyEvents", 0, ""}
Gaudi::Property< int >	m_printDetailedInformation {this, "printDetailedInformation", 0, ""}
Gaudi::Property< std::string >	m_fileName {this, "outputFileName", "TRT_StrawStatusOutput", ""}
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfoKey {this,"EventInfoKey","EventInfo","RHK to retrieve xAOD::EventInfo"}
SG::ReadHandleKey< TRT_RDO_Container >	m_rdoContainerKey {this,"RDO_ContainerKey","TRT_RDOs","RHK to retrieve TRT RDO's"}
SG::ReadHandleKey< DataVector< Trk::Track > >	m_tracksName {this,"tracksCollectionKey","CombinedInDetTracks","RHK to retrieve CombinedInDetTracks"}
SG::ReadHandleKey< xAOD::VertexContainer >	m_vxContainerKey {this,"VxContainerKey","PrimaryVertices","RHK to retrieve VX Primary candidates"}
int	m_nEvents
	returns index of hardware units: board, chip, pad private fix for now, will call TRTStrawNeighbourSvc when available number boards 0-9 barrel, 0-19 endcap (first 12 on A wheels, ordering from smaller to larger |z|) number chips 0-103 barrel, 0-239 endcap number pads: chips x 2
int	m_runNumber
std::unique_ptr< ACCHITS_t >	m_accumulateHits
const TRT_ID *	m_TRTHelper
std::atomic< int >	m_printStatusCount {0}
DataObjIDColl	m_extendedExtraObjects
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

This algorithm finds dead or hot TRT straws / chips / boards based on
occupancy and hits on track.

To be used in calibration stream (to mask off dead regions before processing)

Author

Sasa Fratina sasa..nosp@m.frat.nosp@m.ina@c.nosp@m.ern..nosp@m.ch

Definition at line 56 of file TRT_StrawStatus.h.

Member Typedef Documentation

ACCHITS_t

typedef std::array<std::array<std::array<std::array<int,6>,5482>,32>,2> InDet::TRT_StrawStatus::ACCHITS_t

private

accumulate hits, last index: 0 - all hits, 1 - hits on track, 2 - all HT (TR) hits, 3 - HT (TR) hits on track

Definition at line 113 of file TRT_StrawStatus.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

TRT_StrawStatus()

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

Standard Athena-Algorithm Constructor.

Definition at line 43 of file TRT_StrawStatus.cxx.

:
AthAlgorithm(name,pSvcLocator),
m_nEvents(0), m_runNumber(0),
m_TRTHelper(nullptr)
{}

~TRT_StrawStatus()

InDet::TRT_StrawStatus::~TRT_StrawStatus

(

)

Default Destructor.

Definition at line 52 of file TRT_StrawStatus.cxx.

{}

Member Function Documentation

barrelStrawNumber()

int InDet::TRT_StrawStatus::barrelStrawNumber ( int strawNumber, int strawlayerNumber, int LayerNumber )

private

clear()

void InDet::TRT_StrawStatus::clear ( )

private

Definition at line 284 of file TRT_StrawStatus.cxx.

                                 {
    m_nEvents = 0;
    *m_accumulateHits = {};
    return;
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

endcapStrawNumber()

int InDet::TRT_StrawStatus::endcapStrawNumber ( int strawNumber, int strawLayerNumber, int LayerNumber )

private

evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode InDet::TRT_StrawStatus::execute

(

)

standard Athena-Algorithm method

Definition at line 93 of file TRT_StrawStatus.cxx.

                                        {
 
    SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInfoKey);
    StatusCode sc = StatusCode::SUCCESS;
    if (not eventInfo.isValid()) {
      ATH_MSG_ERROR( "Unable to retrieve Event Info " );
      return StatusCode::FAILURE;
    }
    int runNumber = (int) eventInfo->runNumber();
    if (runNumber != m_runNumber) {
      if (m_nEvents) { ATH_CHECK( reportResults() ); clear(); }
      m_runNumber = runNumber;
    }
    int lumiBlock0 =eventInfo->lumiBlock();
    SG::ReadHandle<TRT_RDO_Container> rdoContainer(m_rdoContainerKey);
 
    if (not rdoContainer.isValid()) {
      ATH_MSG_ERROR( "no TRT_RDO container available " );
      return StatusCode::FAILURE;
    }
    SG::ReadHandle<DataVector<Trk::Track>> trkCollection(m_tracksName);
    if (not trkCollection.isValid()) {
      ATH_MSG_ERROR( "Could not find Tracks Collection: " << m_tracksName );
      return StatusCode::FAILURE;
    }
 
    //================ Event selection
 
    SG::ReadHandle<xAOD::VertexContainer> vertices(m_vxContainerKey);
    if (not vertices.isValid()) {
       ATH_MSG_DEBUG ("Couldn't retrieve VertexContainer with key: PrimaryVertices");
       return StatusCode::SUCCESS;   // just skip to next event in case of no vertexcontainer
    }
 
    int countVertices(0);
    for (const xAOD::Vertex* vx : *(vertices.cptr()) ) {
        if (vx->vertexType() == xAOD::VxType::PriVtx) {
           if ( vx-> nTrackParticles() >= 3) countVertices++;
        }
    }
    if (countVertices < 1) {
       ATH_MSG_INFO( "no vertices with at least 3 tracks found" );
       return StatusCode::SUCCESS;
    }
 
    if (m_skipBusyEvents) { // cosmic running
        int countRDOhitsInEvent(0);
        for (TRT_RDO_Container::const_iterator rdoIt = rdoContainer->begin(); rdoIt != rdoContainer->end(); ++rdoIt) {
            const InDetRawDataCollection<TRT_RDORawData>* TRTCollection(*rdoIt);
            if (not TRTCollection) continue;
            for (DataVector<TRT_RDORawData>::const_iterator trtIt = TRTCollection->begin(); trtIt != TRTCollection->end(); ++trtIt) {
                countRDOhitsInEvent++;
            }
        }
        if (countRDOhitsInEvent>100000) {
            ATH_MSG_INFO( "N RDO hits in event greater than 100000: " << countRDOhitsInEvent << ", exiting" );
            return sc;
        }
 
        if (trkCollection->size() > 10) {
            ATH_MSG_INFO( "N tracks greater than 10: " << trkCollection->size() << ", exiting" );
            return sc;
        }
    }
 
    //================ End event selection
 
    //================ Loop over all tracks, accumulate hits on track, also find holes on track
 
    std::vector<Identifier> holeIdentifiers;
    std::vector<Identifier> holeIdentifiersWithHits; // holes on straws that have hits, it is just that the hit was not associalted to a track
    for ( DataVector<Trk::Track>::const_iterator trackIt = trkCollection->begin(); trackIt != trkCollection->end(); ++trackIt) {
        const Trk::Track *track = *trackIt;
        //=== select track
        const Trk::Perigee* perigee = (*trackIt)->perigeeParameters();
        if ( not perigee  ) { ATH_MSG_ERROR( "Trk::Perigee missing" ); continue; }
        if ( std::abs(perigee->pT())/CLHEP::GeV < 1. ) continue; // 1 GeV pT cut
 
        const DataVector<const Trk::TrackStateOnSurface>* trackStates = (**trackIt).trackStateOnSurfaces();
        if ( not trackStates  ) { ATH_MSG_ERROR( "Trk::TrackStateOnSurface empty" ); continue; }
 
        int n_pixel_hits(0), n_sct_hits(0), n_trt_hits(0);  // count hits, require minimal number of all hits
        for ( DataVector<const Trk::TrackStateOnSurface>::const_iterator trackStatesIt = trackStates->begin(); trackStatesIt != trackStates->end(); ++trackStatesIt ) {
            if ( *trackStatesIt == nullptr ) { ATH_MSG_ERROR( "*trackStatesIt == 0" ); continue; }
 
            if ( !((*trackStatesIt)->type(Trk::TrackStateOnSurface::Measurement)) ) continue; // this skips outliers
 
            if ( dynamic_cast<const InDet::TRT_DriftCircleOnTrack*> ( (*trackStatesIt)->measurementOnTrack() )  ) n_trt_hits++;
            else if ( dynamic_cast<const InDet::SCT_ClusterOnTrack*>     ( (*trackStatesIt)->measurementOnTrack() )  ) n_sct_hits++;
            else if( dynamic_cast<const InDet::PixelClusterOnTrack*>    ( (*trackStatesIt)->measurementOnTrack() )  ) n_pixel_hits++;
        }
        if (n_pixel_hits<2 || n_sct_hits < 6 || n_trt_hits<15) continue; // end count hits
 
        //=== end select track
 
        //=== loop over all hits on track, accumulate them
 
        for ( DataVector<const Trk::TrackStateOnSurface>::const_iterator trackStatesIt = trackStates->begin(); trackStatesIt != trackStates->end(); ++trackStatesIt ) {
 
            if ( *trackStatesIt == nullptr ) { ATH_MSG_ERROR( "*trackStatesIt == 0" ); continue; }
 
            if ( !((*trackStatesIt)->type(Trk::TrackStateOnSurface::Measurement)) ) continue; // this skips outliers
 
            const InDet::TRT_DriftCircleOnTrack *driftCircleOnTrack = dynamic_cast<const InDet::TRT_DriftCircleOnTrack *>( (*trackStatesIt)->measurementOnTrack() );
            if ( not driftCircleOnTrack ) continue; // not TRT measurement - this way, keep both hits and outliers
 
            const Trk::TrackStateOnSurface& hit = **trackStatesIt;
 
            const Trk::TrackParameters* unbiased_track_parameters = m_updator->removeFromState( *(hit.trackParameters()),
                                                                                          hit.measurementOnTrack()->localParameters(),
                                                                                          hit.measurementOnTrack()->localCovariance()).release();
 
            double unbiased_locR = unbiased_track_parameters->parameters()[Trk::locR];
            if ( std::abs(unbiased_locR) > m_locR_cut ) continue; // same cut as the default hole search cut
 
            const InDet::TRT_DriftCircle *driftCircle = driftCircleOnTrack->prepRawData();
            if ( driftCircle == nullptr ) { ATH_MSG_ERROR( "driftCircle == 0" ); continue; }
 
            Identifier id = driftCircle->identify();
            int index[6]; myStrawIndex(id, index); // side, layer, phi, straw_layer, straw_within_layer, straw_index
            (*m_accumulateHits)[(index[0]>0)?0:1][index[2]][index[5]][1]++; // accumulate hits on track
            if (driftCircle->highLevel()) (*m_accumulateHits)[(index[0]>0)?0:1][index[2]][index[5]][3]++; // accumulate hits on track
 
        } // end trackStatesIt loop
 
        // add holeIdentifiers - fill vector
 
        std::unique_ptr<const Trk::TrackStates> holes (m_trt_hole_finder->getHolesOnTrack( *track ));
        if ( holes==nullptr ) continue; // no holes found
        for (const Trk::TrackStateOnSurface* trackStates : *holes) {
 
            if ( !trackStates->type(   Trk::TrackStateOnSurface::Hole  )  ) { ATH_MSG_ERROR( "m_trt_hole_finder returned something that is not a hole" ); continue; }
 
            const Trk::TrackParameters* track_parameters = trackStates->trackParameters();
            if (!track_parameters) { ATH_MSG_WARNING( "m_trt_hole_finder track_parameters missing" ); continue; }
 
            Identifier id = track_parameters->associatedSurface().associatedDetectorElementIdentifier();
            if ( !(m_TRTHelper->is_trt(id)) ) { ATH_MSG_ERROR( "m_trt_hole_finder returned something that is not a TRT hole" ); continue; }
 
            // add se same 1.4 mm locR selection, in case it is not on by default
            if ( std::abs( track_parameters->parameters()[Trk::locR] ) > m_locR_cut ) continue;
 
            holeIdentifiers.push_back( id );
        } // end add holeIdentifiers
 
    } // end trackIt loop
 
    //================ End loop over all tracks
 
    //================ Loop over all hits - it includes hits from dead straws that are masked off in drift circle creation
 
    for (TRT_RDO_Container::const_iterator rdoIt = rdoContainer->begin(); rdoIt != rdoContainer->end(); ++rdoIt) {
        const InDetRawDataCollection<TRT_RDORawData>* TRTCollection(*rdoIt);
        if (TRTCollection==nullptr) continue;
        for (DataVector<TRT_RDORawData>::const_iterator trtIt = TRTCollection->begin(); trtIt != TRTCollection->end(); ++trtIt) {
            Identifier id = (*trtIt)->identify();
            int index[6]; myStrawIndex(id, index); // side, layer, phi, straw_layer, straw_within_layer, straw_index
            (*m_accumulateHits)[(index[0]>0)?0:1][index[2]][index[5]][0]++; // accumulate all hits
            if ((*trtIt)->highLevel()) (*m_accumulateHits)[(index[0]>0)?0:1][index[2]][index[5]][2]++; // accumulate TR hits
 
            if (std::find(holeIdentifiers.begin(), holeIdentifiers.end(), id) != holeIdentifiers.end())  // a hole was found on the same straw, but hits is there
                holeIdentifiersWithHits.push_back( id );
        }
    }
 
    //================ End loop over all hits
 
    //================ End loop over all holes, each time also save whether the straw with a hole had a hit
 
    for (unsigned int i=0; i<holeIdentifiers.size(); i++) {
 
        Identifier id = holeIdentifiers[i];
 
        int index[6]; myStrawIndex(id, index); // side, layer, phi, straw_layer, straw_within_layer, straw_index
 
        (*m_accumulateHits)[(index[0]>0)?0:1][index[2]][index[5]][4]++;
 
        if (std::find(holeIdentifiersWithHits.begin(), holeIdentifiersWithHits.end(), id) != holeIdentifiersWithHits.end())
          (*m_accumulateHits)[(index[0]>0)?0:1][index[2]][index[5]][5]++;
    }
 
    //================ End loop over all hits
 
    m_nEvents++;
    last_lumiBlock0 = lumiBlock0;
    if (m_nEvents%1000==0 && msgLvl(MSG::DEBUG)) ATH_CHECK( reportResults() );
    return sc;
}

extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return Algorithm::extraOutputDeps();
}

finalize()

StatusCode InDet::TRT_StrawStatus::finalize

(

)

standard Athena-Algorithm method

Definition at line 83 of file TRT_StrawStatus.cxx.

                                         {
    ATH_CHECK( reportResults() );
    if (m_printDetailedInformation) printDetailedInformation();
    // Code entered here will be executed once at the end of the program run.
    m_accumulateHits.reset();
    return StatusCode::SUCCESS;
}

initialize()

StatusCode InDet::TRT_StrawStatus::initialize

(

)

standard Athena-Algorithm method

Definition at line 58 of file TRT_StrawStatus.cxx.

{
 
  m_accumulateHits = std::make_unique<ACCHITS_t>();
  assert( (*m_accumulateHits)[0][0].size() == nAllStraws );
  clear();
 
  // Code entered here will be executed once at program start.
  // Initialize ReadHandleKey
  ATH_CHECK(m_eventInfoKey.initialize());
  ATH_CHECK(m_rdoContainerKey.initialize());
  ATH_CHECK(m_tracksName.initialize());
  ATH_CHECK(m_vxContainerKey.initialize());
  ATH_CHECK( detStore()->retrieve(m_TRTHelper, "TRT_ID"));
  ATH_CHECK( m_TRTStrawNeighbourSvc.retrieve()) ;
  ATH_CHECK( m_trt_hole_finder.retrieve() );
  ATH_CHECK(m_updator.retrieve());
  ATH_CHECK( m_mapSvc.retrieve());
  ATH_CHECK( m_DCSSvc.retrieve());
  ATH_MSG_DEBUG( "initialize() successful in " << name() << ", retrieved: ..., locR_cut = " << m_locR_cut << " mm." );
  return StatusCode::SUCCESS;
}

inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

myStrawIndex()

void InDet::TRT_StrawStatus::myStrawIndex ( Identifier id, int * index )

private

function that returns straw index (in range 0-5481; 0-1641 for barrel, the rest for endcap) same convention as for TRT_monitoring (copied from there)

Definition at line 360 of file TRT_StrawStatus.cxx.

                                                            {
    int side = m_TRTHelper->barrel_ec(id);
    int layerNumber = m_TRTHelper->layer_or_wheel(id);
    int strawLayerNumber = m_TRTHelper->straw_layer(id);
    int strawNumber = m_TRTHelper->straw(id);
    int straw(0);
 
    const int numberOfStraws[74] = { 0, 15, 31, 47, 63, 79, 96, 113, 130, 147, 164, 182, 200, 218, 236, 254, 273, 292, 311, 329, // layer 0, 329 straws, strawlayers 0-18
        348, 368, 388, 408, 428, 448, 469, 490, 511, 532, 553, 575, 597, 619, 641, 663, 686, 709, 732, 755, 778, 802, 826, 849, // layer 1, 520 straws, strawLayers 0-23
        872, 896, 920, 944, 968, 993, 1018, 1043, 1068, 1093, 1119, 1145, 1171, 1197, 1223, 1250, 1277, 1304, 1331, 1358, 1386, 1414, 1442, 1470, 1498, 1527, 1556, 1585, 1614, 1642 }; // layer 2
 
    if (abs(side)==1) { // barrel unique straw number
        if (layerNumber==1) strawLayerNumber+= 19;
        else if (layerNumber==2) strawLayerNumber+= 43;
        straw = ( numberOfStraws[strawLayerNumber+1] - strawNumber -1 );
    } else { // end-cap unique straw number
        int board = layerNumber;
        if (board<6) { board *= 2; if (strawLayerNumber>7) board++; }
        else { board += 6; }
        straw = board * 192 + strawNumber * 8 + strawLayerNumber % 8 ;
        straw += 1642;
    }
    index[0] = side;
    index[1] = layerNumber;
    index[2] = m_TRTHelper->phi_module(id);
    index[3] = strawLayerNumber;
    index[4] = strawNumber;
    index[5] = straw;
    return;
}

outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

printDetailedInformation()

void InDet::TRT_StrawStatus::printDetailedInformation ( )

private

Definition at line 312 of file TRT_StrawStatus.cxx.

                                              {
  ATH_MSG_INFO("InDet::TRT_StrawStatus::printDetailedInformation()");
  const std::string fileName =
    std::format("{}.{}{}_printDetailedInformation.txt",
                m_fileName.value(),
                std::string(7 - std::min<std::size_t>(7, std::to_string(m_runNumber).size()), '0'),
                m_runNumber);
  std::ofstream out(fileName, std::ios::out | std::ios::trunc);
  if (!out) {
    ATH_MSG_ERROR(std::format("Failed to open '{}' for writing", fileName));
    return;
  }
 
  // Iterate over straw layers
  for (auto it = m_TRTHelper->straw_layer_begin(); it != m_TRTHelper->straw_layer_end(); ++it) {
    const Identifier layerId = *it;
    const int maxStraw = m_TRTHelper->straw_max(layerId);
 
    for (int i = 0; i <= maxStraw; ++i) {
      const Identifier id = m_TRTHelper->straw_id(layerId, i);
      std::array<int, 6> index{};
      myStrawIndex(id, index.data());
      int chip = 0;
      int HVpad = 0;
      m_TRTStrawNeighbourSvc->getChip(id, chip);
      m_TRTStrawNeighbourSvc->getPad(id, HVpad);
      if (!m_printStatusCount) {
        ATH_MSG_INFO("if the code crashes on the next line, there is a problem with "
                     "m_TRTStrawStatusSummarySvc not being loaded ");
        ATH_MSG_INFO("in that case, running with reco turned on normally solves the problem, "
                     "know of no better solution at the moment");
        ATH_MSG_INFO("if you do not need the detailed print information, you can also just set "
                     "printDetailedInformation to 0 to avoid this crash");
        ++m_printStatusCount;
      }
      const auto& ctx = Gaudi::Hive::currentContext();
      const int status          = m_TRTStrawStatusSummaryTool->get_status(id, ctx);
      const int statusTemporary = m_TRTStrawStatusSummaryTool->getStatus(id, ctx);
      const int statusPermanent = m_TRTStrawStatusSummaryTool->getStatusPermanent(id, ctx);
      // Write: six indices, then chip/HVpad/status trio
      for (int j = 0; j < 6; ++j) out << index[j] << ' ';
      out << chip << ' ' << HVpad << ' ' << status << ' '
          << statusTemporary << ' ' << statusPermanent << '\n';
    }
  }
}

renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

reportResults()

StatusCode InDet::TRT_StrawStatus::reportResults ( )

private

Definition at line 290 of file TRT_StrawStatus.cxx.

                                               {
    ATH_MSG_INFO( "InDet::TRT_StrawStatus::reportResults() for " << m_nEvents << " events." );
    char fileName[300];
    snprintf(fileName, 299,"%s.%07d_newFormat.txt", m_fileName.value().c_str(), m_runNumber);
    FILE *f = fopen(fileName, "w");
    if (!f) {
      ATH_MSG_ERROR( "InDet::TRT_StrawStatus::reportResults: Cannot open " << fileName << " for write" );
      return StatusCode::FAILURE;
    }
    fprintf(f, "%d %d %d %d %d %d %d %d %d \n", 0, 0, 0, 0, 0, 0, 0, 0, m_nEvents);
    for (size_t i=0; i<2; i++) for (size_t j=0; j<32; j++) for (size_t k=0; k<nAllStraws; k++) {
        int side = (i>0)?-1:1;
        if (k>=1642) side *= 2;
        fprintf(f, "%d %zu %zu", side, j, k);
        for (int m=0; m<6; m++) fprintf(f, " %d", (*m_accumulateHits)[i][j][k][m]);
        fprintf(f, "\n");
    }
    fclose(f);
    return StatusCode::SUCCESS;
}

sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

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

Definition at line 66 of file AthAlgorithm.cxx.

                                       {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )

overridevirtualinherited

Handle START transition.

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

updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_accumulateHits

std::unique_ptr<ACCHITS_t> InDet::TRT_StrawStatus::m_accumulateHits

private

Definition at line 114 of file TRT_StrawStatus.h.

m_DCSSvc

ServiceHandle<ITRT_DCS_ConditionsSvc> InDet::TRT_StrawStatus::m_DCSSvc {this,"InDetTRT_DCS_ConditionsSvc","TRT_DCS_ConditionsSvc","" }

private

Definition at line 76 of file TRT_StrawStatus.h.

{this,"InDetTRT_DCS_ConditionsSvc","TRT_DCS_ConditionsSvc","" };

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_eventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> InDet::TRT_StrawStatus::m_eventInfoKey {this,"EventInfoKey","EventInfo","RHK to retrieve xAOD::EventInfo"}

private

Definition at line 88 of file TRT_StrawStatus.h.

{this,"EventInfoKey","EventInfo","RHK to retrieve xAOD::EventInfo"};

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects

privateinherited

Definition at line 79 of file AthAlgorithm.h.

m_fileName

Gaudi::Property<std::string> InDet::TRT_StrawStatus::m_fileName {this, "outputFileName", "TRT_StrawStatusOutput", ""}

private

Definition at line 85 of file TRT_StrawStatus.h.

{this, "outputFileName", "TRT_StrawStatusOutput", ""};

m_locR_cut

Gaudi::Property<double> InDet::TRT_StrawStatus::m_locR_cut {this, "locR_cut", 1.4, ""}

private

Definition at line 82 of file TRT_StrawStatus.h.

{this, "locR_cut", 1.4, ""};

m_mapSvc

ServiceHandle<ITRT_HWMappingSvc> InDet::TRT_StrawStatus::m_mapSvc {this,"HWMapSvc","TRT_HWMappingSvc","" }

private

Definition at line 75 of file TRT_StrawStatus.h.

{this,"HWMapSvc","TRT_HWMappingSvc","" };

m_nEvents

int InDet::TRT_StrawStatus::m_nEvents

private

returns index of hardware units: board, chip, pad private fix for now, will call TRTStrawNeighbourSvc when available number boards 0-9 barrel, 0-19 endcap (first 12 on A wheels, ordering from smaller to larger |z|) number chips 0-103 barrel, 0-239 endcap number pads: chips x 2

Definition at line 109 of file TRT_StrawStatus.h.

m_printDetailedInformation

Gaudi::Property<int> InDet::TRT_StrawStatus::m_printDetailedInformation {this, "printDetailedInformation", 0, ""}

private

Definition at line 84 of file TRT_StrawStatus.h.

{this, "printDetailedInformation", 0, ""};

m_printStatusCount

std::atomic<int> InDet::TRT_StrawStatus::m_printStatusCount {0}

mutableprivate

Definition at line 117 of file TRT_StrawStatus.h.

{0};

m_rdoContainerKey

SG::ReadHandleKey<TRT_RDO_Container> InDet::TRT_StrawStatus::m_rdoContainerKey {this,"RDO_ContainerKey","TRT_RDOs","RHK to retrieve TRT RDO's"}

private

Definition at line 89 of file TRT_StrawStatus.h.

{this,"RDO_ContainerKey","TRT_RDOs","RHK to retrieve TRT RDO's"};

m_runNumber

int InDet::TRT_StrawStatus::m_runNumber

private

Definition at line 110 of file TRT_StrawStatus.h.

m_skipBusyEvents

Gaudi::Property<int> InDet::TRT_StrawStatus::m_skipBusyEvents {this, "skipBusyEvents", 0, ""}

private

Definition at line 83 of file TRT_StrawStatus.h.

{this, "skipBusyEvents", 0, ""};

m_tracksName

SG::ReadHandleKey<DataVector<Trk::Track> > InDet::TRT_StrawStatus::m_tracksName {this,"tracksCollectionKey","CombinedInDetTracks","RHK to retrieve CombinedInDetTracks"}

private

Definition at line 90 of file TRT_StrawStatus.h.

{this,"tracksCollectionKey","CombinedInDetTracks","RHK to retrieve CombinedInDetTracks"};

m_trt_hole_finder

ToolHandle<Trk::ITrackHoleSearchTool> InDet::TRT_StrawStatus::m_trt_hole_finder {this, "trt_hole_finder", "Trk::ITrackHoleSearchTool", ""}

private

Definition at line 79 of file TRT_StrawStatus.h.

{this, "trt_hole_finder", "Trk::ITrackHoleSearchTool", ""};

m_TRTHelper

const TRT_ID* InDet::TRT_StrawStatus::m_TRTHelper

private

Definition at line 116 of file TRT_StrawStatus.h.

m_TRTStrawNeighbourSvc

ServiceHandle<ITRT_StrawNeighbourSvc> InDet::TRT_StrawStatus::m_TRTStrawNeighbourSvc {this,"TRT_StrawNeighbourSvc","TRT_StrawNeighbourSvc","retrieve barrel and end-cap straw number later on, as well as DTMROC" }

private

Definition at line 77 of file TRT_StrawStatus.h.

{this,"TRT_StrawNeighbourSvc","TRT_StrawNeighbourSvc","retrieve barrel and end-cap straw number later on, as well as DTMROC" };

m_TRTStrawStatusSummaryTool

ToolHandle<ITRT_StrawStatusSummaryTool> InDet::TRT_StrawStatus::m_TRTStrawStatusSummaryTool {this, "TRT_StrawStatusSummaryTool", "ITRT_StrawStatusSummaryTool", ""}

private

Definition at line 78 of file TRT_StrawStatus.h.

{this, "TRT_StrawStatusSummaryTool", "ITRT_StrawStatusSummaryTool", ""};

m_updator

PublicToolHandle<Trk::IUpdator> InDet::TRT_StrawStatus::m_updator {this, "KalmanUpdator", "Trk::KalmanUpdator/TrkKalmanUpdator",""}

private

Definition at line 80 of file TRT_StrawStatus.h.

{this, "KalmanUpdator", "Trk::KalmanUpdator/TrkKalmanUpdator",""};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_vxContainerKey

SG::ReadHandleKey<xAOD::VertexContainer> InDet::TRT_StrawStatus::m_vxContainerKey {this,"VxContainerKey","PrimaryVertices","RHK to retrieve VX Primary candidates"}

private

Definition at line 91 of file TRT_StrawStatus.h.

{this,"VxContainerKey","PrimaryVertices","RHK to retrieve VX Primary candidates"};

---

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

• TRT_StrawStatus.h
• TRT_StrawStatus.cxx