CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Types | Private Member Functions | Private Attributes
HLTScoutingMuonProducer Class Reference

#include <HLTScoutingMuonProducer.h>

Inheritance diagram for HLTScoutingMuonProducer:
edm::global::EDProducer<> edm::global::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 HLTScoutingMuonProducer (const edm::ParameterSet &)
 
 ~HLTScoutingMuonProducer () override
 
- Public Member Functions inherited from edm::global::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
 
EDProducer & operator= (const EDProducer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
bool wantsStreamLuminosityBlocks () const final
 
bool wantsStreamRuns () const final
 
- Public Member Functions inherited from edm::global::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducerBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
std::vector< bool > const & recordProvenanceList () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
TypeLabelList const & typeLabelList () const
 used by the fwk to register the list of products of this module More...
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Types

typedef edm::AssociationMap< edm::OneToValue< std::vector< reco::RecoChargedCandidate >, float, unsigned int > > RecoChargedCandMap
 

Private Member Functions

void produce (edm::StreamID sid, edm::Event &iEvent, edm::EventSetup const &setup) const final
 

Private Attributes

const edm::EDGetTokenT< reco::RecoChargedCandidateCollection > ChargedCandidateCollection_
 
const edm::EDGetTokenT< reco::VertexCollection > displacedvertexCollection_
 
const edm::EDGetTokenT< RecoChargedCandMap > EcalPFClusterIsoMap_
 
const edm::EDGetTokenT< RecoChargedCandMap > HcalPFClusterIsoMap_
 
const edm::EDGetTokenT< reco::MuonTrackLinksCollection > linkToken_
 
const double minVtxProbCut
 
const edm::EDGetTokenT< reco::MuonCollection > MuonCollection_
 
const double muonEtaCut
 
const double muonPtCut
 
const edm::EDGetTokenT< reco::TrackCollection > TrackCollection_
 
const edm::EDGetTokenT< edm::ValueMap< double > > TrackIsoMap_
 

Additional Inherited Members

- Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
template<typename T >
using BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >
 
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto produces (std::string instanceName) noexcept
 declare what type of product will make and with which optional label More...
 
template<Transition B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<BranchType B>
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
BranchAliasSetter produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces ()
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType > produces ()
 
template<class ProductType >
BranchAliasSetterT< ProductType > produces (std::string instanceName)
 
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType > produces ()
 
template<Transition Tr = Transition::Event>
auto produces () noexcept
 
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Description: Producer for ScoutingMuon

Description: Producer for Run3ScoutingMuon

Definition at line 50 of file HLTScoutingMuonProducer.h.

Member Typedef Documentation

◆ RecoChargedCandMap

typedef edm::AssociationMap<edm::OneToValue<std::vector<reco::RecoChargedCandidate>, float, unsigned int> > HLTScoutingMuonProducer::RecoChargedCandMap
private

Definition at line 52 of file HLTScoutingMuonProducer.h.

Constructor & Destructor Documentation

◆ HLTScoutingMuonProducer()

HLTScoutingMuonProducer::HLTScoutingMuonProducer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 24 of file HLTScoutingMuonProducer.cc.

25  : ChargedCandidateCollection_(
26  consumes<reco::RecoChargedCandidateCollection>(iConfig.getParameter<edm::InputTag>("ChargedCandidates"))),
27  displacedvertexCollection_(
28  consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("displacedvertexCollection"))),
29  MuonCollection_(consumes<reco::MuonCollection>(iConfig.getParameter<edm::InputTag>("InputMuons"))),
30  linkToken_(consumes<reco::MuonTrackLinksCollection>(iConfig.getParameter<edm::InputTag>("InputLinks"))),
31  TrackCollection_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("Tracks"))),
32  EcalPFClusterIsoMap_(consumes<RecoChargedCandMap>(iConfig.getParameter<edm::InputTag>("EcalPFClusterIsoMap"))),
33  HcalPFClusterIsoMap_(consumes<RecoChargedCandMap>(iConfig.getParameter<edm::InputTag>("HcalPFClusterIsoMap"))),
34  TrackIsoMap_(consumes<edm::ValueMap<double>>(iConfig.getParameter<edm::InputTag>("TrackIsoMap"))),
35  muonPtCut(iConfig.getParameter<double>("muonPtCut")),
36  muonEtaCut(iConfig.getParameter<double>("muonEtaCut")),
37  minVtxProbCut(iConfig.getParameter<double>("minVtxProbCut")) {
38  //register products
39  produces<Run3ScoutingMuonCollection>();
40  produces<Run3ScoutingVertexCollection>("displacedVtx");
41 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
HLTScoutingMuonProducer::MuonCollection_
const edm::EDGetTokenT< reco::MuonCollection > MuonCollection_
Definition: HLTScoutingMuonProducer.h:65
HLTScoutingMuonProducer::EcalPFClusterIsoMap_
const edm::EDGetTokenT< RecoChargedCandMap > EcalPFClusterIsoMap_
Definition: HLTScoutingMuonProducer.h:68
HLTScoutingMuonProducer::TrackCollection_
const edm::EDGetTokenT< reco::TrackCollection > TrackCollection_
Definition: HLTScoutingMuonProducer.h:67
HLTScoutingMuonProducer::HcalPFClusterIsoMap_
const edm::EDGetTokenT< RecoChargedCandMap > HcalPFClusterIsoMap_
Definition: HLTScoutingMuonProducer.h:69
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:161
HLTScoutingMuonProducer::TrackIsoMap_
const edm::EDGetTokenT< edm::ValueMap< double > > TrackIsoMap_
Definition: HLTScoutingMuonProducer.h:70
HLTScoutingMuonProducer::displacedvertexCollection_
const edm::EDGetTokenT< reco::VertexCollection > displacedvertexCollection_
Definition: HLTScoutingMuonProducer.h:64
HLTScoutingMuonProducer::ChargedCandidateCollection_
const edm::EDGetTokenT< reco::RecoChargedCandidateCollection > ChargedCandidateCollection_
Definition: HLTScoutingMuonProducer.h:63
HLTScoutingMuonProducer::linkToken_
const edm::EDGetTokenT< reco::MuonTrackLinksCollection > linkToken_
Definition: HLTScoutingMuonProducer.h:66

◆ ~HLTScoutingMuonProducer()

HLTScoutingMuonProducer::~HLTScoutingMuonProducer ( )
overridedefault

Member Function Documentation

◆ fillDescriptions()

void HLTScoutingMuonProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 298 of file HLTScoutingMuonProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_2022v15_cff::InputTag.

298  {
299  edm::ParameterSetDescription desc;
300  desc.add<edm::InputTag>("ChargedCandidates", edm::InputTag("hltIterL3MuonCandidatesNoVtx"));
301  desc.add<edm::InputTag>("displacedvertexCollection", edm::InputTag("hltDisplacedmumuVtxProducer"));
302  desc.add<edm::InputTag>("InputMuons", edm::InputTag("hltIterL3MuonsNoVtx"));
303  desc.add<edm::InputTag>("InputLinks", edm::InputTag("hltL3MuonsIterL3LinksNoVtx"));
304  desc.add<edm::InputTag>("Tracks", edm::InputTag("hltPixelTracks"));
305  desc.add<edm::InputTag>("EcalPFClusterIsoMap", edm::InputTag("hltMuonEcalMFPFClusterIsoForMuonsNoVtx"));
306  desc.add<edm::InputTag>("HcalPFClusterIsoMap", edm::InputTag("hltMuonHcalPFClusterIsoForMuonsNoVtx"));
307  desc.add<edm::InputTag>("TrackIsoMap",
308  edm::InputTag("hltMuonTkRelIsolationCut0p09MapNoVtx:combinedRelativeIsoDeposits"));
309  desc.add<double>("muonPtCut", 3.0);
310  desc.add<double>("muonEtaCut", 2.4);
311  desc.add<double>("minVtxProbCut", 0.001);
312 
313  descriptions.add("hltScoutingMuonProducer", desc);
314 }
HLT_2022v15_cff.InputTag
InputTag
Definition: HLT_2022v15_cff.py:60417
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57

◆ produce()

void HLTScoutingMuonProducer::produce ( edm::StreamID  sid,
edm::Event &  iEvent,
edm::EventSetup const &  setup 
) const
finalprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 46 of file HLTScoutingMuonProducer.cc.

References funct::abs(), ChargedCandidateCollection_, reco::deltaR2(), HLTMuonOfflineAnalyzer_cfi::deltaR2, HLT_2022v15_cff::displacedvertexCollection, displacedvertexCollection_, MillePedeFileConverter_cfg::e, HLT_2022v15_cff::EcalPFClusterIsoMap, EcalPFClusterIsoMap_, reco::TrackBase::eta(), first, edm::getRef(), HLT_2022v15_cff::HcalPFClusterIsoMap, HcalPFClusterIsoMap_, mps_fire::i, iEvent, edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNull(), MainPageGenerator::link, electronStore::links, linkToken_, minVtxProbCut, eostools::move(), MuonCollection_, muonEtaCut, muonPtCut, TrackCollections2monitor_cff::normalizedChi2, reco::TrackBase::phi(), reco::TrackBase::pt(), edm::second(), HLT_2022v15_cff::track, HLT_2022v15_cff::TrackIsoMap, and TrackIsoMap_.

46  {
47  using namespace edm;
48 
49  std::unique_ptr<Run3ScoutingMuonCollection> outMuons(new Run3ScoutingMuonCollection());
50  std::unique_ptr<Run3ScoutingVertexCollection> dispVertices(new Run3ScoutingVertexCollection());
51 
52  // Get RecoChargedCandidate
53  Handle<reco::RecoChargedCandidateCollection> ChargedCandidateCollection;
54  if (!iEvent.getByToken(ChargedCandidateCollection_, ChargedCandidateCollection)) {
55  iEvent.put(std::move(outMuons));
56  iEvent.put(std::move(dispVertices), "displacedVtx");
57  return;
58  }
59 
60  std::pair<reco::RecoChargedCandidate, reco::RecoChargedCandidate> ivtxMuPair;
61  std::vector<std::pair<reco::RecoChargedCandidate, reco::RecoChargedCandidate>> vtxMuPair;
62 
63  //get displaced vertices
64  Handle<reco::VertexCollection> displacedvertexCollection;
65  if (iEvent.getByToken(displacedvertexCollection_, displacedvertexCollection)) {
66  for (auto& dispvtx : *displacedvertexCollection) {
67  if (!dispvtx.isValid())
68  continue;
69  float vtxProb = 0.0;
70  if ((dispvtx.chi2() >= 0.0) && (dispvtx.ndof() > 0))
71  vtxProb = TMath::Prob(dispvtx.chi2(), dispvtx.ndof());
72  if (vtxProb < minVtxProbCut)
73  continue;
74 
75  // Get the 2 tracks associated to displaced vertex
76  auto trackIt = dispvtx.tracks_begin();
77  reco::TrackRef vertextkRef1 = (*trackIt).castTo<reco::TrackRef>();
78  trackIt++;
79  reco::TrackRef vertextkRef2 = (*trackIt).castTo<reco::TrackRef>();
80 
81  // Get the muons associated with the tracks
82  int iFoundRefs = 0;
83  for (auto const& cand : *ChargedCandidateCollection) {
84  reco::TrackRef tkRef = cand.get<reco::TrackRef>();
85  if (tkRef == vertextkRef1) {
86  ivtxMuPair.first = cand;
87  iFoundRefs++;
88  }
89  if (tkRef == vertextkRef2) {
90  ivtxMuPair.second = cand;
91  iFoundRefs++;
92  }
93  }
94  if (iFoundRefs < 2)
95  continue;
96  vtxMuPair.push_back(ivtxMuPair);
97 
98  dispVertices->emplace_back(dispvtx.x(),
99  dispvtx.y(),
100  dispvtx.z(),
101  dispvtx.zError(),
102  dispvtx.xError(),
103  dispvtx.yError(),
104  dispvtx.tracksSize(),
105  dispvtx.chi2(),
106  dispvtx.ndof(),
107  dispvtx.isValid());
108  }
109  }
110 
111  // Get Muon collection
112  Handle<reco::MuonCollection> MuonCollection;
113  if (!iEvent.getByToken(MuonCollection_, MuonCollection)) {
114  iEvent.put(std::move(outMuons));
115  iEvent.put(std::move(dispVertices), "displacedVtx");
116  return;
117  }
118 
119  // Get Links collection
120  edm::Handle<reco::MuonTrackLinksCollection> links;
121  if (!iEvent.getByToken(linkToken_, links)) {
122  iEvent.put(std::move(outMuons));
123  iEvent.put(std::move(dispVertices), "displacedVtx");
124  return;
125  }
126 
127  // Get EcalPFClusterIsoMap
128  Handle<RecoChargedCandMap> EcalPFClusterIsoMap;
129  iEvent.getByToken(EcalPFClusterIsoMap_, EcalPFClusterIsoMap);
130 
131  // Get HcalPFClusterIsoMap
132  Handle<RecoChargedCandMap> HcalPFClusterIsoMap;
133  iEvent.getByToken(HcalPFClusterIsoMap_, HcalPFClusterIsoMap);
134 
135  // Get TrackIsoMap
136  Handle<ValueMap<double>> TrackIsoMap;
137  iEvent.getByToken(TrackIsoMap_, TrackIsoMap);
138 
139  // Produce muons
140  std::vector<int> vtxInd;
141  float minDR2 = 1e-06;
142  int index = 0;
143 
144  for (auto& muon : *ChargedCandidateCollection) {
145  reco::RecoChargedCandidateRef muonRef = getRef(ChargedCandidateCollection, index);
146  ++index;
147  if (muonRef.isNull() || !muonRef.isAvailable())
148  continue;
149 
150  reco::TrackRef track = muon.track();
151  if (track.isNull() || !track.isAvailable())
152  continue;
153 
154  int validStandAloneMuonHits = 0;
155  int matchedStandAloneMuonStations = 0;
156  for (auto const& link : *links) {
157  const reco::Track& trackerTrack = *link.trackerTrack();
158  float dR2 = deltaR2(track->eta(), track->phi(), trackerTrack.eta(), trackerTrack.phi());
159  float dPt = std::abs(track->pt() - trackerTrack.pt());
160  if (track->pt() != 0)
161  dPt = dPt / track->pt();
162 
163  if (dR2 < 0.02 * 0.02 and dPt < 0.001) {
164  if (link.standAloneTrack().isNonnull()) {
165  validStandAloneMuonHits = link.standAloneTrack()->hitPattern().numberOfValidMuonHits();
166  matchedStandAloneMuonStations = link.standAloneTrack()->hitPattern().muonStationsWithValidHits();
167  }
168  }
169  }
170 
171  unsigned int recoMuonType = 2; // Global muon
172  float normalizedChi2 = 999.0;
173  int nRecoMuonValidMuonHits = 0;
174  int nRecoMuonChambers = 0;
175  int nRecoMuonChambersCSCorDT = 0;
176  int nRecoMuonMatches = 0;
177  int nRecoMuonMatchedStations = 0;
178  unsigned int nRecoMuonExpectedMatchedStations = 0;
179  unsigned int recoMuonStationMask = 0;
180  int nRecoMuonMatchedRPCLayers = 0;
181  unsigned int recoMuonRPClayerMask = 0;
182  for (auto const& recoMu : *MuonCollection) {
183  float dR2 = deltaR2(muon.eta(), muon.phi(), recoMu.eta(), recoMu.phi());
184  float dPt = std::abs(muon.pt() - recoMu.pt());
185  if (muon.pt() != 0)
186  dPt = dPt / muon.pt();
187 
188  if (dR2 < 0.02 * 0.02 and dPt < 0.001) {
189  if (recoMu.globalTrack().isNonnull()) {
190  normalizedChi2 = recoMu.globalTrack()->normalizedChi2();
191  nRecoMuonValidMuonHits = recoMu.globalTrack()->hitPattern().numberOfValidMuonHits();
192  }
193  recoMuonType = recoMu.type();
194  nRecoMuonChambers = recoMu.numberOfChambers();
195  nRecoMuonChambersCSCorDT = recoMu.numberOfChambersCSCorDT();
196  nRecoMuonMatches = recoMu.numberOfMatches();
197  nRecoMuonMatchedStations = recoMu.numberOfMatchedStations();
198  nRecoMuonExpectedMatchedStations = recoMu.expectedNnumberOfMatchedStations();
199  recoMuonStationMask = recoMu.stationMask();
200  nRecoMuonMatchedRPCLayers = recoMu.numberOfMatchedRPCLayers();
201  recoMuonRPClayerMask = recoMu.RPClayerMask();
202  }
203  }
204 
205  if (muon.pt() < muonPtCut)
206  continue;
207  if (fabs(muon.eta()) > muonEtaCut)
208  continue;
209 
210  double ecalisopf = -1.0;
211  if (EcalPFClusterIsoMap.isValid()) {
212  ecalisopf = (*EcalPFClusterIsoMap)[muonRef];
213  }
214 
215  double hcalisopf = -1.0;
216  if (HcalPFClusterIsoMap.isValid()) {
217  hcalisopf = (*HcalPFClusterIsoMap)[muonRef];
218  }
219 
220  double trackiso = -1.0;
221  if (TrackIsoMap.isValid()) {
222  trackiso = (*TrackIsoMap)[muonRef];
223  }
224 
225  vtxInd.reserve(vtxMuPair.size());
226  for (unsigned int i = 0; i < vtxMuPair.size(); i++) {
227  float dr2_1 = reco::deltaR2(((vtxMuPair[i]).first), muon);
228  float dr2_2 = reco::deltaR2(((vtxMuPair[i]).second), muon);
229  if ((dr2_1 < minDR2) || (dr2_2 < minDR2))
230  vtxInd.push_back(i);
231  }
232 
233  outMuons->emplace_back(muon.pt(),
234  muon.eta(),
235  muon.phi(),
236  muon.mass(),
237  recoMuonType,
238  track->charge(),
239  normalizedChi2,
240  ecalisopf,
241  hcalisopf,
242  trackiso,
243  validStandAloneMuonHits,
244  matchedStandAloneMuonStations,
245  nRecoMuonValidMuonHits,
246  nRecoMuonChambers,
247  nRecoMuonChambersCSCorDT,
248  nRecoMuonMatches,
249  nRecoMuonMatchedStations,
250  nRecoMuonExpectedMatchedStations,
251  recoMuonStationMask,
252  nRecoMuonMatchedRPCLayers,
253  recoMuonRPClayerMask,
254  track->hitPattern().numberOfValidPixelHits(),
255  track->hitPattern().numberOfValidStripHits(),
256  track->hitPattern().pixelLayersWithMeasurement(),
257  track->hitPattern().trackerLayersWithMeasurement(),
258  track->chi2(),
259  track->ndof(),
260  track->dxy(),
261  track->dz(),
262  track->qoverp(),
263  track->lambda(),
264  track->pt(),
265  track->phi(),
266  track->eta(),
267  track->dxyError(),
268  track->dzError(),
269  track->qoverpError(),
270  track->lambdaError(),
271  track->phiError(),
272  track->dsz(),
273  track->dszError(),
274  track->covariance(0, 1),
275  track->covariance(0, 2),
276  track->covariance(0, 3),
277  track->covariance(0, 4),
278  track->covariance(1, 2),
279  track->covariance(1, 3),
280  track->covariance(1, 4),
281  track->covariance(2, 3),
282  track->covariance(2, 4),
283  track->covariance(3, 4),
284  track->vx(),
285  track->vy(),
286  track->vz(),
287  vtxInd,
288  track->hitPattern().run3ScoutingHitPatternPOD());
289  vtxInd.clear();
290  }
291 
292  // Put output
293  iEvent.put(std::move(outMuons));
294  iEvent.put(std::move(dispVertices), "displacedVtx");
295 }
HLT_2022v15_cff.TrackIsoMap
TrackIsoMap
Definition: HLT_2022v15_cff.py:60258
edm::getRef
helper::MatcherGetRef< C >::ref_type getRef(const Handle< C > &c, size_t k)
Definition: getRef.h:28
HLT_2022v15_cff.displacedvertexCollection
displacedvertexCollection
Definition: HLT_2022v15_cff.py:60252
mps_fire.i
i
Definition: mps_fire.py:429
Run3ScoutingVertexCollection
std::vector< Run3ScoutingVertex > Run3ScoutingVertexCollection
Definition: Run3ScoutingVertex.h:69
HLT_2022v15_cff.EcalPFClusterIsoMap
EcalPFClusterIsoMap
Definition: HLT_2022v15_cff.py:60256
HLTScoutingMuonProducer::MuonCollection_
const edm::EDGetTokenT< reco::MuonCollection > MuonCollection_
Definition: HLTScoutingMuonProducer.h:65
HLTScoutingMuonProducer::EcalPFClusterIsoMap_
const edm::EDGetTokenT< RecoChargedCandMap > EcalPFClusterIsoMap_
Definition: HLTScoutingMuonProducer.h:68
HLTScoutingMuonProducer::HcalPFClusterIsoMap_
const edm::EDGetTokenT< RecoChargedCandMap > HcalPFClusterIsoMap_
Definition: HLTScoutingMuonProducer.h:69
HLT_2022v15_cff.HcalPFClusterIsoMap
HcalPFClusterIsoMap
Definition: HLT_2022v15_cff.py:60257
TrackCollections2monitor_cff.normalizedChi2
normalizedChi2
Definition: TrackCollections2monitor_cff.py:252
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
reco::MuonCollection
std::vector< Muon > MuonCollection
collection of Muon objects
Definition: MuonFwd.h:9
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637
HLTMuonOfflineAnalyzer_cfi.deltaR2
deltaR2
Definition: HLTMuonOfflineAnalyzer_cfi.py:53
iEvent
int iEvent
Definition: GenABIO.cc:224
HLTScoutingMuonProducer::TrackIsoMap_
const edm::EDGetTokenT< edm::ValueMap< double > > TrackIsoMap_
Definition: HLTScoutingMuonProducer.h:70
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:537
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:649
HLTScoutingMuonProducer::displacedvertexCollection_
const edm::EDGetTokenT< reco::VertexCollection > displacedvertexCollection_
Definition: HLTScoutingMuonProducer.h:64
edm::Ref::isNull
bool isNull() const
Checks for null.
Definition: Ref.h:235
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
reco::deltaR2
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16
HLTScoutingMuonProducer::ChargedCandidateCollection_
const edm::EDGetTokenT< reco::RecoChargedCandidateCollection > ChargedCandidateCollection_
Definition: HLTScoutingMuonProducer.h:63
Run3ScoutingMuonCollection
std::vector< Run3ScoutingMuon > Run3ScoutingMuonCollection
Definition: Run3ScoutingMuon.h:302
HLTScoutingMuonProducer::linkToken_
const edm::EDGetTokenT< reco::MuonTrackLinksCollection > linkToken_
Definition: HLTScoutingMuonProducer.h:66
edm
HLT enums.
Definition: AlignableModifier.h:19
HLT_2022v15_cff.track
track
Definition: HLT_2022v15_cff.py:9644
eostools.move
def move(src, dest)
Definition: eostools.py:511

Member Data Documentation

◆ ChargedCandidateCollection_

const edm::EDGetTokenT<reco::RecoChargedCandidateCollection> HLTScoutingMuonProducer::ChargedCandidateCollection_
private

Definition at line 63 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ displacedvertexCollection_

const edm::EDGetTokenT<reco::VertexCollection> HLTScoutingMuonProducer::displacedvertexCollection_
private

Definition at line 64 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ EcalPFClusterIsoMap_

const edm::EDGetTokenT<RecoChargedCandMap> HLTScoutingMuonProducer::EcalPFClusterIsoMap_
private

Definition at line 68 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ HcalPFClusterIsoMap_

const edm::EDGetTokenT<RecoChargedCandMap> HLTScoutingMuonProducer::HcalPFClusterIsoMap_
private

Definition at line 69 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ linkToken_

const edm::EDGetTokenT<reco::MuonTrackLinksCollection> HLTScoutingMuonProducer::linkToken_
private

Definition at line 66 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ minVtxProbCut

const double HLTScoutingMuonProducer::minVtxProbCut
private

Definition at line 74 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ MuonCollection_

const edm::EDGetTokenT<reco::MuonCollection> HLTScoutingMuonProducer::MuonCollection_
private

Definition at line 65 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ muonEtaCut

const double HLTScoutingMuonProducer::muonEtaCut
private

Definition at line 73 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ muonPtCut

const double HLTScoutingMuonProducer::muonPtCut
private

Definition at line 72 of file HLTScoutingMuonProducer.h.

Referenced by produce().

◆ TrackCollection_

const edm::EDGetTokenT<reco::TrackCollection> HLTScoutingMuonProducer::TrackCollection_
private

Definition at line 67 of file HLTScoutingMuonProducer.h.

◆ TrackIsoMap_

const edm::EDGetTokenT<edm::ValueMap<double> > HLTScoutingMuonProducer::TrackIsoMap_
private

Definition at line 70 of file HLTScoutingMuonProducer.h.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.14