#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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
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<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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.

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

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

◆ ~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.

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

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

References validateAlignments::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;
 }

◆ createSeed() [1/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.

References validateAlignments::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;
 }

◆ createSeed() [2/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.

References validateAlignments::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.

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

◆ createSeeds() [1/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.

References validateAlignments::category, hfClusterShapes_cfi::hits, LogTrace, mysort::results, 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;
 }

◆ createSeeds() [2/2]

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

private

Definition at line 286 of file CosmicMuonSeedGenerator.cc.

References validateAlignments::category, LogTrace, mysort::results, 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;
 }

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

References nano_mu_local_reco_cff::chi2.

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

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

References validateAlignments::category, SiPixelRawToDigiRegional_cfi::deltaPhi, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, and makeGlobalPositionRcd_cfg::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;
 }

◆ produce()

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

override

reconstruct muon's seeds

Definition at line 97 of file CosmicMuonSeedGenerator.cc.

References validateAlignments::category, edm::EventSetup::getHandle(), LogTrace, eostools::move(), fileCollector::seed, HLT_2024v10_cff::seeds, and AlCaHLTBitMon_QueryRunRegistry::string.

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

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

References validateAlignments::category, caHitNtupletGeneratorKernels::good, hfClusterShapes_cfi::hits, 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;
 }

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

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

                                                                                                   {
   PTrajectoryStateOnDet const& seedTSOS = trajectoryStateTransform::persistentState(tsos, id);
   TrajectorySeed seed(seedTSOS, container, alongMomentum);
   return 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.

Classes

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ MuonRecHitPairVector

Constructor & Destructor Documentation

◆ CosmicMuonSeedGenerator()

◆ ~CosmicMuonSeedGenerator()

Member Function Documentation

◆ areCorrelated()

◆ checkQuality()

◆ createSeed() [1/2]

◆ createSeed() [2/2]

◆ createSeeds() [1/2]

◆ createSeeds() [2/2]

◆ leftIsBetter()

◆ makeSegPairs()

◆ produce()

◆ selectSegments()

◆ tsosToSeed() [1/2]

◆ tsosToSeed() [2/2]

Member Data Documentation

◆ magFieldToken

◆ muonLayersToken

◆ muonMeasurements

◆ theCSCRecSegmentLabel

◆ theDTRecSegmentLabel

◆ theEnableCSCFlag

◆ theEnableDTFlag

◆ theField

◆ theForcePointDownFlag

◆ theMaxCSCChi2

◆ theMaxDTChi2

◆ theMaxSeeds

◆ theMuonLayers

◆ theParameters