#include <L2MuonSeedGenerator.h>

Inheritance diagram for L2MuonSeedGenerator:

Public Member Functions
	L2MuonSeedGenerator (const edm::ParameterSet &)
	Constructor. More...

void	produce (edm::Event &, const edm::EventSetup &) override

	~L2MuonSeedGenerator () 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 Member Functions
const TrajectorySeed *	associateOfflineSeedToL1 (edm::Handle< edm::View< TrajectorySeed > > &, std::vector< int > &, TrajectoryStateOnSurface &)

Private Attributes
edm::EDGetTokenT < L1MuGMTReadoutCollection >	gmtToken_

edm::EDGetTokenT < l1extra::L1MuonParticleCollection >	muCollToken_

edm::EDGetTokenT< edm::View < TrajectorySeed > >	offlineSeedToken_

MeasurementEstimator *	theEstimator

edm::InputTag	theL1GMTReadoutCollection

const double	theL1MaxEta

const double	theL1MinPt

const unsigned	theL1MinQuality

edm::InputTag	theOfflineSeedLabel

std::string	thePropagatorName

MuonServiceProxy *	theService
	the event setup proxy, it takes care the services update More...

edm::InputTag	theSource

const bool	useOfflineSeed

const bool	useUnassociatedL1

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

L2 muon seed generator: Transform the L1 informations in seeds for the L2 muon reconstruction

Author: A.Everett, R.Bellan, J. Alcaraz

ORCA's author: N. Neumeister

L2 muon seed generator: Transform the L1 informations in seeds for the L2 muon reconstruction

Author: A.Everett, R.Bellan

ORCA's author: N. Neumeister

Definition at line 54 of file L2MuonSeedGenerator.h.

Constructor & Destructor Documentation

L2MuonSeedGenerator::L2MuonSeedGenerator ( const edm::ParameterSet & iConfig )

explicit

Constructor.

Definition at line 44 of file L2MuonSeedGenerator.cc.

References Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gmtToken_, muCollToken_, MuonServiceProxy_cff::MuonServiceProxy, offlineSeedToken_, theEstimator, theL1GMTReadoutCollection, theOfflineSeedLabel, theService, theSource, and useOfflineSeed.

     : theSource(iConfig.getParameter<InputTag>("InputObjects")),
       theL1GMTReadoutCollection(iConfig.getParameter<InputTag>("GMTReadoutCollection")),
       thePropagatorName(iConfig.getParameter<string>("Propagator")),
       theL1MinPt(iConfig.getParameter<double>("L1MinPt")),
       theL1MaxEta(iConfig.getParameter<double>("L1MaxEta")),
       theL1MinQuality(iConfig.getParameter<unsigned int>("L1MinQuality")),
       useOfflineSeed(iConfig.getUntrackedParameter<bool>("UseOfflineSeed", false)),
       useUnassociatedL1(iConfig.existsAs<bool>("UseUnassociatedL1") ? iConfig.getParameter<bool>("UseUnassociatedL1")
                                                                     : true) {
   gmtToken_ = consumes<L1MuGMTReadoutCollection>(theL1GMTReadoutCollection);
   muCollToken_ = consumes<L1MuonParticleCollection>(theSource);
 
   if (useOfflineSeed) {
     theOfflineSeedLabel = iConfig.getUntrackedParameter<InputTag>("OfflineSeedLabel");
     offlineSeedToken_ = consumes<edm::View<TrajectorySeed> >(theOfflineSeedLabel);
   }
 
   // service parameters
   ParameterSet serviceParameters = iConfig.getParameter<ParameterSet>("ServiceParameters");
 
   // the services
   theService = new MuonServiceProxy(serviceParameters, consumesCollector());
 
   // the estimator
   theEstimator = new Chi2MeasurementEstimator(10000.);
 
   produces<L2MuonTrajectorySeedCollection>();
 }

L2MuonSeedGenerator::~L2MuonSeedGenerator ( )

override

Destructor.

Definition at line 75 of file L2MuonSeedGenerator.cc.

References theEstimator, and theService.

                                           {
   if (theService)
     delete theService;
   if (theEstimator)
     delete theEstimator;
 }

Member Function Documentation

const TrajectorySeed * L2MuonSeedGenerator::associateOfflineSeedToL1	(	edm::Handle< edm::View< TrajectorySeed > > &	offseeds,
		std::vector< int > &	offseedMap,
		TrajectoryStateOnSurface &	newTsos
	)

private

Definition at line 325 of file L2MuonSeedGenerator.cc.

References edm::View< T >::begin(), HLT_FULL_cff::deltaR, MuonPatternRecoDumper::dumpMuonId(), PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), LogDebug, PV3DBase< T, PVType, FrameType >::phi(), Geom::Phi< T1, Range >::phi(), AlCaHLTBitMon_QueryRunRegistry::string, TrajectoryStateOnSurface::surface(), thePropagatorName, and theService.

Referenced by produce().

                                                                                                        {
   const std::string metlabel = "Muon|RecoMuon|L2MuonSeedGenerator";
   MuonPatternRecoDumper debugtmp;
 
   edm::View<TrajectorySeed>::const_iterator offseed, endOffseed = offseeds->end();
   const TrajectorySeed* selOffseed = nullptr;
   double bestDr = 99999.;
   unsigned int nOffseed(0);
   int lastOffseed(-1);
 
   for (offseed = offseeds->begin(); offseed != endOffseed; ++offseed, ++nOffseed) {
     if (offseedMap[nOffseed] != 0)
       continue;
     GlobalPoint glbPos = theService->trackingGeometry()
                              ->idToDet(offseed->startingState().detId())
                              ->surface()
                              .toGlobal(offseed->startingState().parameters().position());
     GlobalVector glbMom = theService->trackingGeometry()
                               ->idToDet(offseed->startingState().detId())
                               ->surface()
                               .toGlobal(offseed->startingState().parameters().momentum());
 
     // Preliminary check
     double preDr = deltaR(newTsos.globalPosition().eta(), newTsos.globalPosition().phi(), glbPos.eta(), glbPos.phi());
     if (preDr > 1.0)
       continue;
 
     const FreeTrajectoryState offseedFTS(
         glbPos, glbMom, offseed->startingState().parameters().charge(), &*theService->magneticField());
     TrajectoryStateOnSurface offseedTsos =
         theService->propagator(thePropagatorName)->propagate(offseedFTS, newTsos.surface());
     LogDebug(metlabel) << "Offline seed info: Det and State" << std::endl;
     LogDebug(metlabel) << debugtmp.dumpMuonId(offseed->startingState().detId()) << std::endl;
     //LogDebug(metlabel) << "(x, y, z) = (" << newTSOS.globalPosition().x() << ", "
     //                   << newTSOS.globalPosition().y() << ", " << newTSOS.globalPosition().z() << ")" << std::endl;
     LogDebug(metlabel) << "pos: (r=" << offseedFTS.position().mag() << ", phi=" << offseedFTS.position().phi()
                        << ", eta=" << offseedFTS.position().eta() << ")" << std::endl;
     LogDebug(metlabel) << "mom: (q*pt=" << offseedFTS.charge() * offseedFTS.momentum().perp()
                        << ", phi=" << offseedFTS.momentum().phi() << ", eta=" << offseedFTS.momentum().eta() << ")"
                        << std::endl
                        << std::endl;
     //LogDebug(metlabel) << debugtmp.dumpFTS(offseedFTS) << std::endl;
 
     if (offseedTsos.isValid()) {
       LogDebug(metlabel) << "Offline seed info after propagation to L1 layer:" << std::endl;
       //LogDebug(metlabel) << "(x, y, z) = (" << offseedTsos.globalPosition().x() << ", "
       //                   << offseedTsos.globalPosition().y() << ", " << offseedTsos.globalPosition().z() << ")" << std::endl;
       LogDebug(metlabel) << "pos: (r=" << offseedTsos.globalPosition().mag()
                          << ", phi=" << offseedTsos.globalPosition().phi()
                          << ", eta=" << offseedTsos.globalPosition().eta() << ")" << std::endl;
       LogDebug(metlabel) << "mom: (q*pt=" << offseedTsos.charge() * offseedTsos.globalMomentum().perp()
                          << ", phi=" << offseedTsos.globalMomentum().phi()
                          << ", eta=" << offseedTsos.globalMomentum().eta() << ")" << std::endl
                          << std::endl;
       //LogDebug(metlabel) << debugtmp.dumpTSOS(offseedTsos) << std::endl;
       double newDr = deltaR(newTsos.globalPosition().eta(),
                             newTsos.globalPosition().phi(),
                             offseedTsos.globalPosition().eta(),
                             offseedTsos.globalPosition().phi());
       LogDebug(metlabel) << "   -- DR = " << newDr << std::endl;
       if (newDr < 0.3 && newDr < bestDr) {
         LogDebug(metlabel) << "          --> OK! " << newDr << std::endl << std::endl;
         selOffseed = &*offseed;
         bestDr = newDr;
         offseedMap[nOffseed] = 1;
         if (lastOffseed > -1)
           offseedMap[lastOffseed] = 0;
         lastOffseed = nOffseed;
       } else {
         LogDebug(metlabel) << "          --> Rejected. " << newDr << std::endl << std::endl;
       }
     } else {
       LogDebug(metlabel) << "Invalid offline seed TSOS after propagation!" << std::endl << std::endl;
     }
   }
 
   return selOffseed;
 }

void L2MuonSeedGenerator::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 82 of file L2MuonSeedGenerator.cc.

References alongMomentum, associateOfflineSeedToL1(), Reference_intrackfit_cff::barrel, RecoTauCleanerPlugins::charge, TrajectoryStateOnSurface::charge(), L1MuGMTCand::charge_valid(), GeometricSearchDet::compatibleDets(), funct::cos(), CSCDetId, debug, DTChamberId, MuonPatternRecoDumper::dumpFTS(), MuonPatternRecoDumper::dumpLayer(), MuonPatternRecoDumper::dumpMuonId(), L1MuGMTCand::empty(), relativeConstraints::error, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), funct::exp(), GeomDet::geographicalId(), L1MuGMTReadoutRecord::getBrlRPCCands(), edm::Event::getByToken(), L1MuGMTReadoutRecord::getCSCCands(), L1MuGMTReadoutRecord::getDTBXCands(), L1MuGMTExtendedCand::getDTCSCIndex(), L1MuGMTReadoutRecord::getFwdRPCCands(), L1MuGMTExtendedCand::getRPCIndex(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), gmtToken_, gpuClustering::id, L1MuGMTExtendedCand::isFwd(), L1MuGMTExtendedCand::isRPC(), TrajectoryStateOnSurface::isValid(), LogDebug, LogTrace, PV3DBase< T, PVType, FrameType >::mag(), metname, eostools::move(), muCollToken_, offlineSeedToken_, convertSQLitetoXML_cfg::output, PV3DBase< T, PVType, FrameType >::perp(), trajectoryStateTransform::persistentState(), phi, PV3DBase< T, PVType, FrameType >::phi(), Geom::pi(), GeometricSearchDet::position(), DiDispStaMuonMonitor_cfi::pt, edm::OwnVector< T, P >::push_back(), edm::Event::put(), L1MuGMTCand::quality(), quality, CosmicsPD_Skims::radius, DetId::rawId(), TrajectorySeed::recHits(), funct::sin(), TrajectorySeed::startingState(), AlCaHLTBitMon_QueryRunRegistry::string, GeometricSearchDet::surface(), theEstimator, theL1MaxEta, theL1MinPt, theL1MinQuality, thePropagatorName, theService, theta(), useOfflineSeed, useUnassociatedL1, and PV3DBase< T, PVType, FrameType >::z().

                                                                                {
   const std::string metname = "Muon|RecoMuon|L2MuonSeedGenerator";
   MuonPatternRecoDumper debug;
 
   auto output = std::make_unique<L2MuonTrajectorySeedCollection>();
 
   // Muon particles and GMT readout collection
   edm::Handle<L1MuGMTReadoutCollection> gmtrc_handle;
   iEvent.getByToken(gmtToken_, gmtrc_handle);
   L1MuGMTReadoutRecord const& gmtrr = gmtrc_handle.product()->getRecord(0);
 
   edm::Handle<L1MuonParticleCollection> muColl;
   iEvent.getByToken(muCollToken_, muColl);
   LogTrace(metname) << "Number of muons " << muColl->size() << endl;
 
   edm::Handle<edm::View<TrajectorySeed> > offlineSeedHandle;
   vector<int> offlineSeedMap;
   if (useOfflineSeed) {
     iEvent.getByToken(offlineSeedToken_, offlineSeedHandle);
     LogTrace(metname) << "Number of offline seeds " << offlineSeedHandle->size() << endl;
     offlineSeedMap = vector<int>(offlineSeedHandle->size(), 0);
   }
 
   L1MuonParticleCollection::const_iterator it;
   L1MuonParticleRef::key_type l1ParticleIndex = 0;
 
   for (it = muColl->begin(); it != muColl->end(); ++it, ++l1ParticleIndex) {
     const L1MuGMTExtendedCand muonCand = (*it).gmtMuonCand();
     unsigned int quality = 0;
     bool valid_charge = false;
     ;
 
     if (muonCand.empty()) {
       LogWarning(metname) << "L2MuonSeedGenerator: WARNING, no L1MuGMTCand! " << endl;
       LogWarning(metname) << "L2MuonSeedGenerator:   this should make sense only within MC tests" << endl;
       // FIXME! Temporary to handle the MC input
       quality = 7;
       valid_charge = true;
     } else {
       quality = muonCand.quality();
       valid_charge = muonCand.charge_valid();
     }
 
     float pt = (*it).pt();
     float eta = (*it).eta();
     float theta = 2 * atan(exp(-eta));
     float phi = (*it).phi();
     int charge = (*it).charge();
     // Set charge=0 for the time being if the valid charge bit is zero
     if (!valid_charge)
       charge = 0;
     bool barrel = !(*it).isForward();
 
     // Get a better eta and charge from regional information
     // Phi has the same resolution in GMT than regionally, is not it?
     if (!(muonCand.empty())) {
       int idx = -1;
       vector<L1MuRegionalCand> rmc;
       if (!muonCand.isRPC()) {
         idx = muonCand.getDTCSCIndex();
         if (muonCand.isFwd())
           rmc = gmtrr.getCSCCands();
         else
           rmc = gmtrr.getDTBXCands();
       } else {
         idx = muonCand.getRPCIndex();
         if (muonCand.isFwd())
           rmc = gmtrr.getFwdRPCCands();
         else
           rmc = gmtrr.getBrlRPCCands();
       }
       if (idx >= 0) {
         eta = rmc[idx].etaValue();
         //phi = rmc[idx].phiValue();
         // Use this charge if the valid charge bit is zero
         if (!valid_charge)
           charge = rmc[idx].chargeValue();
       }
     }
 
     if (pt < theL1MinPt || fabs(eta) > theL1MaxEta)
       continue;
 
     LogTrace(metname) << "New L2 Muon Seed";
     LogTrace(metname) << "Pt = " << pt << " GeV/c";
     LogTrace(metname) << "eta = " << eta;
     LogTrace(metname) << "theta = " << theta << " rad";
     LogTrace(metname) << "phi = " << phi << " rad";
     LogTrace(metname) << "charge = " << charge;
     LogTrace(metname) << "In Barrel? = " << barrel;
 
     if (quality <= theL1MinQuality)
       continue;
     LogTrace(metname) << "quality = " << quality;
 
     // Update the services
     theService->update(iSetup);
 
     const DetLayer* detLayer = nullptr;
     float radius = 0.;
 
     CLHEP::Hep3Vector vec(0., 1., 0.);
     vec.setTheta(theta);
     vec.setPhi(phi);
 
     // Get the det layer on which the state should be put
     if (barrel) {
       LogTrace(metname) << "The seed is in the barrel";
 
       // MB2
       DetId id = DTChamberId(0, 2, 0);
       detLayer = theService->detLayerGeometry()->idToLayer(id);
       LogTrace(metname) << "L2 Layer: " << debug.dumpLayer(detLayer);
 
       const BoundSurface* sur = &(detLayer->surface());
       const BoundCylinder* bc = dynamic_cast<const BoundCylinder*>(sur);
 
       radius = fabs(bc->radius() / sin(theta));
 
       LogTrace(metname) << "radius " << radius;
 
       if (pt < 3.5)
         pt = 3.5;
     } else {
       LogTrace(metname) << "The seed is in the endcap";
 
       DetId id;
       // ME2
       if (theta < Geom::pi() / 2.)
         id = CSCDetId(1, 2, 0, 0, 0);
       else
         id = CSCDetId(2, 2, 0, 0, 0);
 
       detLayer = theService->detLayerGeometry()->idToLayer(id);
       LogTrace(metname) << "L2 Layer: " << debug.dumpLayer(detLayer);
 
       radius = fabs(detLayer->position().z() / cos(theta));
 
       if (pt < 1.0)
         pt = 1.0;
     }
 
     vec.setMag(radius);
 
     GlobalPoint pos(vec.x(), vec.y(), vec.z());
 
     GlobalVector mom(pt * cos(phi), pt * sin(phi), pt * cos(theta) / sin(theta));
 
     GlobalTrajectoryParameters param(pos, mom, charge, &*theService->magneticField());
     AlgebraicSymMatrix55 mat;
 
     mat[0][0] = (0.25 / pt) * (0.25 / pt);  // sigma^2(charge/abs_momentum)
     if (!barrel)
       mat[0][0] = (0.4 / pt) * (0.4 / pt);
 
     //Assign q/pt = 0 +- 1/pt if charge has been declared invalid
     if (!valid_charge)
       mat[0][0] = (1. / pt) * (1. / pt);
 
     mat[1][1] = 0.05 * 0.05;  // sigma^2(lambda)
     mat[2][2] = 0.2 * 0.2;    // sigma^2(phi)
     mat[3][3] = 20. * 20.;    // sigma^2(x_transverse))
     mat[4][4] = 20. * 20.;    // sigma^2(y_transverse))
 
     CurvilinearTrajectoryError error(mat);
 
     const FreeTrajectoryState state(param, error);
 
     LogTrace(metname) << "Free trajectory State from the parameters";
     LogTrace(metname) << debug.dumpFTS(state);
 
     // Propagate the state on the MB2/ME2 surface
     TrajectoryStateOnSurface tsos = theService->propagator(thePropagatorName)->propagate(state, detLayer->surface());
 
     LogTrace(metname) << "State after the propagation on the layer";
     LogTrace(metname) << debug.dumpLayer(detLayer);
     LogTrace(metname) << debug.dumpFTS(state);
 
     if (tsos.isValid()) {
       // Get the compatible dets on the layer
       std::vector<pair<const GeomDet*, TrajectoryStateOnSurface> > detsWithStates =
           detLayer->compatibleDets(tsos, *theService->propagator(thePropagatorName), *theEstimator);
       if (!detsWithStates.empty()) {
         TrajectoryStateOnSurface newTSOS = detsWithStates.front().second;
         const GeomDet* newTSOSDet = detsWithStates.front().first;
 
         LogTrace(metname) << "Most compatible det";
         LogTrace(metname) << debug.dumpMuonId(newTSOSDet->geographicalId());
 
         LogDebug(metname) << "L1 info: Det and State:";
         LogDebug(metname) << debug.dumpMuonId(newTSOSDet->geographicalId());
 
         if (newTSOS.isValid()) {
           //LogDebug(metname) << "(x, y, z) = (" << newTSOS.globalPosition().x() << ", "
           //                  << newTSOS.globalPosition().y() << ", " << newTSOS.globalPosition().z() << ")";
           LogDebug(metname) << "pos: (r=" << newTSOS.globalPosition().mag()
                             << ", phi=" << newTSOS.globalPosition().phi() << ", eta=" << newTSOS.globalPosition().eta()
                             << ")";
           LogDebug(metname) << "mom: (q*pt=" << newTSOS.charge() * newTSOS.globalMomentum().perp()
                             << ", phi=" << newTSOS.globalMomentum().phi() << ", eta=" << newTSOS.globalMomentum().eta()
                             << ")";
 
           //LogDebug(metname) << "State on it";
           //LogDebug(metname) << debug.dumpTSOS(newTSOS);
 
           //PTrajectoryStateOnDet seedTSOS;
           edm::OwnVector<TrackingRecHit> container;
 
           if (useOfflineSeed) {
             const TrajectorySeed* assoOffseed = associateOfflineSeedToL1(offlineSeedHandle, offlineSeedMap, newTSOS);
 
             if (assoOffseed != nullptr) {
               PTrajectoryStateOnDet const& seedTSOS = assoOffseed->startingState();
               for (auto const& tsci : assoOffseed->recHits()) {
                 container.push_back(tsci);
               }
               output->push_back(
                   L2MuonTrajectorySeed(seedTSOS, container, alongMomentum, L1MuonParticleRef(muColl, l1ParticleIndex)));
             } else {
               if (useUnassociatedL1) {
                 // convert the TSOS into a PTSOD
                 PTrajectoryStateOnDet const& seedTSOS =
                     trajectoryStateTransform::persistentState(newTSOS, newTSOSDet->geographicalId().rawId());
                 output->push_back(L2MuonTrajectorySeed(
                     seedTSOS, container, alongMomentum, L1MuonParticleRef(muColl, l1ParticleIndex)));
               }
             }
           } else {
             // convert the TSOS into a PTSOD
             PTrajectoryStateOnDet const& seedTSOS =
                 trajectoryStateTransform::persistentState(newTSOS, newTSOSDet->geographicalId().rawId());
             output->push_back(
                 L2MuonTrajectorySeed(seedTSOS, container, alongMomentum, L1MuonParticleRef(muColl, l1ParticleIndex)));
           }
         }
       }
     }
   }
 
   iEvent.put(std::move(output));
 }