Muon::MuonEDMHelperSvc Class Reference

Helper tool containing functionality needed by multiple tools. More...

#include <MuonEDMHelperSvc.h>

Inheritance diagram for Muon::MuonEDMHelperSvc:

Collaboration diagram for Muon::MuonEDMHelperSvc:

Public Member Functions
	MuonEDMHelperSvc (const std::string &name, ISvcLocator *svc)
	default AlgService constructor
virtual	~MuonEDMHelperSvc () override
	destructor
virtual StatusCode	initialize () override
	initialize method, method taken from bass-class AlgTool
virtual Identifier	getIdentifier (const Trk::MeasurementBase &meas) const override
	tries to get Identifier for measurement, if not possible it will return Identifier()
virtual Identifier	chamberId (const MuonSegment &seg) const override
	returns the first MDT or CSC identifier of the segment
virtual std::set< Identifier >	chamberIds (const MuonSegment &seg) const override
	returns a set containing the chamber Ids of all MDT or CSC chambers on the segment
virtual bool	isEndcap (const MuonSegment &seg) const override
	returns whether the segment is in the endcap
virtual bool	isEndcap (const Trk::Track &seg) const override
	returns whether the track is in the endcap
virtual const Trk::AtaPlane *	createTrackParameters (const MuonSegment &seg, double momentum=1., double charge=0.) const override
	brief create a AtaPlane track parameter for segment
virtual bool	goodTrack (const Trk::Track &track, double chi2Cut) const override
	check whether track satisfies selection chi2 cut
virtual bool	isSLTrack (const Trk::Track &track) const override
	determines whether a track is a straight line track (no momentum fitted).

Private Attributes
ServiceHandle< Muon::IMuonIdHelperSvc >	m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

Detailed Description

Helper tool containing functionality needed by multiple tools.

Definition at line 33 of file MuonEDMHelperSvc.h.

Constructor & Destructor Documentation

MuonEDMHelperSvc()

Muon::MuonEDMHelperSvc::MuonEDMHelperSvc	(	const std::string &	name,
		ISvcLocator *	svc )

default AlgService constructor

Definition at line 24 of file MuonEDMHelperSvc.cxx.

: base_class(name, svc) {}

~MuonEDMHelperSvc()

virtual Muon::MuonEDMHelperSvc::~MuonEDMHelperSvc ( )

inlineoverridevirtual

destructor

Definition at line 39 of file MuonEDMHelperSvc.h.

{};

Member Function Documentation

chamberId()

Identifier Muon::MuonEDMHelperSvc::chamberId ( const MuonSegment & seg ) const

overridevirtual

returns the first MDT or CSC identifier of the segment

Definition at line 50 of file MuonEDMHelperSvc.cxx.

                                                                       {
    Identifier chid;
    for(const Trk::MeasurementBase* meas : seg.containedMeasurements()){
 
      // get Identifier
      Identifier id = getIdentifier(*meas);
      if( !id.is_valid() ) continue;
 
      // create chamber ID
      chid = m_idHelperSvc->chamberId(id);
 
      // stop at first none trigger hit
      if( !m_idHelperSvc->isTrigger(id) ) break;
    }
 
    if( !chid.is_valid() ){
      ATH_MSG_WARNING("Got segment without valid identifiers");
    }
 
    return chid;
  }

chamberIds()

std::set< Identifier > Muon::MuonEDMHelperSvc::chamberIds ( const MuonSegment & seg ) const

overridevirtual

returns a set containing the chamber Ids of all MDT or CSC chambers on the segment

Definition at line 72 of file MuonEDMHelperSvc.cxx.

                                                                                {
 
    std::set<Identifier> chIds;
    Identifier chid;
    Identifier chidTrig;
    for(const Trk::MeasurementBase* meas : seg.containedMeasurements()){
 
      // get Identifier
      Identifier id = getIdentifier(*meas);
      if( !id.is_valid() ) continue;
 
      // create chamber ID
      chid = m_idHelperSvc->chamberId(id);
 
      // stop at first none trigger hit
      if( !m_idHelperSvc->isTrigger(id) ) chIds.insert(chid);
      else chidTrig = chid;      
    }
    if( chIds.empty() ) {
      chIds.insert(chidTrig);
    }
 
    return chIds;
  }

createTrackParameters()

const Trk::AtaPlane * Muon::MuonEDMHelperSvc::createTrackParameters ( const MuonSegment & seg, double momentum = 1., double charge = 0. ) const

overridevirtual

brief create a AtaPlane track parameter for segment

Definition at line 119 of file MuonEDMHelperSvc.cxx.

                                                                                                                           {
 
    // we need a none zero momentum
    if( momentum == 0. ) {
      ATH_MSG_WARNING(" cannot create parameters with zero momentum ");
      return nullptr;
    }
    double locx = seg.localParameters().contains(Trk::locX) ? seg.localParameters()[Trk::locX] : 0.;
    double locy = seg.localParameters().contains(Trk::locY) ? seg.localParameters()[Trk::locY] : 0.;
    double qoverp = charge/momentum;
    return dynamic_cast<const Trk::AtaPlane*>(
      seg.associatedSurface()
        .createUniqueParameters<5, Trk::Charged>(locx,
                                           locy,
                                           seg.globalDirection().phi(),
                                           seg.globalDirection().theta(),
                                           qoverp)
        .release());
  }

getIdentifier()

Identifier Muon::MuonEDMHelperSvc::getIdentifier ( const Trk::MeasurementBase & meas ) const

overridevirtual

tries to get Identifier for measurement, if not possible it will return Identifier()

Definition at line 31 of file MuonEDMHelperSvc.cxx.

                                                                                   {
    const Trk::RIO_OnTrack* rot = dynamic_cast<const Trk::RIO_OnTrack*>(&meas);
    if( rot ) return rot->identify();
    const CompetingMuonClustersOnTrack* crot = dynamic_cast<const CompetingMuonClustersOnTrack*>(&meas);
    if( crot ) {
      if( crot->containedROTs().empty() ){
          ATH_MSG_WARNING(" CompetingMuonClustersOnTrack without contained ROTs ");
          return {};
      }
      if( !crot->containedROTs().front() ) {
        ATH_MSG_WARNING(" CompetingMuonClustersOnTrack contains a ROT pointer that is zero ");
          return {};
      }
      return crot->containedROTs().front()->identify();
    }
    return {};
  }

goodTrack()

bool Muon::MuonEDMHelperSvc::goodTrack ( const Trk::Track & track, double chi2Cut ) const

overridevirtual

check whether track satisfies selection chi2 cut

Definition at line 139 of file MuonEDMHelperSvc.cxx.

                                                                                {
 
    // get reduced chi2
    const Trk::FitQuality* fq = track.fitQuality();
    if( !fq || fq->numberDoF() == 0 ){
      return false;
    }
 
    double reducedChi2 = fq->chiSquared()/fq->numberDoF();
    // reject fit if larger than cut
    return reducedChi2 <= chi2Cut;
  }

initialize()

StatusCode Muon::MuonEDMHelperSvc::initialize ( )

overridevirtual

initialize method, method taken from bass-class AlgTool

Definition at line 25 of file MuonEDMHelperSvc.cxx.

                                          {
    ATH_CHECK(AthService::initialize());
    ATH_CHECK(m_idHelperSvc.retrieve());
    return StatusCode::SUCCESS;
  }

isEndcap() [1/2]

bool Muon::MuonEDMHelperSvc::isEndcap ( const MuonSegment & seg ) const

overridevirtual

returns whether the segment is in the endcap

Definition at line 97 of file MuonEDMHelperSvc.cxx.

                                                                {
    return m_idHelperSvc->isEndcap( chamberId(seg) );
  }

isEndcap() [2/2]

bool Muon::MuonEDMHelperSvc::isEndcap ( const Trk::Track & seg ) const

overridevirtual

returns whether the track is in the endcap

Definition at line 101 of file MuonEDMHelperSvc.cxx.

                                                               {
 
    const DataVector<const Trk::MeasurementBase>* measurements = track.measurementsOnTrack();
    if( !measurements ) return true;
 
    // loop over measurements in reversed order (should be faster as endcap hits tend to be at the end of the track)
    DataVector<const Trk::MeasurementBase>::const_reverse_iterator rit = measurements->rbegin();
    DataVector<const Trk::MeasurementBase>::const_reverse_iterator rit_end = measurements->rend();
    for( ;rit!=rit_end;++rit ){
      Identifier id = getIdentifier(**rit);
      if( !id.is_valid() ) continue;
 
      if( m_idHelperSvc->isEndcap(id) ) return true;
    }
    return false;
  }

isSLTrack()

bool Muon::MuonEDMHelperSvc::isSLTrack ( const Trk::Track & track ) const

overridevirtual

determines whether a track is a straight line track (no momentum fitted).

sum covariance terms of momentum, use it to determine whether fit was SL fit

Definition at line 152 of file MuonEDMHelperSvc.cxx.

                                                                {
    // use track info if set properly
    if( track.info().trackProperties(Trk::TrackInfo::StraightTrack) ) return true;
 
    // else try using error matrix
    const Trk::Perigee* pp = track.perigeeParameters();
    if (!pp || !pp->covariance()) return false;   
      double momCov = 0.;
      for( int i=0;i<4;++i ) momCov += std::abs( (*pp->covariance())(4,i) );
    return momCov < 1.e-10;
  }

Member Data Documentation

m_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> Muon::MuonEDMHelperSvc::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}

private

Definition at line 69 of file MuonEDMHelperSvc.h.

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

---

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

• MuonEDMHelperSvc.h
• MuonEDMHelperSvc.cxx