TevMuonProducer Class Reference

#include <TevMuonProducer.h>

Inheritance diagram for TevMuonProducer:

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override
	reconstruct muons More...
	TevMuonProducer (const edm::ParameterSet &)
	constructor with config More...
	~TevMuonProducer () 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 Member Functions
void	setAlias (std::string alias)

Private Attributes
edm::EDGetTokenT< reco::TrackCollection >	glbMuonsToken
edm::EDGetTokenT< std::vector< Trajectory > >	glbMuonsTrajToken
std::string	theAlias
edm::InputTag	theGLBCollectionLabel
	STA Label. More...
std::vector< int >	theRefitIndex
std::vector< std::string >	theRefits
std::unique_ptr< GlobalMuonRefitter >	theRefitter
std::unique_ptr< MuonServiceProxy >	theService
	the event setup proxy, it takes care the services update More...
std::unique_ptr< MuonTrackLoader >	theTrackLoader

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

TeV muon reconstructor:

Author

Piotr Traczyk (SINS Warsaw)

Global muon reconstructor: reconstructs muons using DT, CSC, RPC and tracker information,
starting from a standalone reonstructed muon.

Author

R.Bellan - INFN TO

Definition at line 43 of file TevMuonProducer.h.

Constructor & Destructor Documentation

TevMuonProducer()

TevMuonProducer::TevMuonProducer

(

const edm::ParameterSet &

parameterSet

)

constructor with config

Definition at line 36 of file TevMuonProducer.cc.

                                                                 {
  LogDebug("Muon|RecoMuon|TevMuonProducer") << "constructor called" << endl;
 
  // GLB Muon Collection Label
  theGLBCollectionLabel = parameterSet.getParameter<InputTag>("MuonCollectionLabel");
  glbMuonsToken = consumes<reco::TrackCollection>(theGLBCollectionLabel);
  glbMuonsTrajToken = consumes<std::vector<Trajectory> >(theGLBCollectionLabel.label());
 
  // service parameters
  ParameterSet serviceParameters = parameterSet.getParameter<ParameterSet>("ServiceParameters");
 
  // the services
  theService = std::make_unique<MuonServiceProxy>(serviceParameters, consumesCollector());
  edm::ConsumesCollector iC = consumesCollector();
 
  // TrackRefitter parameters
  ParameterSet refitterParameters = parameterSet.getParameter<ParameterSet>("RefitterParameters");
  theRefitter = std::make_unique<GlobalMuonRefitter>(refitterParameters, theService.get(), iC);
 
  // TrackLoader parameters
  ParameterSet trackLoaderParameters = parameterSet.getParameter<ParameterSet>("TrackLoaderParameters");
  theTrackLoader = std::make_unique<MuonTrackLoader>(trackLoaderParameters, iC, theService.get());
 
  theRefits = parameterSet.getParameter<std::vector<std::string> >("Refits");
  theRefitIndex = parameterSet.getParameter<std::vector<int> >("RefitIndex");
 
  for (unsigned int ww = 0; ww < theRefits.size(); ww++) {
    LogDebug("Muon|RecoMuon|TevMuonProducer") << "Refit " << theRefits[ww];
    produces<reco::TrackCollection>(theRefits[ww]);
    produces<TrackingRecHitCollection>(theRefits[ww]);
    produces<reco::TrackExtraCollection>(theRefits[ww]);
    produces<vector<Trajectory> >(theRefits[ww]);
    produces<TrajTrackAssociationCollection>(theRefits[ww]);
    produces<reco::TrackToTrackMap>(theRefits[ww]);
  }
  produces<DYTestimators>("dytInfo");
}

References edm::ParameterSet::getParameter(), LogDebug, and edm::parameterSet().

~TevMuonProducer()

TevMuonProducer::~TevMuonProducer ( )

override

destructor

Definition at line 77 of file TevMuonProducer.cc.

{ LogTrace("Muon|RecoMuon|TevMuonProducer") << "destructor called" << endl; }

References LogTrace.

Member Function Documentation

produce()

void TevMuonProducer::produce ( edm::Event &  event, const edm::EventSetup &  eventSetup  )

override

reconstruct muons

Definition at line 82 of file TevMuonProducer.cc.

                                                                        {
  const string metname = "Muon|RecoMuon|TevMuonProducer";
  LogTrace(metname) << endl << endl;
  LogTrace(metname) << "TeV Muon Reconstruction started" << endl;
 
  // Update the services
  theService->update(eventSetup);
 
  theRefitter->setEvent(event);
 
  theRefitter->setServices(theService->eventSetup());
 
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHand;
  eventSetup.get<TrackerTopologyRcd>().get(tTopoHand);
  const TrackerTopology* tTopo = tTopoHand.product();
 
  // Take the GLB muon container(s)
  Handle<reco::TrackCollection> glbMuons;
  event.getByToken(glbMuonsToken, glbMuons);
 
  auto dytInfo = std::make_unique<DYTestimators>();
  DYTestimators::Filler filler(*dytInfo);
  size_t GLBmuonSize = glbMuons->size();
  vector<DYTInfo> dytTmp(GLBmuonSize);
 
  Handle<vector<Trajectory> > glbMuonsTraj;
 
  LogTrace(metname) << "Taking " << glbMuons->size() << " Global Muons " << theGLBCollectionLabel << endl;
 
  vector<MuonTrajectoryBuilder::TrackCand> glbTrackCands;
 
  event.getByToken(glbMuonsTrajToken, glbMuonsTraj);
 
  const reco::TrackCollection* glbTracks = glbMuons.product();
 
  for (unsigned int ww = 0; ww < theRefits.size(); ww++) {
    LogDebug(metname) << "TeVRefit for Refit: " << theRefitIndex[ww];
    std::vector<std::pair<Trajectory*, reco::TrackRef> > miniMap;
    MuonTrackLoader::TrajectoryContainer trajectories;
    reco::TrackRef::key_type trackIndex = 0;
    int glbCounter = 0;
    for (reco::TrackCollection::const_iterator track = glbTracks->begin(); track != glbTracks->end();
         track++, ++trackIndex) {
      reco::TrackRef glbRef(glbMuons, trackIndex);
 
      vector<Trajectory> refitted = theRefitter->refit(*track, theRefitIndex[ww], tTopo);
 
      if (theRefits[ww] == "dyt")
        dytTmp[glbCounter] = *theRefitter->getDYTInfo();
      glbCounter++;
 
      if (!refitted.empty()) {
        auto refit = std::make_unique<Trajectory>(refitted.front());
        LogDebug(metname) << "TeVTrackLoader for Refit: " << theRefits[ww];
        std::pair<Trajectory*, reco::TrackRef> thisPair(refit.get(), glbRef);
        miniMap.push_back(thisPair);
        trajectories.push_back(std::move(refit));
      }
    }
    theTrackLoader->loadTracks(trajectories, event, miniMap, glbMuons, *tTopo, theRefits[ww]);
  }
 
  filler.insert(glbMuons, dytTmp.begin(), dytTmp.end());
  filler.fill();
  event.put(std::move(dytInfo), "dytInfo");
 
  LogTrace(metname) << "Done." << endl;
}

References options_cfi::eventSetup, trigObjTnPSource_cfi::filler, edm::get(), LogDebug, LogTrace, metname, eostools::move(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), HLT_2018_cff::track, and FastTrackerRecHitMaskProducer_cfi::trajectories.

setAlias()

void TevMuonProducer::setAlias ( std::string  alias )

inlineprivate

Definition at line 71 of file TevMuonProducer.h.

                                 {
    alias.erase(alias.size() - 1, alias.size());
    theAlias = alias;
  }

References SiStripOfflineCRack_cfg::alias, and theAlias.

Member Data Documentation

glbMuonsToken

edm::EDGetTokenT<reco::TrackCollection> TevMuonProducer::glbMuonsToken

private

Definition at line 57 of file TevMuonProducer.h.

glbMuonsTrajToken

edm::EDGetTokenT<std::vector<Trajectory> > TevMuonProducer::glbMuonsTrajToken

private

Definition at line 58 of file TevMuonProducer.h.

theAlias

std::string TevMuonProducer::theAlias

private

Definition at line 67 of file TevMuonProducer.h.

Referenced by setAlias().

theGLBCollectionLabel

edm::InputTag TevMuonProducer::theGLBCollectionLabel

private

STA Label.

Definition at line 56 of file TevMuonProducer.h.

theRefitIndex

std::vector<int> TevMuonProducer::theRefitIndex

private

Definition at line 69 of file TevMuonProducer.h.

theRefits

std::vector<std::string> TevMuonProducer::theRefits

private

Definition at line 68 of file TevMuonProducer.h.

theRefitter

std::unique_ptr<GlobalMuonRefitter> TevMuonProducer::theRefitter

private

Definition at line 63 of file TevMuonProducer.h.

theService

std::unique_ptr<MuonServiceProxy> TevMuonProducer::theService

private

the event setup proxy, it takes care the services update

Definition at line 61 of file TevMuonProducer.h.

theTrackLoader

std::unique_ptr<MuonTrackLoader> TevMuonProducer::theTrackLoader

private

Definition at line 65 of file TevMuonProducer.h.