Access to all raw hits with layer/station granularity. More...

#include <MuonCalibExtendedRawHits.h>

Collaboration diagram for MuonCalib::MuonCalibExtendedRawHits:

Classes
struct	ExtendedChannelInfo
struct	ExtendedInfo
struct	MdtRegionInfo

Public Types
using	MdtCalibRawHitPtr = MuonCalibRawHitCollection::MdtCalibRawHitPtr
using	CscCalibRawHitPtr = MuonCalibRawHitCollection::CscCalibRawHitPtr
using	RpcCalibRawHitPtr = MuonCalibRawHitCollection::RpcCalibRawHitPtr
using	TgcCalibRawHitPtr = MuonCalibRawHitCollection::TgcCalibRawHitPtr
using	ExtendedSegPtr = std::shared_ptr<MuonCalibExtendedSegment>
using	ExtendedTrkPtr = std::shared_ptr<MuonCalibExtendedTrack>
typedef std::map< MuonFixedId, ExtendedChannelInfo >	MdtHitInfoMap
typedef MdtHitInfoMap::iterator	MdtHitInfoIt
typedef MdtHitInfoMap::const_iterator	MdtHitInfoCit
typedef std::map< MuonFixedId, MdtRegionInfo >	MdtHitChamberIdMap
typedef std::map< MuonFixedIdManipulator::StationIndex, MdtRegionInfo >	MdtHitStationIndexMap
typedef std::pair< RpcCalibRawHitPtr, ExtendedInfo >	RpcHitInfo
typedef std::vector< RpcHitInfo >	RpcHitInfoVec
typedef std::map< MuonFixedId, RpcHitInfo >	RpcHitChamberIdMap
typedef std::map< MuonFixedIdManipulator::StationIndex, RpcHitInfoVec >	RpcHitStationIndexMap
typedef std::pair< TgcCalibRawHitPtr, ExtendedInfo >	TgcHitInfo
typedef std::vector< TgcHitInfo >	TgcHitInfoVec
typedef std::map< MuonFixedId, TgcHitInfo >	TgcHitChamberIdMap
typedef std::map< MuonFixedIdManipulator::StationIndex, TgcHitInfoVec >	TgcHitStationIndexMap
typedef std::pair< CscCalibRawHitPtr, ExtendedInfo >	CscHitInfo
typedef std::vector< CscHitInfo >	CscHitInfoVec
typedef std::map< MuonFixedId, CscHitInfo >	CscHitChamberIdMap
typedef std::map< MuonFixedIdManipulator::StationIndex, CscHitInfoVec >	CscHitStationIndexMap

Public Member Functions
	MuonCalibExtendedRawHits ()
	Default constructor .
	MuonCalibExtendedRawHits (const MuonCalibRawHitCollection &rawHits, const std::vector< ExtendedSegPtr > &segments, const std::vector< ExtendedTrkPtr > &tracks)
std::string	dumpMdtHits () const
const MdtHitChamberIdMap &	mdtHitChamberIdMap () const
	access to the mdt hits per chamber and station
const MdtHitStationIndexMap &	mdtChamberIndexMap () const
const RpcHitChamberIdMap &	rpcHitChamberIdMap () const
	access to the rpc hits per chamber and station
const RpcHitStationIndexMap &	rpcChamberIndexMap () const
const TgcHitChamberIdMap &	tgcHitChamberIdMap () const
	access to the tgc hits per chamber and station
const TgcHitStationIndexMap &	tgcChamberIndexMap () const
const CscHitChamberIdMap &	cscHitChamberIdMap () const
	access to the csc hits per chamber and station
const CscHitStationIndexMap &	cscChamberIndexMap () const

Private Attributes
double	m_adcThreshold {50.}
MdtHitChamberIdMap	m_mdtChamberIdMap
MdtHitStationIndexMap	m_mdtChamberIndexMap
RpcHitChamberIdMap	m_rpcChamberIdMap
RpcHitStationIndexMap	m_rpcChamberIndexMap
TgcHitChamberIdMap	m_tgcChamberIdMap
TgcHitStationIndexMap	m_tgcChamberIndexMap
CscHitChamberIdMap	m_cscChamberIdMap
CscHitStationIndexMap	m_cscChamberIndexMap

Detailed Description

Access to all raw hits with layer/station granularity.

Also information on whether the hits are associated with segments or tracks

Definition at line 28 of file MuonCalibExtendedRawHits.h.

Member Typedef Documentation

◆ CscCalibRawHitPtr

using MuonCalib::MuonCalibExtendedRawHits::CscCalibRawHitPtr = MuonCalibRawHitCollection::CscCalibRawHitPtr

Definition at line 31 of file MuonCalibExtendedRawHits.h.

◆ CscHitChamberIdMap

typedef std::map<MuonFixedId, CscHitInfo> MuonCalib::MuonCalibExtendedRawHits::CscHitChamberIdMap

Definition at line 76 of file MuonCalibExtendedRawHits.h.

◆ CscHitInfo

typedef std::pair<CscCalibRawHitPtr, ExtendedInfo> MuonCalib::MuonCalibExtendedRawHits::CscHitInfo

Definition at line 73 of file MuonCalibExtendedRawHits.h.

◆ CscHitInfoVec

typedef std::vector<CscHitInfo> MuonCalib::MuonCalibExtendedRawHits::CscHitInfoVec

Definition at line 75 of file MuonCalibExtendedRawHits.h.

◆ CscHitStationIndexMap

typedef std::map<MuonFixedIdManipulator::StationIndex, CscHitInfoVec> MuonCalib::MuonCalibExtendedRawHits::CscHitStationIndexMap

Definition at line 77 of file MuonCalibExtendedRawHits.h.

◆ ExtendedSegPtr

using MuonCalib::MuonCalibExtendedRawHits::ExtendedSegPtr = std::shared_ptr<MuonCalibExtendedSegment>

Definition at line 36 of file MuonCalibExtendedRawHits.h.

◆ ExtendedTrkPtr

using MuonCalib::MuonCalibExtendedRawHits::ExtendedTrkPtr = std::shared_ptr<MuonCalibExtendedTrack>

Definition at line 37 of file MuonCalibExtendedRawHits.h.

◆ MdtCalibRawHitPtr

using MuonCalib::MuonCalibExtendedRawHits::MdtCalibRawHitPtr = MuonCalibRawHitCollection::MdtCalibRawHitPtr

Definition at line 30 of file MuonCalibExtendedRawHits.h.

◆ MdtHitChamberIdMap

typedef std::map<MuonFixedId, MdtRegionInfo> MuonCalib::MuonCalibExtendedRawHits::MdtHitChamberIdMap

Definition at line 55 of file MuonCalibExtendedRawHits.h.

◆ MdtHitInfoCit

typedef MdtHitInfoMap::const_iterator MuonCalib::MuonCalibExtendedRawHits::MdtHitInfoCit

Definition at line 46 of file MuonCalibExtendedRawHits.h.

◆ MdtHitInfoIt

typedef MdtHitInfoMap::iterator MuonCalib::MuonCalibExtendedRawHits::MdtHitInfoIt

Definition at line 45 of file MuonCalibExtendedRawHits.h.

◆ MdtHitInfoMap

typedef std::map<MuonFixedId, ExtendedChannelInfo> MuonCalib::MuonCalibExtendedRawHits::MdtHitInfoMap

Definition at line 44 of file MuonCalibExtendedRawHits.h.

◆ MdtHitStationIndexMap

typedef std::map<MuonFixedIdManipulator::StationIndex, MdtRegionInfo> MuonCalib::MuonCalibExtendedRawHits::MdtHitStationIndexMap

Definition at line 56 of file MuonCalibExtendedRawHits.h.

◆ RpcCalibRawHitPtr

using MuonCalib::MuonCalibExtendedRawHits::RpcCalibRawHitPtr = MuonCalibRawHitCollection::RpcCalibRawHitPtr

Definition at line 33 of file MuonCalibExtendedRawHits.h.

◆ RpcHitChamberIdMap

typedef std::map<MuonFixedId, RpcHitInfo> MuonCalib::MuonCalibExtendedRawHits::RpcHitChamberIdMap

Definition at line 65 of file MuonCalibExtendedRawHits.h.

◆ RpcHitInfo

typedef std::pair<RpcCalibRawHitPtr, ExtendedInfo> MuonCalib::MuonCalibExtendedRawHits::RpcHitInfo

Definition at line 63 of file MuonCalibExtendedRawHits.h.

◆ RpcHitInfoVec

typedef std::vector<RpcHitInfo> MuonCalib::MuonCalibExtendedRawHits::RpcHitInfoVec

Definition at line 64 of file MuonCalibExtendedRawHits.h.

◆ RpcHitStationIndexMap

typedef std::map<MuonFixedIdManipulator::StationIndex, RpcHitInfoVec> MuonCalib::MuonCalibExtendedRawHits::RpcHitStationIndexMap

Definition at line 66 of file MuonCalibExtendedRawHits.h.

◆ TgcCalibRawHitPtr

using MuonCalib::MuonCalibExtendedRawHits::TgcCalibRawHitPtr = MuonCalibRawHitCollection::TgcCalibRawHitPtr

Definition at line 34 of file MuonCalibExtendedRawHits.h.

◆ TgcHitChamberIdMap

typedef std::map<MuonFixedId, TgcHitInfo> MuonCalib::MuonCalibExtendedRawHits::TgcHitChamberIdMap

Definition at line 70 of file MuonCalibExtendedRawHits.h.

◆ TgcHitInfo

typedef std::pair<TgcCalibRawHitPtr, ExtendedInfo> MuonCalib::MuonCalibExtendedRawHits::TgcHitInfo

Definition at line 68 of file MuonCalibExtendedRawHits.h.

◆ TgcHitInfoVec

typedef std::vector<TgcHitInfo> MuonCalib::MuonCalibExtendedRawHits::TgcHitInfoVec

Definition at line 69 of file MuonCalibExtendedRawHits.h.

◆ TgcHitStationIndexMap

typedef std::map<MuonFixedIdManipulator::StationIndex, TgcHitInfoVec> MuonCalib::MuonCalibExtendedRawHits::TgcHitStationIndexMap

Definition at line 71 of file MuonCalibExtendedRawHits.h.

Constructor & Destructor Documentation

◆ MuonCalibExtendedRawHits() [1/2]

MuonCalib::MuonCalibExtendedRawHits::MuonCalibExtendedRawHits ( )

Default constructor .

Definition at line 24 of file MuonCalibExtendedRawHits.cxx.

24: m_adcThreshold(50.) {}

MuonCalib::MuonCalibExtendedRawHits::m_adcThreshold

double m_adcThreshold

Definition MuonCalibExtendedRawHits.h:105

◆ MuonCalibExtendedRawHits() [2/2]

MuonCalib::MuonCalibExtendedRawHits::MuonCalibExtendedRawHits	(	const MuonCalibRawHitCollection &	rawHits,
		const std::vector< ExtendedSegPtr > &	segments,
		const std::vector< ExtendedTrkPtr > &	tracks )

Definition at line 26 of file MuonCalibExtendedRawHits.cxx.

                                                                                                :
        m_adcThreshold(70.) {
        MuonFixedIdManipulator idManip;
        // loop over MDT raw hits and
        for (const MdtCalibRawHitPtr& rawHit : rawHits.MdtContainer()) {
            const MuonFixedId id = rawHit->identify();
            // create ChamberId + StationIndex
            MuonFixedId chId = idManip.chamberIdentifier(id);
            MuonFixedIdManipulator::StationIndex stationIndex = idManip.stationLayerIndex(id);
 
            MdtRegionInfo& chInfo = m_mdtChamberIdMap[chId];
            ExtendedChannelInfo& channelChInfo = chInfo.hits[id];
            MdtRegionInfo& stInfo = m_mdtChamberIndexMap[stationIndex];
            ExtendedChannelInfo& channelStInfo = stInfo.hits[id];
 
            channelChInfo.hits.emplace_back(rawHit);
            channelStInfo.hits.emplace_back(rawHit);
 
            const bool above_thresh = {rawHit->adcCount() > m_adcThreshold};
            chInfo.nhitsAboveAdcThreshold += above_thresh;
            stInfo.nhitsAboveAdcThreshold += above_thresh;
 
            // only increase the hit counters once
            bool isFirst = true;
            // loop over tracks and check whether hit is 'on-track'
            for (const ExtendedTrkPtr& trk : tracks) {
                // look for id in map
                const IdHitMap& idHitMap = trk->idHitMap();
                IdHitMap::const_iterator pos = idHitMap.find(id);
                if (pos == idHitMap.end()) continue;
 
                // if we get here the hit was found, now check that it is a hit on track
                const MuonCalibHit_E& hit = *pos->second;
                if (hit.type() < 4) {
                    if (isFirst) {
                        ++chInfo.nhitsOnTrack;
                        ++stInfo.nhitsOnTrack;
                        isFirst = false;
                    }
                    trk->rawHitAssociationMap().addMdtRawHit(rawHit);
                    channelChInfo.associatedTracks.emplace_back(trk);
                    channelStInfo.associatedTracks.emplace_back(trk);
                }
            }
 
            // only increase the hit counters once
            isFirst = true;
            // loop over segments and check whether hit is 'on-segment'
            for (const ExtendedSegPtr& seg : segments) {
                // look for id in segment
                if (!seg->idSet().count(id)) continue;
 
                if (isFirst) {
                    ++chInfo.nhitsOnSegment;
                    ++stInfo.nhitsOnSegment;
                    isFirst = false;
                }
                seg->rawHitAssociationMap().addMdtRawHit(rawHit);
                channelChInfo.associatedSegments.emplace_back(seg);
                channelStInfo.associatedSegments.emplace_back(seg);
            }
        }
    }

Member Function Documentation

◆ cscChamberIndexMap()

const CscHitStationIndexMap & MuonCalib::MuonCalibExtendedRawHits::cscChamberIndexMap ( ) const

inline

Definition at line 102 of file MuonCalibExtendedRawHits.h.

102{ return m_cscChamberIndexMap; }

MuonCalib::MuonCalibExtendedRawHits::m_cscChamberIndexMap

CscHitStationIndexMap m_cscChamberIndexMap

Definition MuonCalibExtendedRawHits.h:117

◆ cscHitChamberIdMap()

const CscHitChamberIdMap & MuonCalib::MuonCalibExtendedRawHits::cscHitChamberIdMap ( ) const

inline

access to the csc hits per chamber and station

Definition at line 101 of file MuonCalibExtendedRawHits.h.

101{ return m_cscChamberIdMap; }

MuonCalib::MuonCalibExtendedRawHits::m_cscChamberIdMap

CscHitChamberIdMap m_cscChamberIdMap

Definition MuonCalibExtendedRawHits.h:116

◆ dumpMdtHits()

std::string MuonCalib::MuonCalibExtendedRawHits::dumpMdtHits ( ) const

Definition at line 92 of file MuonCalibExtendedRawHits.cxx.

                                                          {
        std::ostringstream sout;
        MuonFixedIdPrinter printer{};
 
        sout << " number of MDT chambers with hits: " << m_mdtChamberIdMap.size() << std::endl;
        for (const auto& chIt : m_mdtChamberIdMap) {
            if (chIt.second.nhitsAboveAdcThreshold < 3) continue;
            const MuonFixedId& chId = chIt.first;
            sout << "    " << std::setw(20) << printer.chamberIdentifier(chId) << "    nhits " << std::setw(3) << chIt.second.hits.size()
                 << "    above ADC cut " << std::setw(3) << chIt.second.nhitsAboveAdcThreshold << "    on segment " << std::setw(3)
                 << chIt.second.nhitsOnSegment << "    on track " << std::setw(3) << chIt.second.nhitsOnTrack;
 
            // look at next iterator to decide whether to add endl or not
            sout << std::endl;
        }
        return sout.str();
    }