#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 MuonTransientTrackingRecHit::MuonRecHitPointer &, const edm::EventSetup &) const
	create TrajectorySeed from MuonTransientTrackingRecHit More...

std::vector< TrajectorySeed >	createSeed (const MuonRecHitPair &, const edm::EventSetup &) const
	create TrajectorySeed from MuonRecHitPair More...

void	createSeeds (TrajectorySeedCollection &results, const MuonTransientTrackingRecHit::MuonRecHitContainer &hits, const edm::EventSetup &eSetup) const
	generate TrajectorySeeds and put them into results More...

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

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

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

private

Definition at line 55 of file CosmicMuonSeedGenerator.h.

Constructor & Destructor Documentation

CosmicMuonSeedGenerator::CosmicMuonSeedGenerator ( const edm::ParameterSet & pset )

Constructor.

Definition at line 44 of file CosmicMuonSeedGenerator.cc.

References edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), and HLT_2018_cff::InputTag.

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

CosmicMuonSeedGenerator::~CosmicMuonSeedGenerator ( )

override

Destructor.

Definition at line 94 of file CosmicMuonSeedGenerator.cc.

                                                   {
   if (muonMeasurements)
     delete muonMeasurements;
 }

Member Function Documentation

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.

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

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

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.

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

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

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.

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.

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

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.

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.

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

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.

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

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

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

private

Definition at line 289 of file CosmicMuonSeedGenerator.cc.

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

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

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.

References hltPixelTracks_cff::chi2.

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

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.

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

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

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

override

reconstruct muon's seeds

Definition at line 100 of file CosmicMuonSeedGenerator.cc.

References taus_updatedMVAIds_cff::category, edm::EventSetup::get(), LogTrace, eostools::move(), convertSQLitetoXML_cfg::output, SurveyInfoScenario_cff::seed, and AlCaHLTBitMon_QueryRunRegistry::string.

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

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

private

select seed candidates from Segments in Event

Definition at line 221 of file CosmicMuonSeedGenerator.cc.

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

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

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

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

private

Definition at line 539 of file CosmicMuonSeedGenerator.cc.

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

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

Member Data Documentation

MuonDetLayerMeasurements* CosmicMuonSeedGenerator::muonMeasurements

private

Definition at line 127 of file CosmicMuonSeedGenerator.h.

edm::InputTag CosmicMuonSeedGenerator::theCSCRecSegmentLabel

private

the name of the CSC rec hits collection

Definition at line 113 of file CosmicMuonSeedGenerator.h.

edm::InputTag CosmicMuonSeedGenerator::theDTRecSegmentLabel

private

the name of the DT rec hits collection

Definition at line 110 of file CosmicMuonSeedGenerator.h.

bool CosmicMuonSeedGenerator::theEnableCSCFlag

private

enable CSCSegment Flag

Definition at line 107 of file CosmicMuonSeedGenerator.h.

bool CosmicMuonSeedGenerator::theEnableDTFlag

private

enable DT Segment Flag

Definition at line 104 of file CosmicMuonSeedGenerator.h.

edm::ESHandle<MagneticField> CosmicMuonSeedGenerator::theField

private

Definition at line 123 of file CosmicMuonSeedGenerator.h.

bool CosmicMuonSeedGenerator::theForcePointDownFlag

private

Definition at line 121 of file CosmicMuonSeedGenerator.h.

double CosmicMuonSeedGenerator::theMaxCSCChi2

private

Definition at line 120 of file CosmicMuonSeedGenerator.h.

double CosmicMuonSeedGenerator::theMaxDTChi2

private

the maximum chi2 required for dt and csc rechits

Definition at line 119 of file CosmicMuonSeedGenerator.h.

unsigned int CosmicMuonSeedGenerator::theMaxSeeds

private

the maximum number of Seeds

Definition at line 116 of file CosmicMuonSeedGenerator.h.

edm::ESHandle<MuonDetLayerGeometry> CosmicMuonSeedGenerator::theMuonLayers

private

Definition at line 122 of file CosmicMuonSeedGenerator.h.

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

private

Definition at line 125 of file CosmicMuonSeedGenerator.h.

Classes

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation