#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

typedef MuonTransientTrackingRecHit::ConstMuonRecHitContainer GlobalTrajectoryBuilderBase::ConstMuonRecHitContainer

Definition at line 53 of file GlobalTrajectoryBuilderBase.h.

typedef MuonTransientTrackingRecHit::ConstMuonRecHitPointer GlobalTrajectoryBuilderBase::ConstMuonRecHitPointer

Definition at line 51 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::ConstRecHitContainer GlobalTrajectoryBuilderBase::ConstRecHitContainer

Definition at line 47 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::ConstRecHitPointer GlobalTrajectoryBuilderBase::ConstRecHitPointer

Definition at line 49 of file GlobalTrajectoryBuilderBase.h.

typedef MuonTransientTrackingRecHit::MuonRecHitContainer GlobalTrajectoryBuilderBase::MuonRecHitContainer

Definition at line 52 of file GlobalTrajectoryBuilderBase.h.

typedef MuonTransientTrackingRecHit::MuonRecHitPointer GlobalTrajectoryBuilderBase::MuonRecHitPointer

Definition at line 50 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::RecHitContainer GlobalTrajectoryBuilderBase::RecHitContainer

Definition at line 46 of file GlobalTrajectoryBuilderBase.h.

typedef TransientTrackingRecHit::RecHitPointer GlobalTrajectoryBuilderBase::RecHitPointer

Definition at line 48 of file GlobalTrajectoryBuilderBase.h.

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

Definition at line 54 of file GlobalTrajectoryBuilderBase.h.

typedef TC::const_iterator GlobalTrajectoryBuilderBase::TI

Definition at line 55 of file GlobalTrajectoryBuilderBase.h.

Member Enumeration Documentation

enum GlobalTrajectoryBuilderBase::RefitDirection

protected

Enumerator
inToOut
outToIn
undetermined

Definition at line 71 of file GlobalTrajectoryBuilderBase.h.

71 {inToOut,outToIn,undetermined};

GlobalTrajectoryBuilderBase::undetermined

Definition: GlobalTrajectoryBuilderBase.h:71

GlobalTrajectoryBuilderBase::outToIn

Definition: GlobalTrajectoryBuilderBase.h:71

GlobalTrajectoryBuilderBase::inToOut

Definition: GlobalTrajectoryBuilderBase.h:71

Constructor & Destructor Documentation

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.

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

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

GlobalTrajectoryBuilderBase::~GlobalTrajectoryBuilderBase ( )

override

destructor

Definition at line 121 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 115 of file GlobalTrajectoryBuilderBase.h.

References checkRecHitsOrdering(), and getTransientRecHits().

115 {}

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

protected

build combined trajectory from sta Track and tracker RecHits

Definition at line 163 of file GlobalTrajectoryBuilderBase.cc.

References funct::abs(), checkRecHitsOrdering(), fixTEC(), getTransientRecHits(), LogDebug, LogTrace, MuonServiceProxy::magneticField(), SiStripPI::max, conversionTrackMerger_cfi::minProb, MuonCandidate::muonTrack(), nullptr, outToIn, TtFullHadEvtBuilder_cfi::prob, GlobalMuonRefitter::refit(), groupFilesInBlocks::reverse, selectTrackerHits(), Trajectory::setSeedRef(), theCategory, theGlbRefitter, theMuonHitsOption, thePCut, thePtCut, theRefitFlag, theService, theTECxScale, theTECyScale, theTopo, theTrackTransformer, MuonCandidate::trackerTrack(), MuonCandidate::trackerTrajectory(), MuonServiceProxy::trackingGeometry(), trackProbability(), MuonCandidate::trajectory(), and TrackTransformer::transform().

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

                                                                                     {
 
   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;
         MuonCandidate* finalTrajectory = nullptr;
         Trajectory *tkTrajectory = nullptr;
     
         // tracker only track
         if ( ! (it->trackerTrajectory() && it->trackerTrajectory()->isValid()) ) {
           refitted0 = theTrackTransformer->transform(it->trackerTrack()) ;
           if (!refitted0.empty()) tkTrajectory = new Trajectory(*(refitted0.begin())); 
           else edm::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, theTopo);
         LogTrace(theCategory)<<"     This track-sta refitted to " << refitted1.size() << " trajectories";
     
         Trajectory *glbTrajectory1 = nullptr;
         if (!refitted1.empty()) glbTrajectory1 = new Trajectory(*(refitted1.begin()));
         else LogDebug(theCategory)<< "     Failed to load global track trajectory 1"; 
         if (glbTrajectory1) glbTrajectory1->setSeedRef(tmpSeed);
         
         finalTrajectory = nullptr;
         if(glbTrajectory1 && tkTrajectory) finalTrajectory = new MuonCandidate(glbTrajectory1, it->muonTrack(), it->trackerTrack(),
                         tkTrajectory? new Trajectory(*tkTrajectory) : nullptr);
     
         if ( finalTrajectory ) refittedResult.push_back(finalTrajectory);
         if(tkTrajectory) delete tkTrajectory;
       }
       else{
         MuonCandidate* finalTrajectory = nullptr;
         edm::RefToBase<TrajectorySeed> tmpSeed;
         if(it->trackerTrack()->seedRef().isAvailable()) tmpSeed = it->trackerTrack()->seedRef();
     
         TC refitted0;
         Trajectory *tkTrajectory = nullptr;
         if ( ! (it->trackerTrajectory() && it->trackerTrajectory()->isValid()) ) {
           refitted0 = theTrackTransformer->transform(it->trackerTrack());
           if (!refitted0.empty()){
               tkTrajectory = new Trajectory(*(refitted0.begin()));
           }
           else edm::LogWarning(theCategory)<< "     Failed to load tracker track trajectory";
         }
         else tkTrajectory = it->trackerTrajectory();
         if (tkTrajectory) tkTrajectory->setSeedRef(tmpSeed);
         // Creating MuonCandidate using only the tracker trajectory:
         finalTrajectory = new MuonCandidate(new Trajectory(*tkTrajectory), it->muonTrack(), it->trackerTrack(),new Trajectory(*tkTrajectory));
         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 = nullptr;
   if ( !refittedResult.empty() ) tmpCand = *(refittedResult.begin());
   double minProb = std::numeric_limits<double>::max();
 
   for (auto&& iter: refittedResult){
     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()) ):nullptr ) );
 
   for (auto&& it: refittedResult){
     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 409 of file GlobalTrajectoryBuilderBase.cc.

References inToOut, outToIn, theCategory, and undetermined.

Referenced by addTraj(), and 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 {
       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;
   }
 }

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.

References defineRegionOfInterest(), deltaR(), TrackingRegion::direction(), PV3DBase< T, PVType, FrameType >::eta(), PV3DBase< T, PVType, FrameType >::phi(), and mps_fire::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();
 
   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; 
 
 }

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

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

Referenced by chooseRegionalTrackerTracks().

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

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

References angle(), MuonTransientTrackingRecHit::build(), SiStripRecHit2D::cluster(), TrackingRecHit::det(), relativeConstraints::error, mps_fire::i, BaseTrackerRecHit::localPosition(), BaseTrackerRecHit::localPositionError(), LocalError::rotate(), digitizers_cfi::strip, 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 528 of file GlobalTrajectoryBuilderBase.cc.

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

Referenced by addTraj(), and 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 383 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().

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

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

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

References Vispa.Plugins.EdmBrowser.EdmDataAccessor::all(), hfClusterShapes_cfi::hits, mps_fire::i, StripSubdetector::TEC, and DetId::Tracker.

Referenced by build().

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

const MuonServiceProxy* GlobalTrajectoryBuilderBase::service ( ) const

inlineprotected

Definition at line 128 of file GlobalTrajectoryBuilderBase.h.

References theService.

128 { return theService; }

GlobalTrajectoryBuilderBase::theService

const MuonServiceProxy * theService

Definition: GlobalTrajectoryBuilderBase.h:160

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

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

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

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

GlobalMuonTrackMatcher* GlobalTrajectoryBuilderBase::trackMatcher ( ) const

inlineprotected

Definition at line 125 of file GlobalTrajectoryBuilderBase.h.

References theTrackMatcher.

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

125 { return theTrackMatcher; }

GlobalTrajectoryBuilderBase::theTrackMatcher

GlobalMuonTrackMatcher * theTrackMatcher

Definition: GlobalTrajectoryBuilderBase.h:156

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

protected

calculate chi2 probability (-ln(P))

Definition at line 369 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;
   }
 
 }