CosmicMuonSeedGenerator Class Reference

#include <CosmicMuonSeedGenerator.h>

Inheritance diagram for CosmicMuonSeedGenerator:

Classes
struct	DecreasingGlobalY
struct	MuonRecHitPair

Public Member Functions
	CosmicMuonSeedGenerator (const edm::ParameterSet &)
	Constructor. More...
void	produce (edm::Event &, const edm::EventSetup &) override
	reconstruct muon's seeds More...
	~CosmicMuonSeedGenerator () override
	Destructor. More...
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Types
typedef std::vector< MuonRecHitPair >	MuonRecHitPairVector

Private Member Functions
bool	areCorrelated (const MuonTransientTrackingRecHit::MuonRecHitPointer &, const MuonTransientTrackingRecHit::MuonRecHitPointer &) const
	check if two rechits are correlated More...
bool	checkQuality (const MuonTransientTrackingRecHit::MuonRecHitPointer &) const
	determine if a MuonTransientTrackingRecHit is qualified to build seed More...
std::vector< TrajectorySeed >	createSeed (const MuonRecHitPair &, const edm::EventSetup &) const
	create TrajectorySeed from MuonRecHitPair More...
std::vector< TrajectorySeed >	createSeed (const MuonTransientTrackingRecHit::MuonRecHitPointer &, const edm::EventSetup &) const
	create TrajectorySeed from MuonTransientTrackingRecHit More...
void	createSeeds (TrajectorySeedCollection &results, const CosmicMuonSeedGenerator::MuonRecHitPairVector &hits, const edm::EventSetup &eSetup) const
void	createSeeds (TrajectorySeedCollection &results, const MuonTransientTrackingRecHit::MuonRecHitContainer &hits, const edm::EventSetup &eSetup) const
	generate TrajectorySeeds and put them into results More...
bool	leftIsBetter (const MuonTransientTrackingRecHit::MuonRecHitPointer &, const MuonTransientTrackingRecHit::MuonRecHitPointer &) const
	compare quality of two rechits More...
std::vector< MuonRecHitPair >	makeSegPairs (const MuonTransientTrackingRecHit::MuonRecHitContainer &, const MuonTransientTrackingRecHit::MuonRecHitContainer &, std::string) const
MuonTransientTrackingRecHit::MuonRecHitContainer	selectSegments (const MuonTransientTrackingRecHit::MuonRecHitContainer &) const
	select seed candidates from Segments in Event More...
TrajectorySeed	tsosToSeed (const TrajectoryStateOnSurface &, uint32_t) const
TrajectorySeed	tsosToSeed (const TrajectoryStateOnSurface &, uint32_t, edm::OwnVector< TrackingRecHit > &) const

Private Attributes
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magFieldToken
edm::ESGetToken< MuonDetLayerGeometry, MuonRecoGeometryRecord >	muonLayersToken
MuonDetLayerMeasurements *	muonMeasurements
edm::InputTag	theCSCRecSegmentLabel
	the name of the CSC rec hits collection More...
edm::InputTag	theDTRecSegmentLabel
	the name of the DT rec hits collection More...
bool	theEnableCSCFlag
	enable CSCSegment Flag More...
bool	theEnableDTFlag
	enable DT Segment Flag More...
edm::ESHandle< MagneticField >	theField
bool	theForcePointDownFlag
double	theMaxCSCChi2
double	theMaxDTChi2
	the maximum chi2 required for dt and csc rechits More...
unsigned int	theMaxSeeds
	the maximum number of Seeds More...
edm::ESHandle< MuonDetLayerGeometry >	theMuonLayers
std::map< std::string, float >	theParameters

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

SeedGenerator for Cosmic Muon

Author

Chang Liu - Purdue University

Definition at line 35 of file CosmicMuonSeedGenerator.h.

Member Typedef Documentation

MuonRecHitPairVector

typedef std::vector<MuonRecHitPair> CosmicMuonSeedGenerator::MuonRecHitPairVector

private

Definition at line 60 of file CosmicMuonSeedGenerator.h.

Constructor & Destructor Documentation

CosmicMuonSeedGenerator()

CosmicMuonSeedGenerator::CosmicMuonSeedGenerator

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 39 of file CosmicMuonSeedGenerator.cc.

                                                                            {
  produces<TrajectorySeedCollection>();
 
  // enable the DT chamber
  theEnableDTFlag = pset.getParameter<bool>("EnableDTMeasurement");
  // enable the CSC chamber
  theEnableCSCFlag = pset.getParameter<bool>("EnableCSCMeasurement");
 
  theDTRecSegmentLabel = pset.getParameter<InputTag>("DTRecSegmentLabel");
 
  theCSCRecSegmentLabel = pset.getParameter<InputTag>("CSCRecSegmentLabel");
 
  // the maximum number of TrajectorySeed
  theMaxSeeds = pset.getParameter<int>("MaxSeeds");
 
  theMaxDTChi2 = pset.getParameter<double>("MaxDTChi2");
  theMaxCSCChi2 = pset.getParameter<double>("MaxCSCChi2");
 
  theForcePointDownFlag = pset.existsAs<bool>("ForcePointDown") ? pset.getParameter<bool>("ForcePointDown") : true;
 
  // pre-determined parameters for seed pt calculation ( pt * dphi )
  theParameters["topmb41"] = 0.87;
  theParameters["bottommb41"] = 1.2;
  theParameters["topmb42"] = 0.67;
  theParameters["bottommb42"] = 0.98;
  theParameters["topmb43"] = 0.34;
  theParameters["bottommb43"] = 0.58;
  theParameters["topmb31"] = 0.54;
  theParameters["bottommb31"] = 0.77;
  theParameters["topmb32"] = 0.35;
  theParameters["bottommb32"] = 0.55;
  theParameters["topmb21"] = 0.21;
  theParameters["bottommb21"] = 0.31;
 
  edm::ConsumesCollector iC = consumesCollector();
  muonMeasurements = new MuonDetLayerMeasurements(theDTRecSegmentLabel,
                                                  theCSCRecSegmentLabel,
 
                                                  InputTag(),
                                                  InputTag(),
                                                  InputTag(),
                                                  iC,
                                                  theEnableDTFlag,
                                                  theEnableCSCFlag,
                                                  false,
                                                  false,
                                                  false);
  muonLayersToken = esConsumes<MuonDetLayerGeometry, MuonRecoGeometryRecord>();
  magFieldToken = esConsumes<MagneticField, IdealMagneticFieldRecord>();
}

References HLT_FULL_cff::InputTag, and muonDTDigis_cfi::pset.

~CosmicMuonSeedGenerator()

CosmicMuonSeedGenerator::~CosmicMuonSeedGenerator ( )

override

Destructor.

Definition at line 91 of file CosmicMuonSeedGenerator.cc.

                                                  {
  if (muonMeasurements)
    delete muonMeasurements;
}

Member Function Documentation

areCorrelated()

bool CosmicMuonSeedGenerator::areCorrelated ( const MuonTransientTrackingRecHit::MuonRecHitPointer &  lhs, const MuonTransientTrackingRecHit::MuonRecHitPointer &  rhs  ) const

private

check if two rechits are correlated

Definition at line 369 of file CosmicMuonSeedGenerator.cc.

                                                                                                            {
  bool result = false;
 
  GlobalVector dir1 = lhs->globalDirection();
  GlobalPoint pos1 = lhs->globalPosition();
  GlobalVector dir2 = rhs->globalDirection();
  GlobalPoint pos2 = rhs->globalPosition();
 
  GlobalVector dis = pos2 - pos1;
 
  if ((deltaR<double>(dir1.eta(), dir1.phi(), dir2.eta(), dir2.phi()) < 0.1 ||
       deltaR<double>(dir1.eta(), dir1.phi(), -dir2.eta(), -dir2.phi()) < 0.1) &&
      dis.mag() < 5.0)
    result = true;
 
  if ((deltaR<double>(dir1.eta(), dir1.phi(), dir2.eta(), dir2.phi()) < 0.1 ||
       deltaR<double>(dir1.eta(), dir1.phi(), -dir2.eta(), -dir2.phi()) < 0.1) &&
      (deltaR<double>(dir1.eta(), dir1.phi(), dis.eta(), dis.phi()) < 0.1 ||
       deltaR<double>(dir2.eta(), dir2.phi(), dis.eta(), dis.phi()) < 0.1))
    result = true;
 
  if (fabs(dir1.eta()) > 4.0 || fabs(dir2.eta()) > 4.0) {
    if ((fabs(dir1.theta() - dir2.theta()) < 0.07 || fabs(dir1.theta() + dir2.theta()) > 3.07) &&
        (fabs(dir1.theta() - dis.theta()) < 0.07 || fabs(dir1.theta() - dis.theta()) < 0.07 ||
         fabs(dir1.theta() + dis.theta()) > 3.07 || fabs(dir1.theta() + dis.theta()) > 3.07))
 
      result = true;
  }
 
  return result;
}

References L1TDiffHarvesting_cfi::dir1, L1TDiffHarvesting_cfi::dir2, PV3DBase< T, PVType, FrameType >::eta(), PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::phi(), mps_fire::result, and PV3DBase< T, PVType, FrameType >::theta().

checkQuality()

bool CosmicMuonSeedGenerator::checkQuality ( const MuonTransientTrackingRecHit::MuonRecHitPointer &  hit ) const

private

determine if a MuonTransientTrackingRecHit is qualified to build seed

Definition at line 196 of file CosmicMuonSeedGenerator.cc.

                                                                             {
  const std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  // only use 4D segments
  if (!hit->isValid())
    return false;
 
  if (hit->dimension() < 4) {
    LogTrace(category) << "dim < 4";
    return false;
  }
 
  if (hit->isDT() && (hit->chi2() > theMaxDTChi2)) {
    LogTrace(category) << "DT chi2 too large";
    return false;
  } else if (hit->isCSC() && (hit->chi2() > theMaxCSCChi2)) {
    LogTrace(category) << "CSC chi2 too large";
    return false;
  }
  return true;
}

References taus_updatedMVAIds_cff::category, LogTrace, and AlCaHLTBitMon_QueryRunRegistry::string.

createSeed() [1/2]

std::vector< TrajectorySeed > CosmicMuonSeedGenerator::createSeed ( const MuonRecHitPair &  hitpair, const edm::EventSetup &  eSetup  ) const

private

create TrajectorySeed from MuonRecHitPair

Definition at line 446 of file CosmicMuonSeedGenerator.cc.

                                                                                                   {
  std::vector<TrajectorySeed> result;
 
  const std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  MuonPatternRecoDumper dumper;
 
  float dphi = deltaPhi((hitpair.first)->globalDirection().barePhi(), (hitpair.second)->globalDirection().barePhi());
 
  LogTrace(category) << "hitpair.type " << hitpair.type;
 
  map<string, float>::const_iterator iterPar = theParameters.find(hitpair.type);
  if (iterPar == theParameters.end()) {
    return result;
  }
 
  // set the pt and charge by dphi
  int charge = (dphi > 0) ? -1 : 1;
 
  double pt = 999.0;
  float paraC = (iterPar->second);
 
  if (fabs(dphi) > 1e-5) {
    pt = paraC / fabs(dphi);
  }
 
  if (pt < 10.0) {
    return result;
  }  //still use the old strategy for low pt
 
  AlgebraicVector t(4);
  AlgebraicSymMatrix mat(5, 0);
 
  MuonTransientTrackingRecHit::MuonRecHitPointer hit = hitpair.first;
  if (hit->dimension() < (hitpair.second)->dimension())
    hit = hitpair.second;
 
  // Fill the LocalTrajectoryParameters
  LocalPoint segPos = hit->localPosition();
 
  GlobalVector polar(GlobalVector::Spherical(hit->globalDirection().theta(), hit->globalDirection().phi(), 1.));
  // Force all track downward for cosmic, not beam-halo
  if (theForcePointDownFlag) {
    if (hit->geographicalId().subdetId() == MuonSubdetId::DT && fabs(hit->globalDirection().eta()) < 4.0 &&
        hit->globalDirection().phi() > 0)
      polar = -polar;
 
    if (hit->geographicalId().subdetId() == MuonSubdetId::CSC && fabs(hit->globalDirection().eta()) > 2.3)
      polar = -polar;
  }
 
  polar *= fabs(pt) / polar.perp();
 
  LocalVector segDir = hit->det()->toLocal(polar);
 
  LocalTrajectoryParameters param(segPos, segDir, charge);
 
  mat = hit->parametersError().similarityT(hit->projectionMatrix());
 
  float p_err = 0.004 / paraC;
  if (pt < 10.01)
    p_err = 0.1;
  mat[0][0] = p_err;
 
  LocalTrajectoryError error(asSMatrix<5>(mat));
 
  // Create the TrajectoryStateOnSurface
  TrajectoryStateOnSurface tsos(param, error, hit->det()->surface(), &*theField);
 
  LogTrace(category) << "Trajectory State on Surface of Seed";
  LogTrace(category) << "mom: " << tsos.globalMomentum() << " phi: " << tsos.globalMomentum().phi();
  LogTrace(category) << "pos: " << tsos.globalPosition();
  LogTrace(category) << "The RecSegment relies on: ";
  LogTrace(category) << dumper.dumpMuonId(hit->geographicalId());
 
  edm::OwnVector<TrackingRecHit> container;
  container.push_back(hitpair.first->hit()->clone());
  container.push_back(hitpair.second->hit()->clone());
 
  result.push_back(tsosToSeed(tsos, hit->geographicalId().rawId(), container));
 
  return result;
}

References taus_updatedMVAIds_cff::category, ALCARECOTkAlJpsiMuMu_cff::charge, MuonSubdetId::CSC, SiPixelRawToDigiRegional_cfi::deltaPhi, pat::helper::ParametrizationHelper::dimension(), MuonSubdetId::DT, MuonPatternRecoDumper::dumpMuonId(), MillePedeFileConverter_cfg::e, relativeConstraints::error, CosmicMuonSeedGenerator::MuonRecHitPair::first, LogTrace, PV3DBase< T, PVType, FrameType >::perp(), DiDispStaMuonMonitor_cfi::pt, edm::OwnVector< T, P >::push_back(), mps_fire::result, CosmicMuonSeedGenerator::MuonRecHitPair::second, AlCaHLTBitMon_QueryRunRegistry::string, submitPVValidationJobs::t, and CosmicMuonSeedGenerator::MuonRecHitPair::type.

createSeed() [2/2]

std::vector< TrajectorySeed > CosmicMuonSeedGenerator::createSeed ( const MuonTransientTrackingRecHit::MuonRecHitPointer &  hit, const edm::EventSetup &  eSetup  ) const

private

create TrajectorySeed from MuonTransientTrackingRecHit

Definition at line 305 of file CosmicMuonSeedGenerator.cc.

                                                                                                   {
  std::vector<TrajectorySeed> result;
 
  const std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  MuonPatternRecoDumper dumper;
 
  // set the pt by hand
  double pt = 10.0;
 
  // AlgebraicVector4 t;
  AlgebraicSymMatrix mat(5, 0);
 
  // Fill the LocalTrajectoryParameters
  LocalPoint segPos = hit->localPosition();
 
  GlobalVector polar(GlobalVector::Spherical(hit->globalDirection().theta(), hit->globalDirection().phi(), 1.));
  // Force all track downward for cosmic, not beam-halo
  if (theForcePointDownFlag) {
    if (hit->geographicalId().subdetId() == MuonSubdetId::DT && fabs(hit->globalDirection().eta()) < 4.0 &&
        hit->globalDirection().phi() > 0)
      polar = -polar;
 
    if (hit->geographicalId().subdetId() == MuonSubdetId::CSC && fabs(hit->globalDirection().eta()) > 2.3)
      polar = -polar;
  }
 
  polar *= fabs(pt) / polar.perp();
 
  LocalVector segDir = hit->det()->toLocal(polar);
 
  int charge = 1;
  LocalTrajectoryParameters param(segPos, segDir, charge);
 
  charge = -1;
  LocalTrajectoryParameters param2(segPos, segDir, charge);
 
  mat = hit->parametersError().similarityT(hit->projectionMatrix());
 
  float p_err = 0.2;
  mat[0][0] = p_err;
 
  LocalTrajectoryError error(asSMatrix<5>(mat));
 
  // Create the TrajectoryStateOnSurface
  TrajectoryStateOnSurface tsos(param, error, hit->det()->surface(), &*theField);
  TrajectoryStateOnSurface tsos2(param2, error, hit->det()->surface(), &*theField);
 
  LogTrace(category) << "Trajectory State on Surface of Seed";
  LogTrace(category) << "mom: " << tsos.globalMomentum() << " phi: " << tsos.globalMomentum().phi();
  LogTrace(category) << "pos: " << tsos.globalPosition();
  LogTrace(category) << "The RecSegment relies on: ";
  LogTrace(category) << dumper.dumpMuonId(hit->geographicalId());
 
  edm::OwnVector<TrackingRecHit> container;
  container.push_back(hit->hit()->clone());
 
  result.push_back(tsosToSeed(tsos, hit->geographicalId().rawId(), container));
  result.push_back(tsosToSeed(tsos2, hit->geographicalId().rawId(), container));
 
  return result;
}

References taus_updatedMVAIds_cff::category, ALCARECOTkAlJpsiMuMu_cff::charge, MuonSubdetId::CSC, MuonSubdetId::DT, MuonPatternRecoDumper::dumpMuonId(), relativeConstraints::error, LogTrace, PV3DBase< T, PVType, FrameType >::perp(), DiDispStaMuonMonitor_cfi::pt, edm::OwnVector< T, P >::push_back(), mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

createSeeds() [1/2]

void CosmicMuonSeedGenerator::createSeeds ( TrajectorySeedCollection &  results, const CosmicMuonSeedGenerator::MuonRecHitPairVector &  hits, const edm::EventSetup &  eSetup  ) const

private

Definition at line 286 of file CosmicMuonSeedGenerator.cc.

                                                                             {
  const std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  if (hitpairs.empty() || results.size() >= theMaxSeeds)
    return;
  for (CosmicMuonSeedGenerator::MuonRecHitPairVector::const_iterator ihitpair = hitpairs.begin();
       ihitpair != hitpairs.end();
       ihitpair++) {
    const std::vector<TrajectorySeed>& sds = createSeed((*ihitpair), eSetup);
    LogTrace(category) << "created seeds from rechit " << sds.size();
    results.insert(results.end(), sds.begin(), sds.end());
    if (results.size() >= theMaxSeeds)
      break;
  }
  return;
}

References taus_updatedMVAIds_cff::category, LogTrace, bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

createSeeds() [2/2]

void CosmicMuonSeedGenerator::createSeeds ( TrajectorySeedCollection &  results, const MuonTransientTrackingRecHit::MuonRecHitContainer &  hits, const edm::EventSetup &  eSetup  ) const

private

generate TrajectorySeeds and put them into results

Definition at line 269 of file CosmicMuonSeedGenerator.cc.

                                                                             {
  const std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  if (hits.empty() || results.size() >= theMaxSeeds)
    return;
  for (MuonRecHitContainer::const_iterator ihit = hits.begin(); ihit != hits.end(); ihit++) {
    const std::vector<TrajectorySeed>& sds = createSeed((*ihit), eSetup);
    LogTrace(category) << "created seeds from rechit " << sds.size();
    results.insert(results.end(), sds.begin(), sds.end());
    if (results.size() >= theMaxSeeds)
      break;
  }
  return;
}

References taus_updatedMVAIds_cff::category, hfClusterShapes_cfi::hits, LogTrace, bookConverter::results, and AlCaHLTBitMon_QueryRunRegistry::string.

leftIsBetter()

bool CosmicMuonSeedGenerator::leftIsBetter ( const MuonTransientTrackingRecHit::MuonRecHitPointer &  lhs, const MuonTransientTrackingRecHit::MuonRecHitPointer &  rhs  ) const

private

compare quality of two rechits

Definition at line 401 of file CosmicMuonSeedGenerator.cc.

                                                                                                          {
  if ((lhs->degreesOfFreedom() > rhs->degreesOfFreedom()) ||
      ((lhs->degreesOfFreedom() == rhs->degreesOfFreedom()) && (lhs)->chi2() < (rhs)->chi2()))
    return true;
  else
    return false;
}

References hltPixelTracks_cff::chi2.

makeSegPairs()

CosmicMuonSeedGenerator::MuonRecHitPairVector CosmicMuonSeedGenerator::makeSegPairs ( const MuonTransientTrackingRecHit::MuonRecHitContainer &  hits1, const MuonTransientTrackingRecHit::MuonRecHitContainer &  hits2, std::string  tag  ) const

private

Definition at line 410 of file CosmicMuonSeedGenerator.cc.

                         {
  MuonRecHitPairVector result;
  const std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  if (hits1.empty() || hits2.empty())
    return result;
 
  for (MuonRecHitContainer::const_iterator ihit1 = hits1.begin(); ihit1 != hits1.end(); ihit1++) {
    if (!checkQuality(*ihit1))
      continue;
 
    for (MuonRecHitContainer::const_iterator ihit2 = hits2.begin(); ihit2 != hits2.end(); ihit2++) {
      if (!checkQuality(*ihit2))
        continue;
 
      float dphi = deltaPhi((*ihit1)->globalPosition().barePhi(), (*ihit2)->globalPosition().barePhi());
      if (dphi < 0.5) {
        if ((*ihit1)->globalPosition().y() > 0.0 && ((*ihit1)->globalPosition().y() > (*ihit2)->globalPosition().y())) {
          std::string tag2 = "top" + tag;
 
          result.push_back(MuonRecHitPair(*ihit1, *ihit2, tag2));
        } else if ((*ihit1)->globalPosition().y() < 0.0 &&
                   ((*ihit1)->globalPosition().y() < (*ihit2)->globalPosition().y())) {
          std::string tag2 = "bottom" + tag;
          result.push_back(MuonRecHitPair(*ihit2, *ihit1, tag2));
        }
      }
    }
  }
 
  return result;
}

References taus_updatedMVAIds_cff::category, SiPixelRawToDigiRegional_cfi::deltaPhi, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, and GlobalPosition_Frontier_DevDB_cff::tag.

produce()

void CosmicMuonSeedGenerator::produce ( edm::Event &  event, const edm::EventSetup &  eSetup  )

override

reconstruct muon's seeds

Definition at line 97 of file CosmicMuonSeedGenerator.cc.

                                                                                  {
  theField = eSetup.getHandle(magFieldToken);
 
  auto output = std::make_unique<TrajectorySeedCollection>();
 
  TrajectorySeedCollection seeds;
 
  std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  // Muon Geometry - DT, CSC and RPC
  theMuonLayers = eSetup.getHandle(muonLayersToken);
 
  // get the DT layers
  vector<const DetLayer*> dtLayers = theMuonLayers->allDTLayers();
 
  // get the CSC layers
  vector<const DetLayer*> cscForwardLayers = theMuonLayers->forwardCSCLayers();
  vector<const DetLayer*> cscBackwardLayers = theMuonLayers->backwardCSCLayers();
 
  muonMeasurements->setEvent(event);
 
  MuonRecHitContainer allHits;
 
  vector<MuonRecHitContainer> RHMBs;
  vector<MuonRecHitContainer> RHMEFs;
  vector<MuonRecHitContainer> RHMEBs;
 
  stable_sort(allHits.begin(), allHits.end(), DecreasingGlobalY());
 
  for (vector<const DetLayer*>::reverse_iterator icsclayer = cscForwardLayers.rbegin();
       icsclayer != cscForwardLayers.rend() - 1;
       ++icsclayer) {
    MuonRecHitContainer RHMF = muonMeasurements->recHits(*icsclayer);
    allHits.insert(allHits.end(), RHMF.begin(), RHMF.end());
  }
 
  for (vector<const DetLayer*>::reverse_iterator icsclayer = cscBackwardLayers.rbegin();
       icsclayer != cscBackwardLayers.rend() - 1;
       ++icsclayer) {
    MuonRecHitContainer RHMF = muonMeasurements->recHits(*icsclayer);
    allHits.insert(allHits.end(), RHMF.begin(), RHMF.end());
  }
 
  for (vector<const DetLayer*>::reverse_iterator idtlayer = dtLayers.rbegin(); idtlayer != dtLayers.rend();
       ++idtlayer) {
    MuonRecHitContainer RHMB = muonMeasurements->recHits(*idtlayer);
    RHMBs.push_back(RHMB);
 
    if (idtlayer != dtLayers.rbegin())
      allHits.insert(allHits.end(), RHMB.begin(), RHMB.end());
  }
 
  //  stable_sort(allHits.begin(),allHits.end(),DecreasingGlobalY());
 
  LogTrace(category) << "all RecHits: " << allHits.size();
 
  //  CosmicMuonSeedGenerator::MuonRecHitPairVector mb41 = makeSegPairs(RHMBs[0], RHMBs[3], "mb41");
  //  createSeeds(seeds, mb41, eSetup);
 
  //  CosmicMuonSeedGenerator::MuonRecHitPairVector mb43 = makeSegPairs(RHMBs[0],RHMBs[1], "mb43");
  //  createSeeds(seeds, mb43, eSetup);
 
  CosmicMuonSeedGenerator::MuonRecHitPairVector mb42 = makeSegPairs(RHMBs[0], RHMBs[2], "mb42");
  createSeeds(seeds, mb42, eSetup);
 
  //  CosmicMuonSeedGenerator::MuonRecHitPairVector mb32 = makeSegPairs(RHMBs[1], RHMBs[2], "mb32");
  //  createSeeds(seeds, mb32, eSetup);
 
  CosmicMuonSeedGenerator::MuonRecHitPairVector mb31 = makeSegPairs(RHMBs[1], RHMBs[3], "mb31");
  createSeeds(seeds, mb31, eSetup);
 
  //  CosmicMuonSeedGenerator::MuonRecHitPairVector mb21 = makeSegPairs(RHMBs[2], RHMBs[3], "mb21");
  //  createSeeds(seeds, mb21, eSetup);
 
  if (!allHits.empty()) {
    MuonRecHitContainer goodhits = selectSegments(allHits);
    LogTrace(category) << "good RecHits: " << goodhits.size();
 
    if (goodhits.empty()) {
      LogTrace(category) << "No qualified Segments in Event! ";
      LogTrace(category) << "Use 2D RecHit";
 
      createSeeds(seeds, allHits, eSetup);
 
    } else {
      createSeeds(seeds, goodhits, eSetup);
    }
  }
 
  LogTrace(category) << "Seeds built: " << seeds.size();
 
  for (std::vector<TrajectorySeed>::iterator seed = seeds.begin(); seed != seeds.end(); ++seed) {
    output->push_back(*seed);
  }
 
  event.put(std::move(output));
  seeds.clear();
}

References taus_updatedMVAIds_cff::category, edm::EventSetup::getHandle(), LogTrace, eostools::move(), convertSQLitetoXML_cfg::output, fileCollector::seed, InitialStep_cff::seeds, and AlCaHLTBitMon_QueryRunRegistry::string.

selectSegments()

MuonRecHitContainer CosmicMuonSeedGenerator::selectSegments ( const MuonTransientTrackingRecHit::MuonRecHitContainer &  hits ) const

private

select seed candidates from Segments in Event

Definition at line 218 of file CosmicMuonSeedGenerator.cc.

                                                                                                 {
  MuonRecHitContainer result;
  const std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  //Only select good quality Segments
  for (MuonRecHitContainer::const_iterator hit = hits.begin(); hit != hits.end(); hit++) {
    if (checkQuality(*hit))
      result.push_back(*hit);
  }
 
  if (result.size() < 2)
    return result;
 
  MuonRecHitContainer result2;
 
  //avoid selecting Segments with similar direction
  for (MuonRecHitContainer::iterator hit = result.begin(); hit != result.end(); hit++) {
    if (*hit == nullptr)
      continue;
    if (!(*hit)->isValid())
      continue;
    bool good = true;
    //UNUSED:    GlobalVector dir1 = (*hit)->globalDirection();
    //UNUSED:    GlobalPoint pos1 = (*hit)->globalPosition();
    for (MuonRecHitContainer::iterator hit2 = hit + 1; hit2 != result.end(); hit2++) {
      if (*hit2 == nullptr)
        continue;
      if (!(*hit2)->isValid())
        continue;
 
      //compare direction and position
      //UNUSED:        GlobalVector dir2 = (*hit2)->globalDirection();
      //UNUSED:        GlobalPoint pos2 = (*hit2)->globalPosition();
      if (!areCorrelated((*hit), (*hit2)))
        continue;
 
      if (!leftIsBetter((*hit), (*hit2))) {
        good = false;
      } else
        (*hit2) = nullptr;
    }
 
    if (good)
      result2.push_back(*hit);
  }
 
  result.clear();
 
  return result2;
}

References taus_updatedMVAIds_cff::category, hfClusterShapes_cfi::hits, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

tsosToSeed() [1/2]

TrajectorySeed CosmicMuonSeedGenerator::tsosToSeed ( const TrajectoryStateOnSurface &  tsos, uint32_t  id  ) const

private

Definition at line 531 of file CosmicMuonSeedGenerator.cc.

                                                                                                          {
  edm::OwnVector<TrackingRecHit> container;
  return tsosToSeed(tsos, id, container);
}

tsosToSeed() [2/2]

TrajectorySeed CosmicMuonSeedGenerator::tsosToSeed ( const TrajectoryStateOnSurface &  tsos, uint32_t  id, edm::OwnVector< TrackingRecHit > &  container  ) const

private

Definition at line 536 of file CosmicMuonSeedGenerator.cc.

                                                                                                  {
  PTrajectoryStateOnDet const& seedTSOS = trajectoryStateTransform::persistentState(tsos, id);
  TrajectorySeed seed(seedTSOS, container, alongMomentum);
  return seed;
}

References alongMomentum, trajectoryStateTransform::persistentState(), and fileCollector::seed.

Member Data Documentation

magFieldToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> CosmicMuonSeedGenerator::magFieldToken

private

Definition at line 130 of file CosmicMuonSeedGenerator.h.

muonLayersToken

edm::ESGetToken<MuonDetLayerGeometry, MuonRecoGeometryRecord> CosmicMuonSeedGenerator::muonLayersToken

private

Definition at line 129 of file CosmicMuonSeedGenerator.h.

muonMeasurements

MuonDetLayerMeasurements* CosmicMuonSeedGenerator::muonMeasurements

private

Definition at line 134 of file CosmicMuonSeedGenerator.h.

theCSCRecSegmentLabel

edm::InputTag CosmicMuonSeedGenerator::theCSCRecSegmentLabel

private

the name of the CSC rec hits collection

Definition at line 118 of file CosmicMuonSeedGenerator.h.

theDTRecSegmentLabel

edm::InputTag CosmicMuonSeedGenerator::theDTRecSegmentLabel

private

the name of the DT rec hits collection

Definition at line 115 of file CosmicMuonSeedGenerator.h.

theEnableCSCFlag

bool CosmicMuonSeedGenerator::theEnableCSCFlag

private

enable CSCSegment Flag

Definition at line 112 of file CosmicMuonSeedGenerator.h.

theEnableDTFlag

bool CosmicMuonSeedGenerator::theEnableDTFlag

private

enable DT Segment Flag

Definition at line 109 of file CosmicMuonSeedGenerator.h.

theField

edm::ESHandle<MagneticField> CosmicMuonSeedGenerator::theField

private

Definition at line 128 of file CosmicMuonSeedGenerator.h.

theForcePointDownFlag

bool CosmicMuonSeedGenerator::theForcePointDownFlag

private

Definition at line 126 of file CosmicMuonSeedGenerator.h.

theMaxCSCChi2

double CosmicMuonSeedGenerator::theMaxCSCChi2

private

Definition at line 125 of file CosmicMuonSeedGenerator.h.

theMaxDTChi2

double CosmicMuonSeedGenerator::theMaxDTChi2

private

the maximum chi2 required for dt and csc rechits

Definition at line 124 of file CosmicMuonSeedGenerator.h.

theMaxSeeds

unsigned int CosmicMuonSeedGenerator::theMaxSeeds

private

the maximum number of Seeds

Definition at line 121 of file CosmicMuonSeedGenerator.h.

theMuonLayers

edm::ESHandle<MuonDetLayerGeometry> CosmicMuonSeedGenerator::theMuonLayers

private

Definition at line 127 of file CosmicMuonSeedGenerator.h.

theParameters

std::map<std::string, float> CosmicMuonSeedGenerator::theParameters

private

Definition at line 132 of file CosmicMuonSeedGenerator.h.