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	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () 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
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 30 of file CosmicMuonSeedGenerator.h.

Member Typedef Documentation

MuonRecHitPairVector

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

private

Definition at line 55 of file CosmicMuonSeedGenerator.h.

Constructor & Destructor Documentation

CosmicMuonSeedGenerator()

CosmicMuonSeedGenerator::CosmicMuonSeedGenerator

(

const edm::ParameterSet &

pset

)

Constructor.

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

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

~CosmicMuonSeedGenerator()

CosmicMuonSeedGenerator::~CosmicMuonSeedGenerator ( )

override

Destructor.

Definition at line 94 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 372 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 199 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 449 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, OrderedSet::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 308 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 289 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 272 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 404 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 413 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 100 of file CosmicMuonSeedGenerator.cc.

                                                                                  {
  eSetup.get<IdealMagneticFieldRecord>().get(theField);
 
  auto output = std::make_unique<TrajectorySeedCollection>();
 
  TrajectorySeedCollection seeds;
 
  std::string category = "Muon|RecoMuon|CosmicMuonSeedGenerator";
 
  // Muon Geometry - DT, CSC and RPC
  eSetup.get<MuonRecoGeometryRecord>().get(theMuonLayers);
 
  // 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::get(), edm::get(), LogTrace, eostools::move(), convertSQLitetoXML_cfg::output, SurveyInfoScenario_cff::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 221 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 534 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 539 of file CosmicMuonSeedGenerator.cc.

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

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

Member Data Documentation

muonMeasurements

MuonDetLayerMeasurements* CosmicMuonSeedGenerator::muonMeasurements

private

Definition at line 127 of file CosmicMuonSeedGenerator.h.

theCSCRecSegmentLabel

edm::InputTag CosmicMuonSeedGenerator::theCSCRecSegmentLabel

private

the name of the CSC rec hits collection

Definition at line 113 of file CosmicMuonSeedGenerator.h.

theDTRecSegmentLabel

edm::InputTag CosmicMuonSeedGenerator::theDTRecSegmentLabel

private

the name of the DT rec hits collection

Definition at line 110 of file CosmicMuonSeedGenerator.h.

theEnableCSCFlag

bool CosmicMuonSeedGenerator::theEnableCSCFlag

private

enable CSCSegment Flag

Definition at line 107 of file CosmicMuonSeedGenerator.h.

theEnableDTFlag

bool CosmicMuonSeedGenerator::theEnableDTFlag

private

enable DT Segment Flag

Definition at line 104 of file CosmicMuonSeedGenerator.h.

theField

edm::ESHandle<MagneticField> CosmicMuonSeedGenerator::theField

private

Definition at line 123 of file CosmicMuonSeedGenerator.h.

theForcePointDownFlag

bool CosmicMuonSeedGenerator::theForcePointDownFlag

private

Definition at line 121 of file CosmicMuonSeedGenerator.h.

theMaxCSCChi2

double CosmicMuonSeedGenerator::theMaxCSCChi2

private

Definition at line 120 of file CosmicMuonSeedGenerator.h.

theMaxDTChi2

double CosmicMuonSeedGenerator::theMaxDTChi2

private

the maximum chi2 required for dt and csc rechits

Definition at line 119 of file CosmicMuonSeedGenerator.h.

theMaxSeeds

unsigned int CosmicMuonSeedGenerator::theMaxSeeds

private

the maximum number of Seeds

Definition at line 116 of file CosmicMuonSeedGenerator.h.

theMuonLayers

edm::ESHandle<MuonDetLayerGeometry> CosmicMuonSeedGenerator::theMuonLayers

private

Definition at line 122 of file CosmicMuonSeedGenerator.h.

theParameters

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

private

Definition at line 125 of file CosmicMuonSeedGenerator.h.