TSGForRoadSearch Class Reference

#include <TSGForRoadSearch.h>

Inheritance diagram for TSGForRoadSearch:

Public Types
typedef std::vector< TrajectorySeed >	BTSeedCollection
typedef std::pair< const Trajectory *, reco::TrackRef >	TrackCand
Public Types inherited from TrackerSeedGenerator
typedef std::vector< TrajectorySeed >	BTSeedCollection
typedef std::pair< const Trajectory *, reco::TrackRef >	TrackCand

Public Member Functions
void	init (const MuonServiceProxy *service) override
	initialize the service More...
void	setEvent (const edm::Event &event) override
	set the event: update the MeasurementTracker More...
void	trackerSeeds (const TrackCand &, const TrackingRegion &, const TrackerTopology *, BTSeedCollection &) override
	generated seed(s) for a track. the tracking region is not used. More...
	TSGForRoadSearch (const edm::ParameterSet &pset, edm::ConsumesCollector &IC)
	~TSGForRoadSearch () override
Public Member Functions inherited from TrackerSeedGenerator
const edm::Event *	getEvent () const
	TrackerSeedGenerator ()
virtual	~TrackerSeedGenerator ()
	destructor More...

Private Member Functions
bool	IPfts (const reco::Track &, FreeTrajectoryState &)
	get the FTS for a Track: adjusting the error matrix if requested More...
void	makeSeeds_0 (const reco::Track &, std::vector< TrajectorySeed > &)
	oseed from inside-out: innermost Strip layer More...
void	makeSeeds_1 (const reco::Track &, std::vector< TrajectorySeed > &)
	not implemented More...
void	makeSeeds_2 (const reco::Track &, std::vector< TrajectorySeed > &)
	not implemented More...
void	makeSeeds_3 (const reco::Track &, std::vector< TrajectorySeed > &)
	outside-in: outermost Strip layer More...
void	makeSeeds_4 (const reco::Track &, std::vector< TrajectorySeed > &)
	inside-out: innermost Pixel/Strip layer More...
bool	notAtIPtsos (TrajectoryStateOnSurface &state)
	make the adjustement away from PCA state if requested More...
void	pushTrajectorySeed (const reco::Track &muon, std::vector< DetLayer::DetWithState > &compatible, PropagationDirection direction, std::vector< TrajectorySeed > &result) const
	add the seed(s) to the collection of seeds More...

Private Attributes
bool	theAdjustAtIp
	adjust the state at IP or where it is defined for the seed More...
std::string	theCategory
Chi2MeasurementEstimator *	theChi2Estimator
edm::ParameterSet	theConfig
bool	theCopyMuonRecHit
MuonErrorMatrix *	theErrorMatrixAdjuster
edm::ESHandle< GeometricSearchTracker >	theGeometricSearchTracker
edm::ESGetToken< GeometricSearchTracker, TrackerRecoGeometryRecord >	theGeometricSearchTrackerToken
bool	theManySeeds
const MeasurementTrackerEvent *	theMeasurementTrackerEvent
edm::InputTag	theMeasurementTrackerEventTag
edm::EDGetTokenT< MeasurementTrackerEvent >	theMeasurementTrackerEventToken
unsigned int	theOption
std::string	thePropagatorCompatibleName
std::string	thePropagatorName
const MuonServiceProxy *	theProxyService
TrajectoryStateUpdator *	theUpdator

Additional Inherited Members
Protected Attributes inherited from TrackerSeedGenerator
const edm::Event *	theEvent
const MuonServiceProxy *	theProxyService

Detailed Description

Description: this class generates hit-less TrajectorySeed from a given Track. the original error matrix of the Track is adjusted (configurable). this class is principally used for muon HLT. for options are available:

• inside-out seeds, on innermost Pixel/Strip layer.
• inside-out seeds, on innermost Strip layer.
• outside-in seeds, on outermost Strip layer.

regular operation is one seed per track, but optionnaly, more than one seed can be madefor one track.

Author

Jean-Roch Vlimant

Definition at line 46 of file TSGForRoadSearch.h.

Member Typedef Documentation

BTSeedCollection

typedef std::vector<TrajectorySeed> TSGForRoadSearch::BTSeedCollection

Definition at line 48 of file TSGForRoadSearch.h.

TrackCand

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

Definition at line 49 of file TSGForRoadSearch.h.

Constructor & Destructor Documentation

TSGForRoadSearch()

TSGForRoadSearch::TSGForRoadSearch	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &	IC
	)

Definition at line 26 of file TSGForRoadSearch.cc.

References Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, edm::ConsumesCollector::consumes(), HLT_2022v12_cff::errorMatrixPset, edm::ParameterSet::getParameter(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, theAdjustAtIp, theCategory, theChi2Estimator, theCopyMuonRecHit, theErrorMatrixAdjuster, theManySeeds, theMeasurementTrackerEvent, theMeasurementTrackerEventTag, theMeasurementTrackerEventToken, theOption, thePropagatorCompatibleName, thePropagatorName, and theUpdator.

    : theGeometricSearchTrackerToken(iC.esConsumes()) {
  theOption = par.getParameter<unsigned int>("option");
  theCopyMuonRecHit = par.getParameter<bool>("copyMuonRecHit");

  double Chi2 = par.getParameter<double>("maxChi2");
  theChi2Estimator = new Chi2MeasurementEstimator(Chi2, sqrt(Chi2));

  thePropagatorName = par.getParameter<std::string>("propagatorName");
  thePropagatorCompatibleName = par.getParameter<std::string>("propagatorCompatibleName");

  theCategory = "TSGForRoadSearch|TrackerSeedGenerator";

  theManySeeds = par.getParameter<bool>("manySeeds");
  if (theManySeeds) {
    theUpdator = new KFUpdator();
  } else {
    theUpdator = nullptr;
  }

  edm::ParameterSet errorMatrixPset = par.getParameter<edm::ParameterSet>("errorMatrixPset");
  if (!errorMatrixPset.empty()) {
    theAdjustAtIp = errorMatrixPset.getParameter<bool>("atIP");
    theErrorMatrixAdjuster = new MuonErrorMatrix(errorMatrixPset);
  } else {
    theAdjustAtIp = false;
    theErrorMatrixAdjuster = nullptr;
  }

  theMeasurementTrackerEventTag = par.getParameter<edm::InputTag>("MeasurementTrackerEvent");
  theMeasurementTrackerEvent = nullptr;

  theMeasurementTrackerEventToken = iC.consumes<MeasurementTrackerEvent>(theMeasurementTrackerEventTag);
}

~TSGForRoadSearch()

TSGForRoadSearch::~TSGForRoadSearch ( )

override

Definition at line 61 of file TSGForRoadSearch.cc.

References theChi2Estimator, theErrorMatrixAdjuster, and theUpdator.

                                    {
  delete theChi2Estimator;
  if (theUpdator)
    delete theUpdator;
  if (theErrorMatrixAdjuster)
    delete theErrorMatrixAdjuster;
}

Member Function Documentation

init()

void TSGForRoadSearch::init ( const MuonServiceProxy *  service )

overridevirtual

initialize the service

Reimplemented from TrackerSeedGenerator.

Definition at line 69 of file TSGForRoadSearch.cc.

References theProxyService.

{ theProxyService = service; }

IPfts()

bool TSGForRoadSearch::IPfts ( const reco::Track &  muon, FreeTrajectoryState &  fts  )

private

get the FTS for a Track: adjusting the error matrix if requested

Definition at line 112 of file TSGForRoadSearch.cc.

References MuonErrorMatrix::adjust(), trajectoryStateTransform::initialFreeState(), LogDebug, MuonServiceProxy::magneticField(), theAdjustAtIp, theCategory, theErrorMatrixAdjuster, and theProxyService.

Referenced by makeSeeds_0(), makeSeeds_3(), and makeSeeds_4().

                                                                            {
  fts = trajectoryStateTransform::initialFreeState(muon, &*theProxyService->magneticField());
  LogDebug(theCategory) << "pure L2 state: " << fts;
  if (fts.position().mag() == 0 && fts.momentum().mag() == 0) {
    edm::LogError(theCategory) << "initial state of muon is (0,0,0)(0,0,0). no seed.";
    return false;
  }

  //rescale the error at IP
  if (theErrorMatrixAdjuster && theAdjustAtIp) {
    theErrorMatrixAdjuster->adjust(fts);
    LogDebug(theCategory) << "after adjusting the error matrix: " << fts;
  }

  return true;
}

makeSeeds_0()

void TSGForRoadSearch::makeSeeds_0 ( const reco::Track &  muon, std::vector< TrajectorySeed > &  result  )

private

oseed from inside-out: innermost Strip layer

Definition at line 132 of file TSGForRoadSearch.cc.

References alongMomentum, GeometricSearchDet::compatibleDetsV(), SurfaceOrientation::inner, IPfts(), LogDebug, GeometricSearchTracker::negTecLayers(), GeometricSearchTracker::negTidLayers(), notAtIPtsos(), GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, GeometricSearchTracker::posTecLayers(), GeometricSearchTracker::posTidLayers(), Propagator::propagate(), MuonServiceProxy::propagator(), pushTrajectorySeed(), mps_fire::result, DetLayer::subDetector(), GeomDetEnumerators::subDetGeom, GeomDetEnumerators::TEC, theCategory, theChi2Estimator, theGeometricSearchTracker, thePropagatorCompatibleName, thePropagatorName, theProxyService, GeomDetEnumerators::TIB, GeometricSearchTracker::tibLayers(), GeomDetEnumerators::TID, GeomDetEnumerators::TOB, and z.

Referenced by makeSeeds_4(), and trackerSeeds().

                                                                                           {
  //get the state at IP
  FreeTrajectoryState cIPFTS;
  if (!IPfts(muon, cIPFTS))
    return;

  //take state at inner surface and check the first part reached
  const std::vector<const BarrelDetLayer *> &blc = theGeometricSearchTracker->tibLayers();
  TrajectoryStateOnSurface inner =
      theProxyService->propagator(thePropagatorName)->propagate(cIPFTS, blc.front()->surface());
  if (!inner.isValid()) {
    LogDebug(theCategory) << "inner state is not valid. no seed.";
    return;
  }

  //rescale the error
  if (!notAtIPtsos(inner))
    return;

  double z = inner.globalPosition().z();

  const std::vector<const ForwardDetLayer *> &ptidc = theGeometricSearchTracker->posTidLayers();
  const std::vector<const ForwardDetLayer *> &ptecc = theGeometricSearchTracker->posTecLayers();
  const std::vector<const ForwardDetLayer *> &ntidc = theGeometricSearchTracker->negTidLayers();
  const std::vector<const ForwardDetLayer *> &ntecc = theGeometricSearchTracker->negTecLayers();

  const DetLayer *inLayer = nullptr;
  if (fabs(z) < ptidc.front()->surface().position().z()) {
    inLayer = blc.front();
  } else if (fabs(z) < ptecc.front()->surface().position().z()) {
    inLayer = (z < 0) ? ntidc.front() : ptidc.front();
  } else {
    inLayer = (z < 0) ? ntecc.front() : ptecc.front();
  }

  //find out at least one compatible detector reached
  std::vector<DetLayer::DetWithState> compatible;
  compatible.reserve(10);
  inLayer->compatibleDetsV(
      inner, *theProxyService->propagator(thePropagatorCompatibleName), *theChi2Estimator, compatible);

  //loop the parts until at least a compatible is found
  while (compatible.empty()) {
    switch (GeomDetEnumerators::subDetGeom[inLayer->subDetector()]) {
      case GeomDetEnumerators::PixelBarrel:
      case GeomDetEnumerators::PixelEndcap:
      case GeomDetEnumerators::TOB:
      case GeomDetEnumerators::TEC:
        LogDebug(theCategory) << "from inside-out, trying TEC or TOB layers. no seed.";
        return;
        break;
      case GeomDetEnumerators::TIB:
        inLayer = (z < 0) ? ntidc.front() : ptidc.front();
        break;
      case GeomDetEnumerators::TID:
        inLayer = (z < 0) ? ntecc.front() : ptecc.front();
        break;
      default:
        LogDebug(theCategory) << "subdetectorid is not a tracker sub-dectector id. skipping.";
        return;
    }
    inLayer->compatibleDetsV(
        inner, *theProxyService->propagator(thePropagatorCompatibleName), *theChi2Estimator, compatible);
  }

  pushTrajectorySeed(muon, compatible, alongMomentum, result);

  return;
}

makeSeeds_1()

void TSGForRoadSearch::makeSeeds_1 ( const reco::Track &  muon, std::vector< TrajectorySeed > &  result  )

private

not implemented

Definition at line 202 of file TSGForRoadSearch.cc.

References theCategory.

Referenced by trackerSeeds().

                                                                                           {
  edm::LogError(theCategory) << "option 1 of TSGForRoadSearch is not implemented yet. Please use 0,3 or 4. no seed.";
  return;
}

makeSeeds_2()

void TSGForRoadSearch::makeSeeds_2 ( const reco::Track &  muon, std::vector< TrajectorySeed > &  result  )

private

not implemented

Definition at line 207 of file TSGForRoadSearch.cc.

References theCategory.

Referenced by trackerSeeds().

                                                                                           {
  edm::LogError(theCategory) << "option 2 of TSGForRoadSearch is not implemented yet. Please use 0,3 or 4. no seed.";
  return;
}

makeSeeds_3()

void TSGForRoadSearch::makeSeeds_3 ( const reco::Track &  muon, std::vector< TrajectorySeed > &  result  )

private

outside-in: outermost Strip layer

Definition at line 215 of file TSGForRoadSearch.cc.

References GeometricSearchDet::compatibleDetsV(), IPfts(), LogDebug, LogTrace, GeometricSearchTracker::negTecLayers(), notAtIPtsos(), oppositeToMomentum, SurfaceOrientation::outer, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, GloballyPositioned< T >::position(), position, GeometricSearchTracker::posTecLayers(), edm::ESHandle< T >::product(), MuonServiceProxy::propagator(), pushTrajectorySeed(), mps_fire::result, DetLayer::subDetector(), GeomDetEnumerators::subDetGeom, GeometricSearchDet::surface(), GeomDetEnumerators::TEC, theCategory, theChi2Estimator, theGeometricSearchTracker, thePropagatorCompatibleName, thePropagatorName, theProxyService, GeomDetEnumerators::TIB, GeomDetEnumerators::TID, GeomDetEnumerators::TOB, GeometricSearchTracker::tobLayers(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by trackerSeeds().

                                                                                           {
  //get the state at IP
  FreeTrajectoryState cIPFTS;
  if (!IPfts(muon, cIPFTS))
    return;

  //take state at outer surface and check the first part reached
  const std::vector<const BarrelDetLayer *> &blc = theGeometricSearchTracker->tobLayers();

  //  TrajectoryStateOnSurface outer = theProxyService->propagator(thePropagatorName)->propagate(cIPFTS,blc.back()->surface());
  StateOnTrackerBound onBounds(theProxyService->propagator(thePropagatorName).product());
  TrajectoryStateOnSurface outer = onBounds(cIPFTS);

  if (!outer.isValid()) {
    LogDebug(theCategory) << "outer state is not valid. no seed.";
    return;
  }

  //rescale the error
  if (!notAtIPtsos(outer))
    return;

  double z = outer.globalPosition().z();

  const std::vector<const ForwardDetLayer *> &ptecc = theGeometricSearchTracker->posTecLayers();
  const std::vector<const ForwardDetLayer *> &ntecc = theGeometricSearchTracker->negTecLayers();

  LogDebug(theCategory) << "starting looking for a compatible layer from: " << outer << "\nz: " << z
                        << "TEC1 z: " << ptecc.front()->surface().position().z();

  unsigned int layerShift = 0;
  const DetLayer *inLayer = nullptr;
  if (fabs(z) < ptecc.front()->surface().position().z()) {
    inLayer = *(blc.rbegin() + layerShift);
    LogTrace(theCategory) << "choosing TOB layer with shift: " << layerShift;
  } else {
    unsigned int tecIt = 1;
    for (; tecIt != ptecc.size(); tecIt++) {
      LogTrace(theCategory) << "checking surface with shift: " << tecIt
                            << "z: " << ptecc[tecIt]->surface().position().z();
      if (fabs(z) < ptecc[tecIt]->surface().position().z()) {
        inLayer = (z < 0) ? ntecc[tecIt - 1] : ptecc[tecIt - 1];
        layerShift = tecIt - 1;
        LogTrace(theCategory) << "choosing TEC layer with shift: " << layerShift
                              << " and z: " << inLayer->surface().position().z();
        break;
      }
    }
    if (!inLayer) {
      inLayer = (z < 0) ? ntecc.back() : ptecc.back();
      LogTrace(theCategory) << "choosing last TEC layer with z: " << inLayer->surface().position().z();
    }
  }

  //find out at least one compatible detector reached
  std::vector<DetLayer::DetWithState> compatible;
  compatible.reserve(10);
  inLayer->compatibleDetsV(
      outer, *theProxyService->propagator(thePropagatorCompatibleName), *theChi2Estimator, compatible);

  //loop the parts until at least a compatible is found
  while (compatible.empty()) {
    switch (GeomDetEnumerators::subDetGeom[inLayer->subDetector()]) {
      case GeomDetEnumerators::PixelBarrel:
      case GeomDetEnumerators::PixelEndcap:
      case GeomDetEnumerators::TIB:
      case GeomDetEnumerators::TID:
      case GeomDetEnumerators::TOB:
        layerShift++;
        if (layerShift >= blc.size()) {
          LogDebug(theCategory) << "all barrel layers are exhausted to find starting state. no seed,";
          return;
        }
        inLayer = *(blc.rbegin() + layerShift);
        break;
      case GeomDetEnumerators::TEC:
        if (layerShift == 0) {
          LogDebug(theCategory) << "failed to get a compatible module on a TEC layer, using the last TOB layer.";
          inLayer = *(blc.rbegin() + layerShift);
        } else {
          layerShift--;
          LogDebug(theCategory) << "reaching more in with layer " << layerShift << " in TEC";
          inLayer = (z < 0) ? ntecc[layerShift] : ptecc[layerShift];
        }
        break;
      default:
        edm::LogError(theCategory) << "subdetectorid is not a tracker sub-dectector id. skipping.";
        return;
    }
    inLayer->compatibleDetsV(
        outer, *theProxyService->propagator(thePropagatorCompatibleName), *theChi2Estimator, compatible);
  }

  pushTrajectorySeed(muon, compatible, oppositeToMomentum, result);

  return;
}

makeSeeds_4()

void TSGForRoadSearch::makeSeeds_4 ( const reco::Track &  muon, std::vector< TrajectorySeed > &  result  )

private

inside-out: innermost Pixel/Strip layer

Definition at line 316 of file TSGForRoadSearch.cc.

References alongMomentum, GeometricSearchDet::compatibleDetsV(), SurfaceOrientation::inner, IPfts(), LogDebug, makeSeeds_0(), GeometricSearchTracker::negPixelForwardLayers(), GeometricSearchTracker::negTecLayers(), GeometricSearchTracker::negTidLayers(), notAtIPtsos(), GeomDetEnumerators::PixelBarrel, GeometricSearchTracker::pixelBarrelLayers(), GeomDetEnumerators::PixelEndcap, position, GeometricSearchTracker::posPixelForwardLayers(), GeometricSearchTracker::posTecLayers(), GeometricSearchTracker::posTidLayers(), Propagator::propagate(), MuonServiceProxy::propagator(), pushTrajectorySeed(), mps_fire::result, GeomDetEnumerators::subDetGeom, GeomDetEnumerators::TEC, theCategory, theChi2Estimator, theGeometricSearchTracker, thePropagatorCompatibleName, thePropagatorName, theProxyService, GeomDetEnumerators::TIB, GeomDetEnumerators::TID, GeomDetEnumerators::TOB, and z.

Referenced by trackerSeeds().

                                                                                           {
  //get the state at IP
  FreeTrajectoryState cIPFTS;
  if (!IPfts(muon, cIPFTS))
    return;

  //take state at inner surface and check the first part reached
  const std::vector<const BarrelDetLayer *> &blc = theGeometricSearchTracker->pixelBarrelLayers();
  if (blc.empty()) {
    edm::LogError(theCategory) << "want to start from pixel layer, but no barrel exists. trying without pixel.";
    makeSeeds_0(muon, result);
    return;
  }

  TrajectoryStateOnSurface inner =
      theProxyService->propagator(thePropagatorName)->propagate(cIPFTS, blc.front()->surface());
  if (!inner.isValid()) {
    LogDebug(theCategory) << "inner state is not valid. no seed.";
    return;
  }

  //rescale the error
  if (!notAtIPtsos(inner))
    return;

  double z = inner.globalPosition().z();

  const std::vector<const ForwardDetLayer *> &ppxlc = theGeometricSearchTracker->posPixelForwardLayers();
  const std::vector<const ForwardDetLayer *> &npxlc = theGeometricSearchTracker->negPixelForwardLayers();
  const std::vector<const ForwardDetLayer *> &ptidc = theGeometricSearchTracker->posTidLayers();
  const std::vector<const ForwardDetLayer *> &ptecc = theGeometricSearchTracker->posTecLayers();
  const std::vector<const ForwardDetLayer *> &ntidc = theGeometricSearchTracker->negTidLayers();
  const std::vector<const ForwardDetLayer *> &ntecc = theGeometricSearchTracker->negTecLayers();

  if ((ppxlc.empty() || npxlc.empty()) && (ptidc.empty() || ptecc.empty())) {
    edm::LogError(theCategory) << "want to start from pixel layer, but no forward layer exists. trying without pixel.";
    makeSeeds_0(muon, result);
    return;
  }

  const DetLayer *inLayer = nullptr;
  std::vector<const ForwardDetLayer *>::const_iterator layerIt;

  double fz = fabs(z);

  //simple way of finding a first layer to try out
  if (fz < fabs(((z > 0) ? ppxlc : npxlc).front()->surface().position().z())) {
    inLayer = blc.front();
  } else if (fz < fabs(((z > 0) ? ppxlc : npxlc).back()->surface().position().z())) {
    layerIt = ((z > 0) ? ppxlc : npxlc).begin();
    inLayer = *layerIt;
  } else if (fz < fabs(((z > 0) ? ptidc : ntidc).front()->surface().position().z())) {
    layerIt = ((z > 0) ? ppxlc : npxlc).end() - 1;
    inLayer = *layerIt;
  } else if (fz < fabs(((z > 0) ? ptecc : ntecc).front()->surface().position().z())) {
    layerIt = ((z > 0) ? ptidc : ntidc).begin();
    inLayer = *layerIt;
  } else if (fz < fabs(((z > 0) ? ptecc : ntecc).back()->surface().position().z())) {
    layerIt = ((z > 0) ? ptecc : ntecc).begin();
    inLayer = *layerIt;
  } else {
    edm::LogWarning(theCategory) << "the state is not consistent with any tracker layer:\n" << inner;
    return;
  }

  //find out at least one compatible detector reached
  std::vector<DetLayer::DetWithState> compatible;
  compatible.reserve(10);
  inLayer->compatibleDetsV(
      inner, *theProxyService->propagator(thePropagatorCompatibleName), *theChi2Estimator, compatible);

  //if none were found. you should do something more.
  if (compatible.empty()) {
    std::vector<const ForwardDetLayer *>::const_iterator pxlEnd = (z > 0) ? ppxlc.end() : npxlc.end();
    std::vector<const ForwardDetLayer *>::const_iterator tidEnd = (z > 0) ? ptidc.end() : ntidc.end();
    std::vector<const ForwardDetLayer *>::const_iterator tecEnd = (z > 0) ? ptecc.end() : ntecc.end();
    std::vector<const ForwardDetLayer *>::const_iterator pxlBegin = (z > 0) ? ppxlc.begin() : npxlc.begin();
    std::vector<const ForwardDetLayer *>::const_iterator tidBegin = (z > 0) ? ptidc.begin() : ntidc.begin();
    std::vector<const ForwardDetLayer *>::const_iterator tecBegin = (z > 0) ? ptecc.begin() : ntecc.begin();

    //go to first disk if not already in a disk situation
    if (!dynamic_cast<const ForwardDetLayer *>(inLayer))
      layerIt = pxlBegin--;

    while (compatible.empty()) {
      switch (GeomDetEnumerators::subDetGeom[(*layerIt)->subDetector()]) {
        case GeomDetEnumerators::PixelEndcap: {
          layerIt++;
          //if end of list reached. go to the first TID
          if (layerIt == pxlEnd)
            layerIt = tidBegin;
          break;
        }
        case GeomDetEnumerators::TID: {
          layerIt++;
          //if end of list reached. go to the first TEC
          if (layerIt == tidEnd)
            layerIt = tecBegin;
          break;
        }
        case GeomDetEnumerators::TEC: {
          layerIt++;
          if (layerIt == tecEnd) {
            edm::LogWarning(theCategory) << "ran out of layers to find a seed: no seed.";
            return;
          }
          break;
        }
        case GeomDetEnumerators::PixelBarrel: {
          edm::LogError(theCategory)
              << "this should not happen... ever. Please report. GeomDetEnumerators::PixelBarrel. no seed.";
          return;
        }
        case GeomDetEnumerators::TIB: {
          edm::LogError(theCategory)
              << "this should not happen... ever. Please report. GeomDetEnumerators::TIB. no seed.";
          return;
        }
        case GeomDetEnumerators::TOB: {
          edm::LogError(theCategory)
              << "this should not happen... ever. Please report. GeomDetEnumerators::TOB. no seed.";
          return;
        }
        default: {
          edm::LogError(theCategory) << "Subdetector id is not a tracker sub-detector id. no seed.";
          return;
        }
      }  //switch

      (*layerIt)->compatibleDetsV(
          inner, *theProxyService->propagator(thePropagatorCompatibleName), *theChi2Estimator, compatible);
    }  //while
  }    //if size==0

  pushTrajectorySeed(muon, compatible, alongMomentum, result);

  return;
}

notAtIPtsos()

bool TSGForRoadSearch::notAtIPtsos ( TrajectoryStateOnSurface &  state )

private

make the adjustement away from PCA state if requested

Definition at line 103 of file TSGForRoadSearch.cc.

References MuonErrorMatrix::adjust(), LogDebug, theAdjustAtIp, theCategory, and theErrorMatrixAdjuster.

Referenced by makeSeeds_0(), makeSeeds_3(), and makeSeeds_4().

                                                                  {
  LogDebug(theCategory) << "outer state: " << state;
  if (theErrorMatrixAdjuster && !theAdjustAtIp) {
    theErrorMatrixAdjuster->adjust(state);
    LogDebug(theCategory) << "outer state after rescale: " << state;
  }
  return true;
}

pushTrajectorySeed()

void TSGForRoadSearch::pushTrajectorySeed ( const reco::Track &  muon, std::vector< DetLayer::DetWithState > &  compatible, PropagationDirection  direction, std::vector< TrajectorySeed > &  result  ) const

private

add the seed(s) to the collection of seeds

Definition at line 458 of file TSGForRoadSearch.cc.

References edm::OwnVector< T, P >::clear(), MeasurementDetWithData::fastMeasurements(), GeomDet::geographicalId(), MeasurementTrackerEvent::idToDet(), LogDebug, trajectoryStateTransform::persistentState(), MuonServiceProxy::propagator(), edm::OwnVector< T, P >::push_back(), DetId::rawId(), mps_fire::result, theCategory, theChi2Estimator, theCopyMuonRecHit, theManySeeds, theMeasurementTrackerEvent, thePropagatorCompatibleName, theProxyService, theUpdator, createJobs::tmp, and TrajectoryStateUpdator::update().

Referenced by makeSeeds_0(), makeSeeds_3(), and makeSeeds_4().

                                                                                   {
  if (compatible.empty()) {
    LogDebug(theCategory) << "pushTrajectorySeed with no compatible module. 0 seed.";
    return;
  }

  if (theManySeeds) {
    //finf out every compatible measurements
    for (std::vector<DetLayer::DetWithState>::iterator DWSit = compatible.begin(); DWSit != compatible.end(); ++DWSit) {
      bool aBareTS = false;
      const GeomDet *gd = DWSit->first;
      if (!gd) {
        edm::LogError(theCategory) << "GeomDet is not valid.";
        continue;
      }
      MeasurementDetWithData md = theMeasurementTrackerEvent->idToDet(gd->geographicalId());
      std::vector<TrajectoryMeasurement> tmp = md.fastMeasurements(
          DWSit->second, DWSit->second, *theProxyService->propagator(thePropagatorCompatibleName), *theChi2Estimator);
      //make a trajectory seed for each of them

      for (std::vector<TrajectoryMeasurement>::iterator Mit = tmp.begin(); Mit != tmp.end(); ++Mit) {
        TrajectoryStateOnSurface predState(Mit->predictedState());
        TrajectoryMeasurement::ConstRecHitPointer hit = Mit->recHit();
        TrajectorySeed::RecHitContainer rhContainer;
        if (theCopyMuonRecHit) {
          LogDebug(theCategory) << "copying (" << muon.recHitsSize() << ") muon recHits";
          //copy the muon rechit into the seed
          for (trackingRecHit_iterator trit = muon.recHitsBegin(); trit != muon.recHitsEnd(); trit++) {
            rhContainer.push_back((*trit)->clone());
          }
        }

        if (hit->isValid()) {
          TrajectoryStateOnSurface upState(theUpdator->update(predState, *hit));

          PTrajectoryStateOnDet const &PTSOD =
              trajectoryStateTransform::persistentState(upState, gd->geographicalId().rawId());
          LogDebug(theCategory) << "state used to build a trajectory seed: \n"
                                << upState << "on detector: " << gd->geographicalId().rawId();
          //add the tracking rechit
          if (theCopyMuonRecHit) {
            edm::LogError(theCategory) << "not a bare seed and muon hits are copied. dumping the muon hits.";
            rhContainer.clear();
          }
          rhContainer.push_back(hit->hit()->clone());

          result.push_back(TrajectorySeed(PTSOD, rhContainer, direction));
        } else {
          //rec hit is not valid. put a bare TrajectorySeed, only once !
          if (!aBareTS) {
            aBareTS = true;

            PTrajectoryStateOnDet const &PTSOD =
                trajectoryStateTransform::persistentState(predState, gd->geographicalId().rawId());
            LogDebug(theCategory) << "state used to build a bare trajectory seed: \n"
                                  << predState << "on detector: " << gd->geographicalId().rawId();

            result.push_back(TrajectorySeed(PTSOD, rhContainer, direction));
          }
        }
      }
    }
  } else {
    //transform it into a PTrajectoryStateOnDet

    PTrajectoryStateOnDet const &PTSOD = trajectoryStateTransform::persistentState(
        compatible.front().second, compatible.front().first->geographicalId().rawId());
    LogDebug(theCategory) << "state used to build a bare trajectory seed: \n"
                          << compatible.front().second
                          << "on detector: " << compatible.front().first->geographicalId().rawId();

    TrajectorySeed::RecHitContainer rhContainer;
    if (theCopyMuonRecHit) {
      LogDebug(theCategory) << "copying (" << muon.recHitsSize() << ") muon recHits";
      //copy the muon rechit into the seed
      for (trackingRecHit_iterator trit = muon.recHitsBegin(); trit != muon.recHitsEnd(); trit++) {
        rhContainer.push_back((*trit)->clone());
      }
    }

    //add this seed to the list and return it
    result.push_back(TrajectorySeed(PTSOD, rhContainer, direction));
  }
  return;
}

setEvent()

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

overridevirtual

set the event: update the MeasurementTracker

Reimplemented from TrackerSeedGenerator.

Definition at line 71 of file TSGForRoadSearch.cc.

References data, MuonServiceProxy::eventSetup(), edm::EventSetup::getHandle(), theGeometricSearchTracker, theGeometricSearchTrackerToken, theMeasurementTrackerEvent, theMeasurementTrackerEventToken, and theProxyService.

                                                     {
  //get the measurementtracker
  theGeometricSearchTracker = theProxyService->eventSetup().getHandle(theGeometricSearchTrackerToken);

  edm::Handle<MeasurementTrackerEvent> data;
  event.getByToken(theMeasurementTrackerEventToken, data);
  theMeasurementTrackerEvent = &*data;
}

trackerSeeds()

void TSGForRoadSearch::trackerSeeds ( const TrackCand &  , const TrackingRegion &  , const TrackerTopology *  , BTSeedCollection &    )

overridevirtual

generated seed(s) for a track. the tracking region is not used.

Reimplemented from TrackerSeedGenerator.

Definition at line 80 of file TSGForRoadSearch.cc.

References makeSeeds_0(), makeSeeds_1(), makeSeeds_2(), makeSeeds_3(), makeSeeds_4(), mps_fire::result, and theOption.

                                                                       {
  switch (theOption) {
    case 0:
      makeSeeds_0(*muonTrackCand.second, result);
      break;
    case 1:
      makeSeeds_1(*muonTrackCand.second, result);
      break;
    case 2:
      makeSeeds_2(*muonTrackCand.second, result);
      break;
    case 3:
      makeSeeds_3(*muonTrackCand.second, result);
      break;
    case 4:
      makeSeeds_4(*muonTrackCand.second, result);
      break;
  }
}

Member Data Documentation

theAdjustAtIp

bool TSGForRoadSearch::theAdjustAtIp

private

adjust the state at IP or where it is defined for the seed

Definition at line 83 of file TSGForRoadSearch.h.

Referenced by IPfts(), notAtIPtsos(), and TSGForRoadSearch().

theCategory

std::string TSGForRoadSearch::theCategory

private

Definition at line 108 of file TSGForRoadSearch.h.

Referenced by IPfts(), makeSeeds_0(), makeSeeds_1(), makeSeeds_2(), makeSeeds_3(), makeSeeds_4(), notAtIPtsos(), pushTrajectorySeed(), and TSGForRoadSearch().

theChi2Estimator

Chi2MeasurementEstimator* TSGForRoadSearch::theChi2Estimator

private

Definition at line 107 of file TSGForRoadSearch.h.

Referenced by makeSeeds_0(), makeSeeds_3(), makeSeeds_4(), pushTrajectorySeed(), TSGForRoadSearch(), and ~TSGForRoadSearch().

theConfig

edm::ParameterSet TSGForRoadSearch::theConfig

private

Definition at line 90 of file TSGForRoadSearch.h.

theCopyMuonRecHit

bool TSGForRoadSearch::theCopyMuonRecHit

private

Definition at line 103 of file TSGForRoadSearch.h.

Referenced by pushTrajectorySeed(), and TSGForRoadSearch().

theErrorMatrixAdjuster

MuonErrorMatrix* TSGForRoadSearch::theErrorMatrixAdjuster

private

Definition at line 110 of file TSGForRoadSearch.h.

Referenced by IPfts(), notAtIPtsos(), TSGForRoadSearch(), and ~TSGForRoadSearch().

theGeometricSearchTracker

edm::ESHandle<GeometricSearchTracker> TSGForRoadSearch::theGeometricSearchTracker

private

Definition at line 92 of file TSGForRoadSearch.h.

Referenced by makeSeeds_0(), makeSeeds_3(), makeSeeds_4(), and setEvent().

theGeometricSearchTrackerToken

edm::ESGetToken<GeometricSearchTracker, TrackerRecoGeometryRecord> TSGForRoadSearch::theGeometricSearchTrackerToken

private

Definition at line 93 of file TSGForRoadSearch.h.

Referenced by setEvent().

theManySeeds

bool TSGForRoadSearch::theManySeeds

private

Definition at line 104 of file TSGForRoadSearch.h.

Referenced by pushTrajectorySeed(), and TSGForRoadSearch().

theMeasurementTrackerEvent

const MeasurementTrackerEvent* TSGForRoadSearch::theMeasurementTrackerEvent

private

Definition at line 97 of file TSGForRoadSearch.h.

Referenced by pushTrajectorySeed(), setEvent(), and TSGForRoadSearch().

theMeasurementTrackerEventTag

edm::InputTag TSGForRoadSearch::theMeasurementTrackerEventTag

private

Definition at line 95 of file TSGForRoadSearch.h.

Referenced by TSGForRoadSearch().

theMeasurementTrackerEventToken

edm::EDGetTokenT<MeasurementTrackerEvent> TSGForRoadSearch::theMeasurementTrackerEventToken

private

Definition at line 96 of file TSGForRoadSearch.h.

Referenced by setEvent(), and TSGForRoadSearch().

theOption

unsigned int TSGForRoadSearch::theOption

private

Definition at line 102 of file TSGForRoadSearch.h.

Referenced by trackerSeeds(), and TSGForRoadSearch().

thePropagatorCompatibleName

std::string TSGForRoadSearch::thePropagatorCompatibleName

private

Definition at line 106 of file TSGForRoadSearch.h.

Referenced by makeSeeds_0(), makeSeeds_3(), makeSeeds_4(), pushTrajectorySeed(), and TSGForRoadSearch().

thePropagatorName

std::string TSGForRoadSearch::thePropagatorName

private

Definition at line 105 of file TSGForRoadSearch.h.

Referenced by makeSeeds_0(), makeSeeds_3(), makeSeeds_4(), and TSGForRoadSearch().

theProxyService

const MuonServiceProxy* TSGForRoadSearch::theProxyService

private

Definition at line 100 of file TSGForRoadSearch.h.

Referenced by init(), IPfts(), makeSeeds_0(), makeSeeds_3(), makeSeeds_4(), pushTrajectorySeed(), and setEvent().

theUpdator

TrajectoryStateUpdator* TSGForRoadSearch::theUpdator

private

Definition at line 99 of file TSGForRoadSearch.h.

Referenced by pushTrajectorySeed(), TSGForRoadSearch(), and ~TSGForRoadSearch().