Muon::MuonEDMPrinterTool Class Reference

Helper tool to print EDM objects to string in a fix format. More...

#include <MuonEDMPrinterTool.h>

Inheritance diagram for Muon::MuonEDMPrinterTool:

Collaboration diagram for Muon::MuonEDMPrinterTool:

Public Member Functions
	MuonEDMPrinterTool (const std::string &, const std::string &, const IInterface *)
	constructor
	~MuonEDMPrinterTool ()=default
	destructor
StatusCode	initialize ()
	AlgTool initilize.
std::string	print (const Trk::Track &track) const
	print track parameters to string
std::string	printStations (const Trk::Track &track) const
	print stations on track to string
std::string	print (const Trk::MuonTrackSummary &summary) const
	print stations on track to string
std::string	print (const MuonSegment &segment) const
	print segment parameters to string
std::string	print (const std::vector< const MuonSegment * > &segs) const
	print vector of segments
std::string	print (std::vector< std::unique_ptr< MuonSegment > > &segs) const
std::string	print (const Trk::MeasurementBase &measurement) const
	print Muon MeasurementBase to string
std::string	print (const Trk::PrepRawData &prd) const
	print Muon PrepRawData to string
std::string	print (const std::vector< const Trk::MeasurementBase * > &measurements) const
	print vector of measurement base to string
std::string	print (const MuonSegmentCombinationCollection &combiCol) const
	print MuonSegmentCombinationCollection
std::string	print (const MuonSegmentCombination &combi) const
	print MuonSegmentCombination
std::string	print (const MuonPattern &pattern) const
	print MuonPattern
std::string	print (const MuonPatternCollection &patCol) const
	print MuonPatternCollection
std::string	print (const MuonPrdPatternCollection &patCol) const
	print MuonPatternCollection
std::string	print (const MuonPatternCombination &pattern) const
	print MuonPatternCombination
std::string	print (const MuonPatternCombinationCollection &combiCol) const
	print MuonPatternCombinationCollection
std::string	print (const MuonPatternChamberIntersect &intersect) const
	print MuonPatternChamberIntersect
std::string	print (const Trk::TrackStateOnSurface &tsos) const
	print the track state on surface to the string
std::string	printMeasurements (const Trk::Track &track) const
	print measurements on track to string
std::string	printData (const Trk::MeasurementBase &measurement) const
	print data part of Muon MeasurementBase to string
std::string	printId (const Trk::MeasurementBase &measurement) const
	print identifier part of Muon MeasurementBase to string
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()
	access to tool interface
static std::string	printPatRec (const Trk::Track &track)
	produce a string indicating who produced the track
static std::string	print (const Trk::TrackParameters &pars)
	print Trk::TrackParameters
static std::string	print (const Trk::ResidualPull &resPull)
	print ResidualPull object to string
static std::string	print (const Trk::MaterialEffectsBase &mat)
	print the material effects object to the string
static std::string	print (const Trk::AlignmentEffectsOnTrack &aeot)
	print the alignment effects on track object to the string

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 ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc
ServiceHandle< IMuonEDMHelperSvc >	m_edmHelperSvc
ToolHandle< Trk::IResidualPullCalculator >	m_pullCalculator
SG::ReadCondHandleKey< MuonGM::MuonDetectorManager >	m_DetectorManagerKey
SG::ReadHandleKey< MdtPrepDataContainer >	m_mdtKey
SG::ReadHandleKey< RpcPrepDataContainer >	m_rpcKey
SG::ReadHandleKey< TgcPrepDataContainer >	m_tgcKey
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

Helper tool to print EDM objects to string in a fix format.

Definition at line 43 of file MuonEDMPrinterTool.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

MuonEDMPrinterTool()

Muon::MuonEDMPrinterTool::MuonEDMPrinterTool	(	const std::string &	ty,
		const std::string &	na,
		const IInterface *	pa )

constructor

Definition at line 35 of file MuonEDMPrinterTool.cxx.

    : AthAlgTool(ty, na, pa)
{
    declareInterface<MuonEDMPrinterTool>(this);
}

~MuonEDMPrinterTool()

Muon::MuonEDMPrinterTool::~MuonEDMPrinterTool ( )

default

destructor

Member Function Documentation

declareGaudiProperty()

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

inlineprivateinherited

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

Definition at line 156 of file AthCommonDataStore.h.

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

declareProperty()

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

evtStore()

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

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

extraDeps_update_handler()

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

protectedinherited

Add StoreName to extra input/output deps as needed.

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

initialize()

StatusCode Muon::MuonEDMPrinterTool::initialize

(

)

AlgTool initilize.

Definition at line 41 of file MuonEDMPrinterTool.cxx.

                                         {
    if (parent() != toolSvc()) {
        ATH_MSG_FATAL("There's no need to make this tool private");
        return StatusCode::FAILURE;
    }
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK(m_edmHelperSvc.retrieve());
    ATH_CHECK(m_pullCalculator.retrieve());
 
    ATH_CHECK(m_DetectorManagerKey.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()));
 
 
    return StatusCode::SUCCESS;
}

inputHandles()

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

overridevirtualinherited

Return this algorithm's input handles.

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

interfaceID()

const InterfaceID & Muon::MuonEDMPrinterTool::interfaceID ( )

inlinestatic

access to tool interface

Definition at line 55 of file MuonEDMPrinterTool.h.

    {
        static const InterfaceID IID_MuonEDMPrinterTool("Muon::MuonEDMPrinterTool", 1, 0);
        return IID_MuonEDMPrinterTool;
    }

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

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

overridevirtualinherited

Return this algorithm's output handles.

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

print() [1/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonPattern &

pattern

)

const

print MuonPattern

Definition at line 585 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    unsigned int nhits = pattern.numberOfContainedPrds();
    sout << "MuonPattern with " << nhits << " hits" << std::setprecision(5) << " theta "
         << pattern.globalDirection().theta() << " phi " << pattern.globalDirection().phi() << std::setprecision(1)
         << " r " << std::setw(5) << (int)pattern.globalPosition().perp() << "  z " << std::setw(5)
         << (int)pattern.globalPosition().z() << std::endl;
    for (unsigned int i = 0; i < nhits; ++i) {
        const Trk::PrepRawData* prd = pattern.prd(i);
        if (!prd) {
            sout << " ERROR found zero pointer on MuonPattern " << std::endl;
            continue;
        }
        sout << "  " << print(*prd);
        if (i < nhits - 1) sout << std::endl;
    }
    return sout.str();
}

print() [2/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonPatternChamberIntersect &

intersect

)

const

print MuonPatternChamberIntersect

&& !detEls.count(detEl)

&& !detEls.count(detEl2)

Definition at line 696 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    std::string chIdString("Unknown id");
    double      chTheta(-99.);
    double      chPhi(-99.);
    // if the intersect is not empty use the Identifier of first hit as chamber ID
    const Trk::PrepRawData* firstPrd = nullptr;
    if (!intersect.prepRawDataVec().empty() && intersect.prepRawDataVec().front()) {
        firstPrd   = intersect.prepRawDataVec().front();
        chIdString = m_idHelperSvc->toStringChamber(firstPrd->identify());
        chTheta    = firstPrd->detectorElement()->center().theta();
        chPhi      = firstPrd->detectorElement()->center().phi();
    } else {
        return chIdString;
    }
    Identifier chId  = m_idHelperSvc->chamberId(firstPrd->identify());
    int        neta  = 0;
    int        nphi  = 0;
    bool       isMdt = false;
 
    unsigned int                                         nchannelsEta = 0;
    unsigned int                                         nchannelsPhi = 0;
    std::set<const Trk::TrkDetElementBase*>              detEls;
    std::vector<const Trk::PrepRawData*>::const_iterator hit     = intersect.prepRawDataVec().begin();
    std::vector<const Trk::PrepRawData*>::const_iterator hit_end = intersect.prepRawDataVec().end();
    for (; hit != hit_end; ++hit) {
        const Trk::PrepRawData* prd = *hit;
        if (!prd) {
            sout << " ERROR found zero pointer on MuonPatternChamberIntersect " << std::endl;
            continue;
        }
        const Identifier& id = prd->identify();
        if (m_idHelperSvc->measuresPhi(id))
            ++nphi;
        else
            ++neta;
        if (!isMdt && m_idHelperSvc->isMdt(id)) isMdt = true;
 
        if (!detEls.count(prd->detectorElement())) {
            detEls.insert(prd->detectorElement());
 
            if (isMdt && detEls.empty()) {
 
                SG::ReadCondHandle<MuonGM::MuonDetectorManager> DetectorManagerHandle{m_DetectorManagerKey};
                const MuonGM::MuonDetectorManager*              MuonDetMgr{*DetectorManagerHandle};
                if (MuonDetMgr == nullptr) {
                    ATH_MSG_DEBUG("Cannot retrieve DetectorManager ");
                } else {
 
                    Identifier                       idml1  = m_idHelperSvc->mdtIdHelper().channelID(id, 1, 1, 1);
                    Identifier                       idml2  = m_idHelperSvc->mdtIdHelper().channelID(id, 2, 1, 1);
                    const MuonGM::MdtReadoutElement* detEl1 = MuonDetMgr->getMdtReadoutElement(idml1);
                    const MuonGM::MdtReadoutElement* detEl2 = nullptr;
                    if (m_idHelperSvc->mdtIdHelper().numberOfMultilayers(id) == 2) {
                        detEl2 = MuonDetMgr->getMdtReadoutElement(idml2);
                    } else {
                        ATH_MSG_DEBUG("A single multilayer for this station " << m_idHelperSvc->toString(id));
                    }
                    if (detEl1 ) {
                        detEls.insert(detEl1);
                        nchannelsEta += detEl1->getNLayers() * detEl1->getNtubesperlayer();
                    } else {
                        ATH_MSG_DEBUG("Couldn't find first detector element ");
                    }
                    if (detEl2 ) {
                        detEls.insert(detEl2);
                        nchannelsEta += detEl2->getNLayers() * detEl2->getNtubesperlayer();
                    }
                }
            } else if (m_idHelperSvc->isTgc(id)) {
                const MuonGM::TgcReadoutElement* detEl =
                    dynamic_cast<const MuonGM::TgcReadoutElement*>(prd->detectorElement());
                if (detEl) {
                    for (int i = 1; i <= detEl->nGasGaps(); ++i) {
                        nchannelsEta += detEl->nWireGangs(i);
                        nchannelsPhi += detEl->nStrips(i);
                    }
                }
            } else if (m_idHelperSvc->isRpc(id)) {
                const MuonGM::RpcReadoutElement* detEl =
                    dynamic_cast<const MuonGM::RpcReadoutElement*>(prd->detectorElement());
                if (detEl) {
                    nchannelsPhi += detEl->nGasGapPerLay()  * detEl->NphiStripPanels() * detEl->NphiStrips();
                    nchannelsEta += detEl->NetaStrips();
                }
            }
        }
    }
 
    unsigned int nchHitsEta = 0;
    unsigned int nchHitsPhi = 0;
 
 
    if (isMdt && !m_mdtKey.empty()) {
        SG::ReadHandle mdtPrdContainer{m_mdtKey};
        const MdtPrepDataCollection *coll = mdtPrdContainer->indexFindPtr(m_idHelperSvc->moduleHash(chId));
        if (coll != nullptr) {
            nchHitsEta = coll->size();
        }
    } else if (m_idHelperSvc->isRpc(chId) && !m_rpcKey.empty()) {
        SG::ReadHandle rpcPrdContainer{m_rpcKey};
        const RpcPrepDataCollection* coll = rpcPrdContainer->indexFindPtr(m_idHelperSvc->moduleHash(chId));
        if (coll != nullptr) {
            for (const RpcPrepData* prd : *coll) {
                if (m_idHelperSvc->measuresPhi(prd->identify()))
                    ++nchHitsPhi;
                else
                    ++nchHitsEta;
            }
        } 
    } else if (m_idHelperSvc->isTgc(chId) && !m_tgcKey.empty()) {
        SG::ReadHandle tgcPrdContainer{m_tgcKey};
        const TgcPrepDataCollection *coll = tgcPrdContainer->indexFindPtr(m_idHelperSvc->moduleHash(chId));
         if (coll != nullptr) {
            for (const TgcPrepData* prd : *coll) {
                if (m_idHelperSvc->measuresPhi(prd->identify()))
                    ++nchHitsPhi;
                else
                    ++nchHitsEta;
            }
        }
    }
 
    sout.setf(std::ios::left);
    sout << std::setw(33) << chIdString << " neta " << std::setw(3) << neta;
    if (!isMdt)
        sout << " nphi " << std::setw(3) << nphi;
    else
        sout << "         ";
 
    double scaleEta = nchannelsEta != 0 ? 1. / nchannelsEta : 0.;
    double scalePhi = nchannelsPhi != 0 ? 1. / nchannelsPhi : 0.;
    sout << " occupancy eta " << std::fixed << std::setprecision(2) << std::setw(5) << scaleEta * nchHitsEta;
    if (!isMdt)
        sout << " phi " << std::fixed << std::setw(5) << scalePhi * nchHitsPhi;
    else
        sout << "          ";
 
    sout << std::setprecision(2) << "  pat theta " << std::setw(5) << intersect.intersectDirection().theta()
         << " ch theta " << std::setw(5) << chTheta << " ch phi " << std::setw(6) << chPhi << std::setprecision(1)
         << " r " << std::setw(5) << (int)intersect.intersectPosition().perp() << "  z " << std::setw(5)
         << (int)intersect.intersectPosition().z();
 
    return sout.str();
}

print() [3/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonPatternCollection &

patCol

)

const

print MuonPatternCollection

Definition at line 607 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    sout << "MuonPatternCollection with " << patCol.size() << " patterns " << std::endl;
    MuonPatternCollection::const_iterator it      = patCol.begin();
    MuonPatternCollection::const_iterator it_end  = patCol.end();
    MuonPatternCollection::const_iterator it_last = it_end;
    --it_last;
    for (; it != it_end; ++it) {
        if (!*it) {
            sout << " WARNING, zero pointer detected in MuonPatternCollection!!! " << std::endl;
            continue;
        }
        sout << " " << print(**it);
        if (it != it_last) sout << std::endl;
    }
    return sout.str();
}

print() [4/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonPatternCombination &

pattern

)

const

print MuonPatternCombination

Definition at line 670 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    unsigned int nchambers = combi.chamberData().size();
    sout << "MuonPatternCombination with " << nchambers << " chambers ";
    const Trk::TrackParameters* pars = combi.trackParameter();
    if (pars) {
        sout << " " << print(*pars);
    } else {
        sout << "    -> combi has no parameters ";
    }
    if (!combi.chamberData().empty()) sout << std::endl;
    std::vector<MuonPatternChamberIntersect>::const_iterator chit      = combi.chamberData().begin();
    std::vector<MuonPatternChamberIntersect>::const_iterator chit_end  = combi.chamberData().end();
    std::vector<MuonPatternChamberIntersect>::const_iterator chit_last = chit_end;
    --chit_last;
    for (; chit != chit_end; ++chit) {
        const MuonPatternChamberIntersect& chamberIntersect = *chit;
        sout << print(chamberIntersect);
        if (chit != chit_last) sout << std::endl;
    }
    return sout.str();
}

print() [5/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonPatternCombinationCollection &

combiCol

)

const

print MuonPatternCombinationCollection

Definition at line 649 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    sout << "MuonPatternCombinationCollection with combis " << combiCol.size() << std::endl;
    MuonPatternCombinationCollection::const_iterator it      = combiCol.begin();
    MuonPatternCombinationCollection::const_iterator it_end  = combiCol.end();
    MuonPatternCombinationCollection::const_iterator it_last = it_end;
    --it_last;
    for (; it != it_end; ++it) {
        if (!*it) {
            sout << " WARNING, zero pointer detected in MuonPatternCombinationCollection!!! " << std::endl;
            continue;
        }
        sout << print(**it);
        if (it != it_last) sout << std::endl;
    }
    return sout.str();
}

print() [6/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonPrdPatternCollection &

patCol

)

const

print MuonPatternCollection

Definition at line 628 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    sout << "MuonPrdPatternCollection with " << patCol.size() << " patterns " << std::endl;
    MuonPrdPatternCollection::const_iterator it      = patCol.begin();
    MuonPrdPatternCollection::const_iterator it_end  = patCol.end();
    MuonPrdPatternCollection::const_iterator it_last = it_end;
    --it_last;
    for (; it != it_end; ++it) {
        if (!*it) {
            sout << " WARNING, zero pointer detected in MuonPatternCollection!!! " << std::endl;
            continue;
        }
        sout << " " << print(**it);
        if (it != it_last) sout << std::endl;
    }
    return sout.str();
}

print() [7/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonSegment &

segment

)

const

print segment parameters to string

Definition at line 315 of file MuonEDMPrinterTool.cxx.

{
 
    std::ostringstream sout;
 
    // get first none-trigger id
    Identifier chid     = m_edmHelperSvc->chamberId(segment);
    int        nphi     = 0;
    int        ntrigEta = 0;
    int        neta     = 0;
 
    const MuonGM::MdtReadoutElement* mdtDetEl = nullptr;
    Identifier                       shortestTubeId;
    double                           shortestTubeLen = 1e9;
 
    std::vector<const Trk::MeasurementBase*>::const_iterator hit     = segment.containedMeasurements().begin();
    std::vector<const Trk::MeasurementBase*>::const_iterator hit_end = segment.containedMeasurements().end();
    for (; hit != hit_end; ++hit) {
        Identifier              id;
        const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(*hit);
        if (rot)
            id = rot->identify();
        else {
            const CompetingMuonClustersOnTrack* crot = dynamic_cast<const CompetingMuonClustersOnTrack*>(*hit);
            if (crot) id = crot->containedROTs().front()->identify();
        }
        if (!id.is_valid()) continue;
        bool measuresPhi = m_idHelperSvc->measuresPhi(id);
        bool isTrigger   = m_idHelperSvc->isTrigger(id);
        if (measuresPhi)
            ++nphi;
        else if (!isTrigger)
            ++neta;
 
        if (!isTrigger) {
            const MdtDriftCircleOnTrack* mdt = dynamic_cast<const MdtDriftCircleOnTrack*>(rot);
            if (mdt && mdt->prepRawData()) {
                double tubelen = 0.5 * mdt->prepRawData()->detectorElement()->tubeLength(id);
                if (tubelen < shortestTubeLen) {
                    shortestTubeId  = id;
                    shortestTubeLen = tubelen;
                    mdtDetEl        = mdt->prepRawData()->detectorElement();
                }
            }
        } else {
            if (!measuresPhi) ++ntrigEta;
        }
    }
 
    sout << m_idHelperSvc->toStringChamber(chid);
 
    const Trk::FitQuality* fq = segment.fitQuality();
    if (fq) {
        sout << std::setprecision(2) << "   chi2 " << std::setw(7) << fq->chiSquared();
    } else {
        sout << " no fit Quality ";
    }
 
    sout << " neta " << std::setw(2) << neta << " nphi " << std::setw(2) << nphi << " nTrigEta " << std::setw(2)
         << ntrigEta;
 
    const MuonSegmentQuality* q = dynamic_cast<const MuonSegmentQuality*>(fq);
    if (q) sout << std::setw(2) << " nholes " << q->numberOfHoles();
 
 
    sout << " theta " << std::fixed << std::setprecision(5) << std::setw(7) << segment.globalDirection().theta()
         << " phi " << std::fixed << std::setprecision(3) << std::setw(6) << segment.globalDirection().phi();
 
    if (segment.hasFittedT0())
        sout << " T0 " << std::fixed << std::setprecision(2) << std::setw(5) << segment.time() << " err "
             << std::setw(5) << segment.errorTime();
 
    if (mdtDetEl) {
        double posAlongTube = std::abs((mdtDetEl->globalToLocalTransf(shortestTubeId)*segment.globalPosition()).z());
        double distFromEdge = posAlongTube - shortestTubeLen;
        if (distFromEdge < -100.)
            sout << " inside chamber";
        else if (distFromEdge < -1.)
            sout << " close to edge";
        else if (distFromEdge < 0.0001)
            sout << " on edge";
        else
            sout << " outside chamber";
    }
    return sout.str();
}

print() [8/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonSegmentCombination &

combi

)

const

print MuonSegmentCombination

Definition at line 559 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    unsigned int nstations = combi.numberOfStations();
    sout << "SegmentCombination with stations " << nstations << std::endl;
 
    // loop over chambers in combi
    for (unsigned int i = 0; i < nstations; ++i) {
 
        // loop over segments in station
        MuonSegmentCombination::SegmentVec* stationSegs = combi.stationSegments(i);
 
        // check if not empty
        if (!stationSegs || stationSegs->empty()) continue;
 
        // get chamber identifier, chamber index and station index
        // Identifier chid = m_edmHelperSvc->chamberId( *stationSegs->front() );
        sout << print(*stationSegs);
        if (i != nstations - 1) sout << std::endl;
    }
 
    return sout.str();
}

print() [9/21]

std::string Muon::MuonEDMPrinterTool::print

(

const MuonSegmentCombinationCollection &

combiCol

)

const

print MuonSegmentCombinationCollection

Definition at line 538 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    sout << "MuonSegmentCombinationCollection with combis " << combiCol.size() << std::endl;
    MuonSegmentCombinationCollection::const_iterator it      = combiCol.begin();
    MuonSegmentCombinationCollection::const_iterator it_end  = combiCol.end();
    MuonSegmentCombinationCollection::const_iterator it_last = it_end;
    --it_last;
    for (; it != it_end; ++it) {
        if (!*it) {
            sout << " WARNING, zero pointer detected in MuonSegmentCombinationCollection!!! " << std::endl;
            continue;
        }
        sout << print(**it);
        if (it != it_last) sout << std::endl;
    }
    return sout.str();
}

print() [10/21]

std::string Muon::MuonEDMPrinterTool::print

(

const std::vector< const MuonSegment * > &

segs

)

const

print vector of segments

Definition at line 494 of file MuonEDMPrinterTool.cxx.

{
 
    std::ostringstream sout;
 
    std::vector<const MuonSegment*>::const_iterator it      = segs.begin();
    std::vector<const MuonSegment*>::const_iterator it_end  = segs.end();
    std::vector<const MuonSegment*>::const_iterator it_last = it_end;
    --it_last;
    for (; it != it_end; ++it) {
        if (!*it) {
            sout << " WARNING, zero pointer detected in MuonSegment vector!!! " << std::endl;
            continue;
        }
        sout << print(**it);
        if (it != it_last) sout << std::endl;
    }
 
    return sout.str();
}

print() [11/21]

std::string Muon::MuonEDMPrinterTool::print

(

const std::vector< const Trk::MeasurementBase * > &

measurements

)

const

print vector of measurement base to string

Definition at line 455 of file MuonEDMPrinterTool.cxx.

{
 
    std::ostringstream sout;
 
    // first loop to get width of Id's for nice alignment
    std::vector<std::string> idStrings;
    std::vector<std::string> dataStrings;
    unsigned int             idWidth = 0;
 
    std::vector<const Trk::MeasurementBase*>::const_iterator hit     = measurements.begin();
    std::vector<const Trk::MeasurementBase*>::const_iterator hit_end = measurements.end();
    // std::vector< const Trk::MeasurementBase* >::const_iterator hit_last = hit_end;
    for (; hit != hit_end; ++hit) {
        if (!*hit) {
            sout << " WARNING, zero pointer detected in MeasurementBase vector!!! " << std::endl;
            continue;
        }
        // Identifier part
        std::string idStr = printId(**hit);
        idStrings.push_back(idStr);
        if (idStr.length() > idWidth) idWidth = idStr.length();
        // Data part
        std::string dataStr = printData(**hit);
        dataStrings.push_back(dataStr);
    }
 
 
    // second loop to print out aligned strings
    unsigned int n = idStrings.size();
    for (unsigned int i = 0; i < n; ++i) {
        sout << std::left << std::setw(idWidth) << idStrings[i] << std::right << "  " << dataStrings[i];
        if (i != n - 1) sout << std::endl;
    }
 
    return sout.str();
}

print() [12/21]

std::string Muon::MuonEDMPrinterTool::print ( const Trk::AlignmentEffectsOnTrack & aeot )

static

print the alignment effects on track object to the string

Definition at line 995 of file MuonEDMPrinterTool.cxx.

                                                                          {
    std::stringstream aeot_string;
    aeot_string<<"dTrans: "<<aeot.deltaTranslation()<<" +- "<<aeot.sigmaDeltaTranslation();
    aeot_string<<" dAngle: "<<aeot.deltaAngle()<<" +- "<<aeot.sigmaDeltaAngle();
    return aeot_string.str();
}

print() [13/21]

std::string Muon::MuonEDMPrinterTool::print ( const Trk::MaterialEffectsBase & mat )

static

print the material effects object to the string

Definition at line 969 of file MuonEDMPrinterTool.cxx.

                                                                     {
    std::stringstream mat_string{};
    const Trk::MaterialEffectsOnTrack* matOnTrk = dynamic_cast<const Trk::MaterialEffectsOnTrack*>(&mat);
    if (matOnTrk) {
        const Trk::ScatteringAngles* scatAng = matOnTrk->scatteringAngles();
        if (scatAng) {
            mat_string<<" dPhi: "<<scatAng->deltaPhi()<<" +- "<<scatAng->sigmaDeltaPhi()<< " / ";
            mat_string<<" dTheta: "<<scatAng->deltaTheta()<<" +- "<<scatAng->sigmaDeltaTheta()<<"   ";
        }
        const Trk::EnergyLoss* eloss = matOnTrk->energyLoss();
        if (eloss) {
            mat_string<<" deltaE: "<<eloss->deltaE();
            if (std::abs(std::abs(eloss->sigmaMinusDeltaE()) - std::abs(eloss->sigmaPlusDeltaE())) > std::numeric_limits<float>::epsilon()) {
                mat_string<<" +/- "<<eloss->sigmaPlusDeltaE()<<"/"<<eloss->sigmaMinusDeltaE();
            } else {
                mat_string<<" +- "<<eloss->sigmaDeltaE();
            }
            mat_string<<" Ion: "<<eloss->meanIoni()<<" +- "<<eloss->sigmaIoni();
            mat_string<<" Rad: "<<eloss->meanRad()<<"  +- "<<eloss->sigmaRad();
            mat_string<<" Length "<<eloss->length();
        }
        mat_string<<" ";
 
    }
    return mat_string.str();
}

print() [14/21]

std::string Muon::MuonEDMPrinterTool::print

(

const Trk::MeasurementBase &

measurement

)

const

print Muon MeasurementBase to string

Definition at line 448 of file MuonEDMPrinterTool.cxx.

{
    return printId(measurement) + "  " + printData(measurement);
}

print() [15/21]

std::string Muon::MuonEDMPrinterTool::print

(

const Trk::MuonTrackSummary &

summary

)

const

print stations on track to string

Definition at line 60 of file MuonEDMPrinterTool.cxx.

{
 
    std::ostringstream sout;
    // create detailed track summary
    typedef Trk::MuonTrackSummary::ChamberHitSummary ChamberHitSummary;
 
    sout << "Hits: eta " << summary.netaHits() << "  phi " << summary.nphiHits() << "  holes " << summary.nholes()
         << "  outliers " << summary.noutliers() << "  pseudo " << summary.npseudoMeasurements() << "  scatterers "
         << summary.nscatterers() << " close Hits " << summary.ncloseHits() << std::endl;
 
 
    const std::vector<ChamberHitSummary>& chamberHitSummary = summary.chamberHitSummary();
    sout.setf(std::ios::left);
    std::vector<ChamberHitSummary>::const_iterator chit      = chamberHitSummary.begin();
    std::vector<ChamberHitSummary>::const_iterator chit_end  = chamberHitSummary.end();
    std::vector<ChamberHitSummary>::const_iterator chit_last = chit_end - 1;
    for (; chit != chit_end; ++chit) {
        const Identifier& chId  = chit->chamberId();
        bool              isMdt = m_idHelperSvc->isMdt(chId);
 
        sout << "  " << std::setw(35) << m_idHelperSvc->toStringChamber(chId);
 
        const ChamberHitSummary::Projection& first  = isMdt ? chit->mdtMl1() : chit->etaProjection();
        const ChamberHitSummary::Projection& second = isMdt ? chit->mdtMl2() : chit->phiProjection();
 
 
        std::string firstString  = isMdt ? "ml1 " : "eta ";
        std::string secondString = isMdt ? "ml2 " : "phi ";
 
        sout << " Hits:  " << firstString << std::setw(3) << first.nhits << "   " << secondString << std::setw(3)
             << second.nhits;
 
        if (first.nholes || second.nholes) {
            sout << "  Holes:  ";
            if (first.nholes != 0) sout << firstString << std::setw(3) << first.nholes;
            if (second.nholes != 0) {
                if (first.nholes != 0) sout << "  ";
                sout << secondString << std::setw(3) << second.nholes;
            }
        }
        if (first.noutliers || second.noutliers) {
            sout << "  Outliers:  ";
            if (first.noutliers != 0) sout << firstString << std::setw(3) << first.noutliers << "  ";
            if (second.noutliers != 0) {
                if (first.noutliers != 0) sout << "  ";
                sout << secondString << std::setw(3) << second.noutliers;
            }
        }
 
        if (first.ndeltas || second.ndeltas) {
            sout << "  Deltas:  ";
            if (first.ndeltas != 0) sout << firstString << std::setw(3) << first.ndeltas << "  ";
            if (second.ndeltas != 0) {
                if (first.ndeltas != 0) sout << "  ";
                sout << secondString << std::setw(3) << second.ndeltas;
            }
        }
 
        if (first.ncloseHits || second.ncloseHits) {
            sout << "  Close Hits:  ";
            if (first.ncloseHits != 0) sout << firstString << std::setw(3) << first.ncloseHits << "  ";
            if (second.ncloseHits != 0) {
                if (first.ncloseHits != 0) sout << "  ";
                sout << secondString << std::setw(3) << second.ncloseHits;
            }
        }
 
        if (chit != chit_last) sout << std::endl;
    }
 
    return sout.str();
}

print() [16/21]

std::string Muon::MuonEDMPrinterTool::print

(

const Trk::PrepRawData &

prd

)

const

print Muon PrepRawData to string

Definition at line 403 of file MuonEDMPrinterTool.cxx.

{
 
    std::ostringstream sout;
 
    Identifier id = prd.identify();
    sout << m_idHelperSvc->toString(id) << "  ";
 
    const Amg::Vector3D* pos = nullptr;
    const MuonCluster*   cl  = dynamic_cast<const MuonCluster*>(&prd);
    if (cl)
        pos = &cl->globalPosition();
    else {
        const MdtPrepData* mdt = dynamic_cast<const MdtPrepData*>(&prd);
        if (mdt)
            pos = &mdt->detectorElement()->surface(id).center();
        else {
            const CscStripPrepData* strip = dynamic_cast<const CscStripPrepData*>(&prd);
            if (strip) pos = &strip->globalPosition();
        }
    }
    if (!pos)
        sout << " unknown type of muon prepdata " << std::endl;
    else {
        double h_r     = pos->perp();
        double h_z     = pos->z();
        double h_phi   = pos->phi();
        double h_theta = pos->theta();
        double error = std::sqrt(prd.localCovariance()(0, 0));
 
        // add time for RPC
        double             rpcTime = 0.0;
        const RpcPrepData* rpc     = dynamic_cast<const RpcPrepData*>(&prd);
        if (rpc) rpcTime = rpc->time();
 
        sout << "r " << std::fixed << std::setprecision(0) << std::setw(5) << h_r << " z " << std::fixed
             << std::setprecision(0) << std::setw(5) << h_z << " theta " << std::fixed << std::setprecision(3)
             << std::setw(4) << h_theta << " phi " << std::fixed << std::setprecision(3) << std::setw(4) << h_phi
              << " error " << std::fixed << std::setprecision(2) << std::setw(5) << error;
        if (rpc) sout << "  time " << std::fixed << std::setprecision(2) << std::setw(5) << rpcTime;
    }
    return sout.str();
}

print() [17/21]

std::string Muon::MuonEDMPrinterTool::print ( const Trk::ResidualPull & resPull )

static

print ResidualPull object to string

Definition at line 861 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    const std::vector<double>& residual = resPull.residual();
    const std::vector<double>& pull     = resPull.pull();
    for (unsigned int i = 0; i < residual.size(); ++i) {
        if (residual[i] != 999. && residual[i] != -999.)
            sout << " residual " << std::setprecision(3) << std::setw(8) << residual[i] << "  ";
    }
    sout << " pull ";
    for (unsigned int i = 0; i < pull.size(); ++i) sout << std::setprecision(3) << std::setw(8) << pull[i] << "  ";
    return sout.str();
}

print() [18/21]

std::string Muon::MuonEDMPrinterTool::print

(

const Trk::Track &

track

)

const

print track parameters to string

Definition at line 279 of file MuonEDMPrinterTool.cxx.

{
 
    std::ostringstream sout;
 
    sout << printPatRec(track);
 
    const Trk::FitQuality* fq = track.fitQuality();
    if (fq) {
        sout << std::setprecision(4) << " : chi2 " << fq->chiSquared() << " ndof " << fq->numberDoF();
    } else {
        sout << " no fit Quality ";
    }
    const Trk::Perigee* pp = track.perigeeParameters();
    if (pp) {
        sout << "  " << print(*pp);
    } else {
        sout << " no perigee ";
    }
    sout <<std::endl;
    sout<<"-----  Track states -------"<<std::endl;
    unsigned int n{1};
    for (const Trk::TrackStateOnSurface* tsos : *track.trackStateOnSurfaces()) {
        if (tsos->measurementOnTrack())      sout<<" **** "<<std::setw(3)<<n<<" Measurement: "<<print(*tsos->measurementOnTrack())<<std::endl;
        if (tsos->materialEffectsOnTrack())  sout<<" **** "<<std::setw(3)<<n<<" Material:    "<<print(*tsos->materialEffectsOnTrack())<<std::endl;
        if (tsos->alignmentEffectsOnTrack()) sout<<" **** "<<std::setw(3)<<n<<" AEOT:        "<<print(*tsos->alignmentEffectsOnTrack())<<std::endl;
        if (tsos->trackParameters())         sout<<" **** "<<std::setw(3)<<n<<" Parameters:  "<<print(*tsos->trackParameters())<<std::endl;
        ++n;
    }
    sout<<"-----------------------------"<<std::endl;
 
    return sout.str();
}

print() [19/21]

std::string Muon::MuonEDMPrinterTool::print ( const Trk::TrackParameters & pars )

static

print Trk::TrackParameters

Definition at line 845 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
    sout << "r " << std::fixed << std::setprecision(0) << std::setw(5)
         << pars.position().perp()
         << " z " << std::fixed << std::setprecision(0) << std::setw(6) << pars.position().z() << " theta "
         << std::fixed << std::setprecision(5) << std::setw(7) << pars.momentum().theta()<< " eta "
         << std::fixed << std::setprecision(5) << std::setw(7) << pars.momentum().eta() << " phi " << std::fixed
         << std::setprecision(3) << std::setw(6) << pars.momentum().phi() << " q*p(GeV) " << std::scientific
         << std::setprecision(3) << std::setw(10) << pars.momentum().mag() * pars.charge() * 1e-3 << " pT(GeV) "
         << std::scientific << std::setprecision(3) << std::setw(9) << pars.momentum().perp() * 1e-3;
 
    return sout.str();
}

print() [20/21]

std::string Muon::MuonEDMPrinterTool::print

(

const Trk::TrackStateOnSurface &

tsos

)

const

print the track state on surface to the string

Definition at line 1002 of file MuonEDMPrinterTool.cxx.

                                                                            {
   std::stringstream tsos_str{};
   if (tsos.measurementOnTrack()) tsos_str<<"Measurement: "<<print(*tsos.measurementOnTrack())<<"\t";
   if (tsos.materialEffectsOnTrack()) tsos_str<<"Material "<<print(*tsos.materialEffectsOnTrack())<<"\t";
   if (tsos.trackParameters()) tsos_str<<"Parameters: "<<print(*tsos.trackParameters())<<"\t";
 
   return tsos_str.str();
}

print() [21/21]

std::string Muon::MuonEDMPrinterTool::print

(

std::vector< std::unique_ptr< MuonSegment > > &

segs

)

const

Definition at line 516 of file MuonEDMPrinterTool.cxx.

{
 
    std::ostringstream sout;
 
    std::vector<std::unique_ptr<MuonSegment> >::iterator it      = segs.begin();
    std::vector<std::unique_ptr<MuonSegment> >::iterator it_end  = segs.end();
    std::vector<std::unique_ptr<MuonSegment> >::iterator it_last = it_end;
    --it_last;
    for (; it != it_end; ++it) {
        if (!*it) {
            sout << " WARNING, zero pointer detected in MuonSegment vector!!! " << std::endl;
            continue;
        }
        sout << print(**it);
        if (it != it_last) sout << std::endl;
    }
 
    return sout.str();
}

printData()

std::string Muon::MuonEDMPrinterTool::printData

(

const Trk::MeasurementBase &

measurement

)

const

print data part of Muon MeasurementBase to string

Definition at line 878 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream sout;
 
    // print position of hit
    double h_r     = measurement.globalPosition().perp();
    double h_z     = measurement.globalPosition().z();
    double h_phi   = measurement.globalPosition().phi();
    double h_theta = measurement.globalPosition().theta();
    double error = std::sqrt(measurement.localCovariance()(0, 0));
    const Trk::LocalParameters& lpars = measurement.localParameters();
    double local_pos =  lpars.contains(Trk::locX) ? lpars[Trk::locX] : (lpars.contains(Trk::locY) ? lpars[Trk::locY] : 0.);
    sout << "r " << std::fixed << std::setprecision(0) << std::setw(5) << h_r << " z " << std::fixed
         << std::setprecision(0) << std::setw(5) << h_z << " phi " << std::fixed << std::setprecision(3) << std::setw(4)
         << h_phi << " theta " << std::fixed << std::setprecision(3) << std::setw(4) << h_theta
         << " lPos " << std::fixed << std::setprecision(2) << std::setw(5) << local_pos
         << " error " << std::fixed << std::setprecision(2) << std::setw(5) << error;
 
 
    // print measurement data
    const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(&measurement);
    if (rot) {
        // add drift time for MDT
        const MdtDriftCircleOnTrack* mdt = dynamic_cast<const MdtDriftCircleOnTrack*>(rot);
        if (mdt) {
            sout << "  r_drift " << std::fixed << std::setprecision(2) << std::setw(5) << mdt->driftRadius();
 
        } else {
            // add time for RPC
            const RpcClusterOnTrack* rpc = dynamic_cast<const RpcClusterOnTrack*>(rot);
            if (rpc) {
                const RpcPrepData* rpcPRD = rpc->prepRawData();
                if (rpcPRD) {
                    sout << "  time " << std::fixed << std::setprecision(2) << std::setw(5) << rpcPRD->time();
                }
            }
        }
    } else {  // !rot
        // if we get here: not a ROT, maybe a CROT
        const CompetingMuonClustersOnTrack* crot = dynamic_cast<const CompetingMuonClustersOnTrack*>(&measurement);
        if (crot) {
            unsigned int                                  nlayers = 0;
            unsigned int                                  nhits   = 0;
            std::set<Identifier>                          layers;
            std::vector<double>                           rpcTimes;
            const std::vector<const MuonClusterOnTrack*>& rots = crot->containedROTs();
            nhits                                              = rots.size();
            rpcTimes.reserve(nhits);
            std::vector<const MuonClusterOnTrack*>::const_iterator itR = rots.begin(), itR_end = rots.end();
            for (; itR != itR_end; ++itR) {
                Identifier layerId = m_idHelperSvc->layerId((*itR)->identify());
                layers.insert(layerId);
                const RpcClusterOnTrack* rpc    = dynamic_cast<const RpcClusterOnTrack*>(*itR);
                const RpcPrepData*       rpcPRD = rpc ? rpc->prepRawData() : nullptr;
                if (rpcPRD) rpcTimes.push_back(rpcPRD->time());
            }
            nlayers = layers.size();
            sout << "  CompRot: hits " << nhits << " layers " << nlayers;
            // add time for RPC
            if (!rpcTimes.empty()) {
                sout << (rpcTimes.size() == 1 ? "  time" : "  times") << std::fixed << std::setprecision(2);
                std::vector<double>::iterator itD = rpcTimes.begin(), itD_end = rpcTimes.end();
                for (; itD != itD_end; ++itD) sout << " " << std::setw(5) << *itD;
            }
        }  // if crot
    }      // else !rot
 
    return sout.str();
}  // printData( Trk::MeasurementBase )

printId()

std::string Muon::MuonEDMPrinterTool::printId

(

const Trk::MeasurementBase &

measurement

)

const

print identifier part of Muon MeasurementBase to string

Definition at line 950 of file MuonEDMPrinterTool.cxx.

{
    std::string idStr;
    Identifier  id = m_edmHelperSvc->getIdentifier(measurement);
    if (!id.is_valid()) {
        const Trk::PseudoMeasurementOnTrack* pseudo = dynamic_cast<const Trk::PseudoMeasurementOnTrack*>(&measurement);
        if (pseudo)
            idStr = "pseudo measurement";
        else
            idStr = "no Identifier";
    } else if (!m_idHelperSvc->isMuon(id)) {
        idStr = "Id hit";
    } else {
        idStr = m_idHelperSvc->toString(id);
    }
 
    return idStr;
}

printMeasurements()

std::string Muon::MuonEDMPrinterTool::printMeasurements

(

const Trk::Track &

track

)

const

print measurements on track to string

Definition at line 154 of file MuonEDMPrinterTool.cxx.

{
    std::ostringstream                                sout;
    const Trk::TrackStates* states = track.trackStateOnSurfaces();
    if (!states) return "";
    Trk::TrackStates::const_iterator it      = states->begin();
    Trk::TrackStates::const_iterator it_end  = states->end();
    Trk::TrackStates::const_iterator it2     = states->begin();
    Trk::TrackStates::const_iterator it2_end = states->end();
    // Build map of AEOTs and the measurements they affect
    std::multimap<const Trk::MeasurementBase*, const Trk::AlignmentEffectsOnTrack*> measAndTheirAlignmentEffects;
    const Trk::MeasurementBase*                                                     m    = nullptr;
    const Trk::AlignmentEffectsOnTrack*                                             aeot = nullptr;
    for (; it != it_end; ++it) {
        aeot = (*it)->alignmentEffectsOnTrack();
        if (aeot) {
            // Now get the list of identifiers which this AEOT impacts.
            const std::vector<Identifier>& identifiers = aeot->vectorOfAffectedTSOS();
            it2                                        = states->begin();
            it2_end                                    = states->end();
            for (; it2 != it2_end; ++it2) {
                m = (*it2)->measurementOnTrack();
                if (m) {
                    Identifier id = m_edmHelperSvc->getIdentifier(*m);
                    if ((id.is_valid() && (std::find(identifiers.begin(), identifiers.end(), id) != identifiers.end()))
                        || (aeot->effectsLastFromNowOn() && it2 > it))
                    {
                        // Either this measurement is explicitly listed, OR it is in a TSOS after an AEOT whose effects
                        // last from now on.
                        measAndTheirAlignmentEffects.insert(
                            std::pair<const Trk::MeasurementBase*, const Trk::AlignmentEffectsOnTrack*>(m, aeot));
                    }
                }
            }
        }
    }
    //    std::cout << " measAndTheirAlignmentEffects " << measAndTheirAlignmentEffects.size() << std::endl;
 
    // Reset
    it     = states->begin();
    it_end = states->end();
    // first loop to get width of Id's for nice alignment
    std::vector<std::string> idStrings;
    std::vector<std::string> dataStrings;
    unsigned int             idWidth = 0;
    for (; it != it_end; ++it) {
        m = (*it)->measurementOnTrack();
        if (m) {
            // Identifier part
            std::string idStr = printId(*m);
            idStrings.push_back(idStr);
            if (idStr.length() > idWidth) idWidth = idStr.length();
            // Data part
            const Trk::TrackParameters* trackParameters = (*it)->trackParameters();
            std::string                 dataStr         = printData(*m);
            if (trackParameters) {
                std::multimap<const Trk::MeasurementBase*, const Trk::AlignmentEffectsOnTrack*>::iterator itMap =
                    measAndTheirAlignmentEffects.begin();
                itMap = measAndTheirAlignmentEffects.find(m);
                if (itMap != measAndTheirAlignmentEffects.end()) {
                    std::vector<const Trk::AlignmentEffectsOnTrack*> aeotos;
                    aeotos.push_back(itMap->second);
                    ++itMap;
                    if (itMap != measAndTheirAlignmentEffects.end() && itMap->first == m)
                        aeotos.push_back(itMap->second);
                    std::optional<Trk::ResidualPull> resPull{m_pullCalculator->residualPull(
                        m, trackParameters, Trk::ResidualPull::Unbiased, Trk::TrackState::unidentified, aeotos)};
                    if (resPull) dataStr += print(*resPull);
                    if (resPull) dataStr += " (AEOT)";
                } else {
                    std::optional<Trk::ResidualPull> resPull{
                        m_pullCalculator->residualPull(m, trackParameters, Trk::ResidualPull::Unbiased)};
                    if (resPull) dataStr += print(*resPull);
                }
            }
            if ((*it)->type(Trk::TrackStateOnSurface::Outlier)) {
                dataStr += " (Outlier)";
            } else if ((*it)->type(Trk::TrackStateOnSurface::Hole)) {
                dataStr += " (Hole)";
            }
            dataStrings.push_back(dataStr);
        }
        aeot = (*it)->alignmentEffectsOnTrack();
        if (aeot) {
            std::string idStr = " AEOT ";
            idStrings.push_back(idStr);
            std::ostringstream souta;
            souta << std::setprecision(3) << " deltaTranslation (mm) " << aeot->deltaTranslation() << " error "
                  << aeot->sigmaDeltaTranslation() << " deltaAngle (mrad) " << 1000 * aeot->deltaAngle() << " error "
                  << 1000 * aeot->sigmaDeltaAngle();
            dataStrings.push_back(souta.str());
        }
    }
 
    // second loop to print out aligned strings
    unsigned int n = idStrings.size();
    for (unsigned int i = 0; i < n; ++i) {
        sout << std::left << std::setw(idWidth) << idStrings[i] << std::right << "  " << dataStrings[i];
        if (i != n - 1) sout << std::endl;
    }
 
    return sout.str();
}

printPatRec()

std::string Muon::MuonEDMPrinterTool::printPatRec ( const Trk::Track & track )

static

produce a string indicating who produced the track

Definition at line 260 of file MuonEDMPrinterTool.cxx.

{
    if (track.info().patternRecoInfo(Trk::TrackInfo::STACO)) return "STACO";
    if (track.info().patternRecoInfo(Trk::TrackInfo::StacoLowPt)) return "StacoLowPt";
    if (track.info().patternRecoInfo(Trk::TrackInfo::Muonboy)) return "Muonboy";
    if (track.info().patternRecoInfo(Trk::TrackInfo::MuTag)) return "MuTag";
    if (track.info().patternRecoInfo(Trk::TrackInfo::MuidComb)
        || track.info().patternRecoInfo(Trk::TrackInfo::MuidCombined))
        return "MuidCombined";
    if (track.info().patternRecoInfo(Trk::TrackInfo::MuidStandAlone)) return "MuidStandAlone";
    if (track.info().patternRecoInfo(Trk::TrackInfo::Moore)) return "Moore";
    if (track.info().patternRecoInfo(Trk::TrackInfo::MuGirl)) return "MuGirl";
    if (track.info().patternRecoInfo(Trk::TrackInfo::MuGirlUnrefitted)) return "MuGirlUnrefitted";
    if (track.info().patternRecoInfo(Trk::TrackInfo::MuidVertexAssociator)) return "MuidVertexAssociator";
    if (track.info().patternRecoInfo(Trk::TrackInfo::MuGirlLowBeta)) return "MuGirlLowBeta";
    return "Unknown";
}

printStations()

std::string Muon::MuonEDMPrinterTool::printStations

(

const Trk::Track &

track

)

const

print stations on track to string

Definition at line 136 of file MuonEDMPrinterTool.cxx.

{
 
    const Trk::MuonTrackSummary* muonSummary = nullptr;
 
    // check if the track already has a MuonTrackSummary, if not calculate it using the helper
    const Trk::TrackSummary* summary = track.trackSummary();
    if (summary) muonSummary = summary->muonTrackSummary();
    if (!muonSummary) {
    ATH_MSG_WARNING("No muon summary is present");
    return "";
    }
 
    return print(*muonSummary);
}

renounce()

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

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

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

overridevirtualinherited

Perform system initialization for an algorithm.

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

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

sysStart()

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

overridevirtualinherited

Handle START transition.

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

updateVHKA()

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

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // 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_DetectorManagerKey

SG::ReadCondHandleKey<MuonGM::MuonDetectorManager> Muon::MuonEDMPrinterTool::m_DetectorManagerKey

private

Initial value:

{
        this,
        "DetectorManagerKey",
        "MuonDetectorManager",
        "Key of input MuonDetectorManager condition data",
    }

Definition at line 155 of file MuonEDMPrinterTool.h.

                                                                       {
        this,
        "DetectorManagerKey",
        "MuonDetectorManager",
        "Key of input MuonDetectorManager condition data",
    };

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_edmHelperSvc

ServiceHandle<IMuonEDMHelperSvc> Muon::MuonEDMPrinterTool::m_edmHelperSvc

private

Initial value:

{
        this,
        "edmHelper",
        "Muon::MuonEDMHelperSvc/MuonEDMHelperSvc",
        "Handle to the service providing the IMuonEDMHelperSvc interface",
    }

Definition at line 142 of file MuonEDMPrinterTool.h.

                                                   {
        this,
        "edmHelper",
        "Muon::MuonEDMHelperSvc/MuonEDMHelperSvc",
        "Handle to the service providing the IMuonEDMHelperSvc interface",
    };

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> Muon::MuonEDMPrinterTool::m_idHelperSvc

private

Initial value:

{
        this,
        "MuonIdHelperSvc",
        "Muon::MuonIdHelperSvc/MuonIdHelperSvc",
    }

Definition at line 137 of file MuonEDMPrinterTool.h.

                                                     {
        this,
        "MuonIdHelperSvc",
        "Muon::MuonIdHelperSvc/MuonIdHelperSvc",
    };

m_mdtKey

SG::ReadHandleKey<MdtPrepDataContainer> Muon::MuonEDMPrinterTool::m_mdtKey

private

Initial value:

{
        this,
        "MdtPrdCollection",
        "MDT_DriftCircles",
        "MDT PRD Container",
    }

Definition at line 161 of file MuonEDMPrinterTool.h.

                                                  {
        this,
        "MdtPrdCollection",
        "MDT_DriftCircles",
        "MDT PRD Container",
    };

m_pullCalculator

ToolHandle<Trk::IResidualPullCalculator> Muon::MuonEDMPrinterTool::m_pullCalculator

private

Initial value:

{
        this,
        "ResidualPullCalculator",
        "Trk::ResidualPullCalculator/ResidualPullCalculator",
    }

Definition at line 149 of file MuonEDMPrinterTool.h.

                                                           {
        this,
        "ResidualPullCalculator",
        "Trk::ResidualPullCalculator/ResidualPullCalculator",
    };

m_rpcKey

SG::ReadHandleKey<RpcPrepDataContainer> Muon::MuonEDMPrinterTool::m_rpcKey

private

Initial value:

{
        this,
        "RpcPrdCollection",
        "RPC_Measurements",
        "RPC PRD Container",
    }

Definition at line 167 of file MuonEDMPrinterTool.h.

                                                  {
        this,
        "RpcPrdCollection",
        "RPC_Measurements",
        "RPC PRD Container",
    };

m_tgcKey

SG::ReadHandleKey<TgcPrepDataContainer> Muon::MuonEDMPrinterTool::m_tgcKey

private

Initial value:

{
        this,
        "TgcPrdCollection",
        "TGC_Measurements",
        "TGC PRD Container",
    }

Definition at line 173 of file MuonEDMPrinterTool.h.

                                                  {
        this,
        "TgcPrdCollection",
        "TGC_Measurements",
        "TGC PRD Container",
    };

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• MuonEDMPrinterTool.h
• MuonEDMPrinterTool.cxx