MuonResidualsFromTrack Class Reference

#include <MuonResidualsFromTrack.h>

Public Types
using	BuilderToken = edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord >

Public Member Functions
const std::vector< DetId >	chamberIds () const
MuonChamberResidual *	chamberResidual (DetId chamberId, int type)
TMatrixD	choleskyCorrMatrix (DetId chamberId)
void	clear ()
bool	contains_TIDTEC () const
TMatrixDSym	corrMatrix (DetId chamberId)
TMatrixDSym	covMatrix (DetId chamberId)
const reco::Muon *	getMuon ()
const reco::Track *	getTrack ()
	MuonResidualsFromTrack (edm::ESHandle< TransientTrackingRecHitBuilder > builder, edm::ESHandle< MagneticField > magneticField, edm::ESHandle< GlobalTrackingGeometry > globalGeometry, edm::ESHandle< DetIdAssociator > muonDetIdAssociator_, edm::ESHandle< Propagator > prop, const Trajectory *traj, const reco::Track *recoTrack, AlignableNavigator *navigator, double maxResidual)
	MuonResidualsFromTrack (edm::ESHandle< GlobalTrackingGeometry > globalGeometry, const reco::Muon *recoMuon, AlignableNavigator *navigator, double maxResidual)
double	normalizedChi2 () const
double	trackerChi2 () const
int	trackerNumHits () const
double	trackerRedChi2 () const
	~MuonResidualsFromTrack ()

Static Public Member Functions
static edm::ESInputTag	builderESInputTag ()

Private Member Functions
void	addTrkCovMatrix (DetId, TrajectoryStateOnSurface &)

Private Attributes
std::vector< DetId >	m_chamberIds
bool	m_contains_TIDTEC
std::map< DetId, MuonChamberResidual * >	m_csc
std::map< DetId, MuonChamberResidual * >	m_dt13
std::map< DetId, MuonChamberResidual * >	m_dt2
const reco::Muon *	m_recoMuon
const reco::Track *	m_recoTrack
double	m_tracker_chi2
int	m_tracker_numHits
std::map< DetId, TMatrixDSym >	m_trkCovMatrix
TrajectoryStateCombiner	m_tsoscomb

Detailed Description

Id

Author

J. Pivarski - Texas A&M University pivar.nosp@m.ski@.nosp@m.physi.nosp@m.cs.t.nosp@m.amu.e.nosp@m.du

Definition at line 52 of file MuonResidualsFromTrack.h.

Member Typedef Documentation

BuilderToken

using MuonResidualsFromTrack::BuilderToken = edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord>

Definition at line 54 of file MuonResidualsFromTrack.h.

Constructor & Destructor Documentation

MuonResidualsFromTrack() [1/2]

MuonResidualsFromTrack::MuonResidualsFromTrack	(	edm::ESHandle< TransientTrackingRecHitBuilder >	builder,
		edm::ESHandle< MagneticField >	magneticField,
		edm::ESHandle< GlobalTrackingGeometry >	globalGeometry,
		edm::ESHandle< DetIdAssociator >	muonDetIdAssociator_,
		edm::ESHandle< Propagator >	prop,
		const Trajectory *	traj,
		const reco::Track *	recoTrack,
		AlignableNavigator *	navigator,
		double	maxResidual
	)

Definition at line 24 of file MuonResidualsFromTrack.cc.

References TrajectoryMeasurement::backwardPredictedState(), Surface::bounds(), TransientTrackingRecHitBuilder::build(), clear(), gather_cfg::cout, MuonSubdetId::CSC, DetId::det(), TrackingRecHit::dimension(), MuonSubdetId::DT, TrajectoryMeasurement::forwardPredictedState(), TrackingRecHit::geographicalId(), DetIdAssociator::getGeomDet(), TrajectoryStateOnSurface::globalPosition(), GlobalTrackingGeometry::idToDet(), TrajectoryStateOnSurface::isValid(), TrackingRecHit::isValid(), Bounds::length(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPosition(), TrackingRecHit::localPositionError(), LogTrace, m_chamberIds, m_contains_TIDTEC, m_csc, m_dt13, m_dt2, m_recoTrack, m_tracker_chi2, m_tracker_numHits, m_tsoscomb, HLT_2022v12_cff::magneticField, Trajectory::measurements(), metname, DetId::Muon, HLT_2022v12_cff::navigator, LocalTrajectoryError::positionError(), Propagator::propagate(), DetId::rawId(), TrajectoryMeasurement::recHit(), reco::Track::recHits(), TrackingRecHit::recHits(), MuonSubdetId::RPC, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), GeomDet::surface(), StripSubdetector::TEC, StripSubdetector::TID, GeomDet::toGlobal(), HLT_2022v12_cff::track, DetId::Tracker, TrajectoryMeasurement::updatedState(), Bounds::width(), x, PV3DBase< T, PVType, FrameType >::x(), hit::x, LocalError::xx(), y, PV3DBase< T, PVType, FrameType >::y(), hit::y, z, PV3DBase< T, PVType, FrameType >::z(), and hit::z.

    : m_recoTrack(recoTrack) {
  bool m_debug = false;

  if (m_debug) {
    std::cout << "BEGIN MuonResidualsFromTrack" << std::endl;
    const std::string metname = " *** MuonResidualsFromTrack *** ";
    LogTrace(metname) << "Tracking Component changed!";
  }

  clear();

  reco::TransientTrack track(*m_recoTrack, &*magneticField, globalGeometry);
  TransientTrackingRecHit::ConstRecHitContainer recHitsForRefit;
  int iT = 0, iM = 0;
  for (auto const& hit : m_recoTrack->recHits()) {
    if (hit->isValid()) {
      DetId hitId = hit->geographicalId();
      if (hitId.det() == DetId::Tracker) {
        iT++;
        if (m_debug)
          std::cout << "Tracker Hit " << iT << " is found. Add to refit. Dimension: " << hit->dimension() << std::endl;

        recHitsForRefit.push_back(trackerRecHitBuilder->build(&*hit));
      } else if (hitId.det() == DetId::Muon) {
        //        if ( hit->geographicalId().subdetId() == 3 && !theRPCInTheFit ) {
        //          LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged";
        //          continue;
        //        }
        iM++;
        if (m_debug)
          std::cout << "Muon Hit " << iM
                    << " is found. We do not add muon hits to refit. Dimension: " << hit->dimension() << std::endl;
        if (hitId.subdetId() == MuonSubdetId::DT) {
          const DTChamberId chamberId(hitId.rawId());
          if (m_debug)
            std::cout << "Muon Hit in DT wheel " << chamberId.wheel() << " station " << chamberId.station()
                      << " sector " << chamberId.sector() << "." << std::endl;
        } else if (hitId.subdetId() == MuonSubdetId::CSC) {
          const CSCDetId cscDetId(hitId.rawId());
          if (m_debug)
            std::cout << "Muon hit in CSC endcap " << cscDetId.endcap() << " station " << cscDetId.station() << " ring "
                      << cscDetId.ring() << " chamber " << cscDetId.chamber() << "." << std::endl;
        } else if (hitId.subdetId() == MuonSubdetId::RPC) {
          if (m_debug)
            std::cout << "Muon Hit in RPC" << std::endl;
        } else {
          if (m_debug)
            std::cout << "Warning! Muon Hit not in DT or CSC or RPC" << std::endl;
        }
        //        recHitsForRefit.push_back(theMuonRecHitBuilder->build(&*hit));
      }
    }
  }

  //  TrackTransformer trackTransformer();
  //  std::vector<Trajectory> vTrackerTrajectory = trackTransformer.transform(track, recHitsForReFit);
  //  std::cout << "Tracker trajectories size " << vTrackerTrajectory.size() << std::endl;

  TrajectoryStateOnSurface lastTrackerTsos;
  double lastTrackerTsosGlobalPositionR = 0.0;

  std::vector<TrajectoryMeasurement> vTrajMeasurement = traj->measurements();
  if (m_debug)
    std::cout << "  Size of vector of TrajectoryMeasurements: " << vTrajMeasurement.size() << std::endl;
  int nTrajMeasurement = 0;
  for (std::vector<TrajectoryMeasurement>::const_iterator iTrajMeasurement = vTrajMeasurement.begin();
       iTrajMeasurement != vTrajMeasurement.end();
       ++iTrajMeasurement) {
    nTrajMeasurement++;
    if (m_debug)
      std::cout << "    TrajectoryMeasurement #" << nTrajMeasurement << std::endl;

    TrajectoryMeasurement trajMeasurement = *iTrajMeasurement;

    TrajectoryStateOnSurface tsos =
        m_tsoscomb(trajMeasurement.forwardPredictedState(), trajMeasurement.backwardPredictedState());
    const TrajectoryStateOnSurface& tsosF = trajMeasurement.forwardPredictedState();
    const TrajectoryStateOnSurface& tsosB = trajMeasurement.backwardPredictedState();
    const TrajectoryStateOnSurface& tsosU = trajMeasurement.updatedState();
    if (m_debug)
      std::cout << "      TrajectoryMeasurement TSOS validity: " << tsos.isValid() << std::endl;
    if (tsos.isValid()) {
      double tsosGlobalPositionR = sqrt(tsos.globalPosition().x() * tsos.globalPosition().x() +
                                        tsos.globalPosition().y() * tsos.globalPosition().y());
      if (m_debug) {
        std::cout << "         TrajectoryMeasurement TSOS localPosition"
                  << " x: " << tsos.localPosition().x() << " y: " << tsos.localPosition().y()
                  << " z: " << tsos.localPosition().z() << std::endl;
        std::cout << "         TrajectoryMeasurement TSOS globalPosition"
                  << " x: " << tsos.globalPosition().x() << " y: " << tsos.globalPosition().y()
                  << " R: " << tsosGlobalPositionR << " z: " << tsos.globalPosition().z() << std::endl;
      }
      if (tsosGlobalPositionR > lastTrackerTsosGlobalPositionR) {
        lastTrackerTsos = tsos;
        lastTrackerTsosGlobalPositionR = tsosGlobalPositionR;
      }
    }

    const TransientTrackingRecHit* trajMeasurementHit = &(*trajMeasurement.recHit());
    if (m_debug)
      std::cout << "      TrajectoryMeasurement hit validity: " << trajMeasurementHit->isValid() << std::endl;
    if (trajMeasurementHit->isValid()) {
      DetId trajMeasurementHitId = trajMeasurementHit->geographicalId();
      int trajMeasurementHitDim = trajMeasurementHit->dimension();
      if (trajMeasurementHitId.det() == DetId::Tracker) {
        if (m_debug)
          std::cout << "      TrajectoryMeasurement hit Det: Tracker" << std::endl;
        if (m_debug)
          std::cout << "      TrajectoryMeasurement hit dimension: " << trajMeasurementHitDim << std::endl;
        m_tracker_numHits++;
        double xresid = tsos.localPosition().x() - trajMeasurementHit->localPosition().x();
        double xresiderr2 = tsos.localError().positionError().xx() + trajMeasurementHit->localPositionError().xx();
        m_tracker_chi2 += xresid * xresid / xresiderr2;

        if (trajMeasurementHitId.subdetId() == StripSubdetector::TID ||
            trajMeasurementHitId.subdetId() == StripSubdetector::TEC) {
          m_contains_TIDTEC = true;
        }
        // YP I add false here. No trajectory measurments in Muon system if we corrected TrackTransformer accordingly
      } else if (false && trajMeasurementHitId.det() == DetId::Muon) {
        //AR: I removed the false criteria for cosmic tests
        // } else if (trajMeasurementHitId.det() == DetId::Muon ) {

        if (m_debug)
          std::cout << "      TrajectoryMeasurement hit Det: Muon" << std::endl;

        if (trajMeasurementHitId.subdetId() == MuonSubdetId::DT) {
          const DTChamberId chamberId(trajMeasurementHitId.rawId());
          if (m_debug)
            std::cout << "        TrajectoryMeasurement hit subDet: DT wheel " << chamberId.wheel() << " station "
                      << chamberId.station() << " sector " << chamberId.sector() << std::endl;

          //          double gChX = globalGeometry->idToDet(chamberId)->position().x();
          //          double gChY = globalGeometry->idToDet(chamberId)->position().y();
          //          double gChZ = globalGeometry->idToDet(chamberId)->position().z();
          //          std::cout << "           The chamber position in global frame x: " << gChX << " y: " << gChY << " z: " << gChZ << std::endl;

          const GeomDet* geomDet = muonDetIdAssociator_->getGeomDet(chamberId);
          double chamber_width = geomDet->surface().bounds().width();
          double chamber_length = geomDet->surface().bounds().length();

          double hitX = trajMeasurementHit->localPosition().x();
          double hitY = trajMeasurementHit->localPosition().y();
          // double hitZ = trajMeasurementHit->localPosition().z();

          double tsosX = tsos.localPosition().x();
          double tsosY = tsos.localPosition().y();
          // double tsosZ = tsos.localPosition().z();

          int residualDT13IsAdded = false;
          int residualDT2IsAdded = false;

          // have we seen this chamber before?
          if (m_dt13.find(chamberId) == m_dt13.end() && m_dt2.find(chamberId) == m_dt2.end()) {
            if (m_debug)
              std::cout << "AR: pushing back chamber: " << chamberId << std::endl;
            m_chamberIds.push_back(chamberId);
            //addTrkCovMatrix(chamberId, tsos); // only for the 1st hit
          }
          if (m_debug)
            std::cout << "AR: size of chamberId: " << m_chamberIds.size() << std::endl;

          if (m_debug)
            std::cout << "        TrajectoryMeasurement hit dimension: " << trajMeasurementHitDim << std::endl;
          if (trajMeasurementHitDim > 1) {
            std::vector<const TrackingRecHit*> vDTSeg2D = trajMeasurementHit->recHits();
            if (m_debug)
              std::cout << "          vDTSeg2D size: " << vDTSeg2D.size() << std::endl;
            for (std::vector<const TrackingRecHit*>::const_iterator itDTSeg2D = vDTSeg2D.begin();
                 itDTSeg2D != vDTSeg2D.end();
                 ++itDTSeg2D) {
              std::vector<const TrackingRecHit*> vDTHits1D = (*itDTSeg2D)->recHits();
              if (m_debug)
                std::cout << "            vDTHits1D size: " << vDTHits1D.size() << std::endl;
              for (std::vector<const TrackingRecHit*>::const_iterator itDTHits1D = vDTHits1D.begin();
                   itDTHits1D != vDTHits1D.end();
                   ++itDTHits1D) {
                const TrackingRecHit* hit = *itDTHits1D;
                if (m_debug)
                  std::cout << "              hit dimension: " << hit->dimension() << std::endl;

                DetId hitId = hit->geographicalId();
                const DTSuperLayerId superLayerId(hitId.rawId());
                const DTLayerId layerId(hitId.rawId());
                if (m_debug)
                  std::cout << "              hit superLayerId: " << superLayerId.superLayer() << std::endl;
                if (m_debug)
                  std::cout << "              hit layerId: " << layerId.layer() << std::endl;

                if (superLayerId.superlayer() == 2 && vDTHits1D.size() >= 3) {
                  if (m_dt2.find(chamberId) == m_dt2.end()) {
                    AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
                    m_dt2[chamberId] =
                        new MuonDT2ChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
                    if (m_debug)
                      std::cout << "              This is first appearance of the DT with hits in superlayer 2"
                                << std::endl;
                  }
                  m_dt2[chamberId]->addResidual(prop, &tsos, hit, chamber_width, chamber_length);
                  residualDT2IsAdded = true;

                } else if ((superLayerId.superlayer() == 1 || superLayerId.superlayer() == 3) &&
                           vDTHits1D.size() >= 6) {
                  if (m_dt13.find(chamberId) == m_dt13.end()) {
                    AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
                    m_dt13[chamberId] =
                        new MuonDT13ChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
                    if (m_debug)
                      std::cout << "              This is first appearance of the DT with hits in superlayers 1 and 3"
                                << std::endl;
                  }
                  m_dt13[chamberId]->addResidual(prop, &tsos, hit, chamber_width, chamber_length);
                  residualDT13IsAdded = true;
                }
              }
            }
          }

          if (residualDT13IsAdded == true && residualDT2IsAdded == true && chamberId.wheel() == 0 &&
              chamberId.station() == 2 && chamberId.sector() == 7) {
            if (m_debug) {
              std::cout << "MYMARK " << tsosX << " " << hitX << " " << tsosX - hitX << " "
                        << m_dt13[chamberId]->trackx() << " " << m_dt13[chamberId]->residual() << " " << tsosY << " "
                        << hitY << " " << tsosY - hitY << " " << m_dt2[chamberId]->tracky() << " "
                        << m_dt2[chamberId]->residual() << " " << tsosF.localPosition().x() << " "
                        << tsosF.localPosition().y() << " " << tsosF.localPosition().z() << " "
                        << tsosB.localPosition().x() << " " << tsosB.localPosition().y() << " "
                        << tsosB.localPosition().z() << " " << tsosU.localPosition().x() << " "
                        << tsosU.localPosition().y() << " " << tsosU.localPosition().z() << std::endl;
            }
          }

          // http://cmslxr.fnal.gov/lxr/source/DataFormats/TrackReco/src/HitPattern.cc#101
          // YP I add false here. No trajectory measurments in Muon system if we corrected TrackTransformer accordingly
        } else if (false && trajMeasurementHitId.subdetId() == MuonSubdetId::CSC) {
          const CSCDetId cscDetId(trajMeasurementHitId.rawId());
          const CSCDetId chamberId2(cscDetId.endcap(), cscDetId.station(), cscDetId.ring(), cscDetId.chamber());
          if (m_debug)
            std::cout << "        TrajectoryMeasurement hit subDet: CSC endcap " << cscDetId.endcap() << " station "
                      << cscDetId.station() << " ring " << cscDetId.ring() << " chamber " << cscDetId.chamber()
                      << std::endl;
          if (m_debug)
            std::cout << "        TrajectoryMeasurement hit dimension: " << trajMeasurementHitDim << std::endl;

          if (trajMeasurementHitDim == 4) {
            std::vector<const TrackingRecHit*> vCSCHits2D = trajMeasurementHit->recHits();
            if (m_debug)
              std::cout << "          vCSCHits2D size: " << vCSCHits2D.size() << std::endl;
            if (vCSCHits2D.size() >= 5) {
              for (std::vector<const TrackingRecHit*>::const_iterator itCSCHits2D = vCSCHits2D.begin();
                   itCSCHits2D != vCSCHits2D.end();
                   ++itCSCHits2D) {
                const TrackingRecHit* cscHit2D = *itCSCHits2D;
                if (m_debug)
                  std::cout << "            cscHit2D dimension: " << cscHit2D->dimension() << std::endl;
                const TrackingRecHit* hit = cscHit2D;
                if (m_debug)
                  std::cout << "              hit dimension: " << hit->dimension() << std::endl;

                DetId hitId = hit->geographicalId();
                const CSCDetId cscDetId(hitId.rawId());
                if (m_debug)
                  std::cout << "              hit layer: " << cscDetId.layer() << std::endl;

                // not sure why we sometimes get layer == 0
                if (cscDetId.layer() == 0)
                  continue;

                // have we seen this chamber before?
                if (m_csc.find(chamberId2) == m_csc.end()) {
                  m_chamberIds.push_back(chamberId2);
                  //addTrkCovMatrix(chamberId, tsos); // only for the 1st hit
                  AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId2);
                  m_csc[chamberId2] =
                      new MuonCSCChamberResidual(globalGeometry, navigator, chamberId2, chamberAlignable);
                  if (m_debug)
                    std::cout << "              This is first appearance of the CSC with hits QQQ" << std::endl;
                }

                m_csc[chamberId2]->addResidual(prop, &tsos, hit, 250.0, 250.0);
              }
            }
          }
        } else {
          if (m_debug)
            std::cout << "        TrajectoryMeasurement hit subDet: UNKNOWN" << std::endl;
          if (m_debug)
            std::cout << "AR: trajMeasurementHitId.det(): " << trajMeasurementHitId.subdetId() << std::endl;
        }
      } else {
        if (m_debug)
          std::cout << "      TrajectoryMeasurement hit det: UNKNOWN" << std::endl;
        if (m_debug)
          std::cout << "AR: trajMeasurementHitId.det(): " << trajMeasurementHitId.det() << std::endl;
        if (m_debug)
          std::cout << "DetId::Tracker: " << DetId::Tracker << std::endl;
      }
    }
  }

  int iT2 = 0, iM2 = 0;
  for (auto const& hit2 : m_recoTrack->recHits()) {
    if (hit2->isValid()) {
      DetId hitId2 = hit2->geographicalId();
      if (hitId2.det() == DetId::Tracker) {
        iT2++;
        if (m_debug)
          std::cout << "Tracker Hit " << iT2 << " is found. We don't calcualte Tsos for it" << std::endl;
      } else if (hitId2.det() == DetId::Muon) {
        //        if ( (*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit ) {
        //          LogTrace("Reco|TrackingTools|TrackTransformer") << "RPC Rec Hit discarged";
        //          continue;
        //        }
        iM2++;
        if (m_debug)
          std::cout << "Muon Hit " << iM2 << " is found. Dimension: " << hit2->dimension() << std::endl;
        if (hitId2.subdetId() == MuonSubdetId::DT) {
          const DTChamberId chamberId(hitId2.rawId());
          if (m_debug)
            std::cout << "Muon Hit in DT wheel " << chamberId.wheel() << " station " << chamberId.station()
                      << " sector " << chamberId.sector() << std::endl;

          const GeomDet* geomDet = muonDetIdAssociator_->getGeomDet(chamberId);
          double chamber_width = geomDet->surface().bounds().width();
          double chamber_length = geomDet->surface().bounds().length();

          if (hit2->dimension() > 1) {
            // std::vector<const TrackingRecHit*> vDTSeg2D = hit2->recHits();
            std::vector<TrackingRecHit*> vDTSeg2D = hit2->recHits();

            if (m_debug)
              std::cout << "          vDTSeg2D size: " << vDTSeg2D.size() << std::endl;

            // for ( std::vector<const TrackingRecHit*>::const_iterator itDTSeg2D =  vDTSeg2D.begin();
            //                                                          itDTSeg2D != vDTSeg2D.end();
            //                                                        ++itDTSeg2D ) {

            for (std::vector<TrackingRecHit*>::const_iterator itDTSeg2D = vDTSeg2D.begin(); itDTSeg2D != vDTSeg2D.end();
                 ++itDTSeg2D) {
              // std::vector<const TrackingRecHit*> vDTHits1D =  (*itDTSeg2D)->recHits();
              std::vector<TrackingRecHit*> vDTHits1D = (*itDTSeg2D)->recHits();
              if (m_debug)
                std::cout << "            vDTHits1D size: " << vDTHits1D.size() << std::endl;
              // for ( std::vector<const TrackingRecHit*>::const_iterator itDTHits1D =  vDTHits1D.begin();
              //                                                          itDTHits1D != vDTHits1D.end();
              //                                                        ++itDTHits1D ) {
              for (std::vector<TrackingRecHit*>::const_iterator itDTHits1D = vDTHits1D.begin();
                   itDTHits1D != vDTHits1D.end();
                   ++itDTHits1D) {
                //const TrackingRecHit* hit = *itDTHits1D;
                TrackingRecHit* hit = *itDTHits1D;
                if (m_debug)
                  std::cout << "              hit dimension: " << hit->dimension() << std::endl;

                DetId hitId = hit->geographicalId();
                const DTSuperLayerId superLayerId(hitId.rawId());
                const DTLayerId layerId(hitId.rawId());
                if (m_debug)
                  std::cout << "              hit superLayerId: " << superLayerId.superLayer() << std::endl;
                if (m_debug)
                  std::cout << "              hit layerId: " << layerId.layer() << std::endl;

                if (superLayerId.superlayer() == 2 && vDTHits1D.size() >= 3) {
                  if (m_dt2.find(chamberId) == m_dt2.end()) {
                    AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
                    m_dt2[chamberId] =
                        new MuonDT2ChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
                    if (m_debug)
                      std::cout << "              This is first appearance of the DT with hits in superlayer 2"
                                << std::endl;

                    // have we seen this chamber before? check if it was in dt13
                    if (m_dt13.find(chamberId) == m_dt13.end()) {
                      m_chamberIds.push_back(chamberId);
                    }
                  }

                  TrajectoryStateOnSurface extrapolation;
                  extrapolation = prop->propagate(lastTrackerTsos, globalGeometry->idToDet(hitId)->surface());

                  if (extrapolation.isValid()) {
                    if (m_debug) {
                      std::cout << " extrapolation localPosition()"
                                << " x: " << extrapolation.localPosition().x()
                                << " y: " << extrapolation.localPosition().y()
                                << " z: " << extrapolation.localPosition().z() << std::endl;
                    }
                    m_dt2[chamberId]->addResidual(prop, &extrapolation, hit, chamber_width, chamber_length);
                  }
                  //                    residualDT2IsAdded = true;

                } else if ((superLayerId.superlayer() == 1 || superLayerId.superlayer() == 3) &&
                           vDTHits1D.size() >= 6) {
                  if (m_dt13.find(chamberId) == m_dt13.end()) {
                    AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
                    m_dt13[chamberId] =
                        new MuonDT13ChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
                    if (m_debug)
                      std::cout << "              This is first appearance of the DT with hits in superlayers 1 and 3"
                                << std::endl;

                    // have we seen this chamber before? check if it was in dt2
                    if (m_dt2.find(chamberId) == m_dt2.end()) {
                      m_chamberIds.push_back(chamberId);
                    }
                  }

                  TrajectoryStateOnSurface extrapolation;
                  extrapolation = prop->propagate(lastTrackerTsos, globalGeometry->idToDet(hitId)->surface());

                  if (extrapolation.isValid()) {
                    if (m_debug) {
                      std::cout << " extrapolation localPosition()"
                                << " x: " << extrapolation.localPosition().x()
                                << " y: " << extrapolation.localPosition().y()
                                << " z: " << extrapolation.localPosition().z() << std::endl;
                    }
                    m_dt13[chamberId]->addResidual(prop, &extrapolation, hit, chamber_width, chamber_length);
                  }
                  //                    residualDT13IsAdded = true;
                }
              }
            }
          }

          //          std::cout << "Extrapolate last Tracker TSOS to muon hit" << std::endl;
          //            TrajectoryStateOnSurface extrapolation;
          //            extrapolation = prop->propagate( lastTrackerTsos, globalGeometry->idToDet(hitId2)->surface() );
          //
          //            if ( chamberId2.wheel() == 0 && chamberId2.station() == 2 && chamberId2.sector() == 7 ) {
          //
          //            double hitX2 = hit2->localPosition().x();
          //          double hitY2 = hit2->localPosition().y();
          //          double hitZ2 = hit2->localPosition().z();
          //
          //          double tsosX2 = extrapolation.localPosition().x();
          //          double tsosY2 = extrapolation.localPosition().y();
          //          double tsosZ2 = extrapolation.localPosition().z();
          //
          //            std::cout << "MYMARK " << tsosX2 << " " << hitX2 << " " << tsosX2 - hitX2 << " " << "0" << " " << "0"
          //                            << " " << tsosY2 << " " << hitY2 << " " << tsosY2 - hitY2 << " " << "0"  << " " << "0"
          //                            << " 0 0 0 0 0 0 0 0 0 " << std::endl;
          //            }

        } else if (hitId2.subdetId() == MuonSubdetId::CSC) {
          const CSCDetId cscDetId2(hitId2.rawId());
          const CSCDetId chamberId(cscDetId2.endcap(), cscDetId2.station(), cscDetId2.ring(), cscDetId2.chamber());
          if (m_debug)
            std::cout << "Muon hit in CSC endcap " << cscDetId2.endcap() << " station " << cscDetId2.station()
                      << " ring " << cscDetId2.ring() << " chamber " << cscDetId2.chamber() << "." << std::endl;

          if (hit2->dimension() == 4) {
            // std::vector<const TrackingRecHit*> vCSCHits2D = hit2->recHits();
            std::vector<TrackingRecHit*> vCSCHits2D = hit2->recHits();
            if (m_debug)
              std::cout << "          vCSCHits2D size: " << vCSCHits2D.size() << std::endl;
            if (vCSCHits2D.size() >= 5) {
              // for ( std::vector<const TrackingRecHit*>::const_iterator itCSCHits2D =  vCSCHits2D.begin();
              //                                                          itCSCHits2D != vCSCHits2D.end();
              //                                                        ++itCSCHits2D ) {

              for (std::vector<TrackingRecHit*>::const_iterator itCSCHits2D = vCSCHits2D.begin();
                   itCSCHits2D != vCSCHits2D.end();
                   ++itCSCHits2D) {
                // const TrackingRecHit* cscHit2D = *itCSCHits2D;
                TrackingRecHit* cscHit2D = *itCSCHits2D;
                if (m_debug)
                  std::cout << "            cscHit2D dimension: " << cscHit2D->dimension() << std::endl;
                // const TrackingRecHit* hit = cscHit2D;
                TrackingRecHit* hit = cscHit2D;
                if (m_debug)
                  std::cout << "              hit dimension: " << hit->dimension() << std::endl;

                DetId hitId = hit->geographicalId();
                const CSCDetId cscDetId(hitId.rawId());

                if (m_debug) {
                  std::cout << "              hit layer: " << cscDetId.layer() << std::endl;

                  std::cout << " hit localPosition"
                            << " x: " << hit->localPosition().x() << " y: " << hit->localPosition().y()
                            << " z: " << hit->localPosition().z() << std::endl;
                  std::cout << " hit globalPosition"
                            << " x: " << globalGeometry->idToDet(hitId)->toGlobal(hit->localPosition()).x()
                            << " y: " << globalGeometry->idToDet(hitId)->toGlobal(hit->localPosition()).y()
                            << " z: " << globalGeometry->idToDet(hitId)->toGlobal(hit->localPosition()).z()
                            << std::endl;
                }

                // not sure why we sometimes get layer == 0
                if (cscDetId.layer() == 0)
                  continue;

                // have we seen this chamber before?
                if (m_debug)
                  std::cout << "Have we seen this chamber before?";
                if (m_csc.find(chamberId) == m_csc.end()) {
                  if (m_debug)
                    std::cout << " NO. m_csc.count() = " << m_csc.count(chamberId) << std::endl;
                  AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
                  m_csc[chamberId] = new MuonCSCChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
                  if (m_debug)
                    std::cout << "              This is first appearance of the CSC with hits m_csc.count() = "
                              << m_csc.count(chamberId) << std::endl;
                  m_chamberIds.push_back(chamberId);
                  //addTrkCovMatrix(chamberId, tsos); // only for the 1st hit
                } else {
                  if (m_debug)
                    std::cout << " YES. m_csc.count() = " << m_csc.count(chamberId) << std::endl;
                }

                if (m_debug) {
                  std::cout << " lastTrackerTsos localPosition"
                            << " x: " << lastTrackerTsos.localPosition().x()
                            << " y: " << lastTrackerTsos.localPosition().y()
                            << " z: " << lastTrackerTsos.localPosition().z() << std::endl;
                  std::cout << " lastTrackerTsos globalPosition"
                            << " x: " << lastTrackerTsos.globalPosition().x()
                            << " y: " << lastTrackerTsos.globalPosition().y()
                            << " z: " << lastTrackerTsos.globalPosition().z() << std::endl;
                  std::cout << " Do extrapolation from lastTrackerTsos to hit surface" << std::endl;
                }
                TrajectoryStateOnSurface extrapolation;
                extrapolation = prop->propagate(lastTrackerTsos, globalGeometry->idToDet(hitId)->surface());
                if (m_debug)
                  std::cout << " extrapolation.isValid() = " << extrapolation.isValid() << std::endl;

                if (extrapolation.isValid()) {
                  if (m_debug) {
                    std::cout << " extrapolation localPosition()"
                              << " x: " << extrapolation.localPosition().x()
                              << " y: " << extrapolation.localPosition().y()
                              << " z: " << extrapolation.localPosition().z() << std::endl;
                  }
                  m_csc[chamberId]->addResidual(prop, &extrapolation, hit, 250.0, 250.0);
                }
              }
            }
          }

        } else if (hitId2.subdetId() == MuonSubdetId::RPC) {
          if (m_debug)
            std::cout << "Muon Hit in RPC" << std::endl;
        } else {
          if (m_debug)
            std::cout << "Warning! Muon Hit not in DT or CSC or RPC" << std::endl;
        }
        //        recHitsForRefit.push_back(theMuonRecHitBuilder->build(&**hit));
        if (hitId2.subdetId() == MuonSubdetId::DT || hitId2.subdetId() == MuonSubdetId::CSC) {
        }
      }
    }
  }

  if (m_debug)
    std::cout << "END MuonResidualsFromTrack" << std::endl << std::endl;
}

MuonResidualsFromTrack() [2/2]

MuonResidualsFromTrack::MuonResidualsFromTrack	(	edm::ESHandle< GlobalTrackingGeometry >	globalGeometry,
		const reco::Muon *	recoMuon,
		AlignableNavigator *	navigator,
		double	maxResidual
	)

Definition at line 594 of file MuonResidualsFromTrack.cc.

References cms::cuda::assert(), reco::MuonSegmentMatch::BelongsToTrackByDR, reco::MuonSegmentMatch::BestInStationByDR, clear(), gather_cfg::cout, MuonSubdetId::CSC, MuonSubdetId::DT, edm::Ref< C, T, F >::get(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), reco::Muon::innerTrack(), edm::Ref< C, T, F >::isNonnull(), reco::Muon::isTrackerMuon(), m_chamberIds, m_csc, m_dt13, m_dt2, m_recoMuon, m_recoTrack, m_tracker_chi2, m_tracker_numHits, reco::Muon::matches(), DetId::Muon, and HLT_2022v12_cff::navigator.

    : m_recoMuon(recoMuon) {
  bool m_debug = false;

  clear();
  assert(m_recoMuon->isTrackerMuon() && m_recoMuon->innerTrack().isNonnull());
  m_recoTrack = m_recoMuon->innerTrack().get();

  m_tracker_chi2 = m_recoMuon->innerTrack()->chi2();
  m_tracker_numHits = m_recoMuon->innerTrack()->ndof() + 5;
  m_tracker_numHits = m_tracker_numHits > 0 ? m_tracker_numHits : 0;

  /*
                   for(auto const& hit : m_recoMuon->innerTrack()->recHits())
                   {
                   DetId id = hit->geographicalId();
                   if (id.det() == DetId::Tracker)
                   {
                   m_tracker_numHits++;
                   if (id.subdetId() == StripSubdetector::TID  ||  id.subdetId() == StripSubdetector::TEC) m_contains_TIDTEC = true;
                   }
                   }
                   */

  for (std::vector<reco::MuonChamberMatch>::const_iterator chamberMatch = m_recoMuon->matches().begin();
       chamberMatch != m_recoMuon->matches().end();
       chamberMatch++) {
    if (chamberMatch->id.det() != DetId::Muon)
      continue;

    for (std::vector<reco::MuonSegmentMatch>::const_iterator segMatch = chamberMatch->segmentMatches.begin();
         segMatch != chamberMatch->segmentMatches.end();
         ++segMatch) {
      // select the only segment that belongs to track and is the best in station by dR
      if (!(segMatch->isMask(reco::MuonSegmentMatch::BestInStationByDR) &&
            segMatch->isMask(reco::MuonSegmentMatch::BelongsToTrackByDR)))
        continue;

      if (chamberMatch->id.subdetId() == MuonSubdetId::DT) {
        const DTChamberId chamberId(chamberMatch->id.rawId());

        DTRecSegment4DRef segmentDT = segMatch->dtSegmentRef;
        const DTRecSegment4D* segment = segmentDT.get();
        if (segment == nullptr)
          continue;

        if (segment->hasPhi() && fabs(chamberMatch->x - segMatch->x) > maxResidual)
          continue;
        if (segment->hasZed() && fabs(chamberMatch->y - segMatch->y) > maxResidual)
          continue;

        // have we seen this chamber before?
        if (m_dt13.find(chamberId) == m_dt13.end() && m_dt2.find(chamberId) == m_dt2.end()) {
          m_chamberIds.push_back(chamberId);
        }

        if (segment->hasZed()) {
          if (m_dt2.find(chamberId) == m_dt2.end()) {
            AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
            // YP
            //            m_dt2[chamberId] = new MuonTrackDT2ChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
          } else if (m_debug)
            std::cout << "multi segment match to tmuon: dt2  -- should not happen!" << std::endl;
          m_dt2[chamberId]->setSegmentResidual(&(*chamberMatch), &(*segMatch));
        }
        if (segment->hasPhi()) {
          if (m_dt13.find(chamberId) == m_dt13.end()) {
            AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
            // YP
            //            m_dt13[chamberId] = new MuonTrackDT13ChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
          } else if (m_debug)
            std::cout << "multi segment match to tmuon: dt13  -- should not happen!" << std::endl;
          m_dt13[chamberId]->setSegmentResidual(&(*chamberMatch), &(*segMatch));
        }
      }

      else if (chamberMatch->id.subdetId() == MuonSubdetId::CSC) {
        const CSCDetId cscDetId(chamberMatch->id.rawId());
        const CSCDetId chamberId(cscDetId.chamberId());

        if (fabs(chamberMatch->x - segMatch->x) > maxResidual)
          continue;

        // have we seen this chamber before?
        if (m_csc.find(chamberId) == m_csc.end()) {
          m_chamberIds.push_back(chamberId);
          AlignableDetOrUnitPtr chamberAlignable = navigator->alignableFromDetId(chamberId);
          // YP
          //          m_csc[chamberId] = new MuonTrackCSCChamberResidual(globalGeometry, navigator, chamberId, chamberAlignable);
        } else if (m_debug)
          std::cout << "multi segment match to tmuon: csc  -- should not happen!" << std::endl;
        m_csc[chamberId]->setSegmentResidual(&(*chamberMatch), &(*segMatch));
      }
    }
  }
}

~MuonResidualsFromTrack()

MuonResidualsFromTrack::~MuonResidualsFromTrack

(

)

Definition at line 696 of file MuonResidualsFromTrack.cc.

References m_csc, m_dt13, and m_dt2.

                                                {
  for (std::map<DetId, MuonChamberResidual*>::const_iterator residual = m_dt13.begin(); residual != m_dt13.end();
       ++residual) {
    delete residual->second;
  }
  for (std::map<DetId, MuonChamberResidual*>::const_iterator residual = m_dt2.begin(); residual != m_dt2.end();
       ++residual) {
    delete residual->second;
  }
  for (std::map<DetId, MuonChamberResidual*>::const_iterator residual = m_csc.begin(); residual != m_csc.end();
       ++residual) {
    delete residual->second;
  }
}

Member Function Documentation

addTrkCovMatrix()

void MuonResidualsFromTrack::addTrkCovMatrix ( DetId  chamberId, TrajectoryStateOnSurface &  tsos  )

private

Definition at line 752 of file MuonResidualsFromTrack.cc.

References mps_fire::i, dqmiolumiharvest::j, TrajectoryStateOnSurface::localError(), m_trkCovMatrix, and LocalTrajectoryError::matrix().

                                                                                            {
  const AlgebraicSymMatrix55 cov55 = tsos.localError().matrix();
  TMatrixDSym cov44(4);
  // change indices from q/p,dxdz,dydz,x,y   to   x,y,dxdz,dydz
  int subs[4] = {3, 4, 1, 2};
  for (int i = 0; i < 4; i++)
    for (int j = 0; j < 4; j++)
      cov44(i, j) = cov55(subs[i], subs[j]);
  m_trkCovMatrix[chamberId] = cov44;
}

builderESInputTag()

edm::ESInputTag MuonResidualsFromTrack::builderESInputTag ( )

static

Definition at line 22 of file MuonResidualsFromTrack.cc.

{ return edm::ESInputTag("", "WithTrackAngle"); }

chamberIds()

const std::vector<DetId> MuonResidualsFromTrack::chamberIds ( ) const

inline

Definition at line 89 of file MuonResidualsFromTrack.h.

References m_chamberIds.

Referenced by AlignmentMonitorMuonSystemMap1D::processMuonResidualsFromTrack(), AlignmentMonitorMuonVsCurvature::processMuonResidualsFromTrack(), AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack(), and MuonAlignmentFromReference::processMuonResidualsFromTrack().

{ return m_chamberIds; }

chamberResidual()

MuonChamberResidual * MuonResidualsFromTrack::chamberResidual	(	DetId	chamberId,
		int	type
	)

Definition at line 735 of file MuonResidualsFromTrack.cc.

References MuonChamberResidual::kCSC, MuonChamberResidual::kDT13, MuonChamberResidual::kDT2, m_csc, m_dt13, and m_dt2.

Referenced by AlignmentMonitorMuonSystemMap1D::processMuonResidualsFromTrack(), AlignmentMonitorMuonVsCurvature::processMuonResidualsFromTrack(), AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack(), and MuonAlignmentFromReference::processMuonResidualsFromTrack().

                                                                                      {
  if (type == MuonChamberResidual::kDT13) {
    if (m_dt13.find(chamberId) == m_dt13.end())
      return nullptr;
    return m_dt13[chamberId];
  } else if (type == MuonChamberResidual::kDT2) {
    if (m_dt2.find(chamberId) == m_dt2.end())
      return nullptr;
    return m_dt2[chamberId];
  } else if (type == MuonChamberResidual::kCSC) {
    if (m_csc.find(chamberId) == m_csc.end())
      return nullptr;
    return m_csc[chamberId];
  } else
    return nullptr;
}

choleskyCorrMatrix()

TMatrixD MuonResidualsFromTrack::choleskyCorrMatrix

(

DetId

chamberId

)

Definition at line 826 of file MuonResidualsFromTrack.cc.

References corrMatrix(), gather_cfg::cout, convertSQLiteXML::ok, and mps_fire::result.

                                                                   {
  bool m_debug = false;

  TMatrixD result(4, 4);
  TMatrixDSym corr44 = corrMatrix(chamberId);

  // get an upper triangular matrix U such that corr = U^T * U
  TDecompChol decomp(corr44);
  bool ok = decomp.Decompose();
  result = decomp.GetU();

  if (m_debug)
    std::cout << "MuonResidualsFromTrack:: corr cholesky after:" << std::endl;
  result.Print();

  if (!ok && m_debug)
    std::cout << "MuonResidualsFromTrack:: corr decomposition failed!" << std::endl;
  return result;
}

clear()

void MuonResidualsFromTrack::clear

(

)

Definition at line 711 of file MuonResidualsFromTrack.cc.

References m_chamberIds, m_contains_TIDTEC, m_csc, m_dt13, m_dt2, m_tracker_chi2, m_tracker_numHits, and m_trkCovMatrix.

Referenced by MuonResidualsFromTrack(), and BeautifulSoup.Tag::setString().

                                   {
  m_tracker_numHits = 0;
  m_tracker_chi2 = 0.;
  m_contains_TIDTEC = false;
  m_chamberIds.clear();
  m_dt13.clear();
  m_dt2.clear();
  m_csc.clear();
  m_trkCovMatrix.clear();
}

contains_TIDTEC()

bool MuonResidualsFromTrack::contains_TIDTEC ( ) const

inline

Definition at line 87 of file MuonResidualsFromTrack.h.

References m_contains_TIDTEC.

Referenced by AlignmentMonitorMuonSystemMap1D::processMuonResidualsFromTrack(), AlignmentMonitorMuonVsCurvature::processMuonResidualsFromTrack(), AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack(), and MuonAlignmentFromReference::processMuonResidualsFromTrack().

{ return m_contains_TIDTEC; }

corrMatrix()

TMatrixDSym MuonResidualsFromTrack::corrMatrix

(

DetId

chamberId

)

Definition at line 808 of file MuonResidualsFromTrack.cc.

References gather_cfg::cout, covMatrix(), convertSQLiteXML::ok, and mps_fire::result.

Referenced by choleskyCorrMatrix().

                                                              {
  bool m_debug = false;

  TMatrixDSym result(4);
  TMatrixDSym cov44 = covMatrix(chamberId);

  // invert it using cholesky decomposition
  TDecompChol decomp(cov44);
  bool ok = decomp.Invert(result);
  if (m_debug)
    std::cout << "MuonResidualsFromTrack:: corr after:" << std::endl;
  result.Print();

  if (!ok && m_debug)
    std::cout << "MuonResidualsFromTrack:: cov inversion failed!" << std::endl;
  return result;
}

covMatrix()

TMatrixDSym MuonResidualsFromTrack::covMatrix

(

DetId

chamberId

)

Definition at line 763 of file MuonResidualsFromTrack.cc.

References gather_cfg::cout, m_csc, m_dt13, m_dt2, m_trkCovMatrix, and mps_fire::result.

Referenced by corrMatrix().

                                                             {
  bool m_debug = false;

  TMatrixDSym result(4);
  if (m_debug)
    std::cout << "MuonResidualsFromTrack:: cov initial:" << std::endl;
  result.Print();
  if (m_trkCovMatrix.find(chamberId) == m_trkCovMatrix.end()) {
    if (m_debug)
      std::cout << "MuonResidualsFromTrack:: cov does not exist!" << std::endl;
    return result;
  }
  result = m_trkCovMatrix[chamberId];

  if (m_debug)
    std::cout << "MuonResidualsFromTrack:: cov before:" << std::endl;
  result.Print();

  // add segment's errors in quadratures to track's covariance matrix
  double r_err;
  if (m_csc.find(chamberId) == m_csc.end()) {
    r_err = m_csc[chamberId]->residual_error();
    result(0, 0) += r_err * r_err;
    r_err = m_csc[chamberId]->resslope_error();
    result(2, 2) += r_err * r_err;
  }
  if (m_dt13.find(chamberId) == m_dt13.end()) {
    r_err = m_dt13[chamberId]->residual_error();
    result(0, 0) += r_err * r_err;
    r_err = m_dt13[chamberId]->resslope_error();
    result(2, 2) += r_err * r_err;
  }
  if (m_dt2.find(chamberId) == m_dt2.end()) {
    r_err = m_dt2[chamberId]->residual_error();
    result(1, 1) += r_err * r_err;
    r_err = m_dt2[chamberId]->resslope_error();
    result(3, 3) += r_err * r_err;
  }
  if (m_debug)
    std::cout << "MuonResidualsFromTrack:: cov after:" << std::endl;
  result.Print();

  return result;
}

getMuon()

const reco::Muon* MuonResidualsFromTrack::getMuon ( )

inline

Definition at line 79 of file MuonResidualsFromTrack.h.

References m_recoMuon.

{ return m_recoMuon; }

getTrack()

const reco::Track* MuonResidualsFromTrack::getTrack ( )

inline

Definition at line 78 of file MuonResidualsFromTrack.h.

References m_recoTrack.

Referenced by AlignmentMonitorMuonSystemMap1D::processMuonResidualsFromTrack(), AlignmentMonitorMuonVsCurvature::processMuonResidualsFromTrack(), AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack(), and MuonAlignmentFromReference::processMuonResidualsFromTrack().

{ return m_recoTrack; }

normalizedChi2()

double MuonResidualsFromTrack::normalizedChi2

(

)

const

Definition at line 729 of file MuonResidualsFromTrack.cc.

References m_recoMuon, m_recoTrack, reco::TrackBase::normalizedChi2(), and trackerRedChi2().

Referenced by AlignmentMonitorMuonSystemMap1D::processMuonResidualsFromTrack(), AlignmentMonitorMuonVsCurvature::processMuonResidualsFromTrack(), AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack(), and MuonAlignmentFromReference::processMuonResidualsFromTrack().

                                                    {
  if (m_recoMuon)
    return m_recoTrack->normalizedChi2();
  return trackerRedChi2();
}

trackerChi2()

double MuonResidualsFromTrack::trackerChi2 ( ) const

inline

Definition at line 83 of file MuonResidualsFromTrack.h.

References m_tracker_chi2.

{ return m_tracker_chi2; }

trackerNumHits()

int MuonResidualsFromTrack::trackerNumHits ( ) const

inline

Definition at line 81 of file MuonResidualsFromTrack.h.

References m_tracker_numHits.

Referenced by AlignmentMonitorMuonSystemMap1D::processMuonResidualsFromTrack(), AlignmentMonitorMuonVsCurvature::processMuonResidualsFromTrack(), AlignmentMonitorSegmentDifferences::processMuonResidualsFromTrack(), and MuonAlignmentFromReference::processMuonResidualsFromTrack().

{ return m_tracker_numHits; }

trackerRedChi2()

double MuonResidualsFromTrack::trackerRedChi2

(

)

const

Definition at line 722 of file MuonResidualsFromTrack.cc.

References m_tracker_chi2, and m_tracker_numHits.

Referenced by normalizedChi2(), and AlignmentMonitorMuonVsCurvature::processMuonResidualsFromTrack().

                                                    {
  if (m_tracker_numHits > 5)
    return m_tracker_chi2 / double(m_tracker_numHits - 5);
  else
    return -1.;
}

Member Data Documentation

m_chamberIds

std::vector<DetId> MuonResidualsFromTrack::m_chamberIds

private

Definition at line 104 of file MuonResidualsFromTrack.h.

Referenced by chamberIds(), clear(), and MuonResidualsFromTrack().

m_contains_TIDTEC

bool MuonResidualsFromTrack::m_contains_TIDTEC

private

Definition at line 102 of file MuonResidualsFromTrack.h.

Referenced by clear(), contains_TIDTEC(), and MuonResidualsFromTrack().

m_csc

std::map<DetId, MuonChamberResidual *> MuonResidualsFromTrack::m_csc

private

Definition at line 105 of file MuonResidualsFromTrack.h.

Referenced by chamberResidual(), clear(), covMatrix(), MuonResidualsFromTrack(), and ~MuonResidualsFromTrack().

m_dt13

std::map<DetId, MuonChamberResidual *> MuonResidualsFromTrack::m_dt13

private

Definition at line 105 of file MuonResidualsFromTrack.h.

Referenced by chamberResidual(), clear(), covMatrix(), MuonResidualsFromTrack(), and ~MuonResidualsFromTrack().

m_dt2

std::map<DetId, MuonChamberResidual *> MuonResidualsFromTrack::m_dt2

private

Definition at line 105 of file MuonResidualsFromTrack.h.

Referenced by chamberResidual(), clear(), covMatrix(), MuonResidualsFromTrack(), and ~MuonResidualsFromTrack().

m_recoMuon

const reco::Muon* MuonResidualsFromTrack::m_recoMuon

private

Definition at line 114 of file MuonResidualsFromTrack.h.

Referenced by getMuon(), MuonResidualsFromTrack(), and normalizedChi2().

m_recoTrack

const reco::Track* MuonResidualsFromTrack::m_recoTrack

private

Definition at line 111 of file MuonResidualsFromTrack.h.

Referenced by getTrack(), MuonResidualsFromTrack(), and normalizedChi2().

m_tracker_chi2

double MuonResidualsFromTrack::m_tracker_chi2

private

Definition at line 101 of file MuonResidualsFromTrack.h.

Referenced by clear(), MuonResidualsFromTrack(), trackerChi2(), and trackerRedChi2().

m_tracker_numHits

int MuonResidualsFromTrack::m_tracker_numHits

private

Definition at line 100 of file MuonResidualsFromTrack.h.

Referenced by clear(), MuonResidualsFromTrack(), trackerNumHits(), and trackerRedChi2().

m_trkCovMatrix

std::map<DetId, TMatrixDSym> MuonResidualsFromTrack::m_trkCovMatrix

private

Definition at line 106 of file MuonResidualsFromTrack.h.

Referenced by addTrkCovMatrix(), clear(), and covMatrix().

m_tsoscomb

TrajectoryStateCombiner MuonResidualsFromTrack::m_tsoscomb

private

Definition at line 98 of file MuonResidualsFromTrack.h.

Referenced by MuonResidualsFromTrack().