MuonR4::SpacePointMakerAlg Class Reference

#include <SpacePointMakerAlg.h>

Inheritance diagram for MuonR4::SpacePointMakerAlg:

Collaboration diagram for MuonR4::SpacePointMakerAlg:

Classes
struct	SpacePointsPerChamber
	: Helper struct to collect the space point per muon chamber, which are later sorted into the space point buckets. More...
class	SpacePointStatistics
	Helper class to keep track of how many eta+phi, eta and phi only space points are built in various detector regions. More...

Public Member Functions
	~SpacePointMakerAlg ()=default
StatusCode	execute (const EventContext &ctx) const override
StatusCode	initialize () override
StatusCode	finalize () override
template<>
bool	passOccupancy2D (const std::vector< const xAOD::TgcStrip * > &etaHits, const std::vector< const xAOD::TgcStrip * > &phiHits) const
template<>
bool	passOccupancy2D (const std::vector< const xAOD::RpcMeasurement * > &etaHits, const std::vector< const xAOD::RpcMeasurement * > &phiHits) const
template<>
bool	passOccupancy2D (const std::vector< const xAOD::sTgcMeasurement * > &etaHits, const std::vector< const xAOD::sTgcMeasurement * > &phiHits) const
virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...
virtual bool	isClonable () const override
	Specify if the algorithm is clonable. More...
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm. More...
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...
virtual bool	filterPassed (const EventContext &ctx) const
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
virtual StatusCode	sysStart () override
	Handle START transition. More...
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...
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. More...

Private Types
using	PreSortedSpacePointMap = std::unordered_map< const MuonGMR4::SpectrometerSector *, SpacePointsPerChamber >
	Container abrivation of the presorted space point container per MuonChambers. More...
using	SpacePointBucketVec = std::vector< SpacePointBucket >
	Abrivation of a MuonSapcePoint bucket vector. More...
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
template<class MeasType >
Amg::Transform3D	toChamberTransform (const ActsGeometryContext &gctx, const Amg::Transform3D &sectorTrans, const MeasType *meas) const
	Returns the transform from the attached Muon chamber frame to the sector frame. More...
template<class MeasType >
Amg::Vector3D	positionInChamber (const MeasType *meas, const Amg::Transform3D &toChamberTrans) const
	Returns the position of the uncalibrated muon measurement in the sector frame. More...
template<class MeasType >
Amg::Vector3D	channelDirInChamber (const MeasType *meas, const Amg::Transform3D &toChamberTrans) const
	Returns the direction of the measurement channel in the sector frame. More...
template<class MeasType >
Amg::Vector3D	channelNormalInChamber (const MeasType *meas, const Amg::Transform3D &toChamberTrans) const
	Returns the direction, in the sector frame, of the precision axis of the measurement, i.e. More...
template<class MeasType >
	AmgSymMatrix (2) computeCov(const MeasType *primaryMeas
	Helper function to calculate the covariance for spacepoints having only one measurement. More...
template<class MeasType >
void	fillSpacePoint (std::vector< SpacePoint > &pointColl, const MeasType *primaryMeas, const MeasType *secondaryMeas, const Amg::Transform3D &toChamberTrans_eta, const Amg::Transform3D &toChamberTrans_phi) const
	Helper function that creates the spacepoint and populates it, when having two measurement. More...
template<class ContType >
StatusCode	loadContainerAndSort (const EventContext &ctx, const SG::ReadHandleKey< ContType > &key, PreSortedSpacePointMap &fillContainer) const
	Retrieve an uncalibrated measurement container <ContType> and fill the hits into the presorted space point map. More...
void	distributePointsAndStore (SpacePointsPerChamber &&hitsPerChamber, SpacePointContainer &finalContainer) const
	Distribute the premade spacepoints per chamber into their individual SpacePoint buckets. More...
void	distributePrimaryPoints (std::vector< SpacePoint > &&spacePoints, SpacePointBucketVec &splittedContainer) const
	Distributes the vector of primary eta or eta + phi space points and fills them into the buckets. More...
void	distributePhiPoints (std::vector< SpacePoint > &&spacePoints, SpacePointBucketVec &splittedContainer) const
	Distributs the vector phi space points into the buckets. More...
bool	splitBucket (const SpacePoint &spacePoint, const double firstSpPos, const SpacePointBucketVec &sortedPoints) const
	Returns whether the space point is beyond the bucket boundary. More...
void	newBucket (const SpacePoint &refSp, SpacePointBucketVec &sortedPoints) const
	Closes the current processed bucket and creates a new one. More...
template<class PrdType >
bool	passOccupancy2D (const std::vector< const PrdType * > &etaHits, const std::vector< const PrdType * > &phiHits) const
	: Check whether the occupancy cuts of hits in a gasGap are surpassed. More...
template<class MeasType >
AmgSymMatrix(2) computeCov(const xAOD	AmgSymMatrix (std::vector< SpacePoint > &pointColl, const MeasType *primaryMeas, const Amg::Transform3D &toChamberTrans) constvoid
	Helper function to calculate the covariance for spacepoints having two measurements. More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
const Amg::Vector3D &	dir
const Amg::Vector3D const Amg::Vector3D &nor	const
SG::ReadHandleKey< xAOD::MdtDriftCircleContainer >	m_mdtKey
SG::ReadHandleKey< xAOD::RpcMeasurementContainer >	m_rpcKey
SG::ReadHandleKey< xAOD::TgcStripContainer >	m_tgcKey
SG::ReadHandleKey< xAOD::MMClusterContainer >	m_mmKey
SG::ReadHandleKey< xAOD::sTgcMeasContainer >	m_stgcKey {this, "sTgcKey", "xAODsTgcMeasurements"}
SG::ReadHandleKey< ActsGeometryContext >	m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
SG::WriteHandleKey< SpacePointContainer >	m_writeKey {this, "WriteKey", "MuonSpacePoints"}
Gaudi::Property< double >	m_spacePointWindow
Gaudi::Property< double >	m_maxBucketLength
Gaudi::Property< double >	m_spacePointOverlap
Gaudi::Property< bool >	m_doStat
Gaudi::Property< unsigned int >	m_capacityBucket {this,"CapacityBucket" , 50}
std::unique_ptr< SpacePointStatistics > m_statCounter	ATLAS_THREAD_SAFE {}
Gaudi::Property< double >	m_maxOccRpcEta
Gaudi::Property< double >	m_maxOccRpcPhi
Gaudi::Property< double >	m_maxOccTgcEta
Gaudi::Property< double >	m_maxOccTgcPhi
Gaudi::Property< double >	m_maxOccStgcEta
Gaudi::Property< double >	m_maxOccStgcPhi
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 27 of file SpacePointMakerAlg.h.

Member Typedef Documentation

PreSortedSpacePointMap

using MuonR4::SpacePointMakerAlg::PreSortedSpacePointMap = std::unordered_map<const MuonGMR4::SpectrometerSector*, SpacePointsPerChamber>

private

Container abrivation of the presorted space point container per MuonChambers.

Definition at line 153 of file SpacePointMakerAlg.h.

SpacePointBucketVec

using MuonR4::SpacePointMakerAlg::SpacePointBucketVec = std::vector<SpacePointBucket>

private

Abrivation of a MuonSapcePoint bucket vector.

Definition at line 171 of file SpacePointMakerAlg.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

~SpacePointMakerAlg()

MuonR4::SpacePointMakerAlg::~SpacePointMakerAlg ( )

default

Member Function Documentation

AmgSymMatrix() [1/2]

template<class MeasType >

MuonR4::SpacePointMakerAlg::AmgSymMatrix ( 2  ) const

private

Helper function to calculate the covariance for spacepoints having only one measurement.

Parameters

primaryMeas	Uncalibrated muon measurement
dir	direction of the measurement channel in the sector frame
nor	direction of the precision axis of the measurement in the sector frame

AmgSymMatrix() [2/2]

template<class MeasType >

AmgSymMatrix (2) computeCov( const xAOD MuonR4::SpacePointMakerAlg::AmgSymMatrix ( std::vector< SpacePoint > &  pointColl, const MeasType *  primaryMeas, const Amg::Transform3D &  toChamberTrans  ) const

private

Helper function to calculate the covariance for spacepoints having two measurements.

Parameters

primaryMeas	Primary uncalibrated muon measurement
secondaryMeas	Secondary uncalibrated muon measurement
nor1	direction of the precision axis of the primary measurement in the sector frame
nor2	direction of the precision axis of the secondary measurement in the sector frame

Helper function that creates the spacepoint and populates it, when having one measurement

Parameters

pointColl	Collection of spacepoints where the new one is being push
primaryMeas	Uncalibrated muon measurement
toChamberTrans	transform from the attached Muon chamber frame to the sector frame

cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality

overridevirtualinherited

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

Override this to return 0 for reentrant algorithms.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

channelDirInChamber()

template<class MeasType >

Amg::Vector3D MuonR4::SpacePointMakerAlg::channelDirInChamber ( const MeasType *  meas, const Amg::Transform3D &  toChamberTrans  ) const

private

Returns the direction of the measurement channel in the sector frame.

Parameters

meas	Uncalibrated muon measurement
toChamberTrans	transform from the attached Muon chamber frame to the sector frame

Definition at line 75 of file SpacePointMakerAlg.cxx.

                                                                               {       
    if constexpr (std::is_same_v<MeasType, xAOD::MdtDriftCircle>){
        return toChamberTrans.linear() * Amg::Vector3D::UnitZ();
    }
    else if constexpr (std::is_same_v<MeasType, xAOD::RpcMeasurement> or 
                       std::is_same_v<MeasType, xAOD::MMCluster> or 
                       std::is_same_v<MeasType, xAOD::sTgcMeasurement>){
        return toChamberTrans.linear() * Amg::Vector3D::UnitY();
    }
    else if constexpr (std::is_same_v<MeasType, xAOD::TgcStrip>) {           
        const auto& stripLay = meas->readoutElement()->sensorLayout(meas->layerHash());
        if (meas->measuresPhi()) {
            const auto& sDesign = static_cast<const MuonGMR4::RadialStripDesign&>(stripLay->design(true));
            return toChamberTrans.linear() * stripLay->to3D(sDesign.stripDir(meas->channelNumber()), true);
        }
        return toChamberTrans.linear() * stripLay->to3D(stripLay->design().stripDir(), false);
    }
    else static_assert(std::false_type::value, "Unsupported measurement type.");
    return Amg::Vector3D::Zero();      
}

channelNormalInChamber()

template<class MeasType >

Amg::Vector3D MuonR4::SpacePointMakerAlg::channelNormalInChamber ( const MeasType *  meas, const Amg::Transform3D &  toChamberTrans  ) const

private

Returns the direction, in the sector frame, of the precision axis of the measurement, i.e.

the vector pointing to the next strip or tube

Parameters

meas	Uncalibrated muon measurement
toChamberTrans	transform from the attached Muon chamber frame to the sector frame

Definition at line 98 of file SpacePointMakerAlg.cxx.

                                                                                                    {
    if constexpr (std::is_same_v<MeasType, xAOD::MdtDriftCircle>){
        return toChamberTrans.linear() * Amg::Vector3D::UnitY();
    }
    else if constexpr (std::is_same_v<MeasType, xAOD::RpcMeasurement> or 
                       std::is_same_v<MeasType, xAOD::MMCluster> or 
                       std::is_same_v<MeasType, xAOD::sTgcMeasurement>){
        return toChamberTrans.linear() * Amg::Vector3D::UnitX();
    }
    else if constexpr (std::is_same_v<MeasType, xAOD::TgcStrip>) {                     
        const auto& stripLay = meas->readoutElement()->sensorLayout(meas->layerHash());
        if (meas->measuresPhi()) {
            const auto& sDesign = static_cast<const MuonGMR4::RadialStripDesign&>(stripLay->design(true));
            return toChamberTrans.linear() * stripLay->to3D(sDesign.stripNormal(meas->channelNumber()), true);
        }
        return toChamberTrans.linear() * stripLay->to3D(stripLay->design().stripNormal(), false);
    }
    else static_assert(std::false_type::value, "Unsupported measurement type.");
    return Amg::Vector3D::Zero();     
}

declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

declareGaudiProperty() [2/4]

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

declareGaudiProperty() [4/4]

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

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(t);
  }

declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleBase &  hndl, const std::string &  doc, const SG::VarHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKey &  hndl, const std::string &  doc, const SG::VarHandleKeyType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name, SG::VarHandleKeyArray &  hndArr, const std::string &  doc, const SG::VarHandleKeyArrayType &    )

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc, const SG::NotHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(name, property, doc);
  }

declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc = "none"  )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

  {
    typedef typename SG::HandleClassifier<T>::type htype;
    return declareProperty (name, property, doc, htype());
  }

declareProperty() [6/6]

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

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

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

detStore()

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

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

distributePhiPoints()

void MuonR4::SpacePointMakerAlg::distributePhiPoints ( std::vector< SpacePoint > &&  spacePoints, SpacePointBucketVec &  splittedContainer  ) const

private

Distributs the vector phi space points into the buckets.

In contrast to the primary distribution no new buckets are created and the points are distributed into the existing ones instead.

Parameters

spacePoint	Vecotr of phi space points to sort into the buckets
splittedContainer	Output vector containing all defined bucket

If maxY is smaller than the lower cov boundary or minY is bigger than the other boundary, there's definetely no overlap

Definition at line 535 of file SpacePointMakerAlg.cxx.

                                                                                          {
    for (SpacePoint& sp : spacePoints) {
        auto phiPoint = std::make_shared<SpacePoint>(std::move(sp));
        const double minY = phiPoint->positionInChamber().y() - phiPoint->uncertainty().y();
        const double maxY = phiPoint->positionInChamber().y() + phiPoint->uncertainty().y();
        for (SpacePointBucket& bucket : splittedContainer){
            if (! (maxY < bucket.coveredMin() || bucket.coveredMax() < minY) ) {
                bucket.emplace_back(phiPoint);
            }
        }
    }
}

distributePointsAndStore()

void MuonR4::SpacePointMakerAlg::distributePointsAndStore ( SpacePointsPerChamber &&  hitsPerChamber, SpacePointContainer &  finalContainer  ) const

private

Distribute the premade spacepoints per chamber into their individual SpacePoint buckets.

A new bucket is created everytime if the hit to fill is along the z-axis farther away from the first point in the bucket than the <spacePointWindowSize>. Hit in the previous bucket which are <spacePointOverlap> away from the first hit in the new bucket are also mirrored. The bucket formation starts with the eta Muon space points and then consumes the phi hits.

Parameters

ctx	Event context of the current event
hitsPerChamber	List of all premade space points which have to be sorted
finalContainer	Output SpacePoint bucket container.

Definition at line 500 of file SpacePointMakerAlg.cxx.

                                                                                             {
    SpacePointBucketVec splittedHits{};
    splittedHits.emplace_back();
    if (m_statCounter){
        m_statCounter->addToStat(hitsPerChamber.etaHits);
        m_statCounter->addToStat(hitsPerChamber.phiHits);
 
    }
    distributePrimaryPoints(std::move(hitsPerChamber.etaHits), splittedHits);
    splittedHits.erase(std::remove_if(splittedHits.begin(), splittedHits.end(),
                       [](const SpacePointBucket& bucket) {
                           return bucket.size() <= 1;
                       }), splittedHits.end());
    distributePhiPoints(std::move(hitsPerChamber.phiHits), splittedHits);
    
    for (SpacePointBucket& bucket : splittedHits) {
 
        std::ranges::sort(bucket, MuonR4::SpacePointPerLayerSorter{});
 
        if (msgLvl(MSG::VERBOSE)){
            std::stringstream spStr{};
            for (const std::shared_ptr<SpacePoint>& sp : bucket){
                spStr<<"SpacePoint: PrimaryMeas: " << m_idHelperSvc->toString(sp->identify()) << " SecondaryMeas: " 
                 <<(sp->secondaryMeasurement() ? m_idHelperSvc->toString(xAOD::identify(sp->secondaryMeasurement())) : " None ") 
                 << " Pos: " <<  Amg::toString(sp->positionInChamber())<<std::endl;
            }
            ATH_MSG_VERBOSE("Created a bucket, printing all spacepoints..."<<std::endl<<spStr.str());
        }
 
        bucket.populateChamberLocations();
        finalContainer.push_back(std::make_unique<SpacePointBucket>(std::move(bucket)));
    }
 
}

distributePrimaryPoints()

void MuonR4::SpacePointMakerAlg::distributePrimaryPoints ( std::vector< SpacePoint > &&  spacePoints, SpacePointBucketVec &  splittedContainer  ) const

private

Distributes the vector of primary eta or eta + phi space points and fills them into the buckets.

The buckets are dynamically created based on the distance of the new space point to sort to the previous or the first space point in the bucket.

Parameters

spacePoints	Vector of space points to sort into the buckets
splittedContainer	Output vector containing all defined bucket

Order the space points by local chamber y which is along the tube-plane.

Definition at line 588 of file SpacePointMakerAlg.cxx.

                                                                                          {
 
    if (spacePoints.empty()) return;
   
    std::ranges::sort(spacePoints, 
                      [] (const SpacePoint& a, const SpacePoint& b) {
                        return a.positionInChamber().y() < b.positionInChamber().y();
                      });
 
    double firstPointPos = spacePoints.front().positionInChamber().y();
    
    for (SpacePoint& toSort : spacePoints) {        
        ATH_MSG_VERBOSE("Add new primary space point "<<m_idHelperSvc->toString(toSort.identify())<<", "
                     <<" @ "<<Amg::toString(toSort.positionInChamber()));
 
        if (splitBucket(toSort, firstPointPos, splittedHits)){
            newBucket(toSort, splittedHits);
            firstPointPos = splittedHits.back().empty() ? toSort.positionInChamber().y() : splittedHits.back().front()->positionInChamber().y();
            ATH_MSG_VERBOSE("New bucket: id " << splittedHits.back().bucketId() << " Coverage: " << firstPointPos);
        }
        std::shared_ptr<SpacePoint> spacePoint = std::make_shared<SpacePoint>(std::move(toSort));
        splittedHits.back().emplace_back(spacePoint);
    }
    SpacePointBucket& lastBucket{splittedHits.back()};
    lastBucket.setCoveredRange(lastBucket.front()->positionInChamber().y(), 
                               lastBucket.back()->positionInChamber().y());
}

evtStore() [1/2]

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode MuonR4::SpacePointMakerAlg::execute ( const EventContext &  ctx ) const

override

Definition at line 482 of file SpacePointMakerAlg.cxx.

                                                                    {
    PreSortedSpacePointMap preSortedContainer{};
    ATH_CHECK(loadContainerAndSort(ctx, m_mdtKey, preSortedContainer));
    ATH_CHECK(loadContainerAndSort(ctx, m_rpcKey, preSortedContainer));
    ATH_CHECK(loadContainerAndSort(ctx, m_tgcKey, preSortedContainer));
    ATH_CHECK(loadContainerAndSort(ctx, m_mmKey, preSortedContainer));
    ATH_CHECK(loadContainerAndSort(ctx, m_stgcKey, preSortedContainer));
    std::unique_ptr<SpacePointContainer> outContainer = std::make_unique<SpacePointContainer>();
    
    for (auto &[chamber, hitsPerChamber] : preSortedContainer){
        ATH_MSG_DEBUG("Fill space points for chamber "<<chamber->identString());
        distributePointsAndStore(std::move(hitsPerChamber), *outContainer);
    }
    SG::WriteHandle<SpacePointContainer> writeHandle{m_writeKey, ctx};
    ATH_CHECK(writeHandle.record(std::move(outContainer)));
    return StatusCode::SUCCESS;
}

extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

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

fillSpacePoint()

template<class MeasType >

void MuonR4::SpacePointMakerAlg::fillSpacePoint ( std::vector< SpacePoint > &  pointColl, const MeasType *  primaryMeas, const MeasType *  secondaryMeas, const Amg::Transform3D &  toChamberTrans_eta, const Amg::Transform3D &  toChamberTrans_phi  ) const

private

Helper function that creates the spacepoint and populates it, when having two measurement.

Parameters

pointColl	Collection of spacepoints where the new one is being push
primaryMeas	Primary uncalibrated muon measurement
secondaryMeas	Secondary uncalibrated muon measurement
toChamberTrans_eta	transform from the attached Muon chamber frame to the sector frame for the primary meas
toChamberTrans_phi	transform from the attached Muon chamber frame to the sector frame for the secondary meas

filterPassed()

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

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

finalize()

StatusCode MuonR4::SpacePointMakerAlg::finalize ( )

override

Definition at line 292 of file SpacePointMakerAlg.cxx.

                                        {
    if (m_statCounter) {
        m_statCounter->dumpStatisics(msgStream());
    }
    return StatusCode::SUCCESS;
}

initialize()

StatusCode MuonR4::SpacePointMakerAlg::initialize ( )

override

Definition at line 298 of file SpacePointMakerAlg.cxx.

                                          {
    ATH_CHECK(m_geoCtxKey.initialize());
    ATH_CHECK(m_mdtKey.initialize(!m_mdtKey.empty()));
    ATH_CHECK(m_rpcKey.initialize(!m_rpcKey.empty()));
    ATH_CHECK(m_tgcKey.initialize(!m_tgcKey.empty()));
    ATH_CHECK(m_mmKey.initialize(!m_mmKey.empty()));
    ATH_CHECK(m_stgcKey.initialize(!m_stgcKey.empty()));
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK(m_writeKey.initialize());
    if (m_doStat) m_statCounter = std::make_unique<SpacePointStatistics>(m_idHelperSvc.get());
    return StatusCode::SUCCESS;
}

inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

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

loadContainerAndSort()

template<class ContType >

StatusCode MuonR4::SpacePointMakerAlg::loadContainerAndSort ( const EventContext &  ctx, const SG::ReadHandleKey< ContType > &  key, PreSortedSpacePointMap &  fillContainer  ) const

private

Retrieve an uncalibrated measurement container <ContType> and fill the hits into the presorted space point map.

Per associated MuonChamber, hits from Tgc, Rpc, sTgcs are grouped by their gasGap location and then divided into eta & phi measurements. If both are found, each eta measurement is combined with phi measurement into a SpacePoint. In any other case, the measurements are just transformed into a SpacePoint.

Parameters

ctx	Event context of the current event
key	ReadHandleKey to access the container of data type <ContType>
fillContainer	Global container into which all space points are filled.

Loop over the chamber hits to split the hits per gasGap

Make directly to a space point

Wires measure the phi coordinate

Tgc & Rpcs have the measuresPhi property

Sort in the hit

Create the space points

Simple combination by taking the cross-product

There's no valid combination with another phi hit

Tgc measurements with different bunch crossing tags cannot be combined

If there's a phi measuremnt which cannot be combined with the others or no eta measurement is suitable, then manually push_back the phi hits

Definition at line 346 of file SpacePointMakerAlg.cxx.

                                                                                                     {
    const ContType* measurementCont{nullptr};
    ATH_CHECK(SG::get(measurementCont, key, ctx));
    if (!measurementCont || measurementCont->empty()){
        ATH_MSG_DEBUG("nothing to do"); 
        return StatusCode::SUCCESS;
    }
    const ActsGeometryContext* gctx{nullptr};
    ATH_CHECK(SG::get(gctx, m_geoCtxKey, ctx));
    
    using PrdType = typename ContType::const_value_type;
    using PrdVec = std::vector<PrdType>;
    xAOD::ChamberViewer viewer{*measurementCont};
    do {
 
      SpacePointsPerChamber& pointsInChamb = fillContainer[viewer.at(0)->readoutElement()->msSector()];
      const Amg::Transform3D& sectorTrans = viewer.at(0)->readoutElement()->msSector()->globalToLocalTrans(*gctx);
      ATH_MSG_DEBUG("Fill space points for chamber "<<m_idHelperSvc->toStringDetEl(viewer.at(0)->identify()));
      if constexpr( std::is_same_v<ContType, xAOD::MdtDriftCircleContainer> ||
                    std::is_same_v<ContType, xAOD::MMClusterContainer>) {
 
            pointsInChamb.etaHits.reserve(pointsInChamb.etaHits.capacity() + viewer.size());
            for (const PrdType prd: viewer) {
                ATH_MSG_VERBOSE("Create space point from "<<m_idHelperSvc->toString(prd->identify())
                              <<", hash: "<<prd->identifierHash());
 
                Amg::Transform3D toChamberTrans{ toChamberTransform(*gctx, sectorTrans,prd)};
                fillSpacePoint(pointsInChamb.etaHits, prd, toChamberTrans);
            }
       } else {
            using EtaPhiHits = std::array<PrdVec, 2>;
            std::vector<EtaPhiHits> hitsPerGasGap{};
            for (const PrdType prd : viewer) {
                ATH_MSG_VERBOSE("Create space point from "<<m_idHelperSvc->toString(prd->identify())<<", hash: "<<prd->identifierHash());
                
                unsigned int gapIdx = prd->gasGap() -1;
                if constexpr (std::is_same_v<ContType, xAOD::RpcMeasurementContainer>) {
                    gapIdx = prd->readoutElement()->createHash(0, prd->gasGap(), prd->doubletPhi(), false);
                }
 
                bool measPhi{false};
                if constexpr(std::is_same_v<ContType, xAOD::sTgcMeasContainer>) {
                    if (prd->channelType() == sTgcIdHelper::sTgcChannelTypes::Pad) {
                        Amg::Transform3D toChamberTrans{ toChamberTransform(*gctx, sectorTrans, prd)};
                        fillSpacePoint(pointsInChamb.etaHits, prd, toChamberTrans);
                        continue;
                    }
                    measPhi = prd->channelType() == sTgcIdHelper::sTgcChannelTypes::Wire;
                } else {
                    measPhi = prd->measuresPhi();
                }
 
                if (hitsPerGasGap.size() <= gapIdx) {
                    hitsPerGasGap.resize(gapIdx + 1);
                }
                PrdVec& toPush = hitsPerGasGap[gapIdx][measPhi];
                if (toPush.capacity() == toPush.size()) {
                    toPush.reserve(toPush.size() + m_capacityBucket);
                }
                toPush.push_back(prd);
            }
            for (auto& [etaHits, phiHits] : hitsPerGasGap) {
 
                Amg::Transform3D toChamberTrans_eta {}, toChamberTrans_phi {};
                if (!etaHits.empty()){
                    toChamberTrans_eta = toChamberTransform(*gctx, sectorTrans, etaHits.at(0));
                }
                if (!phiHits.empty()){
                    toChamberTrans_phi = toChamberTransform(*gctx, sectorTrans, phiHits.at(0));
                }
 
                if (!passOccupancy2D(etaHits, phiHits)) {
                    ATH_MSG_VERBOSE("Occupancy cut not passed "<<etaHits.size()<<", "<<phiHits.size());
                    pointsInChamb.etaHits.reserve(pointsInChamb.etaHits.size() + etaHits.size());
                    pointsInChamb.phiHits.reserve(pointsInChamb.phiHits.size() + phiHits.size());
                    for (const PrdType etaPrd : etaHits) {
 
                        fillSpacePoint(pointsInChamb.etaHits, etaPrd, toChamberTrans_eta);
                        ATH_MSG_VERBOSE("Add new eta hit "<<m_idHelperSvc->toString(pointsInChamb.etaHits.back().identify())
                                <<" "<<Amg::toString(pointsInChamb.etaHits.back().positionInChamber()));
 
                    }
                    for (const PrdType phiPrd : phiHits) {    
                        fillSpacePoint(pointsInChamb.phiHits, phiPrd, toChamberTrans_phi);
                        ATH_MSG_VERBOSE("Add new phi hit "<<m_idHelperSvc->toString(pointsInChamb.phiHits.back().identify())
                                <<" "<<Amg::toString(pointsInChamb.phiHits.back().positionInChamber()));
                    }
                    continue;
                }
                std::vector<std::shared_ptr<unsigned>> etaCounts{matchCountVec(etaHits.size())}, 
                                                       phiCounts{matchCountVec(phiHits.size())};
                pointsInChamb.etaHits.reserve(etaHits.size()*phiHits.size());
                for (unsigned int etaP = 0; etaP < etaHits.size(); ++etaP) {
                    for (unsigned int phiP = 0; phiP < phiHits.size(); ++ phiP) {
                        if constexpr(std::is_same<xAOD::TgcStripContainer, ContType>::value) {
                            if (!(etaHits[etaP]->bcBitMap() & phiHits[phiP]->bcBitMap())){
                                continue;
                            }
                        }
                        fillSpacePoint(pointsInChamb.etaHits, etaHits[etaP], phiHits[phiP], toChamberTrans_eta, toChamberTrans_phi);
 
                        SpacePoint& sp {pointsInChamb.etaHits.back()};
                        ATH_MSG_VERBOSE("Create new spacepoint from "<<m_idHelperSvc->toString(etaHits[etaP]->identify())
                        <<" & "<<m_idHelperSvc->toString(phiHits[phiP]->identify())<<" at "<<Amg::toString(sp.positionInChamber()));
                        sp.setInstanceCounts(etaCounts[etaP], phiCounts[phiP]);
                    }
                    if (!(*etaCounts[etaP])) {
                        fillSpacePoint(pointsInChamb.etaHits, etaHits[etaP], toChamberTrans_eta);
                        continue;
                    }
                }
                for (unsigned int phiP = 0; phiP < phiHits.size(); ++ phiP){
                    if (!(*phiCounts[phiP])) {
                        fillSpacePoint(pointsInChamb.phiHits, phiHits[phiP], toChamberTrans_phi);
                    }
                }
            }
       }
    } while (viewer.next());
    return StatusCode::SUCCESS;
}

msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

newBucket()

void MuonR4::SpacePointMakerAlg::newBucket ( const SpacePoint &  refSp, SpacePointBucketVec &  sortedPoints  ) const

private

Closes the current processed bucket and creates a new one.

Space points of the previous bucket within the overlap region to the first space point of the new bucket are copied over

Parameters

refSp	First new space point which will be added to the new bucket.
sortedPoints	List of all processed buckets in the chamber. The list is augmented by 1 element

Set the boundaries from the previous bucket

Copy space points that could be within the overlap region to the next bucket

Definition at line 564 of file SpacePointMakerAlg.cxx.

                                                                            {
    SpacePointBucket& newContainer = sortedPoints.emplace_back();
    newContainer.setBucketId(sortedPoints.size() -1);
 
    SpacePointBucket& overlap{sortedPoints[sortedPoints.size() - 2]};
    overlap.setCoveredRange(overlap.front()->positionInChamber().y(), 
                            overlap.back()->positionInChamber().y());
    
    const double refBound = refSpacePoint.positionInChamber().y();
                                
    for (const std::shared_ptr<SpacePoint>& pointInBucket : overlap | std::views::reverse) {
        const double overlapPos = pointInBucket->positionInChamber().y() + pointInBucket->uncertainty()[Amg::y];
        if (refBound - overlapPos < m_spacePointOverlap) {    
            newContainer.insert(newContainer.begin(), pointInBucket);
        } else {
            break;
        }
    }    
 
}

outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

passOccupancy2D() [1/4]

template<class PrdType >

bool MuonR4::SpacePointMakerAlg::passOccupancy2D ( const std::vector< const PrdType * > &  etaHits, const std::vector< const PrdType * > &  phiHits  ) const

private

: Check whether the occupancy cuts of hits in a gasGap are surpassed.

The method is specified for each of the 3 strip technologies, Rpc, Tgc, sTgc and applies a technology-dependent upper bound on the number of phi & eta hits. If the threshold is surpassed, only 1D space points are built intsead of 2D ones

Parameters

etaHits	List of all presorted eta measurements in a gas gap
phiHits	List of all presorted phi measurements in a gas gap

passOccupancy2D() [2/4]

template<>

bool MuonR4::SpacePointMakerAlg::passOccupancy2D	(	const std::vector< const xAOD::RpcMeasurement * > &	etaHits,
		const std::vector< const xAOD::RpcMeasurement * > &	phiHits
	)		const

Definition at line 323 of file SpacePointMakerAlg.cxx.

                                                                                                        {
        if (etaHits.empty() || phiHits.empty()) {
            return false;
        }
        const MuonGMR4::RpcReadoutElement* re = etaHits[0]->readoutElement();
        ATH_MSG_VERBOSE("Collected "<<etaHits.size()<<"/"<<phiHits.size()<<" hits in "<<m_idHelperSvc->toStringGasGap(etaHits[0]->identify()));
        return ((1.*etaHits.size()) / (1.*re->nEtaStrips())) < m_maxOccRpcEta &&
               ((1.*phiHits.size()) / (1.*re->nPhiStrips())) < m_maxOccRpcPhi;
    }

passOccupancy2D() [3/4]

template<>

bool MuonR4::SpacePointMakerAlg::passOccupancy2D	(	const std::vector< const xAOD::sTgcMeasurement * > &	etaHits,
		const std::vector< const xAOD::sTgcMeasurement * > &	phiHits
	)		const

Definition at line 335 of file SpacePointMakerAlg.cxx.

                                                                                                         {
        if (etaHits.empty() || phiHits.empty()) {
            return false;
        }
        const MuonGMR4::sTgcReadoutElement* re = etaHits[0]->readoutElement();
        return ((1.*etaHits.size()) / (1.*re->numChannels(etaHits[0]->measurementHash()))) < m_maxOccStgcEta &&
               ((1.*phiHits.size()) / (1.*re->numChannels(phiHits[0]->measurementHash()))) < m_maxOccStgcPhi;
    }

passOccupancy2D() [4/4]

template<>

bool MuonR4::SpacePointMakerAlg::passOccupancy2D	(	const std::vector< const xAOD::TgcStrip * > &	etaHits,
		const std::vector< const xAOD::TgcStrip * > &	phiHits
	)		const

Definition at line 312 of file SpacePointMakerAlg.cxx.

                                                                                                  {
        if (etaHits.empty() || phiHits.empty()) {
            return false;
        }
        const MuonGMR4::TgcReadoutElement* re = etaHits[0]->readoutElement();
        ATH_MSG_VERBOSE("Collected "<<etaHits.size()<<"/"<<phiHits.size()<<" hits in "<<m_idHelperSvc->toStringGasGap(etaHits[0]->identify()));
        return ((1.*etaHits.size()) / ((1.*re->numChannels(etaHits[0]->measurementHash())))) < m_maxOccTgcEta &&
               ((1.*phiHits.size()) / ((1.*re->numChannels(phiHits[0]->measurementHash())))) < m_maxOccTgcPhi;
    }

positionInChamber()

template<class MeasType >

Amg::Vector3D MuonR4::SpacePointMakerAlg::positionInChamber ( const MeasType *  meas, const Amg::Transform3D &  toChamberTrans  ) const

private

Returns the position of the uncalibrated muon measurement in the sector frame.

Parameters

meas	Uncalibrated muon measurement
toChamberTrans	transform from the attached Muon chamber frame to the sector frame

Definition at line 48 of file SpacePointMakerAlg.cxx.

                                                                                               {
    if constexpr (std::is_same_v<MeasType, xAOD::MdtDriftCircle>){
        return toChamberTrans * meas->localCirclePosition();
    } else if constexpr (std::is_same_v<MeasType, xAOD::RpcMeasurement>){
        return toChamberTrans * meas->localMeasurementPos();
    } else if constexpr (std::is_same_v<MeasType, xAOD::TgcStrip>){
        const auto& stripLay = meas->readoutElement()->sensorLayout(meas->layerHash());
        return toChamberTrans * stripLay->to3D(meas->template localPosition<1>()[0]*Amg::Vector2D::UnitX(),
                                               meas->measuresPhi());
    } else if constexpr (std::is_same_v<MeasType, xAOD::MMCluster>){
        return toChamberTrans * (meas->template localPosition<1>()[Trk::locX] * Amg::Vector3D::UnitX());
    }
    else if constexpr (std::is_same_v<MeasType, xAOD::sTgcMeasurement>){
        if (meas->channelType() == sTgcIdHelper::sTgcChannelTypes::Strip ||
            meas->channelType() == sTgcIdHelper::sTgcChannelTypes::Wire) {
                return toChamberTrans * (meas->template localPosition<1>()[Trk::locX] * Amg::Vector3D::UnitX());
        }
        Amg::Vector3D locPos{Amg::Vector3D::Zero()};
        locPos.block<2,1>(0,0) = xAOD::toEigen(meas->template localPosition<2>());
        return toChamberTrans * locPos;
    }
    else static_assert(std::false_type::value, "Unsupported measurement type.");
    return Amg::Vector3D::Zero();
}

renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

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

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

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

splitBucket()

bool MuonR4::SpacePointMakerAlg::splitBucket ( const SpacePoint &  spacePoint, const double  firstSpPos, const SpacePointBucketVec &  sortedPoints  ) const

private

Returns whether the space point is beyond the bucket boundary.

Parameters

spacePoint	Space point candidate to add to the bucket
sortedPoints	Container of all defined buckets in the chamber

Distance between this point and the first one exceeds the maximum length

Definition at line 550 of file SpacePointMakerAlg.cxx.

                                                                                    {
    const double spY = spacePoint.positionInChamber().y();
    if (spY - firstSpPos > m_maxBucketLength){
        return true;
    }
    
    if (sortedPoints.empty() || sortedPoints.back().empty()) {
        return false;
    }
    return spY - sortedPoints.back().back()->positionInChamber().y() > m_spacePointWindow;
}

sysExecute()

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

overridevirtualinherited

Execute an algorithm.

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

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

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

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

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in InputMakerBase, and HypoBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

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

sysStart()

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

overridevirtualinherited

Handle START transition.

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

toChamberTransform()

template<class MeasType >

Amg::Transform3D MuonR4::SpacePointMakerAlg::toChamberTransform ( const ActsGeometryContext &  gctx, const Amg::Transform3D &  sectorTrans, const MeasType *  meas  ) const

private

Returns the transform from the attached Muon chamber frame to the sector frame.

Parameters

gctx	Geometry context to calculate the relative alignment between chamber & measurement
meas	Uncalibrated muon measurement

Definition at line 34 of file SpacePointMakerAlg.cxx.

                                                                                   {
    IdentifierHash hash = {};
    if constexpr(std::is_same_v<MeasType, xAOD::MdtDriftCircle>) {
        hash = meas->measurementHash();
    } else {
        hash = meas->layerHash();
    }
    const MuonGMR4::MuonReadoutElement* reEle{meas->readoutElement()};
    return sectorTrans * reEle->localToGlobalTrans(gctx, hash);
}

updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

ATLAS_THREAD_SAFE

std::unique_ptr<SpacePointStatistics> m_statCounter MuonR4::SpacePointMakerAlg::ATLAS_THREAD_SAFE {}

private

Definition at line 256 of file SpacePointMakerAlg.h.

const

const Amg::Vector3D const Amg::Vector3D& nor MuonR4::SpacePointMakerAlg::const

private

Definition at line 66 of file SpacePointMakerAlg.h.

dir

const Amg::Vector3D& MuonR4::SpacePointMakerAlg::dir

private

Definition at line 65 of file SpacePointMakerAlg.h.

m_capacityBucket

Gaudi::Property<unsigned int> MuonR4::SpacePointMakerAlg::m_capacityBucket {this,"CapacityBucket" , 50}

private

Definition at line 255 of file SpacePointMakerAlg.h.

m_detStore

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

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_doStat

Gaudi::Property<bool> MuonR4::SpacePointMakerAlg::m_doStat

private

Initial value:

{this, "doStats", false, 
                                           "If enabled the algorithm keeps track how many hits have been made" }

Definition at line 252 of file SpacePointMakerAlg.h.

m_evtStore

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

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_geoCtxKey

SG::ReadHandleKey<ActsGeometryContext> MuonR4::SpacePointMakerAlg::m_geoCtxKey {this, "AlignmentKey", "ActsAlignment", "cond handle key"}

private

Definition at line 236 of file SpacePointMakerAlg.h.

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> MuonR4::SpacePointMakerAlg::m_idHelperSvc {this, "IdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 238 of file SpacePointMakerAlg.h.

m_maxBucketLength

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_maxBucketLength

private

Initial value:

{this, "maxBucketLength", 2.*Gaudi::Units::m,
                                                       "Maximal size of a space point bucket"}

Definition at line 245 of file SpacePointMakerAlg.h.

m_maxOccRpcEta

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_maxOccRpcEta

private

Initial value:

{this, "maxRpcEtaOccupancy", 0.1, 
                                                   "Maximum occpancy of Rpc eta hits in a gasGap"}

Definition at line 258 of file SpacePointMakerAlg.h.

m_maxOccRpcPhi

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_maxOccRpcPhi

private

Initial value:

{this, "maxRpcPhiOccupancy", 0.1, 
                                                   "Maximum occpancy of Rpc phi hits in a gasGap"}

Definition at line 260 of file SpacePointMakerAlg.h.

m_maxOccStgcEta

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_maxOccStgcEta

private

Initial value:

{this, "maxSTGCEtaOccupancy", 0.1, 
                                                   "Maximum occpancy of sTgc eta hits in a gasGap"}

Definition at line 268 of file SpacePointMakerAlg.h.

m_maxOccStgcPhi

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_maxOccStgcPhi

private

Initial value:

{this, "maxSTGCPhiOccupancy", 0.1, 
                                                    "Maximum occpancy of sTgc phi hits in a gasGap"}

Definition at line 270 of file SpacePointMakerAlg.h.

m_maxOccTgcEta

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_maxOccTgcEta

private

Initial value:

{this, "maxTgcEtaOccupancy", 0.1, 
                                                   "Maximum occpancy of Tgc eta hits in a gasGap"}

Definition at line 263 of file SpacePointMakerAlg.h.

m_maxOccTgcPhi

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_maxOccTgcPhi

private

Initial value:

{this, "maxTgcPhiOccupancy", 0.1, 
                                                   "Maximum occpancy of Tgc phi hits in a gasGap"}

Definition at line 265 of file SpacePointMakerAlg.h.

m_mdtKey

SG::ReadHandleKey<xAOD::MdtDriftCircleContainer> MuonR4::SpacePointMakerAlg::m_mdtKey

private

Initial value:

{this, "MdtKey", "xMdtMeasurements",
                                                                      "Key to the uncalibrated Drift circle measurements"}

Definition at line 221 of file SpacePointMakerAlg.h.

m_mmKey

SG::ReadHandleKey<xAOD::MMClusterContainer> MuonR4::SpacePointMakerAlg::m_mmKey

private

Initial value:

{this, "MmKey", "xAODMMClusters",
                                                                "Key to the uncalibrated 1D Mm hits"}

Definition at line 230 of file SpacePointMakerAlg.h.

m_rpcKey

SG::ReadHandleKey<xAOD::RpcMeasurementContainer> MuonR4::SpacePointMakerAlg::m_rpcKey

private

Initial value:

{this, "RpcKey", "xRpcMeasurements",
                                                                "Key to the uncalibrated 1D rpc hits"}

Definition at line 224 of file SpacePointMakerAlg.h.

m_spacePointOverlap

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_spacePointOverlap

private

Initial value:

{this, "spacePointOverlap", 25.*Gaudi::Units::cm,
                                                        "Hits that are within <spacePointOverlap> of the bucket margin. "
                                                        "Are copied to the next bucket"}

Definition at line 248 of file SpacePointMakerAlg.h.

m_spacePointWindow

Gaudi::Property<double> MuonR4::SpacePointMakerAlg::m_spacePointWindow

private

Initial value:

{this, "spacePointWindowSize", 0.8*Gaudi::Units::m,
                                                       "Maximal distance between consecutive hits in a bucket"}

Definition at line 242 of file SpacePointMakerAlg.h.

m_stgcKey

SG::ReadHandleKey<xAOD::sTgcMeasContainer> MuonR4::SpacePointMakerAlg::m_stgcKey {this, "sTgcKey", "xAODsTgcMeasurements"}

private

Definition at line 233 of file SpacePointMakerAlg.h.

m_tgcKey

SG::ReadHandleKey<xAOD::TgcStripContainer> MuonR4::SpacePointMakerAlg::m_tgcKey

private

Initial value:

{this, "TgcKey", "xTgcStrips",
                                                                "Key to the uncalibrated 1D tgc hits"}

Definition at line 227 of file SpacePointMakerAlg.h.

m_varHandleArraysDeclared

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

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_writeKey

SG::WriteHandleKey<SpacePointContainer> MuonR4::SpacePointMakerAlg::m_writeKey {this, "WriteKey", "MuonSpacePoints"}

private

Definition at line 240 of file SpacePointMakerAlg.h.

---

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

• SpacePointMakerAlg.h
• SpacePointMakerAlg.cxx