#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...

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

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

virtual	~GlobalTrajectoryBuilderBase ()
	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

MuonTrackingRegionBuilder *	theRegionBuilder

bool	theRPCInTheFit

const MuonServiceProxy *	theService

float	theTECxScale

float	theTECyScale

std::string	theTrackerPropagatorName

edm::ESHandle < TransientTrackingRecHitBuilder >	theTrackerRecHitBuilder

std::string	theTrackerRecHitBuilderName

GlobalMuonTrackMatcher *	theTrackMatcher

TrackTransformer *	theTrackTransformer

const TrackerTopology *	tTopo_

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 49 of file GlobalTrajectoryBuilderBase.h.

Member Typedef Documentation

typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer GlobalTrajectoryBuilderBase::ConstMuonRecHitContainer

Definition at line 61 of file GlobalTrajectoryBuilderBase.h.

typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer GlobalTrajectoryBuilderBase::ConstMuonRecHitPointer

Definition at line 59 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::ConstRecHitContainer GlobalTrajectoryBuilderBase::ConstRecHitContainer

Definition at line 54 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::ConstRecHitPointer GlobalTrajectoryBuilderBase::ConstRecHitPointer

Definition at line 56 of file GlobalTrajectoryBuilderBase.h.

typedef MuonTransientTrackingRecHit::MuonRecHitContainer GlobalTrajectoryBuilderBase::MuonRecHitContainer

Definition at line 60 of file GlobalTrajectoryBuilderBase.h.

typedef MuonTransientTrackingRecHit::MuonRecHitPointer GlobalTrajectoryBuilderBase::MuonRecHitPointer

Definition at line 58 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::RecHitContainer GlobalTrajectoryBuilderBase::RecHitContainer

Definition at line 53 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::RecHitPointer GlobalTrajectoryBuilderBase::RecHitPointer

Definition at line 55 of file GlobalTrajectoryBuilderBase.h.

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

Definition at line 63 of file GlobalTrajectoryBuilderBase.h.

typedef TC::const_iterator GlobalTrajectoryBuilderBase::TI

Definition at line 64 of file GlobalTrajectoryBuilderBase.h.

Member Enumeration Documentation

enum GlobalTrajectoryBuilderBase::RefitDirection

protected

Enumerator
inToOut
outToIn
undetermined

Definition at line 82 of file GlobalTrajectoryBuilderBase.h.

82 {inToOut,outToIn,undetermined};

GlobalTrajectoryBuilderBase::undetermined

Definition: GlobalTrajectoryBuilderBase.h:82

GlobalTrajectoryBuilderBase::outToIn

Definition: GlobalTrajectoryBuilderBase.h:82

GlobalTrajectoryBuilderBase::inToOut

Definition: GlobalTrajectoryBuilderBase.h:82

Constructor & Destructor Documentation

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

constructor with Parameter Set and MuonServiceProxy

Definition at line 86 of file GlobalTrajectoryBuilderBase.cc.

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

                                                        : 
   theTrackMatcher(0),theLayerMeasurements(0),theTrackTransformer(0),theRegionBuilder(0), theService(service),theGlbRefitter(0) {
 
   theCategory = par.getUntrackedParameter<string>("Category", "Muon|RecoMuon|GlobalMuon|GlobalTrajectoryBuilderBase");
 
   
   ParameterSet trackMatcherPSet = par.getParameter<ParameterSet>("GlobalMuonTrackMatcher");
   theTrackMatcher = new GlobalMuonTrackMatcher(trackMatcherPSet,theService);
   
   theTrackerPropagatorName = par.getParameter<string>("TrackerPropagator");
 
   ParameterSet trackTransformerPSet = par.getParameter<ParameterSet>("TrackTransformer");
   theTrackTransformer = new TrackTransformer(trackTransformerPSet);
 
   ParameterSet regionBuilderPSet = par.getParameter<ParameterSet>("MuonTrackingRegionBuilder");
 
   theRegionBuilder = new MuonTrackingRegionBuilder(regionBuilderPSet,theService,iC);
 
   // TrackRefitter parameters
   ParameterSet refitterParameters = par.getParameter<ParameterSet>("GlbRefitterParameters");
   theGlbRefitter = new GlobalMuonRefitter(refitterParameters, theService, iC);
 
   theMuonHitsOption = refitterParameters.getParameter<int>("MuonHitsOption");
 
   theTrackerRecHitBuilderName = par.getParameter<string>("TrackerRecHitBuilder");
   theMuonRecHitBuilderName = par.getParameter<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;
 
 }

GlobalTrajectoryBuilderBase::~GlobalTrajectoryBuilderBase ( )

virtual

destructor

Definition at line 130 of file GlobalTrajectoryBuilderBase.cc.

References theGlbRefitter, theRegionBuilder, theTrackMatcher, and theTrackTransformer.

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

Member Function Documentation

void GlobalTrajectoryBuilderBase::addTraj ( TrackCand & )

inlineprotected

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

Definition at line 126 of file GlobalTrajectoryBuilderBase.h.

126 {}

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

protected

build combined trajectory from sta Track and tracker RecHits

Definition at line 172 of file GlobalTrajectoryBuilderBase.cc.

References checkRecHitsOrdering(), fixTEC(), getTransientRecHits(), getDQMSummary::iter, LogDebug, LogTrace, MuonCandidate::muonTrack(), outToIn, GlobalMuonRefitter::refit(), selectTrackerHits(), Trajectory::setSeedRef(), theCategory, theGlbRefitter, theMuonHitsOption, thePCut, thePtCut, theService, theTECxScale, theTECyScale, theTrackTransformer, MuonCandidate::trackerTrack(), MuonCandidate::trackerTrajectory(), trackProbability(), MuonCandidate::trajectory(), TrackTransformer::transform(), and tTopo_.

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

                                                                                     {
 
   LogTrace(theCategory) << " Begin Build" << 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 ( CandidateContainer::const_iterator it = tkTrajs.begin(); it != tkTrajs.end(); it++ ) {
 
     // 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;
 
     ConstRecHitContainer trackerRecHits;
     if ((*it)->trackerTrack().isNonnull()) {
       trackerRecHits = getTransientRecHits(*(*it)->trackerTrack());
     } else {
       LogDebug(theCategory)<<"     NEED HITS FROM TRAJ";
       //trackerRecHits = (*it)->trackerTrajectory()->recHits();
     }
 
     // check for single TEC RecHits in trajectories in the overalp region
     if ( fabs((*it)->trackerTrack()->eta()) > 0.95 && fabs((*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;
     MuonCandidate* finalTrajectory = 0;
     Trajectory *tkTrajectory = 0;
 
     // tracker only track
     if ( ! ((*it)->trackerTrajectory() && (*it)->trackerTrajectory()->isValid()) ) { 
       refitted0 = theTrackTransformer->transform((*it)->trackerTrack()) ;
       if (!refitted0.empty()) tkTrajectory = new Trajectory(*(refitted0.begin())); 
       else LogWarning(theCategory)<< "     Failed to load tracker track trajectory";
     } else tkTrajectory = (*it)->trackerTrajectory();
     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, tTopo_);
     LogTrace(theCategory)<<"     This track-sta refitted to " << refitted1.size() << " trajectories";
 
     Trajectory *glbTrajectory1 = 0;
     if (!refitted1.empty()) glbTrajectory1 = new Trajectory(*(refitted1.begin()));
     else LogDebug(theCategory)<< "     Failed to load global track trajectory 1"; 
     if (glbTrajectory1) glbTrajectory1->setSeedRef(tmpSeed);
     
     finalTrajectory = 0;
     if(glbTrajectory1 && tkTrajectory) finalTrajectory = new MuonCandidate(glbTrajectory1, (*it)->muonTrack(), (*it)->trackerTrack(), 
                     tkTrajectory? new Trajectory(*tkTrajectory) : 0);
 
     if ( finalTrajectory ) 
       refittedResult.push_back(finalTrajectory);
      
     if(tkTrajectory) delete tkTrajectory;
   }
 
   // choose the best global fit for this Standalone Muon based on the track probability
   CandidateContainer selectedResult;
   MuonCandidate* tmpCand = 0;
   if ( refittedResult.size() > 0 ) tmpCand = *(refittedResult.begin());
   double minProb = 9999;
 
   for (CandidateContainer::const_iterator iter=refittedResult.begin(); iter != refittedResult.end(); iter++) {
     double prob = trackProbability(*(*iter)->trajectory());
     LogTrace(theCategory)<<"   refitted-track-sta with pT " << (*iter)->trackerTrack()->pt() << " has probability " << prob;
 
     if (prob < minProb) {
       minProb = prob;
       tmpCand = (*iter);
     }
   }
 
   if ( tmpCand )  selectedResult.push_back(new MuonCandidate(new Trajectory(*(tmpCand->trajectory())), tmpCand->muonTrack(), tmpCand->trackerTrack(), 
                                  (tmpCand->trackerTrajectory())? new Trajectory( *(tmpCand->trackerTrajectory()) ):0 ) );
 
   for (CandidateContainer::const_iterator it = refittedResult.begin(); it != refittedResult.end(); ++it) {
     if ( (*it)->trajectory() ) delete (*it)->trajectory();
     if ( (*it)->trackerTrajectory() ) delete (*it)->trackerTrajectory();
     if ( *it ) delete (*it);
   }
   refittedResult.clear();
 
   return selectedResult;
 
 }

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

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

protected

This does nothing now.

check order of RechIts on a trajectory

Definition at line 402 of file GlobalTrajectoryBuilderBase.cc.

References inToOut, outToIn, theCategory, and undetermined.

Referenced by build().

                                                                                                                   {
 
   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 {
       LogError(theCategory) << "Impossible to determine the rechits order" << endl;
       return undetermined;
     }
   }
   else {
     LogError(theCategory) << "Impossible to determine the rechits order" << endl;
     return undetermined;
   }
 }

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

protected

choose tracker tracks within region of interest

Definition at line 297 of file GlobalTrajectoryBuilderBase.cc.

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

Referenced by GlobalMuonTrajectoryBuilder::makeTkCandCollection().

                                                                                         {
   
   // 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();
 
   vector<TrackCand> result;
 
   double deltaR_max = 1.0;
 
   for ( vector<TrackCand>::const_iterator is = tkTs.begin(); is != tkTs.end(); ++is ) {
     // check if each trackCand is in region of interest
 //    bool inEtaRange = etaRange.inside(is->second->eta());
 //    bool inPhiRange = (fabs(Geom::Phi<float>(is->second->phi()) - Geom::Phi<float>(regionOfInterest.direction().phi())) < phiMargin.right() ) ? true : false ;
 
     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; 
 
 }

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

protected

choose final trajectory

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

protected

define region of interest with tracker

Definition at line 337 of file GlobalTrajectoryBuilderBase.cc.

References TrackingRegion::direction(), RectangularEtaPhiTrackingRegion::etaRange(), PixelRecoRange< T >::max(), PixelRecoRange< T >::mean(), PixelRecoRange< T >::min(), TrackingRegion::origin(), TrackingRegion::originRBound(), TrackingRegion::originZBound(), RectangularEtaPhiTrackingRegion::phiMargin(), TrackingRegion::ptMin(), MuonTrackingRegionBuilder::region(), and theRegionBuilder.

Referenced by chooseRegionalTrackerTracks().

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

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 455 of file GlobalTrajectoryBuilderBase.cc.

References angle(), MuonTransientTrackingRecHit::build(), SiStripRecHit2D::cluster(), TrackingRecHit::det(), relativeConstraints::error, i, BaseTrackerRecHit::localPosition(), BaseTrackerRecHit::localPositionError(), LocalError::rotate(), StripTopology::strip(), StripTopology::stripAngle(), StripSubdetector::TEC, ecaldqm::topology(), DetId::Tracker, LocalError::xx(), and LocalError::yy().

Referenced by build().

                                                              {
 
   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);
  
         }
       }
     }
   }
 
 }

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

protected

get transient RecHits of a Track

Definition at line 521 of file GlobalTrajectoryBuilderBase.cc.

References TkTransientTrackingRecHitBuilder::cloner(), DetId::det(), h, trajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), LogDebug, DetId::Muon, edm::ESHandle< class >::product(), reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), query::result, theCategory, theMuonRecHitBuilder, theRPCInTheFit, theService, theTrackerPropagatorName, theTrackerRecHitBuilder, and DetId::Tracker.

Referenced by build().

                                                                              {
 
   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;
 }

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

protectedpure virtual

make a TrackCand collection using tracker Track, Trajectory information

Implemented in L3MuonTrajectoryBuilder, and GlobalMuonTrajectoryBuilder.

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

protected

print all RecHits of a trajectory

Definition at line 376 of file GlobalTrajectoryBuilderBase.cc.

References LogTrace, mathSSE::sqrt(), theCategory, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                   {
 
   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()
       << "  " << (*ir)->det()->geographicalId().det()
       << "  " << (*ir)->det()->subDetector();
 
   }
 
 }

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

protected

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

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

protected

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

Definition at line 431 of file GlobalTrajectoryBuilderBase.cc.

References cond::ecalcond::all, i, StripSubdetector::TEC, and DetId::Tracker.

Referenced by build().

                                                                                     {
  
   int nTEC(0);
 
   ConstRecHitContainer hits;
   for (ConstRecHitContainer::const_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) {
       nTEC++;
     } else {
       hits.push_back(*i);
     }
     if ( nTEC > 1 ) return all;
   }
 
   return hits;
 
 }

const MuonServiceProxy* GlobalTrajectoryBuilderBase::service ( ) const

inlineprotected

Definition at line 139 of file GlobalTrajectoryBuilderBase.h.

References theService.

139 { return theService; }

GlobalTrajectoryBuilderBase::theService

const MuonServiceProxy * theService

Definition: GlobalTrajectoryBuilderBase.h:173

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

virtual

pass the Event to the algo at each event

Implements MuonTrajectoryBuilder.

Reimplemented in L3MuonTrajectoryBuilder, and GlobalMuonTrajectoryBuilder.

Definition at line 142 of file GlobalTrajectoryBuilderBase.cc.

References event(), LogDebug, edm::ESHandle< class >::product(), MuonTrackingRegionBuilder::setEvent(), GlobalMuonRefitter::setEvent(), TrackTransformer::setServices(), GlobalMuonRefitter::setServices(), theCacheId_TRH, theCategory, theEvent, theGlbRefitter, theMuonRecHitBuilder, theMuonRecHitBuilderName, theRegionBuilder, theService, theTrackerRecHitBuilder, theTrackerRecHitBuilderName, theTrackTransformer, and tTopo_.

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

                                                                 {
   
   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<IdealGeometryRecord>().get(tTopoHand);
   tTopo_=tTopoHand.product();
 
 }

GlobalMuonTrackMatcher* GlobalTrajectoryBuilderBase::trackMatcher ( ) const

inlineprotected

Definition at line 136 of file GlobalTrajectoryBuilderBase.h.

References theTrackMatcher.

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

136 { return theTrackMatcher; }

GlobalTrajectoryBuilderBase::theTrackMatcher

GlobalMuonTrackMatcher * theTrackMatcher

Definition: GlobalTrajectoryBuilderBase.h:169

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

protected

calculate chi2 probability (-ln(P))

Definition at line 362 of file GlobalTrajectoryBuilderBase.cc.

References Trajectory::chiSquared(), LnChiSquaredProbability(), and Trajectory::ndof().

Referenced by build().

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