#include <GlobalTrajectoryBuilderBase.h>

Inheritance diagram for GlobalTrajectoryBuilderBase:

Classes
struct	ComparatorInOut

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

typedef std::vector< Trajectory >	TC

typedef TC::const_iterator	TI

Public Types inherited from MuonTrajectoryBuilder
typedef MuonCandidate::CandidateContainer	CandidateContainer

typedef std::pair< const Trajectory *, reco::TrackRef >	TrackCand

typedef MuonCandidate::TrajectoryContainer	TrajectoryContainer

Public Member Functions
	GlobalTrajectoryBuilderBase (const edm::ParameterSet &, const MuonServiceProxy *, edm::ConsumesCollector &)
	constructor with Parameter Set and MuonServiceProxy More...

void	setEvent (const edm::Event &) override
	pass the Event to the algo at each event More...

MuonTrajectoryBuilder::TrajectoryContainer	trajectories (const TrajectorySeed &) override
	dummy implementation, unused in this class More...

	~GlobalTrajectoryBuilderBase () override
	destructor More...

Public Member Functions inherited from MuonTrajectoryBuilder
	MuonTrajectoryBuilder ()
	constructor More...

virtual CandidateContainer	trajectories (const TrackCand &)=0
	return a container reconstructed muons starting from a given track More...

virtual	~MuonTrajectoryBuilder ()
	destructor More...

Protected Types
enum	RefitDirection { inToOut, outToIn, undetermined }

Protected Member Functions
void	addTraj (TrackCand &)
	if TrackCand has only a TrackRef, attempt to add Trajectory* More...

MuonTrajectoryBuilder::CandidateContainer	build (const TrackCand &, MuonTrajectoryBuilder::CandidateContainer &) const
	build combined trajectory from sta Track and tracker RecHits More...

void	checkMuonHits (const reco::Track &, ConstRecHitContainer &, ConstRecHitContainer &, std::vector< int > &) const
	check muon RecHits, calculate chamber occupancy and select hits to be used in the final fit More...

RefitDirection	checkRecHitsOrdering (const ConstRecHitContainer &) const
	This does nothing now. More...

std::vector< TrackCand >	chooseRegionalTrackerTracks (const TrackCand &, const std::vector< TrackCand > &)
	choose tracker tracks within region of interest More...

const Trajectory *	chooseTrajectory (const std::vector< Trajectory * > &, int) const
	choose final trajectory More...

RectangularEtaPhiTrackingRegion	defineRegionOfInterest (const reco::TrackRef &) const
	define region of interest with tracker More...

void	fixTEC (ConstRecHitContainer &all, double scl_x, double scl_y) const
	rescale errors of outermost TEC RecHit More...

TransientTrackingRecHit::ConstRecHitContainer	getTransientRecHits (const reco::Track &) const
	get transient RecHits of a Track More...

virtual std::vector< TrackCand >	makeTkCandCollection (const TrackCand &)=0
	make a TrackCand collection using tracker Track, Trajectory information More...

void	printHits (const ConstRecHitContainer &) const
	print all RecHits of a trajectory More...

ConstRecHitContainer	selectMuonHits (const Trajectory &, const std::vector< int > &) const
	select muon hits compatible with trajectory; check hits in chambers with showers More...

ConstRecHitContainer	selectTrackerHits (const ConstRecHitContainer &) const
	select tracker hits; exclude some tracker hits in the global trajectory More...

const MuonServiceProxy *	service () const

GlobalMuonTrackMatcher *	trackMatcher () const

double	trackProbability (const Trajectory &) const
	calculate chi2 probability (-ln(P)) More...

Protected Attributes
std::string	theCategory

float	thePCut

float	thePtCut

Private Attributes
unsigned long long	theCacheId_TRH

const edm::Event *	theEvent

GlobalMuonRefitter *	theGlbRefitter

MuonDetLayerMeasurements *	theLayerMeasurements

int	theMuonHitsOption

edm::ESHandle< TransientTrackingRecHitBuilder >	theMuonRecHitBuilder

std::string	theMuonRecHitBuilderName

bool	theRefitFlag

MuonTrackingRegionBuilder *	theRegionBuilder

bool	theRPCInTheFit

const MuonServiceProxy *	theService

float	theTECxScale

float	theTECyScale

const TrackerTopology *	theTopo

std::string	theTrackerPropagatorName

edm::ESHandle< TransientTrackingRecHitBuilder >	theTrackerRecHitBuilder

std::string	theTrackerRecHitBuilderName

GlobalMuonTrackMatcher *	theTrackMatcher

TrackTransformer *	theTrackTransformer

Detailed Description

Base class for GlobalMuonTrajectoryBuilder and L3MuonTrajectoryBuilder Provide common tools and interface to reconstruct muons starting from a muon track reconstructed in the standalone muon system (with DT, CSC and RPC information). It tries to reconstruct the corresponding track in the tracker and performs matching between the reconstructed tracks in the muon system and the tracker.

Author: N. Neumeister Purdue University; C. Liu Purdue University; A. Everett Purdue University

Definition at line 42 of file GlobalTrajectoryBuilderBase.h.

Member Typedef Documentation

◆ ConstMuonRecHitContainer

typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer GlobalTrajectoryBuilderBase::ConstMuonRecHitContainer

Definition at line 51 of file GlobalTrajectoryBuilderBase.h.

◆ ConstMuonRecHitPointer

typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer GlobalTrajectoryBuilderBase::ConstMuonRecHitPointer

Definition at line 49 of file GlobalTrajectoryBuilderBase.h.

◆ ConstRecHitContainer

typedef TransientTrackingRecHit::ConstRecHitContainer GlobalTrajectoryBuilderBase::ConstRecHitContainer

Definition at line 45 of file GlobalTrajectoryBuilderBase.h.

◆ ConstRecHitPointer

typedef TransientTrackingRecHit::ConstRecHitPointer GlobalTrajectoryBuilderBase::ConstRecHitPointer

Definition at line 47 of file GlobalTrajectoryBuilderBase.h.

◆ MuonRecHitContainer

typedef MuonTransientTrackingRecHit::MuonRecHitContainer GlobalTrajectoryBuilderBase::MuonRecHitContainer

Definition at line 50 of file GlobalTrajectoryBuilderBase.h.

◆ MuonRecHitPointer

typedef MuonTransientTrackingRecHit::MuonRecHitPointer GlobalTrajectoryBuilderBase::MuonRecHitPointer

Definition at line 48 of file GlobalTrajectoryBuilderBase.h.

◆ RecHitContainer

typedef TransientTrackingRecHit::RecHitContainer GlobalTrajectoryBuilderBase::RecHitContainer

Definition at line 44 of file GlobalTrajectoryBuilderBase.h.

◆ RecHitPointer

typedef TransientTrackingRecHit::RecHitPointer GlobalTrajectoryBuilderBase::RecHitPointer

Definition at line 46 of file GlobalTrajectoryBuilderBase.h.

◆ TC

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

Definition at line 52 of file GlobalTrajectoryBuilderBase.h.

◆ TI

typedef TC::const_iterator GlobalTrajectoryBuilderBase::TI

Definition at line 53 of file GlobalTrajectoryBuilderBase.h.

Member Enumeration Documentation

◆ RefitDirection

enum GlobalTrajectoryBuilderBase::RefitDirection

protected

Enumerator
inToOut
outToIn
undetermined

Definition at line 70 of file GlobalTrajectoryBuilderBase.h.

70 { inToOut, outToIn, undetermined };

Constructor & Destructor Documentation

◆ GlobalTrajectoryBuilderBase()

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

constructor with Parameter Set and MuonServiceProxy

Class: GlobalTrajectoryBuilderBase

Base class for GlobalMuonTrajectoryBuilder and L3MuonTrajectoryBuilder Provide common tools and interface to reconstruct muons starting from a muon track reconstructed in the standalone muon system (with DT, CSC and RPC information). It tries to reconstruct the corresponding track in the tracker and performs matching between the reconstructed tracks in the muon system and the tracker.

Author: N. Neumeister Purdue University; C. Liu Purdue University; A. Everett Purdue University

Definition at line 76 of file GlobalTrajectoryBuilderBase.cc.

     : theLayerMeasurements(nullptr),
       theTrackTransformer(nullptr),
       theRegionBuilder(nullptr),
       theService(service),
       theGlbRefitter(nullptr) {
   theCategory =
       par.getUntrackedParameter<std::string>("Category", "Muon|RecoMuon|GlobalMuon|GlobalTrajectoryBuilderBase");
  
   edm::ParameterSet trackMatcherPSet = par.getParameter<edm::ParameterSet>("GlobalMuonTrackMatcher");
   theTrackMatcher = new GlobalMuonTrackMatcher(trackMatcherPSet, theService);
  
   theTrackerPropagatorName = par.getParameter<std::string>("TrackerPropagator");
  
   edm::ParameterSet trackTransformerPSet = par.getParameter<edm::ParameterSet>("TrackTransformer");
   theTrackTransformer = new TrackTransformer(trackTransformerPSet);
  
   edm::ParameterSet regionBuilderPSet = par.getParameter<edm::ParameterSet>("MuonTrackingRegionBuilder");
  
   theRegionBuilder = new MuonTrackingRegionBuilder(regionBuilderPSet, iC);
  
   // TrackRefitter parameters
   edm::ParameterSet refitterParameters = par.getParameter<edm::ParameterSet>("GlbRefitterParameters");
   theGlbRefitter = new GlobalMuonRefitter(refitterParameters, theService, iC);
  
   theMuonHitsOption = refitterParameters.getParameter<int>("MuonHitsOption");
   theRefitFlag = refitterParameters.getParameter<bool>("RefitFlag");
  
   theTrackerRecHitBuilderName = par.getParameter<std::string>("TrackerRecHitBuilder");
   theMuonRecHitBuilderName = par.getParameter<std::string>("MuonRecHitBuilder");
  
   theRPCInTheFit = par.getParameter<bool>("RefitRPCHits");
  
   theTECxScale = par.getParameter<double>("ScaleTECxFactor");
   theTECyScale = par.getParameter<double>("ScaleTECyFactor");
   thePtCut = par.getParameter<double>("PtCut");
   thePCut = par.getParameter<double>("PCut");
  
   theCacheId_TRH = 0;
 }

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), GlobalMuonRefitter_cff::GlobalMuonRefitter, HLT_2018_cff::GlobalMuonTrackMatcher, HLT_2018_cff::MuonTrackingRegionBuilder, AlCaHLTBitMon_QueryRunRegistry::string, theCacheId_TRH, theCategory, theGlbRefitter, theMuonHitsOption, theMuonRecHitBuilderName, thePCut, thePtCut, theRefitFlag, theRegionBuilder, theRPCInTheFit, theService, theTECxScale, theTECyScale, theTrackerPropagatorName, theTrackerRecHitBuilderName, theTrackMatcher, theTrackTransformer, and HLT_2018_cff::TrackTransformer.

◆ ~GlobalTrajectoryBuilderBase()

GlobalTrajectoryBuilderBase::~GlobalTrajectoryBuilderBase ( )

override

destructor

Definition at line 122 of file GlobalTrajectoryBuilderBase.cc.

                                                           {
   if (theTrackMatcher)
     delete theTrackMatcher;
   if (theRegionBuilder)
     delete theRegionBuilder;
   if (theTrackTransformer)
     delete theTrackTransformer;
   if (theGlbRefitter)
     delete theGlbRefitter;
 }

References theGlbRefitter, theRegionBuilder, theTrackMatcher, and theTrackTransformer.

Member Function Documentation

◆ addTraj()

void GlobalTrajectoryBuilderBase::addTraj ( TrackCand & )

inlineprotected

if TrackCand has only a TrackRef, attempt to add Trajectory*

Definition at line 106 of file GlobalTrajectoryBuilderBase.h.

106 {}

◆ build()

MuonCandidate::CandidateContainer GlobalTrajectoryBuilderBase::build	(	const TrackCand &	staCand,
		MuonTrajectoryBuilder::CandidateContainer &	tkTrajs
	)		const

protected

build combined trajectory from sta Track and tracker RecHits

Definition at line 162 of file GlobalTrajectoryBuilderBase.cc.

                                                                                                                      {
   LogTrace(theCategory) << " Begin Build" << std::endl;
  
   // tracker trajectory should be built and refit before this point
   if (tkTrajs.empty())
     return CandidateContainer();
  
   // add muon hits and refit/smooth trajectories
   CandidateContainer refittedResult;
   ConstRecHitContainer muonRecHits = getTransientRecHits(*(staCand.second));
  
   // check order of muon measurements
   if ((muonRecHits.size() > 1) &&
       (muonRecHits.front()->globalPosition().mag() > muonRecHits.back()->globalPosition().mag())) {
     LogTrace(theCategory) << "   reverse order: ";
   }
  
   for (auto&& it : tkTrajs) {
     // cut on tracks with low momenta
     LogTrace(theCategory) << "   Track p and pT " << it->trackerTrack()->p() << " " << it->trackerTrack()->pt();
     if (it->trackerTrack()->p() < thePCut || it->trackerTrack()->pt() < thePtCut)
       continue;
  
     // If true we will run theGlbRefitter->refit from all hits
     if (theRefitFlag) {
       ConstRecHitContainer trackerRecHits;
       if (it->trackerTrack().isNonnull()) {
         trackerRecHits = getTransientRecHits(*it->trackerTrack());
       } else {
         LogDebug(theCategory) << "     NEED HITS FROM TRAJ";
       }
  
       // ToDo: Do we need the following ?:
       // check for single TEC RecHits in trajectories in the overalp region
       if (std::abs(it->trackerTrack()->eta()) > 0.95 && std::abs(it->trackerTrack()->eta()) < 1.15 &&
           it->trackerTrack()->pt() < 60) {
         if (theTECxScale < 0 || theTECyScale < 0)
           trackerRecHits = selectTrackerHits(trackerRecHits);
         else
           fixTEC(trackerRecHits, theTECxScale, theTECyScale);
       }
  
       RefitDirection recHitDir = checkRecHitsOrdering(trackerRecHits);
       if (recHitDir == outToIn)
         reverse(trackerRecHits.begin(), trackerRecHits.end());
  
       reco::TransientTrack tTT(it->trackerTrack(), &*theService->magneticField(), theService->trackingGeometry());
       TrajectoryStateOnSurface innerTsos = tTT.innermostMeasurementState();
  
       edm::RefToBase<TrajectorySeed> tmpSeed;
       if (it->trackerTrack()->seedRef().isAvailable())
         tmpSeed = it->trackerTrack()->seedRef();
  
       if (!innerTsos.isValid()) {
         LogTrace(theCategory) << "     inner Trajectory State is invalid. ";
         continue;
       }
  
       innerTsos.rescaleError(100.);
  
       TC refitted0, refitted1;
       std::unique_ptr<Trajectory> tkTrajectory;
  
       // tracker only track
       if (!(it->trackerTrajectory() && it->trackerTrajectory()->isValid())) {
         refitted0 = theTrackTransformer->transform(it->trackerTrack());
         if (!refitted0.empty())
           tkTrajectory = std::make_unique<Trajectory>(*(refitted0.begin()));
         else
           edm::LogWarning(theCategory) << "     Failed to load tracker track trajectory";
       } else
         tkTrajectory = it->releaseTrackerTrajectory();
       if (tkTrajectory)
         tkTrajectory->setSeedRef(tmpSeed);
  
       // full track with all muon hits using theGlbRefitter
       ConstRecHitContainer allRecHits = trackerRecHits;
       allRecHits.insert(allRecHits.end(), muonRecHits.begin(), muonRecHits.end());
       refitted1 = theGlbRefitter->refit(*it->trackerTrack(), tTT, allRecHits, theMuonHitsOption, theTopo);
       LogTrace(theCategory) << "     This track-sta refitted to " << refitted1.size() << " trajectories";
  
       std::unique_ptr<Trajectory> glbTrajectory1;
       if (!refitted1.empty())
         glbTrajectory1 = std::make_unique<Trajectory>(*(refitted1.begin()));
       else
         LogDebug(theCategory) << "     Failed to load global track trajectory 1";
       if (glbTrajectory1)
         glbTrajectory1->setSeedRef(tmpSeed);
  
       if (glbTrajectory1 && tkTrajectory) {
         refittedResult.emplace_back(std::make_unique<MuonCandidate>(
             std::move(glbTrajectory1), it->muonTrack(), it->trackerTrack(), std::move(tkTrajectory)));
       }
     } else {
       edm::RefToBase<TrajectorySeed> tmpSeed;
       if (it->trackerTrack()->seedRef().isAvailable())
         tmpSeed = it->trackerTrack()->seedRef();
  
       TC refitted0;
       std::unique_ptr<Trajectory> tkTrajectory;
       if (!(it->trackerTrajectory() && it->trackerTrajectory()->isValid())) {
         refitted0 = theTrackTransformer->transform(it->trackerTrack());
         if (!refitted0.empty()) {
           tkTrajectory = std::make_unique<Trajectory>(*(refitted0.begin()));
         } else
           edm::LogWarning(theCategory) << "     Failed to load tracker track trajectory";
       } else
         tkTrajectory = it->releaseTrackerTrajectory();
       std::unique_ptr<Trajectory> cpy;
       if (tkTrajectory) {
         tkTrajectory->setSeedRef(tmpSeed);
         cpy = std::make_unique<Trajectory>(*tkTrajectory);
       }
       // Creating MuonCandidate using only the tracker trajectory:
       refittedResult.emplace_back(std::make_unique<MuonCandidate>(
           std::move(tkTrajectory), it->muonTrack(), it->trackerTrack(), std::move(cpy)));
     }
   }
  
   // choose the best global fit for this Standalone Muon based on the track probability
   CandidateContainer selectedResult;
   std::unique_ptr<MuonCandidate> tmpCand;
   double minProb = std::numeric_limits<double>::max();
  
   for (auto&& cand : refittedResult) {
     double prob = trackProbability(*cand->trajectory());
     LogTrace(theCategory) << "   refitted-track-sta with pT " << cand->trackerTrack()->pt() << " has probability "
                           << prob;
  
     if (prob < minProb or not tmpCand) {
       minProb = prob;
       tmpCand = std::move(cand);
     }
   }
  
   if (tmpCand)
     selectedResult.push_back(std::move(tmpCand));
  
   refittedResult.clear();
  
   return selectedResult;
 }

References funct::abs(), checkRecHitsOrdering(), fixTEC(), getTransientRecHits(), LogDebug, LogTrace, MuonServiceProxy::magneticField(), SiStripPI::max, conversionTrackMerger_cfi::minProb, eostools::move(), or, outToIn, TtFullHadEvtBuilder_cfi::prob, GlobalMuonRefitter::refit(), groupFilesInBlocks::reverse, selectTrackerHits(), theCategory, theGlbRefitter, theMuonHitsOption, thePCut, thePtCut, theRefitFlag, theService, theTECxScale, theTECyScale, theTopo, theTrackTransformer, MuonServiceProxy::trackingGeometry(), trackProbability(), and TrackTransformer::transform().

Referenced by GlobalMuonTrajectoryBuilder::trajectories(), and L3MuonTrajectoryBuilder::trajectories().

◆ checkMuonHits()

void GlobalTrajectoryBuilderBase::checkMuonHits	(	const reco::Track &	,
		ConstRecHitContainer &	,
		ConstRecHitContainer &	,
		std::vector< int > &
	)		const

protected

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

◆ checkRecHitsOrdering()

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

protected

This does nothing now.

check order of RechIts on a trajectory

Definition at line 393 of file GlobalTrajectoryBuilderBase.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 inToOut;
     else if (rFirst > rLast)
       return outToIn;
     else {
       edm::LogError(theCategory) << "Impossible to determine the rechits order" << std::endl;
       return undetermined;
     }
   } else {
     edm::LogError(theCategory) << "Impossible to determine the rechits order" << std::endl;
     return undetermined;
   }
 }

References inToOut, outToIn, FastTrackerRecHitMaskProducer_cfi::recHits, theCategory, and undetermined.

Referenced by build().

◆ chooseRegionalTrackerTracks()

std::vector< GlobalTrajectoryBuilderBase::TrackCand > GlobalTrajectoryBuilderBase::chooseRegionalTrackerTracks	(	const TrackCand &	staCand,
		const std::vector< TrackCand > &	tkTs
	)

protected

choose tracker tracks within region of interest

Definition at line 309 of file GlobalTrajectoryBuilderBase.cc.

                                                                 {
   // define eta-phi region
   RectangularEtaPhiTrackingRegion regionOfInterest = defineRegionOfInterest(staCand.second);
  
   // get region's etaRange and phiMargin
   //UNUSED:  PixelRecoRange<float> etaRange = regionOfInterest.etaRange();
   //UNUSED:  TkTrackingRegionsMargin<float> phiMargin = regionOfInterest.phiMargin();
  
   std::vector<TrackCand> result;
  
   double deltaR_max = 1.0;
  
   for (auto&& is : tkTs) {
     double deltaR_tmp = deltaR(static_cast<double>(regionOfInterest.direction().eta()),
                                static_cast<double>(regionOfInterest.direction().phi()),
                                is.second->eta(),
                                is.second->phi());
  
     // for each trackCand in region, add trajectory and add to result
     //if ( inEtaRange && inPhiRange ) {
     if (deltaR_tmp < deltaR_max) {
       TrackCand tmpCand = TrackCand(is);
       result.push_back(tmpCand);
     }
   }
  
   return result;
 }

References defineRegionOfInterest(), PbPb_ZMuSkimMuonDPG_cff::deltaR, TrackingRegion::direction(), PV3DBase< T, PVType, FrameType >::eta(), PV3DBase< T, PVType, FrameType >::phi(), and mps_fire::result.

Referenced by GlobalMuonTrajectoryBuilder::makeTkCandCollection().

◆ chooseTrajectory()

const Trajectory* GlobalTrajectoryBuilderBase::chooseTrajectory	(	const std::vector< Trajectory * > &	,
		int
	)		const

protected

choose final trajectory

◆ defineRegionOfInterest()

RectangularEtaPhiTrackingRegion GlobalTrajectoryBuilderBase::defineRegionOfInterest ( const reco::TrackRef & staTrack ) const

protected

define region of interest with tracker

Definition at line 342 of file GlobalTrajectoryBuilderBase.cc.

                                         {
   std::unique_ptr<RectangularEtaPhiTrackingRegion> region1 = theRegionBuilder->region(staTrack);
  
   TkTrackingRegionsMargin<float> etaMargin(std::abs(region1->etaRange().min() - region1->etaRange().mean()),
                                            std::abs(region1->etaRange().max() - region1->etaRange().mean()));
  
   RectangularEtaPhiTrackingRegion region2(region1->direction(),
                                           region1->origin(),
                                           region1->ptMin(),
                                           region1->originRBound(),
                                           region1->originZBound(),
                                           etaMargin,
                                           region1->phiMargin());
  
   return region2;
 }

References funct::abs(), MuonTrackingRegionBuilder::region(), and theRegionBuilder.

Referenced by chooseRegionalTrackerTracks().

◆ fixTEC()

void GlobalTrajectoryBuilderBase::fixTEC	(	ConstRecHitContainer &	all,
		double	scl_x,
		double	scl_y
	)		const

protected

rescale errors of outermost TEC RecHit

freeze this hit, make sure it will not be recomputed during fitting / the implemetantion below works with cloning / to get a RecHitPointer to SiStripRecHit2D, the only method that works is / RecHitPointer MuonTransientTrackingRecHit::build(const GeomDet*,const TrackingRecHit*)

Definition at line 449 of file GlobalTrajectoryBuilderBase.cc.

                                                                                                     {
   int nTEC(0);
   ConstRecHitContainer::iterator lone_tec;
  
   for (ConstRecHitContainer::iterator i = all.begin(); i != all.end(); i++) {
     if (!(*i)->isValid())
       continue;
  
     if ((*i)->det()->geographicalId().det() == DetId::Tracker &&
         (*i)->det()->geographicalId().subdetId() == StripSubdetector::TEC) {
       lone_tec = i;
       nTEC++;
  
       if ((i + 1) != all.end() && (*(i + 1))->isValid() &&
           (*(i + 1))->det()->geographicalId().det() == DetId::Tracker &&
           (*(i + 1))->det()->geographicalId().subdetId() == StripSubdetector::TEC) {
         nTEC++;
         break;
       }
     }
  
     if (nTEC > 1)
       break;
   }
  
   int hitDet = (*lone_tec)->hit()->geographicalId().det();
   int hitSubDet = (*lone_tec)->hit()->geographicalId().subdetId();
   if (nTEC == 1 && (*lone_tec)->hit()->isValid() && hitDet == DetId::Tracker && hitSubDet == StripSubdetector::TEC) {
     // rescale the TEC rechit error matrix in its rotated frame
     const SiStripRecHit2D* strip = dynamic_cast<const SiStripRecHit2D*>((*lone_tec)->hit());
     if (strip && strip->det()) {
       LocalPoint pos = strip->localPosition();
       if ((*lone_tec)->detUnit()) {
         const StripTopology* topology = dynamic_cast<const StripTopology*>(&(*lone_tec)->detUnit()->topology());
         if (topology) {
           // rescale the local error along/perp the strip by a factor
           float angle = topology->stripAngle(topology->strip((*lone_tec)->hit()->localPosition()));
           LocalError error = strip->localPositionError();
           LocalError rotError = error.rotate(angle);
           LocalError scaledError(rotError.xx() * scl_x * scl_x, 0, rotError.yy() * scl_y * scl_y);
           error = scaledError.rotate(-angle);
           SiStripRecHit2D* st = new SiStripRecHit2D(pos, error, *strip->det(), strip->cluster());
           *lone_tec = MuonTransientTrackingRecHit::build((*lone_tec)->det(), st);
         }
       }
     }
   }
 }

References python.cmstools::all(), angle(), MuonTransientTrackingRecHit::build(), relativeConstraints::error, mps_fire::i, digitizers_cfi::strip, StripSubdetector::TEC, ecaldqm::topology(), DetId::Tracker, LocalError::xx(), and LocalError::yy().

Referenced by build().

◆ getTransientRecHits()

TransientTrackingRecHit::ConstRecHitContainer GlobalTrajectoryBuilderBase::getTransientRecHits ( const reco::Track & track ) const

protected

get transient RecHits of a Track

Definition at line 506 of file GlobalTrajectoryBuilderBase.cc.

                                   {
   TransientTrackingRecHit::ConstRecHitContainer result;
  
   TrajectoryStateOnSurface currTsos = trajectoryStateTransform::innerStateOnSurface(
       track, *theService->trackingGeometry(), &*theService->magneticField());
  
   auto tkbuilder = static_cast<TkTransientTrackingRecHitBuilder const*>(theTrackerRecHitBuilder.product());
   auto hitCloner = tkbuilder->cloner();
   for (trackingRecHit_iterator hit = track.recHitsBegin(); hit != track.recHitsEnd(); ++hit) {
     if ((*hit)->isValid()) {
       DetId recoid = (*hit)->geographicalId();
       if (recoid.det() == DetId::Tracker) {
         if (!(*hit)->hasPositionAndError()) {
           TrajectoryStateOnSurface predTsos =
               theService->propagator(theTrackerPropagatorName)
                   ->propagate(currTsos, theService->trackingGeometry()->idToDet(recoid)->surface());
  
           if (!predTsos.isValid()) {
             edm::LogError("MissingTransientHit")
                 << "Could not get a tsos on the hit surface. We will miss a tracking hit.";
             continue;
           }
           currTsos = predTsos;
           auto h = (**hit).cloneForFit(*tkbuilder->geometry()->idToDet((**hit).geographicalId()));
           result.emplace_back(hitCloner.makeShared(h, predTsos));
         } else {
           result.push_back((*hit)->cloneSH());
         }
       } else if (recoid.det() == DetId::Muon) {
         if ((*hit)->geographicalId().subdetId() == 3 && !theRPCInTheFit) {
           LogDebug(theCategory) << "RPC Rec Hit discarded";
           continue;
         }
         result.push_back(theMuonRecHitBuilder->build(&**hit));
       }
     }
   }
  
   return result;
 }

References TransientTrackingRecHitBuilder::build(), DetId::det(), GlobalTrackingGeometry::idToDet(), trajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), LogDebug, MuonServiceProxy::magneticField(), DetId::Muon, edm::ESHandle< T >::product(), Propagator::propagate(), MuonServiceProxy::propagator(), mps_fire::result, GeomDet::surface(), theCategory, theMuonRecHitBuilder, theRPCInTheFit, theService, theTrackerPropagatorName, theTrackerRecHitBuilder, HLT_2018_cff::track, DetId::Tracker, and MuonServiceProxy::trackingGeometry().

Referenced by build().

◆ makeTkCandCollection()

virtual std::vector<TrackCand> GlobalTrajectoryBuilderBase::makeTkCandCollection ( const TrackCand & )

protectedpure virtual

make a TrackCand collection using tracker Track, Trajectory information

Implemented in L3MuonTrajectoryBuilder, and GlobalMuonTrajectoryBuilder.

◆ printHits()

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

protected

print all RecHits of a trajectory

Definition at line 374 of file GlobalTrajectoryBuilderBase.cc.

                                                                                   {
   LogTrace(theCategory) << "Used RecHits: " << hits.size();
   for (auto&& ir : hits) {
     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() << "  " << ir->det()->geographicalId().det() << "  "
                           << ir->det()->subDetector();
   }
 }

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(), and ntuplePrintersDiff.TrackPrinter::printTrack().

◆ selectMuonHits()

ConstRecHitContainer GlobalTrajectoryBuilderBase::selectMuonHits	(	const Trajectory &	,
		const std::vector< int > &
	)		const

protected

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

◆ selectTrackerHits()

GlobalTrajectoryBuilderBase::ConstRecHitContainer GlobalTrajectoryBuilderBase::selectTrackerHits ( const ConstRecHitContainer & all ) const

protected

select tracker hits; exclude some tracker hits in the global trajectory

Definition at line 425 of file GlobalTrajectoryBuilderBase.cc.

                                            {
   int nTEC(0);
  
   ConstRecHitContainer hits;
   for (auto&& i : all) {
     if (!i->isValid())
       continue;
     if (i->det()->geographicalId().det() == DetId::Tracker &&
         i->det()->geographicalId().subdetId() == StripSubdetector::TEC) {
       nTEC++;
     } else {
       hits.push_back(i);
     }
     if (nTEC > 1)
       return all;
   }
  
   return hits;
 }

References python.cmstools::all(), hfClusterShapes_cfi::hits, mps_fire::i, StripSubdetector::TEC, and DetId::Tracker.

Referenced by build().

◆ service()

const MuonServiceProxy* GlobalTrajectoryBuilderBase::service ( ) const

inlineprotected

Definition at line 118 of file GlobalTrajectoryBuilderBase.h.

118 { return theService; }

References theService.

◆ setEvent()

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

overridevirtual

pass the Event to the algo at each event

Implements MuonTrajectoryBuilder.

Reimplemented in L3MuonTrajectoryBuilder.

Definition at line 136 of file GlobalTrajectoryBuilderBase.cc.

                                                                 {
   theEvent = &event;
  
   theTrackTransformer->setServices(theService->eventSetup());
   theRegionBuilder->setEvent(event);
  
   theGlbRefitter->setEvent(event);
   theGlbRefitter->setServices(theService->eventSetup());
  
   unsigned long long newCacheId_TRH = theService->eventSetup().get<TransientRecHitRecord>().cacheIdentifier();
   if (newCacheId_TRH != theCacheId_TRH) {
     LogDebug(theCategory) << "TransientRecHitRecord changed!";
     theCacheId_TRH = newCacheId_TRH;
     theService->eventSetup().get<TransientRecHitRecord>().get(theTrackerRecHitBuilderName, theTrackerRecHitBuilder);
     theService->eventSetup().get<TransientRecHitRecord>().get(theMuonRecHitBuilderName, theMuonRecHitBuilder);
   }
  
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHand;
   theService->eventSetup().get<TrackerTopologyRcd>().get(tTopoHand);
   theTopo = tTopoHand.product();
 }

References event(), MuonServiceProxy::eventSetup(), edm::EventSetup::get(), get, LogDebug, edm::ESHandle< T >::product(), MuonTrackingRegionBuilder::setEvent(), GlobalMuonRefitter::setEvent(), GlobalMuonRefitter::setServices(), TrackTransformer::setServices(), theCacheId_TRH, theCategory, theEvent, theGlbRefitter, theMuonRecHitBuilder, theMuonRecHitBuilderName, theRegionBuilder, theService, theTopo, theTrackerRecHitBuilder, theTrackerRecHitBuilderName, and theTrackTransformer.

Referenced by GlobalMuonTrajectoryBuilder::setEvent(), and L3MuonTrajectoryBuilder::setEvent().

◆ trackMatcher()

GlobalMuonTrackMatcher* GlobalTrajectoryBuilderBase::trackMatcher ( ) const

inlineprotected

Definition at line 115 of file GlobalTrajectoryBuilderBase.h.

115 { return theTrackMatcher; }

References theTrackMatcher.

Referenced by GlobalMuonTrajectoryBuilder::trajectories(), and L3MuonTrajectoryBuilder::trajectories().

◆ trackProbability()

double GlobalTrajectoryBuilderBase::trackProbability ( const Trajectory & track ) const

protected

calculate chi2 probability (-ln(P))

Definition at line 363 of file GlobalTrajectoryBuilderBase.cc.

                                                                                   {
   if (track.ndof() > 0 && track.chiSquared() > 0) {
     return -LnChiSquaredProbability(track.chiSquared(), track.ndof());
   } else {
     return 0.0;
   }
 }

References LnChiSquaredProbability(), and HLT_2018_cff::track.

Referenced by build().

◆ trajectories()

MuonTrajectoryBuilder::TrajectoryContainer GlobalTrajectoryBuilderBase::trajectories ( const TrajectorySeed & )

inlineoverridevirtual

dummy implementation, unused in this class

Implements MuonTrajectoryBuilder.

Definition at line 62 of file GlobalTrajectoryBuilderBase.h.

                                                                                         {
     return MuonTrajectoryBuilder::TrajectoryContainer();
   }