GlobalTrajectoryBuilderBase Class Reference

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

{inToOut,outToIn,undetermined};

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().

{} 

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(), 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.

{ return theService; }

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().

{ return theTrackMatcher; }

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

}

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

inlineoverridevirtual

dummy implementation, unused in this class

Implements MuonTrajectoryBuilder.

Definition at line 64 of file GlobalTrajectoryBuilderBase.h.

References setEvent().

{ return MuonTrajectoryBuilder::TrajectoryContainer(); }

Member Data Documentation

unsigned long long GlobalTrajectoryBuilderBase::theCacheId_TRH

private

Definition at line 162 of file GlobalTrajectoryBuilderBase.h.

Referenced by GlobalTrajectoryBuilderBase(), and setEvent().

std::string GlobalTrajectoryBuilderBase::theCategory

protected

Definition at line 150 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), checkRecHitsOrdering(), getTransientRecHits(), GlobalTrajectoryBuilderBase(), printHits(), and setEvent().

const edm::Event* GlobalTrajectoryBuilderBase::theEvent

private

Definition at line 169 of file GlobalTrajectoryBuilderBase.h.

Referenced by setEvent().

GlobalMuonRefitter* GlobalTrajectoryBuilderBase::theGlbRefitter

private

Definition at line 161 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), GlobalTrajectoryBuilderBase(), setEvent(), and ~GlobalTrajectoryBuilderBase().

MuonDetLayerMeasurements* GlobalTrajectoryBuilderBase::theLayerMeasurements

private

Definition at line 157 of file GlobalTrajectoryBuilderBase.h.

int GlobalTrajectoryBuilderBase::theMuonHitsOption

private

Definition at line 165 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), and GlobalTrajectoryBuilderBase().

edm::ESHandle<TransientTrackingRecHitBuilder> GlobalTrajectoryBuilderBase::theMuonRecHitBuilder

private

Definition at line 173 of file GlobalTrajectoryBuilderBase.h.

Referenced by getTransientRecHits(), and setEvent().

std::string GlobalTrajectoryBuilderBase::theMuonRecHitBuilderName

private

Definition at line 172 of file GlobalTrajectoryBuilderBase.h.

Referenced by GlobalTrajectoryBuilderBase(), and setEvent().

float GlobalTrajectoryBuilderBase::thePCut

protected

Definition at line 152 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), and GlobalTrajectoryBuilderBase().

float GlobalTrajectoryBuilderBase::thePtCut

protected

Definition at line 151 of file GlobalTrajectoryBuilderBase.h.

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

bool GlobalTrajectoryBuilderBase::theRefitFlag

private

Definition at line 164 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), and GlobalTrajectoryBuilderBase().

MuonTrackingRegionBuilder* GlobalTrajectoryBuilderBase::theRegionBuilder

private

Definition at line 159 of file GlobalTrajectoryBuilderBase.h.

Referenced by defineRegionOfInterest(), GlobalTrajectoryBuilderBase(), setEvent(), and ~GlobalTrajectoryBuilderBase().

bool GlobalTrajectoryBuilderBase::theRPCInTheFit

private

Definition at line 163 of file GlobalTrajectoryBuilderBase.h.

Referenced by getTransientRecHits(), and GlobalTrajectoryBuilderBase().

const MuonServiceProxy* GlobalTrajectoryBuilderBase::theService

private

Definition at line 160 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), getTransientRecHits(), GlobalTrajectoryBuilderBase(), service(), and setEvent().

float GlobalTrajectoryBuilderBase::theTECxScale

private

Definition at line 166 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), and GlobalTrajectoryBuilderBase().

float GlobalTrajectoryBuilderBase::theTECyScale

private

Definition at line 167 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), and GlobalTrajectoryBuilderBase().

const TrackerTopology* GlobalTrajectoryBuilderBase::theTopo

private

Definition at line 174 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), and setEvent().

std::string GlobalTrajectoryBuilderBase::theTrackerPropagatorName

private

Definition at line 168 of file GlobalTrajectoryBuilderBase.h.

Referenced by getTransientRecHits(), and GlobalTrajectoryBuilderBase().

edm::ESHandle<TransientTrackingRecHitBuilder> GlobalTrajectoryBuilderBase::theTrackerRecHitBuilder

private

Definition at line 171 of file GlobalTrajectoryBuilderBase.h.

Referenced by getTransientRecHits(), and setEvent().

std::string GlobalTrajectoryBuilderBase::theTrackerRecHitBuilderName

private

Definition at line 170 of file GlobalTrajectoryBuilderBase.h.

Referenced by GlobalTrajectoryBuilderBase(), and setEvent().

GlobalMuonTrackMatcher* GlobalTrajectoryBuilderBase::theTrackMatcher

private

Definition at line 156 of file GlobalTrajectoryBuilderBase.h.

Referenced by GlobalTrajectoryBuilderBase(), trackMatcher(), and ~GlobalTrajectoryBuilderBase().

TrackTransformer* GlobalTrajectoryBuilderBase::theTrackTransformer

private

Definition at line 158 of file GlobalTrajectoryBuilderBase.h.

Referenced by build(), GlobalTrajectoryBuilderBase(), setEvent(), and ~GlobalTrajectoryBuilderBase().