Csc2dSegmentMaker Class Reference

#include <Csc2dSegmentMaker.h>

Inheritance diagram for Csc2dSegmentMaker:

Collaboration diagram for Csc2dSegmentMaker:

Public Types
typedef std::vector< Cluster >	TrkClusters
typedef std::vector< Cluster >	ChamberTrkClusters[4]
typedef std::vector< const Trk::RIO_OnTrack * >	RioList
typedef DataVector< const Trk::MeasurementBase >	MbaseList
typedef Muon::MuonSegmentCombination::SegmentVec	SegmentVec
typedef std::vector< Segment >	Segments

Public Member Functions
	Csc2dSegmentMaker (const std::string &, const std::string &, const IInterface *)
	~Csc2dSegmentMaker ()=default
StatusCode	initialize ()
std::unique_ptr< MuonSegmentCombinationCollection >	find (const std::vector< const Muon::CscPrepDataCollection * > &pcols, const EventContext &ctx) const
std::unique_ptr< MuonSegmentCombinationCollection >	find (const MuonSegmentCombinationCollection &, const EventContext &ctx) const
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 ()

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
Muon::MuonSegmentCombination *	findSegmentCombination (const Muon::CscPrepDataCollection &pcol, const EventContext &ctx) const
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 {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
Gaudi::Property< std::string >	m_cscdig_sg_inkey {this, "scdig_sg_inkey", "CSC_Measurements"}
ToolHandle< ICscSegmentUtilTool >	m_segmentTool {this, "segmentTool", "CscSegmentUtilTool/CscSegmentUtilTool"}
ToolHandle< Muon::IMuonClusterOnTrackCreator >	m_cscClusterOnTrackCreator
PublicToolHandle< Muon::MuonEDMPrinterTool >	m_printer {this, "printerTool", "Muon::MuonEDMPrinterTool/MuonEDMPrinterTool"}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 24 of file Csc2dSegmentMaker.h.

Member Typedef Documentation

ChamberTrkClusters

typedef std::vector<Cluster> ICscSegmentFinder::ChamberTrkClusters[4]

inherited

Definition at line 83 of file ICscSegmentFinder.h.

MbaseList

typedef DataVector<const Trk::MeasurementBase> ICscSegmentFinder::MbaseList

inherited

Definition at line 86 of file ICscSegmentFinder.h.

RioList

typedef std::vector<const Trk::RIO_OnTrack*> ICscSegmentFinder::RioList

inherited

Definition at line 85 of file ICscSegmentFinder.h.

Segments

typedef std::vector<Segment> ICscSegmentFinder::Segments

inherited

Definition at line 132 of file ICscSegmentFinder.h.

SegmentVec

typedef Muon::MuonSegmentCombination::SegmentVec ICscSegmentFinder::SegmentVec

inherited

Definition at line 88 of file ICscSegmentFinder.h.

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

TrkClusters

typedef std::vector<Cluster> ICscSegmentFinder::TrkClusters

inherited

Definition at line 82 of file ICscSegmentFinder.h.

Constructor & Destructor Documentation

Csc2dSegmentMaker()

Csc2dSegmentMaker::Csc2dSegmentMaker	(	const std::string &	type,
		const std::string &	aname,
		const IInterface *	parent )

Definition at line 65 of file Csc2dSegmentMaker.cxx.

                                                                                                            :
    AthAlgTool(type, aname, parent) {
    declareInterface<ICscSegmentFinder>(this);
}

~Csc2dSegmentMaker()

Csc2dSegmentMaker::~Csc2dSegmentMaker ( )

default

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

find() [1/2]

std::unique_ptr< MuonSegmentCombinationCollection > Csc2dSegmentMaker::find ( const MuonSegmentCombinationCollection & , const EventContext & ctx ) const

virtual

Implements ICscSegmentFinder.

Definition at line 249 of file Csc2dSegmentMaker.cxx.

                                                                                                     {
    return {};
}

find() [2/2]

std::unique_ptr< MuonSegmentCombinationCollection > Csc2dSegmentMaker::find ( const std::vector< const Muon::CscPrepDataCollection * > & pcols, const EventContext & ctx ) const

virtual

Implements ICscSegmentFinder.

Definition at line 85 of file Csc2dSegmentMaker.cxx.

                                                                                                         {
    // Construct output segment collection.
    std::unique_ptr<MuonSegmentCombinationCollection> mpsegs(new MuonSegmentCombinationCollection);
 
    for (unsigned int icol = 0; icol < pcols.size(); ++icol) {
        const CscPrepDataCollection* clus = pcols[icol];
        if (clus->empty()) continue;  // skip zero cluster collection
        break;
    }
 
    unsigned int col_count = 0;
    // Loop over collections in the container.
    for (unsigned int icol = 0; icol < pcols.size(); ++icol) {
        const CscPrepDataCollection* clus = pcols[icol];
        ++col_count;
        ATH_MSG_VERBOSE("** Collection " << col_count << " has " << clus->size() << " clusters");
 
        if (clus->empty()) continue;  // skip zero cluster collection
        MuonSegmentCombination* pcol = findSegmentCombination(*clus, ctx);
        if (pcol) {
            mpsegs->push_back(pcol);
            ATH_MSG_DEBUG("Found 2d CSC segment " << m_printer->print(*pcol));
        }
    }
 
    if (!mpsegs->empty())
        ATH_MSG_DEBUG("found " << mpsegs->size() << " 2d csc segments");
    else
        mpsegs.reset();
 
    return mpsegs;
}

findSegmentCombination()

MuonSegmentCombination * Csc2dSegmentMaker::findSegmentCombination ( const Muon::CscPrepDataCollection & pcol, const EventContext & ctx ) const

private

Definition at line 121 of file Csc2dSegmentMaker.cxx.

                                                                                                                                  {
    // check whether input not empty
    if (clus.empty()) return nullptr;
 
    int col_station = 0;
    int col_eta = 0;
    int col_phisec = 0;
    Identifier eta_id, phi_id;
 
    ChamberTrkClusters eta_clus, phi_clus;
 
    const MuonGM::CscReadoutElement* detEl = clus.front()->detectorElement();
    if (!detEl) {
        ATH_MSG_WARNING("Failed to obtain CscReadoutElement for cluster");
        return nullptr;
    }
    Amg::Transform3D gToLocal = detEl->GlobalToAmdbLRSTransform();
    Amg::Vector3D lpos000 = gToLocal * Amg::Vector3D(0.0, 0.0, 0.0);
 
    int stationEta = detEl->getStationEta();
    int stationPhi = detEl->getStationPhi();
    int chamberLayer = detEl->ChamberLayer();
    std::string stationName = detEl->getStationName();
    std::string redName = stationName.substr(0, 3);  // CSS or CSL, all that the CscIdHelper wants
    std::string isPhi;
    ATH_MSG_DEBUG("in station " << stationName << " eta " << stationEta << ", phi " << stationPhi << ", chamber layer " << chamberLayer);
    // new method: we no longer assume automatically that eta and phi go together
    // this is the usual reason for problems due to broken wires, but it need not always be the case
    // instead, pass the status of every layer to the CscSegmentUtilTool
    // if there are only 2 layers for eta or phi, do 2-layer segment finding for eta or phi
    // if there are fewer than two good layers for eta or phi, so no eta or phi segment can be made, the final segment need not include
    // information from the missing side
 
    int layStatus[2] = {
        0, 0};  // eta=0,phi=1; each decimal place will be 0 for good, 1 for bad, so 1001 means first and 4th layers are bad, etc.
    for (int iLay = 0; iLay < 4; iLay++) {      // loop over the layers in this chamber
        for (int iPhi = 0; iPhi < 2; iPhi++) {  // eta or phi
            int nbad = 0;
            if (iPhi)
                isPhi = "phi";
            else
                isPhi = "eta";
            // ATH_MSG_DEBUG("get hashes for "<<detEl->maxNumberOfStrips(iPhi)<<" strips ");
            for (int iStrip = 0; iStrip < detEl->maxNumberOfStrips(iPhi); iStrip++) {
                ATH_MSG_DEBUG("get strip quality for " << isPhi << " layer " << iLay << " strip " << iStrip);
                Identifier stripId =
                    m_idHelperSvc->cscIdHelper().channelID(redName, stationEta, stationPhi, chamberLayer, iLay + 1, iPhi, iStrip + 1);
                IdentifierHash hashID;
                m_idHelperSvc->cscIdHelper().get_channel_hash(stripId, hashID);
                if (!m_segmentTool->isGood(hashID, ctx)) {
                    ATH_MSG_DEBUG("bad strip");
                    nbad++;
                }
            }
            ATH_MSG_DEBUG("found " << nbad << " bad strips out of " << detEl->maxNumberOfStrips(iPhi) << " in " << isPhi << " layer "
                                   << iLay);
            if (nbad == detEl->maxNumberOfStrips(iPhi)) { layStatus[iPhi] += pow(10, iLay); }
        }
    }
    ATH_MSG_DEBUG("final status for eta is " << layStatus[0] << " and phi " << layStatus[1]);
 
    for (CscPrepDataCollection::const_iterator iclu = clus.begin(); iclu != clus.end(); ++iclu) {
        const CscPrepData* pclu = *iclu;
        Identifier id = pclu->identify();
        int station = m_idHelperSvc->cscIdHelper().stationName(id) - 49;
        int eta = m_idHelperSvc->cscIdHelper().stationEta(id);
        int phisec = m_idHelperSvc->cscIdHelper().stationPhi(id);
        int iwlay = m_idHelperSvc->cscIdHelper().wireLayer(id);
        bool measphi = m_idHelperSvc->cscIdHelper().measuresPhi(id);
 
        if (iclu == clus.begin()) {
            col_station = station;
            col_eta = eta;
            col_phisec = phisec;
        }
        ATH_MSG_DEBUG("****** " << chamber(station, eta, phisec) << measphi_name(measphi) << iwlay);
 
        if (station != col_station || eta != col_eta || phisec != col_phisec) {
            ATH_MSG_WARNING("Inconsistent collection contents");
            return nullptr;
        }
 
        const Muon::MuonClusterOnTrack* clu = m_cscClusterOnTrackCreator->createRIO_OnTrack(*pclu, pclu->globalPosition());
        const Muon::CscClusterOnTrack* ptclu = dynamic_cast<const Muon::CscClusterOnTrack*>(clu);
        if (!ptclu) {
            ATH_MSG_WARNING("Failed to create CscClusterOnTrack from PRD");
            delete clu;
            continue;
        }
        Amg::Vector3D lpos = gToLocal * ptclu->globalPosition();
        if (measphi) {
            phi_id = id;
            phi_clus[iwlay - 1].emplace_back(lpos, ptclu, measphi);
        } else {
            eta_id = id;
            eta_clus[iwlay - 1].emplace_back(lpos, ptclu, measphi);
        }
    }
 
    // Insert r-phi-segments.
 
    int nHitLayer_eta = 0;
    int nHitLayer_phi = 0;
    for (int i = 0; i < 4; ++i) {
        if (!eta_clus[i].empty()) ++nHitLayer_eta;
        if (!phi_clus[i].empty()) ++nHitLayer_phi;
    }
    ATH_MSG_DEBUG("there are " << nHitLayer_eta << " eta hit layers and " << nHitLayer_phi << " phi hit layers");
 
    //  MuonSegmentCombination* pcol;
    MuonSegmentCombination* pcol = nullptr;
    if (nHitLayer_eta >= 2 || nHitLayer_phi >= 2) {
        ATH_MSG_DEBUG("Csc2dSegment calls get2dMuonSegmentCombination !!!");
        pcol = m_segmentTool->get2dMuonSegmentCombination(eta_id, phi_id, eta_clus, phi_clus, lpos000, ctx, layStatus[0], layStatus[1]);
    }
 
    // to avoid memory leak
    for (unsigned int i = 0; i < 4; ++i) {
        for (TrkClusters::iterator itclus = eta_clus[i].begin(); itclus != eta_clus[i].end(); ++itclus) { delete itclus->cl; }
        eta_clus[i].clear();
        for (TrkClusters::iterator itclus = phi_clus[i].begin(); itclus != phi_clus[i].end(); ++itclus) { delete itclus->cl; }
        phi_clus[i].clear();
    }
 
    return pcol;
}

initialize()

StatusCode Csc2dSegmentMaker::initialize

(

)

Definition at line 72 of file Csc2dSegmentMaker.cxx.

                                         {
    ATH_MSG_DEBUG("Initializing " << name());
    // Show keys.
    ATH_MSG_DEBUG("  SegmentTool: " << m_segmentTool.typeAndName());
    ATH_MSG_DEBUG("  Input cscdig key: " << m_cscdig_sg_inkey);
    ATH_CHECK(m_idHelperSvc.retrieve());
    ATH_CHECK(m_segmentTool.retrieve());
    ATH_CHECK(m_cscClusterOnTrackCreator.retrieve());
    ATH_CHECK(m_printer.retrieve());
    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 & ICscSegmentFinder::interfaceID ( )

inlinestaticinherited

Definition at line 31 of file ICscSegmentFinder.h.

{ return IID_ICscSegmentFinder; }

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.

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_cscClusterOnTrackCreator

ToolHandle<Muon::IMuonClusterOnTrackCreator> Csc2dSegmentMaker::m_cscClusterOnTrackCreator

private

Initial value:

{this, "cscRotCreator",
                                                                            "Muon::CscClusterOnTrackCreator/CscClusterOnTrackCreator"}

Definition at line 44 of file Csc2dSegmentMaker.h.

                                                                         {this, "cscRotCreator",
                                                                            "Muon::CscClusterOnTrackCreator/CscClusterOnTrackCreator"};

m_cscdig_sg_inkey

Gaudi::Property<std::string> Csc2dSegmentMaker::m_cscdig_sg_inkey {this, "scdig_sg_inkey", "CSC_Measurements"}

private

Definition at line 41 of file Csc2dSegmentMaker.h.

{this, "scdig_sg_inkey", "CSC_Measurements"};

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_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> Csc2dSegmentMaker::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 39 of file Csc2dSegmentMaker.h.

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

m_printer

PublicToolHandle<Muon::MuonEDMPrinterTool> Csc2dSegmentMaker::m_printer {this, "printerTool", "Muon::MuonEDMPrinterTool/MuonEDMPrinterTool"}

private

Definition at line 46 of file Csc2dSegmentMaker.h.

{this, "printerTool", "Muon::MuonEDMPrinterTool/MuonEDMPrinterTool"};

m_segmentTool

ToolHandle<ICscSegmentUtilTool> Csc2dSegmentMaker::m_segmentTool {this, "segmentTool", "CscSegmentUtilTool/CscSegmentUtilTool"}

private

Definition at line 43 of file Csc2dSegmentMaker.h.

{this, "segmentTool", "CscSegmentUtilTool/CscSegmentUtilTool"};

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:

• Csc2dSegmentMaker.h
• Csc2dSegmentMaker.cxx