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
	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
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
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	tTopoToken

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

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 39 of file TevMuonProducer.h.

Constructor & Destructor Documentation

TevMuonProducer()

TevMuonProducer::TevMuonProducer

(

const edm::ParameterSet &

parameterSet

)

constructor with config

Definition at line 35 of file TevMuonProducer.cc.

References edm::ParameterSet::getParameter(), glbMuonsToken, glbMuonsTrajToken, edm::InputTag::label(), LogDebug, edm::parameterSet(), theGLBCollectionLabel, theRefitIndex, theRefits, theRefitter, theService, and theTrackLoader.

                                                                 : tTopoToken(esConsumes()) {
  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");
}

~TevMuonProducer()

TevMuonProducer::~TevMuonProducer ( )

override

destructor

Definition at line 76 of file TevMuonProducer.cc.

References LogTrace.

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

Member Function Documentation

produce()

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

override

reconstruct muons

Definition at line 81 of file TevMuonProducer.cc.

References options_cfi::eventSetup, trigObjTnPSource_cfi::filler, glbMuonsToken, glbMuonsTrajToken, LogDebug, LogTrace, metname, eostools::move(), edm::Handle< T >::product(), theGLBCollectionLabel, theRefitIndex, theRefits, theRefitter, theService, theTrackLoader, HLT_2023v12_cff::track, FastTrackerRecHitMaskProducer_cfi::trajectories, and tTopoToken.

                                                                        {
  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
  const TrackerTopology* tTopo = &eventSetup.getData(tTopoToken);

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

setAlias()

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

inlineprivate

Definition at line 68 of file TevMuonProducer.h.

References SiStripOfflineCRack_cfg::alias, and theAlias.

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

Member Data Documentation

glbMuonsToken

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

private

Definition at line 53 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

glbMuonsTrajToken

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

private

Definition at line 54 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

theAlias

std::string TevMuonProducer::theAlias

private

Definition at line 64 of file TevMuonProducer.h.

Referenced by setAlias().

theGLBCollectionLabel

edm::InputTag TevMuonProducer::theGLBCollectionLabel

private

STA Label.

Definition at line 52 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

theRefitIndex

std::vector<int> TevMuonProducer::theRefitIndex

private

Definition at line 66 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

theRefits

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

private

Definition at line 65 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

theRefitter

std::unique_ptr<GlobalMuonRefitter> TevMuonProducer::theRefitter

private

Definition at line 60 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

theService

std::unique_ptr<MuonServiceProxy> TevMuonProducer::theService

private

the event setup proxy, it takes care the services update

Definition at line 58 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

theTrackLoader

std::unique_ptr<MuonTrackLoader> TevMuonProducer::theTrackLoader

private

Definition at line 62 of file TevMuonProducer.h.

Referenced by produce(), and TevMuonProducer().

tTopoToken

const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> TevMuonProducer::tTopoToken

private

Definition at line 55 of file TevMuonProducer.h.

Referenced by produce().