GlobalMuonRefitter Class Reference

#include <GlobalMuonRefitter.h>

Public Types
typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer	ConstMuonRecHitContainer
typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer	ConstMuonRecHitPointer
typedef TransientTrackingRecHit::ConstRecHitContainer	ConstRecHitContainer
typedef TransientTrackingRecHit::ConstRecHitPointer	ConstRecHitPointer
typedef MuonTransientTrackingRecHit::MuonRecHitContainer	MuonRecHitContainer
typedef MuonTransientTrackingRecHit::MuonRecHitPointer	MuonRecHitPointer
typedef TransientTrackingRecHit::RecHitContainer	RecHitContainer
typedef TransientTrackingRecHit::RecHitPointer	RecHitPointer
enum	subDetector {   PXB = 1, PXF = 2, TIB = 3, TID = 4,   TOB = 5, TEC = 6 }
typedef std::vector< Trajectory >	TC
typedef TC::const_iterator	TI

Public Member Functions
const reco::DYTInfo *	getDYTInfo ()
ConstRecHitContainer	getRidOfSelectStationHits (const ConstRecHitContainer &hits, const TrackerTopology *tTopo) const
	GlobalMuonRefitter (const edm::ParameterSet &, const MuonServiceProxy *, edm::ConsumesCollector &)
	constructor with Parameter Set and MuonServiceProxy More...
std::vector< Trajectory >	refit (const reco::Track &globalTrack, const int theMuonHitsOption, const TrackerTopology *tTopo) const
	build combined trajectory from sta Track and tracker RecHits More...
std::vector< Trajectory >	refit (const reco::Track &globalTrack, const reco::TransientTrack track, const TransientTrackingRecHit::ConstRecHitContainer &allRecHitsTemp, const int theMuonHitsOption, const TrackerTopology *tTopo) const
	build combined trajectory from subset of sta Track and tracker RecHits More...
virtual void	setEvent (const edm::Event &)
	pass the Event to the algo at each event More...
void	setServices (const edm::EventSetup &)
	set the services needed by the TrackTransformer More...
std::vector< Trajectory >	transform (const reco::Track &newTrack, const reco::TransientTrack track, const TransientTrackingRecHit::ConstRecHitContainer &recHitsForReFit) const
	refit the track with a new set of RecHits More...
virtual	~GlobalMuonRefitter ()
	destructor More...

Protected Types
enum	RefitDirection { insideOut, outsideIn, undetermined }

Protected Member Functions
void	checkMuonHits (const reco::Track &, ConstRecHitContainer &, std::map< DetId, int > &) const
	check muon RecHits, calculate chamber occupancy and select hits to be used in the final fit More...
RefitDirection	checkRecHitsOrdering (const ConstRecHitContainer &) const
void	getFirstHits (const reco::Track &, ConstRecHitContainer &, ConstRecHitContainer &) const
	get the RecHits in the tracker and the first muon chamber with hits More...
void	printHits (const ConstRecHitContainer &) const
	print all RecHits of a trajectory More...
ConstRecHitContainer	selectMuonHits (const Trajectory &, const std::map< DetId, int > &) const
	select muon hits compatible with trajectory; check hits in chambers with showers More...
const MuonServiceProxy *	service () const

Protected Attributes
std::string	theCategory
float	thePtCut
bool	theTkTrajsAvailableFlag

Private Attributes
edm::Handle< DTRecSegment4DCollection >	all4DSegments
edm::EDGetTokenT< DTRecSegment4DCollection >	all4DSegmentsToken
edm::Handle< CSCSegmentCollection >	CSCSegments
edm::EDGetTokenT< CSCSegmentCollection >	CSCSegmentsToken
reco::DYTInfo *	dytInfo
TkClonerImpl	hitCloner
unsigned long long	theCacheId_TRH
bool	theCosmicFlag
float	theCSCChi2Cut
edm::InputTag	theCSCRecHitLabel
edm::Handle< CSCRecHit2DCollection >	theCSCRecHits
edm::EDGetTokenT< CSCRecHit2DCollection >	theCSCRecHitToken
float	theDTChi2Cut
edm::InputTag	theDTRecHitLabel
edm::Handle< DTRecHitCollection >	theDTRecHits
edm::EDGetTokenT< DTRecHitCollection >	theDTRecHitToken
bool	theDYTParThrsMode
int	theDYTselector
std::vector< int >	theDYTthrs
edm::ParameterSet	theDYTthrsParameters
bool	theDYTupdator
bool	theDYTuseAPE
const edm::Event *	theEvent
std::unique_ptr< TrajectoryFitter >	theFitter
std::string	theFitterName
float	theGEMChi2Cut
edm::InputTag	theGEMRecHitLabel
edm::Handle< GEMRecHitCollection >	theGEMRecHits
edm::EDGetTokenT< GEMRecHitCollection >	theGEMRecHitToken
int	theHitThreshold
float	theME0Chi2Cut
edm::InputTag	theME0RecHitLabel
edm::Handle< ME0SegmentCollection >	theME0RecHits
edm::EDGetTokenT< ME0SegmentCollection >	theME0RecHitToken
int	theMuonHitsOption
edm::ESHandle< TransientTrackingRecHitBuilder >	theMuonRecHitBuilder
std::string	theMuonRecHitBuilderName
float	theProbCut
std::string	thePropagatorName
RefitDirection	theRefitDirection
double	theRescaleErrorFactor
float	theRPCChi2Cut
bool	theRPCInTheFit
const MuonServiceProxy *	theService
int	theSkipStation
edm::ESHandle< TransientTrackingRecHitBuilder >	theTrackerRecHitBuilder
std::string	theTrackerRecHitBuilderName
int	theTrackerSkipSection
int	theTrackerSkipSystem

Detailed Description

class to build muon trajectory

Author

N. Neumeister Purdue University

C. Liu Purdue University

A. Everett Purdue University

\modified by C. Calabria INFN & Universita Bari \modified by D. Nash Northeastern University \modified by C. Caputo UCLouvain

Definition at line 52 of file GlobalMuonRefitter.h.

Member Typedef Documentation

ConstMuonRecHitContainer

typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer GlobalMuonRefitter::ConstMuonRecHitContainer

Definition at line 62 of file GlobalMuonRefitter.h.

ConstMuonRecHitPointer

typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer GlobalMuonRefitter::ConstMuonRecHitPointer

Definition at line 60 of file GlobalMuonRefitter.h.

ConstRecHitContainer

typedef TransientTrackingRecHit::ConstRecHitContainer GlobalMuonRefitter::ConstRecHitContainer

Definition at line 55 of file GlobalMuonRefitter.h.

ConstRecHitPointer

typedef TransientTrackingRecHit::ConstRecHitPointer GlobalMuonRefitter::ConstRecHitPointer

Definition at line 57 of file GlobalMuonRefitter.h.

MuonRecHitContainer

typedef MuonTransientTrackingRecHit::MuonRecHitContainer GlobalMuonRefitter::MuonRecHitContainer

Definition at line 61 of file GlobalMuonRefitter.h.

MuonRecHitPointer

typedef MuonTransientTrackingRecHit::MuonRecHitPointer GlobalMuonRefitter::MuonRecHitPointer

Definition at line 59 of file GlobalMuonRefitter.h.

RecHitContainer

typedef TransientTrackingRecHit::RecHitContainer GlobalMuonRefitter::RecHitContainer

Definition at line 54 of file GlobalMuonRefitter.h.

RecHitPointer

typedef TransientTrackingRecHit::RecHitPointer GlobalMuonRefitter::RecHitPointer

Definition at line 56 of file GlobalMuonRefitter.h.

TC

typedef std::vector<Trajectory> GlobalMuonRefitter::TC

Definition at line 64 of file GlobalMuonRefitter.h.

TI

typedef TC::const_iterator GlobalMuonRefitter::TI

Definition at line 65 of file GlobalMuonRefitter.h.

Member Enumeration Documentation

RefitDirection

enum GlobalMuonRefitter::RefitDirection

protected

Enumerator
insideOut
outsideIn
undetermined

Definition at line 106 of file GlobalMuonRefitter.h.

{ insideOut, outsideIn, undetermined };

subDetector

enum GlobalMuonRefitter::subDetector

Enumerator
PXB
PXF
TIB
TID
TOB
TEC

Definition at line 67 of file GlobalMuonRefitter.h.

{ PXB = 1, PXF = 2, TIB = 3, TID = 4, TOB = 5, TEC = 6 };

Constructor & Destructor Documentation

GlobalMuonRefitter()

GlobalMuonRefitter::GlobalMuonRefitter	(	const edm::ParameterSet &	par,
		const MuonServiceProxy *	service,
		edm::ConsumesCollector &	iC
	)

constructor with Parameter Set and MuonServiceProxy

Definition at line 75 of file GlobalMuonRefitter.cc.

    : theCosmicFlag(par.getParameter<bool>("PropDirForCosmics")),
      theDTRecHitLabel(par.getParameter<InputTag>("DTRecSegmentLabel")),
      theCSCRecHitLabel(par.getParameter<InputTag>("CSCRecSegmentLabel")),
      theGEMRecHitLabel(par.getParameter<InputTag>("GEMRecHitLabel")),
      theME0RecHitLabel(par.getParameter<InputTag>("ME0RecHitLabel")),
      theService(service) {
  theCategory = par.getUntrackedParameter<string>("Category", "Muon|RecoMuon|GlobalMuon|GlobalMuonRefitter");
 
  theHitThreshold = par.getParameter<int>("HitThreshold");
  theDTChi2Cut = par.getParameter<double>("Chi2CutDT");
  theCSCChi2Cut = par.getParameter<double>("Chi2CutCSC");
  theRPCChi2Cut = par.getParameter<double>("Chi2CutRPC");
  theGEMChi2Cut = par.getParameter<double>("Chi2CutGEM");
  theME0Chi2Cut = par.getParameter<double>("Chi2CutME0");
 
  // Refit direction
  string refitDirectionName = par.getParameter<string>("RefitDirection");
 
  if (refitDirectionName == "insideOut")
    theRefitDirection = insideOut;
  else if (refitDirectionName == "outsideIn")
    theRefitDirection = outsideIn;
  else
    throw cms::Exception("TrackTransformer constructor")
        << "Wrong refit direction chosen in TrackTransformer ParameterSet"
        << "\n"
        << "Possible choices are:"
        << "\n"
        << "RefitDirection = insideOut or RefitDirection = outsideIn";
 
  theFitterName = par.getParameter<string>("Fitter");
  thePropagatorName = par.getParameter<string>("Propagator");
 
  theSkipStation = par.getParameter<int>("SkipStation");
  theTrackerSkipSystem = par.getParameter<int>("TrackerSkipSystem");
  theTrackerSkipSection = par.getParameter<int>("TrackerSkipSection");  //layer, wheel, or disk depending on the system
 
  theTrackerRecHitBuilderName = par.getParameter<string>("TrackerRecHitBuilder");
  theMuonRecHitBuilderName = par.getParameter<string>("MuonRecHitBuilder");
 
  theRPCInTheFit = par.getParameter<bool>("RefitRPCHits");
 
  theDYTthrs = par.getParameter<std::vector<int> >("DYTthrs");
  theDYTselector = par.getParameter<int>("DYTselector");
  theDYTupdator = par.getParameter<bool>("DYTupdator");
  theDYTuseAPE = par.getParameter<bool>("DYTuseAPE");
  theDYTParThrsMode = par.getParameter<bool>("DYTuseThrsParametrization");
  if (theDYTParThrsMode)
    theDYTthrsParameters = par.getParameter<edm::ParameterSet>("DYTthrsParameters");
  dytInfo = new reco::DYTInfo();
 
  if (par.existsAs<double>("RescaleErrorFactor")) {
    theRescaleErrorFactor = par.getParameter<double>("RescaleErrorFactor");
    edm::LogWarning("GlobalMuonRefitter") << "using error rescale factor " << theRescaleErrorFactor;
  } else
    theRescaleErrorFactor = 1000.;
 
  theCacheId_TRH = 0;
  theDTRecHitToken = iC.consumes<DTRecHitCollection>(theDTRecHitLabel);
  theCSCRecHitToken = iC.consumes<CSCRecHit2DCollection>(theCSCRecHitLabel);
  theGEMRecHitToken = iC.consumes<GEMRecHitCollection>(theGEMRecHitLabel);
  theME0RecHitToken = iC.consumes<ME0SegmentCollection>(theME0RecHitLabel);
  CSCSegmentsToken = iC.consumes<CSCSegmentCollection>(InputTag("cscSegments"));
  all4DSegmentsToken = iC.consumes<DTRecSegment4DCollection>(InputTag("dt4DSegments"));
}

References all4DSegmentsToken, edm::ConsumesCollector::consumes(), CSCSegmentsToken, dytInfo, Exception, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2018_cff::InputTag, insideOut, outsideIn, theCacheId_TRH, theCategory, theCSCChi2Cut, theCSCRecHitLabel, theCSCRecHitToken, theDTChi2Cut, theDTRecHitLabel, theDTRecHitToken, theDYTParThrsMode, theDYTselector, theDYTthrs, theDYTthrsParameters, theDYTupdator, theDYTuseAPE, theFitterName, theGEMChi2Cut, theGEMRecHitLabel, theGEMRecHitToken, theHitThreshold, theME0Chi2Cut, theME0RecHitLabel, theME0RecHitToken, theMuonRecHitBuilderName, thePropagatorName, theRefitDirection, theRescaleErrorFactor, theRPCChi2Cut, theRPCInTheFit, theSkipStation, theTrackerRecHitBuilderName, theTrackerSkipSection, and theTrackerSkipSystem.

~GlobalMuonRefitter()

GlobalMuonRefitter::~GlobalMuonRefitter ( )

virtual

destructor

Definition at line 148 of file GlobalMuonRefitter.cc.

{ delete dytInfo; }

References dytInfo.

Member Function Documentation

checkMuonHits()

void GlobalMuonRefitter::checkMuonHits ( const reco::Track &  muon, ConstRecHitContainer &  all, std::map< DetId, int > &  hitMap  ) const

protected

check muon RecHits, calculate chamber occupancy and select hits to be used in the final fit

Definition at line 306 of file GlobalMuonRefitter.cc.

                                                                      {
  LogTrace(theCategory) << " GlobalMuonRefitter::checkMuonHits " << endl;
 
  float coneSize = 20.0;
 
  // loop through all muon hits and calculate the maximum # of hits in each chamber
  for (ConstRecHitContainer::const_iterator imrh = all.begin(); imrh != all.end(); imrh++) {
    if ((*imrh != nullptr) && !(*imrh)->isValid())
      continue;
 
    int detRecHits = 0;
    MuonRecHitContainer dRecHits;
 
    DetId id = (*imrh)->geographicalId();
    DetId chamberId;
 
    // Skip tracker hits
    if (id.det() != DetId::Muon)
      continue;
 
    if (id.subdetId() == MuonSubdetId::DT) {
      DTChamberId did(id.rawId());
      chamberId = did;
 
      if ((*imrh)->dimension() > 1) {
        std::vector<const TrackingRecHit*> hits2d = (*imrh)->recHits();
        for (std::vector<const TrackingRecHit*>::const_iterator hit2d = hits2d.begin(); hit2d != hits2d.end();
             hit2d++) {
          if ((*hit2d)->dimension() > 1) {
            std::vector<const TrackingRecHit*> hits1d = (*hit2d)->recHits();
            for (std::vector<const TrackingRecHit*>::const_iterator hit1d = hits1d.begin(); hit1d != hits1d.end();
                 hit1d++) {
              DetId id1 = (*hit1d)->geographicalId();
              DTLayerId lid(id1.rawId());
              // Get the 1d DT RechHits from this layer
              DTRecHitCollection::range dRecHits = theDTRecHits->get(lid);
              int layerHits = 0;
              for (DTRecHitCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
                double rhitDistance = fabs(ir->localPosition().x() - (**hit1d).localPosition().x());
                if (rhitDistance < coneSize)
                  layerHits++;
                LogTrace(theCategory) << "       " << (ir)->localPosition() << "  " << (**hit1d).localPosition()
                                      << " Distance: " << rhitDistance << " recHits: " << layerHits
                                      << "  SL: " << lid.superLayer() << endl;
              }
              if (layerHits > detRecHits)
                detRecHits = layerHits;
            }
          } else {
            DTLayerId lid(id.rawId());
            // Get the 1d DT RechHits from this layer
            DTRecHitCollection::range dRecHits = theDTRecHits->get(lid);
            for (DTRecHitCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
              double rhitDistance = fabs(ir->localPosition().x() - (**imrh).localPosition().x());
              if (rhitDistance < coneSize)
                detRecHits++;
              LogTrace(theCategory) << "       " << (ir)->localPosition() << "  " << (**imrh).localPosition()
                                    << " Distance: " << rhitDistance << " recHits: " << detRecHits << endl;
            }
          }
        }
 
      } else {
        DTLayerId lid(id.rawId());
 
        // Get the 1d DT RechHits from this layer
        DTRecHitCollection::range dRecHits = theDTRecHits->get(lid);
 
        for (DTRecHitCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
          double rhitDistance = fabs(ir->localPosition().x() - (**imrh).localPosition().x());
          if (rhitDistance < coneSize)
            detRecHits++;
          LogTrace(theCategory) << "       " << (ir)->localPosition() << "  " << (**imrh).localPosition()
                                << " Distance: " << rhitDistance << " recHits: " << detRecHits << endl;
        }
      }
    }  // end of if DT
    else if (id.subdetId() == MuonSubdetId::CSC) {
      CSCDetId did(id.rawId());
      chamberId = did.chamberId();
 
      if ((*imrh)->recHits().size() > 1) {
        std::vector<const TrackingRecHit*> hits2d = (*imrh)->recHits();
        for (std::vector<const TrackingRecHit*>::const_iterator hit2d = hits2d.begin(); hit2d != hits2d.end();
             hit2d++) {
          DetId id1 = (*hit2d)->geographicalId();
          CSCDetId lid(id1.rawId());
 
          // Get the CSC Rechits from this layer
          CSCRecHit2DCollection::range dRecHits = theCSCRecHits->get(lid);
          int layerHits = 0;
 
          for (CSCRecHit2DCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
            double rhitDistance = (ir->localPosition() - (**hit2d).localPosition()).mag();
            if (rhitDistance < coneSize)
              layerHits++;
            LogTrace(theCategory) << ir->localPosition() << "  " << (**hit2d).localPosition()
                                  << " Distance: " << rhitDistance << " recHits: " << layerHits << endl;
          }
          if (layerHits > detRecHits)
            detRecHits = layerHits;
        }
      } else {
        // Get the CSC Rechits from this layer
        CSCRecHit2DCollection::range dRecHits = theCSCRecHits->get(did);
 
        for (CSCRecHit2DCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
          double rhitDistance = (ir->localPosition() - (**imrh).localPosition()).mag();
          if (rhitDistance < coneSize)
            detRecHits++;
          LogTrace(theCategory) << ir->localPosition() << "  " << (**imrh).localPosition()
                                << " Distance: " << rhitDistance << " recHits: " << detRecHits << endl;
        }
      }
    }  //end of CSC if
    else if (id.subdetId() == MuonSubdetId::GEM) {
      GEMDetId did(id.rawId());
      chamberId = did.chamberId();
 
      if ((*imrh)->recHits().size() > 1) {
        std::vector<const TrackingRecHit*> hits2d = (*imrh)->recHits();
        for (std::vector<const TrackingRecHit*>::const_iterator hit2d = hits2d.begin(); hit2d != hits2d.end();
             hit2d++) {
          DetId id1 = (*hit2d)->geographicalId();
          GEMDetId lid(id1.rawId());
 
          // Get the GEM Rechits from this layer
          GEMRecHitCollection::range dRecHits = theGEMRecHits->get(lid);
          int layerHits = 0;
 
          for (GEMRecHitCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
            double rhitDistance = (ir->localPosition() - (**hit2d).localPosition()).mag();
            if (rhitDistance < coneSize)
              layerHits++;
            LogTrace(theCategory) << ir->localPosition() << "  " << (**hit2d).localPosition()
                                  << " Distance: " << rhitDistance << " recHits: " << layerHits << endl;
          }
          if (layerHits > detRecHits)
            detRecHits = layerHits;
        }
      } else {
        // Get the GEM Rechits from this layer
        GEMRecHitCollection::range dRecHits = theGEMRecHits->get(did);
 
        for (GEMRecHitCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
          double rhitDistance = (ir->localPosition() - (**imrh).localPosition()).mag();
          if (rhitDistance < coneSize)
            detRecHits++;
          LogTrace(theCategory) << ir->localPosition() << "  " << (**imrh).localPosition()
                                << " Distance: " << rhitDistance << " recHits: " << detRecHits << endl;
        }
      }
    }  //end of GEM if
    else if (id.subdetId() == MuonSubdetId::ME0) {
      ME0DetId did(id.rawId());
      chamberId = did.chamberId();
 
      if ((*imrh)->recHits().size() > 1) {
        std::vector<const TrackingRecHit*> hits2d = (*imrh)->recHits();
        for (std::vector<const TrackingRecHit*>::const_iterator hit2d = hits2d.begin(); hit2d != hits2d.end();
             hit2d++) {
          DetId id1 = (*hit2d)->geographicalId();
          ME0DetId lid(id1.rawId());
 
          // Get the ME0 Rechits from this layer
          ME0SegmentCollection::range dRecHits = theME0RecHits->get(lid);
          int layerHits = 0;
 
          for (ME0SegmentCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
            double rhitDistance = (ir->localPosition() - (**hit2d).localPosition()).mag();
            if (rhitDistance < coneSize)
              layerHits++;
            LogTrace(theCategory) << ir->localPosition() << "  " << (**hit2d).localPosition()
                                  << " Distance: " << rhitDistance << " recHits: " << layerHits << endl;
          }
          if (layerHits > detRecHits)
            detRecHits = layerHits;
        }
      } else {
        // Get the ME0 Rechits from this layer
        ME0SegmentCollection::range dRecHits = theME0RecHits->get(did);
 
        for (ME0SegmentCollection::const_iterator ir = dRecHits.first; ir != dRecHits.second; ir++) {
          double rhitDistance = (ir->localPosition() - (**imrh).localPosition()).mag();
          if (rhitDistance < coneSize)
            detRecHits++;
          LogTrace(theCategory) << ir->localPosition() << "  " << (**imrh).localPosition()
                                << " Distance: " << rhitDistance << " recHits: " << detRecHits << endl;
        }
      }
    }  //end of ME0 if
    else {
      if (id.subdetId() != MuonSubdetId::RPC)
        LogError(theCategory) << " Wrong Hit Type ";
      continue;
    }
 
    map<DetId, int>::iterator imap = hitMap.find(chamberId);
    if (imap != hitMap.end()) {
      if (detRecHits > imap->second)
        imap->second = detRecHits;
    } else
      hitMap[chamberId] = detRecHits;
 
  }  // end of loop over muon rechits
 
  for (map<DetId, int>::iterator imap = hitMap.begin(); imap != hitMap.end(); imap++)
    LogTrace(theCategory) << " Station " << imap->first.rawId() << ": " << imap->second << endl;
 
  LogTrace(theCategory) << "CheckMuonHits: " << all.size();
 
  // check order of muon measurements
  if ((all.size() > 1) && (all.front()->globalPosition().mag() > all.back()->globalPosition().mag())) {
    LogTrace(theCategory) << "reverse order: ";
    stable_sort(all.begin(), all.end(), RecHitLessByDet(alongMomentum));
  }
}

References python.cmstools::all(), alongMomentum, CSCDetId::chamberId(), ME0DetId::chamberId(), GEMDetId::chamberId(), HLT_2018_cff::coneSize, MuonSubdetId::CSC, MuonSubdetId::DT, MuonSubdetId::GEM, globals_cff::id1, LogTrace, mag(), MuonSubdetId::ME0, DetId::Muon, MuonSubdetId::RPC, theCategory, theCSCRecHits, theDTRecHits, theGEMRecHits, and theME0RecHits.

Referenced by refit().

checkRecHitsOrdering()

GlobalMuonRefitter::RefitDirection GlobalMuonRefitter::checkRecHitsOrdering ( const ConstRecHitContainer &  recHits ) const

protected

Definition at line 695 of file GlobalMuonRefitter.cc.

                                                                      {
  if (!recHits.empty()) {
    ConstRecHitContainer::const_iterator frontHit = recHits.begin();
    ConstRecHitContainer::const_iterator backHit = recHits.end() - 1;
    while (!(*frontHit)->isValid() && frontHit != backHit) {
      frontHit++;
    }
    while (!(*backHit)->isValid() && backHit != frontHit) {
      backHit--;
    }
 
    double rFirst = (*frontHit)->globalPosition().mag();
    double rLast = (*backHit)->globalPosition().mag();
 
    if (rFirst < rLast)
      return insideOut;
    else if (rFirst > rLast)
      return outsideIn;
    else {
      LogError(theCategory) << "Impossible determine the rechits order" << endl;
      return undetermined;
    }
  } else {
    LogError(theCategory) << "Impossible determine the rechits order" << endl;
    return undetermined;
  }
}

References insideOut, outsideIn, FastTrackerRecHitMaskProducer_cfi::recHits, theCategory, and undetermined.

Referenced by transform().

getDYTInfo()

const reco::DYTInfo* GlobalMuonRefitter::getDYTInfo ( )

inline

Definition at line 103 of file GlobalMuonRefitter.h.

{ return dytInfo; }

References dytInfo.

getFirstHits()

void GlobalMuonRefitter::getFirstHits ( const reco::Track &  muon, ConstRecHitContainer &  all, ConstRecHitContainer &  first  ) const

protected

get the RecHits in the tracker and the first muon chamber with hits

Definition at line 529 of file GlobalMuonRefitter.cc.

                                                                         {
  LogTrace(theCategory) << " GlobalMuonRefitter::getFirstHits\nall rechits length:" << all.size() << endl;
  first.clear();
 
  int station_to_keep = 999;
  vector<int> stations;
  for (ConstRecHitContainer::const_iterator ihit = all.begin(); ihit != all.end(); ++ihit) {
    int station = 0;
    bool use_it = true;
    DetId id = (*ihit)->geographicalId();
    unsigned raw_id = id.rawId();
    if (!(*ihit)->isValid())
      station = -1;
    else {
      if (id.det() == DetId::Muon) {
        switch (id.subdetId()) {
          case MuonSubdetId::DT:
            station = DTChamberId(raw_id).station();
            break;
          case MuonSubdetId::CSC:
            station = CSCDetId(raw_id).station();
            break;
          case MuonSubdetId::GEM:
            station = GEMDetId(raw_id).station();
            break;
          case MuonSubdetId::ME0:
            station = ME0DetId(raw_id).station();
            break;
          case MuonSubdetId::RPC:
            station = RPCDetId(raw_id).station();
            use_it = false;
            break;
        }
      }
    }
 
    if (use_it && station > 0 && station < station_to_keep)
      station_to_keep = station;
    stations.push_back(station);
    LogTrace(theCategory) << "rawId: " << raw_id << " station = " << station << " station_to_keep is now "
                          << station_to_keep;
  }
 
  if (station_to_keep <= 0 || station_to_keep > 4 || stations.size() != all.size())
    LogInfo(theCategory) << "failed to getFirstHits (all muon hits are outliers/bad ?)! station_to_keep = "
                         << station_to_keep << " stations.size " << stations.size() << " all.size " << all.size();
 
  for (unsigned i = 0; i < stations.size(); ++i)
    if (stations[i] >= 0 && stations[i] <= station_to_keep)
      first.push_back(all[i]);
 
  return;
}

References python.cmstools::all(), MuonSubdetId::CSC, MuonSubdetId::DT, dqmdumpme::first, MuonSubdetId::GEM, mps_fire::i, LogTrace, MuonSubdetId::ME0, DetId::Muon, DetId::rawId(), MuonSubdetId::RPC, DTChamberId::station(), ME0DetId::station(), relativeConstraints::station, RPCDetId::station(), CSCDetId::station(), GEMDetId::station(), and theCategory.

Referenced by refit().

getRidOfSelectStationHits()

GlobalMuonRefitter::ConstRecHitContainer GlobalMuonRefitter::getRidOfSelectStationHits	(	const ConstRecHitContainer &	hits,
		const TrackerTopology *	tTopo
	)		const

Definition at line 883 of file GlobalMuonRefitter.cc.

                                                                          {
  ConstRecHitContainer results;
  ConstRecHitContainer::const_iterator it = hits.begin();
  for (; it != hits.end(); it++) {
    DetId id = (*it)->geographicalId();
 
    //Check that this is a Muon hit that we're toying with -- else pass on this because the hacker is a moron / not careful
 
    if (id.det() == DetId::Tracker && theTrackerSkipSystem > 0) {
      int layer = -999;
      int disk = -999;
      int wheel = -999;
      if (id.subdetId() == theTrackerSkipSystem) {
        //                              continue;  //caveat that just removes the whole system from refitting
 
        if (theTrackerSkipSystem == PXB) {
          layer = tTopo->pxbLayer(id);
        }
        if (theTrackerSkipSystem == TIB) {
          layer = tTopo->tibLayer(id);
        }
 
        if (theTrackerSkipSystem == TOB) {
          layer = tTopo->tobLayer(id);
        }
        if (theTrackerSkipSystem == PXF) {
          disk = tTopo->pxfDisk(id);
        }
        if (theTrackerSkipSystem == TID) {
          wheel = tTopo->tidWheel(id);
        }
        if (theTrackerSkipSystem == TEC) {
          wheel = tTopo->tecWheel(id);
        }
        if (theTrackerSkipSection >= 0 && layer == theTrackerSkipSection)
          continue;
        if (theTrackerSkipSection >= 0 && disk == theTrackerSkipSection)
          continue;
        if (theTrackerSkipSection >= 0 && wheel == theTrackerSkipSection)
          continue;
      }
    }
 
    if (id.det() == DetId::Muon && theSkipStation) {
      int station = -999;
      //UNUSED:      int wheel = -999;
      if (id.subdetId() == MuonSubdetId::DT) {
        DTChamberId did(id.rawId());
        station = did.station();
        //UNUSED:   wheel = did.wheel();
      } else if (id.subdetId() == MuonSubdetId::CSC) {
        CSCDetId did(id.rawId());
        station = did.station();
      } else if (id.subdetId() == MuonSubdetId::GEM) {
        GEMDetId did(id.rawId());
        station = did.station();
      } else if (id.subdetId() == MuonSubdetId::ME0) {
        ME0DetId did(id.rawId());
        station = did.station();
      } else if (id.subdetId() == MuonSubdetId::RPC) {
        RPCDetId rpcid(id.rawId());
        station = rpcid.station();
      }
      if (station == theSkipStation)
        continue;
    }
    results.push_back(*it);
  }
  return results;
}

References MuonSubdetId::CSC, MuonSubdetId::DT, MuonSubdetId::GEM, hfClusterShapes_cfi::hits, MuonSubdetId::ME0, DetId::Muon, PXB, TrackerTopology::pxbLayer(), PXF, TrackerTopology::pxfDisk(), bookConverter::results, MuonSubdetId::RPC, DTChamberId::station(), ME0DetId::station(), relativeConstraints::station, RPCDetId::station(), CSCDetId::station(), GEMDetId::station(), TEC, TrackerTopology::tecWheel(), theSkipStation, theTrackerSkipSection, theTrackerSkipSystem, TIB, TrackerTopology::tibLayer(), TID, TrackerTopology::tidWheel(), TOB, TrackerTopology::tobLayer(), DetId::Tracker, and makeMuonMisalignmentScenario::wheel.

Referenced by refit().

printHits()

void GlobalMuonRefitter::printHits ( const ConstRecHitContainer &  hits ) const

protected

print all RecHits of a trajectory

Definition at line 674 of file GlobalMuonRefitter.cc.

                                                                         {
  LogTrace(theCategory) << "Used RecHits: " << hits.size();
  for (ConstRecHitContainer::const_iterator ir = hits.begin(); ir != hits.end(); ir++) {
    if (!(*ir)->isValid()) {
      LogTrace(theCategory) << "invalid RecHit";
      continue;
    }
 
    const GlobalPoint& pos = (*ir)->globalPosition();
 
    LogTrace(theCategory) << "r = " << sqrt(pos.x() * pos.x() + pos.y() * pos.y()) << "  z = " << pos.z()
                          << "  dimension = " << (*ir)->dimension()
                          << "  det = " << (*ir)->det()->geographicalId().det()
                          << "  subdet = " << (*ir)->det()->subDetector()
                          << "  raw id = " << (*ir)->det()->geographicalId().rawId();
  }
}

References hfClusterShapes_cfi::hits, LogTrace, mathSSE::sqrt(), and theCategory.

Referenced by ntuplePrintersDiff.TrackingParticlePrinter::__call__(), ntuplePrintersDiff.SeedPrinter::diff(), ntuplePrintersDiff.TrackPrinter::diff(), ntuplePrintersDiff.TrackingParticlePrinter::diff(), ntuplePrintersDiff.SeedPrinter::printSeed(), ntuplePrintersDiff.TrackPrinter::printTrack(), and transform().

refit() [1/2]

vector< Trajectory > GlobalMuonRefitter::refit	(	const reco::Track &	globalTrack,
		const int	theMuonHitsOption,
		const TrackerTopology *	tTopo
	)		const

build combined trajectory from sta Track and tracker RecHits

Definition at line 182 of file GlobalMuonRefitter.cc.

                                                                                 {
  LogTrace(theCategory) << " *** GlobalMuonRefitter *** option " << theMuonHitsOption << endl;
 
  ConstRecHitContainer allRecHitsTemp;  // all muon rechits temp
 
  reco::TransientTrack track(globalTrack, &*(theService->magneticField()), theService->trackingGeometry());
 
  auto tkbuilder = static_cast<TkTransientTrackingRecHitBuilder const*>(theTrackerRecHitBuilder.product());
 
  for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit)
    if ((*hit)->isValid()) {
      if ((*hit)->geographicalId().det() == DetId::Tracker)
        allRecHitsTemp.push_back((**hit).cloneForFit(*tkbuilder->geometry()->idToDet((**hit).geographicalId())));
      else if ((*hit)->geographicalId().det() == DetId::Muon) {
        if ((*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit) {
          LogTrace(theCategory) << "RPC Rec Hit discarged";
          continue;
        }
        allRecHitsTemp.push_back(theMuonRecHitBuilder->build(&**hit));
      }
    }
  vector<Trajectory> refitted = refit(globalTrack, track, allRecHitsTemp, theMuonHitsOption, tTopo);
  return refitted;
}

References TransientTrackingRecHitBuilder::build(), cutBasedMuonId_MuonPOG_V0_cff::globalTrack, LogTrace, DetId::Muon, edm::ESHandle< T >::product(), theCategory, theMuonHitsOption, theMuonRecHitBuilder, theRPCInTheFit, theService, theTrackerRecHitBuilder, HLT_2018_cff::track, and DetId::Tracker.

Referenced by GlobalTrajectoryBuilderBase::build(), and GlobalTrackQualityProducer::produce().

refit() [2/2]

vector< Trajectory > GlobalMuonRefitter::refit	(	const reco::Track &	globalTrack,
		const reco::TransientTrack	track,
		const TransientTrackingRecHit::ConstRecHitContainer &	allRecHitsTemp,
		const int	theMuonHitsOption,
		const TrackerTopology *	tTopo
	)		const

build combined trajectory from subset of sta Track and tracker RecHits

Definition at line 212 of file GlobalMuonRefitter.cc.

                                                                                 {
  // MuonHitsOption: 0 - tracker only
  //                 1 - include all muon hits
  //                 2 - include only first muon hit(s)
  //                 3 - include only selected muon hits
  //                 4 - redo pattern recognition with dynamic truncation
 
  vector<int> stationHits(4, 0);
  map<DetId, int> hitMap;
 
  ConstRecHitContainer allRecHits;       // all muon rechits
  ConstRecHitContainer fmsRecHits;       // only first muon rechits
  ConstRecHitContainer selectedRecHits;  // selected muon rechits
  ConstRecHitContainer DYTRecHits;       // rec hits from dynamic truncation algorithm
 
  LogTrace(theCategory) << " *** GlobalMuonRefitter *** option " << theMuonHitsOption << endl;
  LogTrace(theCategory) << " Track momentum before refit: " << globalTrack.pt() << endl;
  LogTrace(theCategory) << " Hits size before : " << allRecHitsTemp.size() << endl;
 
  allRecHits = getRidOfSelectStationHits(allRecHitsTemp, tTopo);
  //    printHits(allRecHits);
  LogTrace(theCategory) << " Hits size: " << allRecHits.size() << endl;
 
  vector<Trajectory> outputTraj;
 
  if ((theMuonHitsOption == 1) || (theMuonHitsOption == 3) || (theMuonHitsOption == 4)) {
    // refit the full track with all muon hits
    vector<Trajectory> globalTraj = transform(globalTrack, track, allRecHits);
 
    if (globalTraj.empty()) {
      LogTrace(theCategory) << "No trajectory from the TrackTransformer!" << endl;
      return vector<Trajectory>();
    }
 
    LogTrace(theCategory) << " Initial trajectory state: "
                          << globalTraj.front().lastMeasurement().updatedState().freeState()->parameters() << endl;
 
    if (theMuonHitsOption == 1)
      outputTraj.push_back(globalTraj.front());
 
    if (theMuonHitsOption == 3) {
      checkMuonHits(globalTrack, allRecHits, hitMap);
      selectedRecHits = selectMuonHits(globalTraj.front(), hitMap);
      LogTrace(theCategory) << " Selected hits size: " << selectedRecHits.size() << endl;
      outputTraj = transform(globalTrack, track, selectedRecHits);
    }
 
    if (theMuonHitsOption == 4) {
      //
      // DYT 2.0
      //
      DynamicTruncation dytRefit(*theEvent, *theService);
      dytRefit.setProd(all4DSegments, CSCSegments);
      dytRefit.setSelector(theDYTselector);
      dytRefit.setThr(theDYTthrs);
      dytRefit.setUpdateState(theDYTupdator);
      dytRefit.setUseAPE(theDYTuseAPE);
      if (theDYTParThrsMode) {
        dytRefit.setParThrsMode(theDYTParThrsMode);
        dytRefit.setThrsMap(theDYTthrsParameters);
        dytRefit.setRecoP(globalTrack.p());
        dytRefit.setRecoEta(globalTrack.eta());
      }
      DYTRecHits = dytRefit.filter(globalTraj.front());
      dytInfo->CopyFrom(dytRefit.getDYTInfo());
      if ((DYTRecHits.size() > 1) &&
          (DYTRecHits.front()->globalPosition().mag() > DYTRecHits.back()->globalPosition().mag()))
        stable_sort(DYTRecHits.begin(), DYTRecHits.end(), RecHitLessByDet(alongMomentum));
      outputTraj = transform(globalTrack, track, DYTRecHits);
    }
 
  } else if (theMuonHitsOption == 2) {
    getFirstHits(globalTrack, allRecHits, fmsRecHits);
    outputTraj = transform(globalTrack, track, fmsRecHits);
  }
 
  if (!outputTraj.empty()) {
    LogTrace(theCategory) << "Refitted pt: "
                          << outputTraj.front().firstMeasurement().updatedState().globalParameters().momentum().perp()
                          << endl;
    return outputTraj;
  } else {
    LogTrace(theCategory) << "No refitted Tracks... " << endl;
    return vector<Trajectory>();
  }
}

References all4DSegments, alongMomentum, checkMuonHits(), reco::DYTInfo::CopyFrom(), CSCSegments, dytInfo, DynamicTruncation::filter(), DynamicTruncation::getDYTInfo(), getFirstHits(), getRidOfSelectStationHits(), cutBasedMuonId_MuonPOG_V0_cff::globalTrack, LogTrace, selectMuonHits(), DynamicTruncation::setParThrsMode(), DynamicTruncation::setProd(), DynamicTruncation::setRecoEta(), DynamicTruncation::setRecoP(), DynamicTruncation::setSelector(), DynamicTruncation::setThr(), DynamicTruncation::setThrsMap(), DynamicTruncation::setUpdateState(), DynamicTruncation::setUseAPE(), theCategory, theDYTParThrsMode, theDYTselector, theDYTthrs, theDYTthrsParameters, theDYTupdator, theDYTuseAPE, theEvent, theMuonHitsOption, theService, HLT_2018_cff::track, and transform().

selectMuonHits()

GlobalMuonRefitter::ConstRecHitContainer GlobalMuonRefitter::selectMuonHits ( const Trajectory &  traj, const std::map< DetId, int > &  hitMap  ) const

protected

select muon hits compatible with trajectory; check hits in chambers with showers

Definition at line 589 of file GlobalMuonRefitter.cc.

                                                                                                                 {
  ConstRecHitContainer muonRecHits;
  const double globalChi2Cut = 200.0;
 
  vector<TrajectoryMeasurement> muonMeasurements = traj.measurements();
 
  // loop through all muon hits and skip hits with bad chi2 in chambers with high occupancy
  for (std::vector<TrajectoryMeasurement>::const_iterator im = muonMeasurements.begin(); im != muonMeasurements.end();
       im++) {
    if (!(*im).recHit()->isValid())
      continue;
    if ((*im).recHit()->det()->geographicalId().det() != DetId::Muon) {
      //      if ( ( chi2ndf < globalChi2Cut ) )
      muonRecHits.push_back((*im).recHit());
      continue;
    }
    const MuonTransientTrackingRecHit* immrh = dynamic_cast<const MuonTransientTrackingRecHit*>((*im).recHit().get());
 
    DetId id = immrh->geographicalId();
    DetId chamberId;
    int threshold = 0;
    double chi2Cut = 0.0;
 
    // get station of hit if it is in DT
    if ((*immrh).isDT()) {
      DTChamberId did(id.rawId());
      chamberId = did;
      threshold = theHitThreshold;
      chi2Cut = theDTChi2Cut;
    }
    // get station of hit if it is in CSC
    else if ((*immrh).isCSC()) {
      CSCDetId did(id.rawId());
      chamberId = did.chamberId();
      threshold = theHitThreshold;
      chi2Cut = theCSCChi2Cut;
    }
    // get station of hit if it is in GEM
    else if ((*immrh).isGEM()) {
      GEMDetId did(id.rawId());
      chamberId = did.chamberId();
      threshold = theHitThreshold;
      chi2Cut = theGEMChi2Cut;
    }
    // get station of hit if it is in ME0
    else if ((*immrh).isME0()) {
      ME0DetId did(id.rawId());
      chamberId = did.chamberId();
      threshold = theHitThreshold;
      chi2Cut = theME0Chi2Cut;
    }
    // get station of hit if it is in RPC
    else if ((*immrh).isRPC()) {
      RPCDetId rpcid(id.rawId());
      chamberId = rpcid;
      threshold = theHitThreshold;
      chi2Cut = theRPCChi2Cut;
    } else
      continue;
 
    double chi2ndf = (*im).estimate() / (*im).recHit()->dimension();
 
    bool keep = true;
    map<DetId, int>::const_iterator imap = hitMap.find(chamberId);
    if (imap != hitMap.end())
      if (imap->second > threshold)
        keep = false;
 
    if ((keep || (chi2ndf < chi2Cut)) && (chi2ndf < globalChi2Cut)) {
      muonRecHits.push_back((*im).recHit());
    } else {
      LogTrace(theCategory) << "Skip hit: " << id.rawId() << " chi2=" << chi2ndf << " ( threshold: " << chi2Cut
                            << ") Det: " << imap->second << endl;
    }
  }
 
  // check order of rechits
  reverse(muonRecHits.begin(), muonRecHits.end());
  return muonRecHits;
}

References CSCDetId::chamberId(), ME0DetId::chamberId(), GEMDetId::chamberId(), ZMuMuCategoriesSequences_cff::chi2Cut, TrackingRecHit::geographicalId(), keep, LogTrace, Trajectory::measurements(), DetId::Muon, groupFilesInBlocks::reverse, theCategory, theCSCChi2Cut, theDTChi2Cut, theGEMChi2Cut, theHitThreshold, theME0Chi2Cut, theRPCChi2Cut, and remoteMonitoring_LED_IterMethod_cfg::threshold.

Referenced by refit().

service()

const MuonServiceProxy* GlobalMuonRefitter::service ( ) const

inlineprotected

Definition at line 122 of file GlobalMuonRefitter.h.

{ return theService; }

References theService.

setEvent()

void GlobalMuonRefitter::setEvent ( const edm::Event &  event )

virtual

pass the Event to the algo at each event

Definition at line 153 of file GlobalMuonRefitter.cc.

                                                       {
  theEvent = &event;
  event.getByToken(theDTRecHitToken, theDTRecHits);
  event.getByToken(theCSCRecHitToken, theCSCRecHits);
  event.getByToken(theGEMRecHitToken, theGEMRecHits);
  event.getByToken(theME0RecHitToken, theME0RecHits);
  event.getByToken(CSCSegmentsToken, CSCSegments);
  event.getByToken(all4DSegmentsToken, all4DSegments);
}

References all4DSegments, all4DSegmentsToken, CSCSegments, CSCSegmentsToken, event(), theCSCRecHits, theCSCRecHitToken, theDTRecHits, theDTRecHitToken, theEvent, theGEMRecHits, theGEMRecHitToken, theME0RecHits, and theME0RecHitToken.

Referenced by GlobalTrackQualityProducer::produce(), and GlobalTrajectoryBuilderBase::setEvent().

setServices()

void GlobalMuonRefitter::setServices

(

const edm::EventSetup &

setup

)

set the services needed by the TrackTransformer

Definition at line 163 of file GlobalMuonRefitter.cc.

                                                            {
  edm::ESHandle<TrajectoryFitter> aFitter;
  theService->eventSetup().get<TrajectoryFitter::Record>().get(theFitterName, aFitter);
  theFitter = aFitter->clone();
 
  // Transient Rechit Builders
  unsigned long long newCacheId_TRH = setup.get<TransientRecHitRecord>().cacheIdentifier();
  if (newCacheId_TRH != theCacheId_TRH) {
    LogDebug(theCategory) << "TransientRecHitRecord changed!";
    setup.get<TransientRecHitRecord>().get(theTrackerRecHitBuilderName, theTrackerRecHitBuilder);
    setup.get<TransientRecHitRecord>().get(theMuonRecHitBuilderName, theMuonRecHitBuilder);
    hitCloner = static_cast<TkTransientTrackingRecHitBuilder const*>(theTrackerRecHitBuilder.product())->cloner();
  }
  theFitter->setHitCloner(&hitCloner);
}

References TrajectoryFitter::clone(), get, hitCloner, LogDebug, edm::ESHandle< T >::product(), singleTopDQM_cfi::setup, theCacheId_TRH, theCategory, theFitter, theFitterName, theMuonRecHitBuilder, theMuonRecHitBuilderName, theService, theTrackerRecHitBuilder, and theTrackerRecHitBuilderName.

Referenced by GlobalTrackQualityProducer::produce(), and GlobalTrajectoryBuilderBase::setEvent().

transform()

vector< Trajectory > GlobalMuonRefitter::transform	(	const reco::Track &	newTrack,
		const reco::TransientTrack	track,
		const TransientTrackingRecHit::ConstRecHitContainer &	recHitsForReFit
	)		const

refit the track with a new set of RecHits

Definition at line 727 of file GlobalMuonRefitter.cc.

                                                                               {
  TransientTrackingRecHit::ConstRecHitContainer recHitsForReFit = urecHitsForReFit;
  LogTrace(theCategory) << "GlobalMuonRefitter::transform: " << recHitsForReFit.size() << " hits:";
  printHits(recHitsForReFit);
 
  if (recHitsForReFit.size() < 2)
    return vector<Trajectory>();
 
  // Check the order of the rechits
  RefitDirection recHitsOrder = checkRecHitsOrdering(recHitsForReFit);
 
  LogTrace(theCategory) << "checkRecHitsOrdering() returned " << recHitsOrder << ", theRefitDirection is "
                        << theRefitDirection << " (insideOut == " << insideOut << ", outsideIn == " << outsideIn << ")";
 
  // Reverse the order in the case of inconsistency between the fit direction and the rechit order
  if (theRefitDirection != recHitsOrder)
    reverse(recHitsForReFit.begin(), recHitsForReFit.end());
 
  // Even though we checked the rechits' ordering above, we may have
  // already flipped them elsewhere (getFirstHits() is such a
  // culprit). Use the global positions of the states and the desired
  // refit direction to find the starting TSOS.
  TrajectoryStateOnSurface firstTSOS, lastTSOS;
  unsigned int innerId;  //UNUSED: outerId;
  bool order_swapped = track.outermostMeasurementState().globalPosition().mag() <
                       track.innermostMeasurementState().globalPosition().mag();
  bool inner_is_first;
  LogTrace(theCategory) << "order swapped? " << order_swapped;
 
  // Fill the starting state, depending on the ordering above.
  if ((theRefitDirection == insideOut && !order_swapped) || (theRefitDirection == outsideIn && order_swapped)) {
    innerId = newTrack.innerDetId();
    //UNUSED:    outerId   = newTrack.outerDetId();
    firstTSOS = track.innermostMeasurementState();
    lastTSOS = track.outermostMeasurementState();
    inner_is_first = true;
  } else {
    innerId = newTrack.outerDetId();
    //UNUSED:    outerId   = newTrack.innerDetId();
    firstTSOS = track.outermostMeasurementState();
    lastTSOS = track.innermostMeasurementState();
    inner_is_first = false;
  }
 
  LogTrace(theCategory) << "firstTSOS: inner_is_first? " << inner_is_first << " globalPosition is "
                        << firstTSOS.globalPosition() << " innerId is " << innerId;
 
  if (!firstTSOS.isValid()) {
    LogWarning(theCategory) << "Error wrong initial state!" << endl;
    return vector<Trajectory>();
  }
 
  firstTSOS.rescaleError(theRescaleErrorFactor);
 
  // This is the only way to get a TrajectorySeed with settable propagation direction
  PTrajectoryStateOnDet garbage1;
  edm::OwnVector<TrackingRecHit> garbage2;
 
  // These lines cause the code to ignore completely what was set
  // above, and force propDir for tracks from collisions!
  //  if(propDir == alongMomentum && theRefitDirection == outsideIn)  propDir=oppositeToMomentum;
  //  if(propDir == oppositeToMomentum && theRefitDirection == insideOut) propDir=alongMomentum;
 
  const TrajectoryStateOnSurface& tsosForDir = inner_is_first ? lastTSOS : firstTSOS;
  PropagationDirection propDir =
      (tsosForDir.globalPosition().basicVector().dot(tsosForDir.globalMomentum().basicVector()) > 0)
          ? alongMomentum
          : oppositeToMomentum;
  LogTrace(theCategory) << "propDir based on firstTSOS x dot p is " << propDir << " (alongMomentum == " << alongMomentum
                        << ", oppositeToMomentum == " << oppositeToMomentum << ")";
 
  // Additional propagation diretcion determination logic for cosmic muons
  if (theCosmicFlag) {
    PropagationDirection propDir_first =
        (firstTSOS.globalPosition().basicVector().dot(firstTSOS.globalMomentum().basicVector()) > 0)
            ? alongMomentum
            : oppositeToMomentum;
    PropagationDirection propDir_last =
        (lastTSOS.globalPosition().basicVector().dot(lastTSOS.globalMomentum().basicVector()) > 0) ? alongMomentum
                                                                                                   : oppositeToMomentum;
    LogTrace(theCategory) << "propDir_first " << propDir_first << ", propdir_last " << propDir_last << " : they "
                          << (propDir_first == propDir_last ? "agree" : "disagree");
 
    int y_count = 0;
    for (TransientTrackingRecHit::ConstRecHitContainer::const_iterator it = recHitsForReFit.begin();
         it != recHitsForReFit.end();
         ++it) {
      if ((*it)->globalPosition().y() > 0)
        ++y_count;
      else
        --y_count;
    }
 
    PropagationDirection propDir_ycount = alongMomentum;
    if (y_count > 0) {
      if (theRefitDirection == insideOut)
        propDir_ycount = oppositeToMomentum;
      else if (theRefitDirection == outsideIn)
        propDir_ycount = alongMomentum;
    } else {
      if (theRefitDirection == insideOut)
        propDir_ycount = alongMomentum;
      else if (theRefitDirection == outsideIn)
        propDir_ycount = oppositeToMomentum;
    }
 
    LogTrace(theCategory) << "y_count = " << y_count << "; based on geometrically-outermost TSOS, propDir is "
                          << propDir << ": " << (propDir == propDir_ycount ? "agrees" : "disagrees")
                          << " with ycount determination";
 
    if (propDir_first != propDir_last) {
      LogTrace(theCategory) << "since first/last disagreed, using y_count propDir";
      propDir = propDir_ycount;
    }
  }
 
  TrajectorySeed seed(garbage1, garbage2, propDir);
 
  if (recHitsForReFit.front()->geographicalId() != DetId(innerId)) {
    LogDebug(theCategory) << "Propagation occured" << endl;
    LogTrace(theCategory) << "propagating firstTSOS at " << firstTSOS.globalPosition()
                          << " to first rechit with surface pos "
                          << recHitsForReFit.front()->det()->surface().toGlobal(LocalPoint(0, 0, 0));
    firstTSOS =
        theService->propagator(thePropagatorName)->propagate(firstTSOS, recHitsForReFit.front()->det()->surface());
    if (!firstTSOS.isValid()) {
      LogDebug(theCategory) << "Propagation error!" << endl;
      return vector<Trajectory>();
    }
  }
 
  LogDebug(theCategory) << " GlobalMuonRefitter : theFitter " << propDir << endl;
  LogDebug(theCategory) << "                      First TSOS: " << firstTSOS.globalPosition()
                        << "  p=" << firstTSOS.globalMomentum() << " = " << firstTSOS.globalMomentum().mag() << endl;
 
  LogDebug(theCategory) << "                      Starting seed: "
                        << " nHits= " << seed.nHits() << " tsos: " << seed.startingState().parameters().position()
                        << "  p=" << seed.startingState().parameters().momentum() << endl;
 
  LogDebug(theCategory) << "                      RecHits: " << recHitsForReFit.size() << endl;
 
  vector<Trajectory> trajectories = theFitter->fit(seed, recHitsForReFit, firstTSOS);
 
  if (trajectories.empty()) {
    LogDebug(theCategory) << "No Track refitted!" << endl;
    return vector<Trajectory>();
  }
  return trajectories;
}

References alongMomentum, PV3DBase< T, PVType, FrameType >::basicVector(), checkRecHitsOrdering(), Basic3DVector< T >::dot(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), reco::Track::innerDetId(), insideOut, LogDebug, LogTrace, oppositeToMomentum, reco::Track::outerDetId(), outsideIn, printHits(), groupFilesInBlocks::reverse, SurveyInfoScenario_cff::seed, theCategory, theCosmicFlag, theFitter, thePropagatorName, theRefitDirection, theRescaleErrorFactor, theService, HLT_2018_cff::track, and FastTrackerRecHitMaskProducer_cfi::trajectories.

Referenced by refit().

Member Data Documentation

all4DSegments

edm::Handle<DTRecSegment4DCollection> GlobalMuonRefitter::all4DSegments

private

Definition at line 191 of file GlobalMuonRefitter.h.

Referenced by refit(), and setEvent().

all4DSegmentsToken

edm::EDGetTokenT<DTRecSegment4DCollection> GlobalMuonRefitter::all4DSegmentsToken

private

Definition at line 189 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setEvent().

CSCSegments

edm::Handle<CSCSegmentCollection> GlobalMuonRefitter::CSCSegments

private

Definition at line 190 of file GlobalMuonRefitter.h.

Referenced by refit(), and setEvent().

CSCSegmentsToken

edm::EDGetTokenT<CSCSegmentCollection> GlobalMuonRefitter::CSCSegmentsToken

private

Definition at line 188 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setEvent().

dytInfo

reco::DYTInfo* GlobalMuonRefitter::dytInfo

private

Definition at line 173 of file GlobalMuonRefitter.h.

Referenced by getDYTInfo(), GlobalMuonRefitter(), refit(), and ~GlobalMuonRefitter().

hitCloner

TkClonerImpl GlobalMuonRefitter::hitCloner

private

Definition at line 180 of file GlobalMuonRefitter.h.

Referenced by setServices().

theCacheId_TRH

unsigned long long GlobalMuonRefitter::theCacheId_TRH

private

Definition at line 157 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setServices().

theCategory

std::string GlobalMuonRefitter::theCategory

protected

Definition at line 125 of file GlobalMuonRefitter.h.

Referenced by checkMuonHits(), checkRecHitsOrdering(), getFirstHits(), GlobalMuonRefitter(), printHits(), refit(), selectMuonHits(), setServices(), and transform().

theCosmicFlag

bool GlobalMuonRefitter::theCosmicFlag

private

Definition at line 138 of file GlobalMuonRefitter.h.

Referenced by transform().

theCSCChi2Cut

float GlobalMuonRefitter::theCSCChi2Cut

private

Definition at line 134 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and selectMuonHits().

theCSCRecHitLabel

edm::InputTag GlobalMuonRefitter::theCSCRecHitLabel

private

Definition at line 141 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter().

theCSCRecHits

edm::Handle<CSCRecHit2DCollection> GlobalMuonRefitter::theCSCRecHits

private

Definition at line 145 of file GlobalMuonRefitter.h.

Referenced by checkMuonHits(), and setEvent().

theCSCRecHitToken

edm::EDGetTokenT<CSCRecHit2DCollection> GlobalMuonRefitter::theCSCRecHitToken

private

Definition at line 149 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setEvent().

theDTChi2Cut

float GlobalMuonRefitter::theDTChi2Cut

private

Definition at line 133 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and selectMuonHits().

theDTRecHitLabel

edm::InputTag GlobalMuonRefitter::theDTRecHitLabel

private

Definition at line 140 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter().

theDTRecHits

edm::Handle<DTRecHitCollection> GlobalMuonRefitter::theDTRecHits

private

Definition at line 144 of file GlobalMuonRefitter.h.

Referenced by checkMuonHits(), and setEvent().

theDTRecHitToken

edm::EDGetTokenT<DTRecHitCollection> GlobalMuonRefitter::theDTRecHitToken

private

Definition at line 148 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setEvent().

theDYTParThrsMode

bool GlobalMuonRefitter::theDYTParThrsMode

private

Definition at line 171 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and refit().

theDYTselector

int GlobalMuonRefitter::theDYTselector

private

Definition at line 168 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and refit().

theDYTthrs

std::vector<int> GlobalMuonRefitter::theDYTthrs

private

Definition at line 167 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and refit().

theDYTthrsParameters

edm::ParameterSet GlobalMuonRefitter::theDYTthrsParameters

private

Definition at line 172 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and refit().

theDYTupdator

bool GlobalMuonRefitter::theDYTupdator

private

Definition at line 169 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and refit().

theDYTuseAPE

bool GlobalMuonRefitter::theDYTuseAPE

private

Definition at line 170 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and refit().

theEvent

const edm::Event* GlobalMuonRefitter::theEvent

private

Definition at line 186 of file GlobalMuonRefitter.h.

Referenced by refit(), and setEvent().

theFitter

std::unique_ptr<TrajectoryFitter> GlobalMuonRefitter::theFitter

private

Definition at line 176 of file GlobalMuonRefitter.h.

Referenced by setServices(), and transform().

theFitterName

std::string GlobalMuonRefitter::theFitterName

private

Definition at line 175 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setServices().

theGEMChi2Cut

float GlobalMuonRefitter::theGEMChi2Cut

private

Definition at line 136 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and selectMuonHits().

theGEMRecHitLabel

edm::InputTag GlobalMuonRefitter::theGEMRecHitLabel

private

Definition at line 142 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter().

theGEMRecHits

edm::Handle<GEMRecHitCollection> GlobalMuonRefitter::theGEMRecHits

private

Definition at line 146 of file GlobalMuonRefitter.h.

Referenced by checkMuonHits(), and setEvent().

theGEMRecHitToken

edm::EDGetTokenT<GEMRecHitCollection> GlobalMuonRefitter::theGEMRecHitToken

private

Definition at line 150 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setEvent().

theHitThreshold

int GlobalMuonRefitter::theHitThreshold

private

Definition at line 132 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and selectMuonHits().

theME0Chi2Cut

float GlobalMuonRefitter::theME0Chi2Cut

private

Definition at line 137 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and selectMuonHits().

theME0RecHitLabel

edm::InputTag GlobalMuonRefitter::theME0RecHitLabel

private

Definition at line 143 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter().

theME0RecHits

edm::Handle<ME0SegmentCollection> GlobalMuonRefitter::theME0RecHits

private

Definition at line 147 of file GlobalMuonRefitter.h.

Referenced by checkMuonHits(), and setEvent().

theME0RecHitToken

edm::EDGetTokenT<ME0SegmentCollection> GlobalMuonRefitter::theME0RecHitToken

private

Definition at line 151 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setEvent().

theMuonHitsOption

int GlobalMuonRefitter::theMuonHitsOption

private

Definition at line 130 of file GlobalMuonRefitter.h.

Referenced by refit().

theMuonRecHitBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> GlobalMuonRefitter::theMuonRecHitBuilder

private

Definition at line 183 of file GlobalMuonRefitter.h.

Referenced by refit(), and setServices().

theMuonRecHitBuilderName

std::string GlobalMuonRefitter::theMuonRecHitBuilderName

private

Definition at line 182 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setServices().

theProbCut

float GlobalMuonRefitter::theProbCut

private

Definition at line 131 of file GlobalMuonRefitter.h.

thePropagatorName

std::string GlobalMuonRefitter::thePropagatorName

private

Definition at line 159 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and transform().

thePtCut

float GlobalMuonRefitter::thePtCut

protected

Definition at line 127 of file GlobalMuonRefitter.h.

theRefitDirection

RefitDirection GlobalMuonRefitter::theRefitDirection

private

Definition at line 165 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and transform().

theRescaleErrorFactor

double GlobalMuonRefitter::theRescaleErrorFactor

private

Definition at line 163 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and transform().

theRPCChi2Cut

float GlobalMuonRefitter::theRPCChi2Cut

private

Definition at line 135 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and selectMuonHits().

theRPCInTheFit

bool GlobalMuonRefitter::theRPCInTheFit

private

Definition at line 161 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and refit().

theService

const MuonServiceProxy* GlobalMuonRefitter::theService

private

Definition at line 185 of file GlobalMuonRefitter.h.

Referenced by refit(), service(), setServices(), and transform().

theSkipStation

int GlobalMuonRefitter::theSkipStation

private

Definition at line 153 of file GlobalMuonRefitter.h.

Referenced by getRidOfSelectStationHits(), and GlobalMuonRefitter().

theTkTrajsAvailableFlag

bool GlobalMuonRefitter::theTkTrajsAvailableFlag

protected

Definition at line 126 of file GlobalMuonRefitter.h.

theTrackerRecHitBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> GlobalMuonRefitter::theTrackerRecHitBuilder

private

Definition at line 179 of file GlobalMuonRefitter.h.

Referenced by refit(), and setServices().

theTrackerRecHitBuilderName

std::string GlobalMuonRefitter::theTrackerRecHitBuilderName

private

Definition at line 178 of file GlobalMuonRefitter.h.

Referenced by GlobalMuonRefitter(), and setServices().

theTrackerSkipSection

int GlobalMuonRefitter::theTrackerSkipSection

private

Definition at line 155 of file GlobalMuonRefitter.h.

Referenced by getRidOfSelectStationHits(), and GlobalMuonRefitter().

theTrackerSkipSystem

int GlobalMuonRefitter::theTrackerSkipSystem

private

Definition at line 154 of file GlobalMuonRefitter.h.

Referenced by getRidOfSelectStationHits(), and GlobalMuonRefitter().