CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
MuonsFromRefitTracksProducer Class Reference
Inheritance diagram for MuonsFromRefitTracksProducer:
edm::stream::EDProducer<>

Public Member Functions

 MuonsFromRefitTracksProducer (const edm::ParameterSet &)
 
 ~MuonsFromRefitTracksProducer () override
 
- 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

reco::MuoncloneAndSwitchTrack (const reco::Muon &muon, const reco::Muon::MuonTrackTypePair &newTrack) const
 
void produce (edm::Event &, const edm::EventSetup &) override
 
bool storeMatchMaps (const edm::Event &event)
 

Private Attributes

bool fromCocktail
 
bool fromGlobalTrack
 
bool fromSigmaSwitch
 
bool fromTeVRefit
 
bool fromTMR
 
bool fromTrackerTrack
 
double nSigmaSwitch
 
double ptThreshold
 
edm::InputTag src
 
edm::EDGetTokenT< edm::View< reco::Muon > > srcToken_
 
std::string tevMuonTracks
 
double TMRcut
 
edm::Handle< reco::TrackToTrackMaptrackMap
 
edm::Handle< reco::TrackToTrackMaptrackMapDefault
 
edm::EDGetTokenT< reco::TrackToTrackMaptrackMapDefaultToken_
 
edm::Handle< reco::TrackToTrackMaptrackMapFirstHit
 
edm::EDGetTokenT< reco::TrackToTrackMaptrackMapFirstHitToken_
 
edm::Handle< reco::TrackToTrackMaptrackMapPicky
 
edm::EDGetTokenT< reco::TrackToTrackMaptrackMapPickyToken_
 
edm::EDGetTokenT< reco::TrackToTrackMaptrackMapToken_
 

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, GlobalCacheLuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCacheRunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 

Detailed Description

Definition at line 67 of file MuonsFromRefitTracksProducer.cc.

Constructor & Destructor Documentation

MuonsFromRefitTracksProducer::MuonsFromRefitTracksProducer ( const edm::ParameterSet cfg)
explicit

Definition at line 147 of file MuonsFromRefitTracksProducer.cc.

References fromTeVRefit, HLT_2018_cff::InputTag, src, srcToken_, tevMuonTracks, trackMapDefaultToken_, trackMapFirstHitToken_, trackMapPickyToken_, and trackMapToken_.

148  : src(cfg.getParameter<edm::InputTag>("src")),
149  fromTrackerTrack(cfg.getParameter<bool>("fromTrackerTrack")),
150  fromGlobalTrack(cfg.getParameter<bool>("fromGlobalTrack")),
151  tevMuonTracks(cfg.getParameter<std::string>("tevMuonTracks")),
152  fromCocktail(cfg.getParameter<bool>("fromCocktail")),
153  fromTMR(cfg.getParameter<bool>("fromTMR")),
154  TMRcut(cfg.getParameter<double>("TMRcut")),
155  fromSigmaSwitch(cfg.getParameter<bool>("fromSigmaSwitch")),
156  nSigmaSwitch(cfg.getParameter<double>("nSigmaSwitch")),
157  ptThreshold(cfg.getParameter<double>("ptThreshold")) {
158  fromTeVRefit = tevMuonTracks != "none";
159 
160  srcToken_ = consumes<edm::View<reco::Muon> >(src);
161  trackMapToken_ = consumes<reco::TrackToTrackMap>(edm::InputTag(tevMuonTracks, "default"));
162  trackMapDefaultToken_ = consumes<reco::TrackToTrackMap>(edm::InputTag(tevMuonTracks));
163  trackMapFirstHitToken_ = consumes<reco::TrackToTrackMap>(edm::InputTag(tevMuonTracks, "firstHit"));
164  trackMapPickyToken_ = consumes<reco::TrackToTrackMap>(edm::InputTag(tevMuonTracks, "picky"));
165 
166  produces<reco::MuonCollection>();
167 }
T getParameter(std::string const &) const
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapPickyToken_
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapToken_
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapDefaultToken_
edm::EDGetTokenT< edm::View< reco::Muon > > srcToken_
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapFirstHitToken_
MuonsFromRefitTracksProducer::~MuonsFromRefitTracksProducer ( )
inlineoverride

Definition at line 70 of file MuonsFromRefitTracksProducer.cc.

References cloneAndSwitchTrack(), produce(), and storeMatchMaps().

70 {}

Member Function Documentation

reco::Muon * MuonsFromRefitTracksProducer::cloneAndSwitchTrack ( const reco::Muon muon,
const reco::Muon::MuonTrackTypePair newTrack 
) const
private

Definition at line 181 of file MuonsFromRefitTracksProducer.cc.

References reco::Muon::clone(), reco::Muon::innerTrack(), amptDefaultParameters_cff::mu, fastsim::Constants::muMass, reco::Muon::outerTrack(), AlCaHLTBitMon_ParallelJobs::p, p4, reco::Muon::setBestTrack(), reco::LeafCandidate::setCharge(), reco::Muon::setGlobalTrack(), reco::Muon::setInnerTrack(), reco::Muon::setOuterTrack(), reco::LeafCandidate::setP4(), reco::LeafCandidate::setVertex(), mathSSE::sqrt(), and badGlobalMuonTaggersAOD_cff::vtx.

Referenced by produce(), and ~MuonsFromRefitTracksProducer().

182  {
183  // Muon mass to make a four-vector out of the new track.
184  static const double muMass = 0.10566;
185 
186  reco::TrackRef tkTrack = muon.innerTrack();
187  reco::TrackRef muTrack = muon.outerTrack();
188 
189  // Make up a real Muon from the tracker track.
190  reco::Particle::Point vtx(newTrack.first->vx(), newTrack.first->vy(), newTrack.first->vz());
192  double p = newTrack.first->p();
193  p4.SetXYZT(newTrack.first->px(), newTrack.first->py(), newTrack.first->pz(), sqrt(p * p + muMass * muMass));
194 
195  reco::Muon* mu = muon.clone();
196  mu->setCharge(newTrack.first->charge());
197  mu->setP4(p4);
198  mu->setVertex(vtx);
199  mu->setGlobalTrack(newTrack.first);
200  mu->setInnerTrack(tkTrack);
201  mu->setOuterTrack(muTrack);
202  mu->setBestTrack(newTrack.second);
203  return mu;
204 }
virtual TrackRef innerTrack() const
Definition: Muon.h:45
static double constexpr muMass
Muon mass [GeV].
Definition: Constants.h:14
void setVertex(const Point &vertex) override
set vertex
void setCharge(Charge q) final
set electric charge
virtual void setInnerTrack(const TrackRef &t)
set reference to Track
T sqrt(T t)
Definition: SSEVec.h:19
double p4[4]
Definition: TauolaWrapper.h:92
math::XYZPoint Point
point in the space
Definition: Particle.h:25
virtual TrackRef outerTrack() const
reference to Track reconstructed in the muon detector only
Definition: Muon.h:48
virtual void setOuterTrack(const TrackRef &t)
set reference to Track
virtual void setBestTrack(MuonTrackType muonType)
Definition: Muon.h:86
math::XYZTLorentzVector LorentzVector
Lorentz vector.
Definition: Particle.h:21
void setP4(const LorentzVector &p4) final
set 4-momentum
Muon * clone() const override
create a clone
virtual void setGlobalTrack(const TrackRef &t)
set reference to Track
void MuonsFromRefitTracksProducer::produce ( edm::Event event,
const edm::EventSetup eSetup 
)
overrideprivate

Definition at line 206 of file MuonsFromRefitTracksProducer.cc.

References HLT_2018_cff::cands, edm::ViewBase::clone(), cloneAndSwitchTrack(), reco::Muon::CombinedTrack, DEFINE_FWK_MODULE, edm::AssociationMap< Tag >::end(), edm::HandleBase::failedToGet(), edm::AssociationMap< Tag >::find(), fromCocktail, fromGlobalTrack, fromSigmaSwitch, fromTeVRefit, fromTMR, fromTrackerTrack, reco::Muon::InnerTrack, dqmdumpme::last, eostools::move(), HLT_2018_cff::muon, PDWG_BPHSkim_cff::muons, nSigmaSwitch, convertSQLiteXML::ok, ptThreshold, reco::Muon::setGlobalTrack(), reco::Muon::setInnerTrack(), reco::Muon::setOuterTrack(), sigmaSwitch(), src, srcToken_, storeMatchMaps(), tevMuonTracks, muon::tevOptimized(), tevOptimizedTMR(), TMRcut, trackMap, trackMapFirstHit, and edm::helpers::KeyVal< K, V >::val.

Referenced by ~MuonsFromRefitTracksProducer().

206  {
207  // Get the global muons from the event.
209  event.getByToken(srcToken_, muons);
210 
211  // If we can't get the global muon collection, or below the
212  // track-to-track maps needed, still produce an empty collection of
213  // muons so consumers don't throw an exception.
214  bool ok = !muons.failedToGet();
215 
216  // If we're instructed to use the TeV refit tracks in some way, we
217  // need the track-to-track maps. If we're making a cocktail muon,
218  // get all three track maps (the cocktail ingredients); else just
219  // get the map which takes the above global tracks to the desired
220  // TeV-muon refitted tracks (firstHit or picky).
221  if (ok && fromTeVRefit)
222  ok = storeMatchMaps(event);
223 
224  // Make the output collection.
225  auto cands = std::make_unique<reco::MuonCollection>();
226 
227  if (ok) {
229  for (muon = muons->begin(); muon != muons->end(); muon++) {
230  // Filter out the so-called trackerMuons and stand-alone muons
231  // (and caloMuons, if they were ever to get into the input muons
232  // collection).
233  if (!muon->isGlobalMuon())
234  continue;
235 
236  if (fromTeVRefit || fromSigmaSwitch) {
237  // Start out with a null TrackRef.
239 
240  // If making a cocktail muon, use tevOptimized() to get the track
241  // desired. Otherwise, get the refit track from the desired track
242  // map.
243  if (fromTMR)
244  tevTk = tevOptimizedTMR(*muon, *trackMapFirstHit, TMRcut);
245  else if (fromCocktail)
246  tevTk = muon::tevOptimized(*muon);
247  else if (fromSigmaSwitch)
248  tevTk = sigmaSwitch(*muon, nSigmaSwitch, ptThreshold);
249  else {
250  reco::TrackToTrackMap::const_iterator tevTkRef = trackMap->find(muon->combinedMuon());
251  if (tevTkRef != trackMap->end())
252  tevTk = make_pair(tevTkRef->val, reco::Muon::CombinedTrack);
253  }
254 
255  // If the TrackRef is valid, make a new Muon that has the same
256  // tracker and stand-alone tracks, but has the refit track as
257  // its global track.
258  if (tevTk.first.isNonnull())
259  cands->push_back(*cloneAndSwitchTrack(*muon, tevTk));
260  } else if (fromTrackerTrack)
261  cands->push_back(*cloneAndSwitchTrack(*muon, make_pair(muon->innerTrack(), reco::Muon::InnerTrack)));
262  else if (fromGlobalTrack)
263  cands->push_back(*cloneAndSwitchTrack(*muon, make_pair(muon->globalTrack(), reco::Muon::CombinedTrack)));
264  else {
265  cands->push_back(*muon->clone());
266 
267  // Just cloning does not work in the case of the source being
268  // a pat::Muon with embedded track references -- these do not
269  // get copied. Explicitly set them.
270  reco::Muon& last = cands->at(cands->size() - 1);
271  if (muon->globalTrack().isTransient())
272  last.setGlobalTrack(muon->globalTrack());
273  if (muon->innerTrack().isTransient())
274  last.setInnerTrack(muon->innerTrack());
275  if (muon->outerTrack().isTransient())
276  last.setOuterTrack(muon->outerTrack());
277  }
278  }
279  } else
280  edm::LogWarning("MuonsFromRefitTracksProducer") << "either " << src << " or the track map(s) " << tevMuonTracks
281  << " not present in the event; producing empty collection";
282 
283  event.put(std::move(cands));
284 }
reco::Muon::MuonTrackTypePair sigmaSwitch(const reco::Muon &muon, const double nSigma, const double ptThreshold)
const_iterator end() const
last iterator over the map (read only)
reco::Muon * cloneAndSwitchTrack(const reco::Muon &muon, const reco::Muon::MuonTrackTypePair &newTrack) const
const_iterator find(const key_type &k) const
find element with specified reference key
edm::Handle< reco::TrackToTrackMap > trackMap
virtual void setInnerTrack(const TrackRef &t)
set reference to Track
bool storeMatchMaps(const edm::Event &event)
reco::Muon::MuonTrackTypePair tevOptimized(const reco::TrackRef &combinedTrack, const reco::TrackRef &trackerTrack, const reco::TrackRef &tpfmsTrack, const reco::TrackRef &pickyTrack, const reco::TrackRef &dytTrack, const double ptThreshold=200., const double tune1=17., const double tune2=40., const double dptcut=0.25)
Definition: MuonCocktails.cc:9
std::unique_ptr< ViewBase > clone() const
Definition: View.cc:11
edm::Handle< reco::TrackToTrackMap > trackMapFirstHit
reco::Muon::MuonTrackTypePair tevOptimizedTMR(const reco::Muon &muon, const reco::TrackToTrackMap &fmsMap, const double cut)
bool failedToGet() const
Definition: HandleBase.h:72
edm::EDGetTokenT< edm::View< reco::Muon > > srcToken_
virtual void setOuterTrack(const TrackRef &t)
set reference to Track
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
std::pair< TrackRef, Muon::MuonTrackType > MuonTrackTypePair
Definition: Muon.h:38
virtual void setGlobalTrack(const TrackRef &t)
set reference to Track
def move(src, dest)
Definition: eostools.py:511
bool MuonsFromRefitTracksProducer::storeMatchMaps ( const edm::Event event)
private

Definition at line 169 of file MuonsFromRefitTracksProducer.cc.

References edm::HandleBase::failedToGet(), fromCocktail, fromTMR, trackMap, trackMapDefault, trackMapDefaultToken_, trackMapFirstHit, trackMapFirstHitToken_, trackMapPicky, trackMapPickyToken_, and trackMapToken_.

Referenced by produce(), and ~MuonsFromRefitTracksProducer().

169  {
170  if (fromCocktail || fromTMR) {
171  event.getByToken(trackMapDefaultToken_, trackMapDefault);
172  event.getByToken(trackMapFirstHitToken_, trackMapFirstHit);
173  event.getByToken(trackMapPickyToken_, trackMapPicky);
175  } else {
176  event.getByToken(trackMapToken_, trackMap);
177  return !trackMap.failedToGet();
178  }
179 }
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapPickyToken_
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapToken_
edm::Handle< reco::TrackToTrackMap > trackMapDefault
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapDefaultToken_
edm::Handle< reco::TrackToTrackMap > trackMap
edm::Handle< reco::TrackToTrackMap > trackMapFirstHit
bool failedToGet() const
Definition: HandleBase.h:72
edm::Handle< reco::TrackToTrackMap > trackMapPicky
edm::EDGetTokenT< reco::TrackToTrackMap > trackMapFirstHitToken_

Member Data Documentation

bool MuonsFromRefitTracksProducer::fromCocktail
private

Definition at line 110 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce(), and storeMatchMaps().

bool MuonsFromRefitTracksProducer::fromGlobalTrack
private

Definition at line 96 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce().

bool MuonsFromRefitTracksProducer::fromSigmaSwitch
private

Definition at line 121 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce().

bool MuonsFromRefitTracksProducer::fromTeVRefit
private

Definition at line 100 of file MuonsFromRefitTracksProducer.cc.

Referenced by MuonsFromRefitTracksProducer(), and produce().

bool MuonsFromRefitTracksProducer::fromTMR
private

Definition at line 114 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce(), and storeMatchMaps().

bool MuonsFromRefitTracksProducer::fromTrackerTrack
private

Definition at line 90 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce().

double MuonsFromRefitTracksProducer::nSigmaSwitch
private

Definition at line 124 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce().

double MuonsFromRefitTracksProducer::ptThreshold
private

Definition at line 127 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce().

edm::InputTag MuonsFromRefitTracksProducer::src
private

Definition at line 84 of file MuonsFromRefitTracksProducer.cc.

Referenced by MuonsFromRefitTracksProducer(), and produce().

edm::EDGetTokenT<edm::View<reco::Muon> > MuonsFromRefitTracksProducer::srcToken_
private

Definition at line 140 of file MuonsFromRefitTracksProducer.cc.

Referenced by MuonsFromRefitTracksProducer(), and produce().

std::string MuonsFromRefitTracksProducer::tevMuonTracks
private

Definition at line 105 of file MuonsFromRefitTracksProducer.cc.

Referenced by MuonsFromRefitTracksProducer(), and produce().

double MuonsFromRefitTracksProducer::TMRcut
private

Definition at line 117 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce().

edm::Handle<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMap
private

Definition at line 132 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce(), and storeMatchMaps().

edm::Handle<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMapDefault
private

Definition at line 135 of file MuonsFromRefitTracksProducer.cc.

Referenced by storeMatchMaps().

edm::EDGetTokenT<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMapDefaultToken_
private
edm::Handle<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMapFirstHit
private

Definition at line 136 of file MuonsFromRefitTracksProducer.cc.

Referenced by produce(), and storeMatchMaps().

edm::EDGetTokenT<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMapFirstHitToken_
private
edm::Handle<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMapPicky
private

Definition at line 137 of file MuonsFromRefitTracksProducer.cc.

Referenced by storeMatchMaps().

edm::EDGetTokenT<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMapPickyToken_
private
edm::EDGetTokenT<reco::TrackToTrackMap> MuonsFromRefitTracksProducer::trackMapToken_
private