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
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > 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
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices 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 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:307
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:159
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 301 of file HLTScoutingMuonProducer.cc.

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

301  {
302  edm::ParameterSetDescription desc;
303  desc.add<edm::InputTag>("ChargedCandidates", edm::InputTag("hltIterL3MuonCandidatesNoVtx"));
304  desc.add<edm::InputTag>("displacedvertexCollection", edm::InputTag("hltDisplacedmumuVtxProducer"));
305  desc.add<edm::InputTag>("InputMuons", edm::InputTag("hltIterL3MuonsNoVtx"));
306  desc.add<edm::InputTag>("InputLinks", edm::InputTag("hltL3MuonsIterL3LinksNoVtx"));
307  desc.add<edm::InputTag>("Tracks", edm::InputTag("hltPixelTracks"));
308  desc.add<edm::InputTag>("EcalPFClusterIsoMap", edm::InputTag("hltMuonEcalMFPFClusterIsoForMuonsNoVtx"));
309  desc.add<edm::InputTag>("HcalPFClusterIsoMap", edm::InputTag("hltMuonHcalPFClusterIsoForMuonsNoVtx"));
310  desc.add<edm::InputTag>("TrackIsoMap",
311  edm::InputTag("hltMuonTkRelIsolationCut0p09MapNoVtx:combinedRelativeIsoDeposits"));
312  desc.add<double>("muonPtCut", 3.0);
313  desc.add<double>("muonEtaCut", 2.4);
314  desc.add<double>("minVtxProbCut", 0.001);
315 
316  descriptions.add("hltScoutingMuonProducer", desc);
317 }
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
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_2023v12_cff::displacedvertexCollection, displacedvertexCollection_, ALPAKA_ACCELERATOR_NAMESPACE::dR2(), MillePedeFileConverter_cfg::e, HLT_2023v12_cff::EcalPFClusterIsoMap, EcalPFClusterIsoMap_, reco::TrackBase::eta(), dqmdumpme::first, edm::getRef(), HLT_2023v12_cff::HcalPFClusterIsoMap, HcalPFClusterIsoMap_, mps_fire::i, iEvent, edm::Ref< C, T, F >::isAvailable(), edm::Ref< C, T, F >::isNull(), electronStore::links, linkToken_, minVtxProbCut, eostools::move(), MuonCollection_, muonEtaCut, muonPtCut, TrackCollections2monitor_cff::normalizedChi2, run3scouting_cff::nRecoMuonChambers, run3scouting_cff::nRecoMuonChambersCSCorDT, run3scouting_cff::nRecoMuonExpectedMatchedStations, run3scouting_cff::nRecoMuonMatchedRPCLayers, run3scouting_cff::nRecoMuonMatchedStations, run3scouting_cff::nRecoMuonMatches, reco::TrackBase::phi(), reco::TrackBase::pt(), run3scouting_cff::recoMuonRPClayerMask, run3scouting_cff::recoMuonStationMask, edm::second(), HLT_2023v12_cff::track, HLT_2023v12_cff::TrackIsoMap, TrackIsoMap_, and run3scouting_cff::vtxInd.

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  dispvtx.covariance(0, 1),
109  dispvtx.covariance(0, 2),
110  dispvtx.covariance(1, 2));
111  }
112  }
113 
114  // Get Muon collection
115  Handle<reco::MuonCollection> MuonCollection;
116  if (!iEvent.getByToken(MuonCollection_, MuonCollection)) {
117  iEvent.put(std::move(outMuons));
118  iEvent.put(std::move(dispVertices), "displacedVtx");
119  return;
120  }
121 
122  // Get Links collection
123  edm::Handle<reco::MuonTrackLinksCollection> links;
124  if (!iEvent.getByToken(linkToken_, links)) {
125  iEvent.put(std::move(outMuons));
126  iEvent.put(std::move(dispVertices), "displacedVtx");
127  return;
128  }
129 
130  // Get EcalPFClusterIsoMap
131  Handle<RecoChargedCandMap> EcalPFClusterIsoMap;
132  iEvent.getByToken(EcalPFClusterIsoMap_, EcalPFClusterIsoMap);
133 
134  // Get HcalPFClusterIsoMap
135  Handle<RecoChargedCandMap> HcalPFClusterIsoMap;
136  iEvent.getByToken(HcalPFClusterIsoMap_, HcalPFClusterIsoMap);
137 
138  // Get TrackIsoMap
139  Handle<ValueMap<double>> TrackIsoMap;
140  iEvent.getByToken(TrackIsoMap_, TrackIsoMap);
141 
142  // Produce muons
143  std::vector<int> vtxInd;
144  float minDR2 = 1e-06;
145  int index = 0;
146 
147  for (auto& muon : *ChargedCandidateCollection) {
148  reco::RecoChargedCandidateRef muonRef = getRef(ChargedCandidateCollection, index);
149  ++index;
150  if (muonRef.isNull() || !muonRef.isAvailable())
151  continue;
152 
153  reco::TrackRef track = muon.track();
154  if (track.isNull() || !track.isAvailable())
155  continue;
156 
157  int validStandAloneMuonHits = 0;
158  int matchedStandAloneMuonStations = 0;
159  for (auto const& link : *links) {
160  const reco::Track& trackerTrack = *link.trackerTrack();
161  float dR2 = deltaR2(track->eta(), track->phi(), trackerTrack.eta(), trackerTrack.phi());
162  float dPt = std::abs(track->pt() - trackerTrack.pt());
163  if (track->pt() != 0)
164  dPt = dPt / track->pt();
165 
166  if (dR2 < 0.02 * 0.02 and dPt < 0.001) {
167  if (link.standAloneTrack().isNonnull()) {
168  validStandAloneMuonHits = link.standAloneTrack()->hitPattern().numberOfValidMuonHits();
169  matchedStandAloneMuonStations = link.standAloneTrack()->hitPattern().muonStationsWithValidHits();
170  }
171  }
172  }
173 
174  unsigned int recoMuonType = 2; // Global muon
175  float normalizedChi2 = 999.0;
176  int nRecoMuonValidMuonHits = 0;
177  int nRecoMuonChambers = 0;
178  int nRecoMuonChambersCSCorDT = 0;
179  int nRecoMuonMatches = 0;
180  int nRecoMuonMatchedStations = 0;
181  unsigned int nRecoMuonExpectedMatchedStations = 0;
182  unsigned int recoMuonStationMask = 0;
183  int nRecoMuonMatchedRPCLayers = 0;
184  unsigned int recoMuonRPClayerMask = 0;
185  for (auto const& recoMu : *MuonCollection) {
186  float dR2 = deltaR2(muon.eta(), muon.phi(), recoMu.eta(), recoMu.phi());
187  float dPt = std::abs(muon.pt() - recoMu.pt());
188  if (muon.pt() != 0)
189  dPt = dPt / muon.pt();
190 
191  if (dR2 < 0.02 * 0.02 and dPt < 0.001) {
192  if (recoMu.globalTrack().isNonnull()) {
193  normalizedChi2 = recoMu.globalTrack()->normalizedChi2();
194  nRecoMuonValidMuonHits = recoMu.globalTrack()->hitPattern().numberOfValidMuonHits();
195  }
196  recoMuonType = recoMu.type();
197  nRecoMuonChambers = recoMu.numberOfChambers();
198  nRecoMuonChambersCSCorDT = recoMu.numberOfChambersCSCorDT();
199  nRecoMuonMatches = recoMu.numberOfMatches();
200  nRecoMuonMatchedStations = recoMu.numberOfMatchedStations();
201  nRecoMuonExpectedMatchedStations = recoMu.expectedNnumberOfMatchedStations();
202  recoMuonStationMask = recoMu.stationMask();
203  nRecoMuonMatchedRPCLayers = recoMu.numberOfMatchedRPCLayers();
204  recoMuonRPClayerMask = recoMu.RPClayerMask();
205  }
206  }
207 
208  if (muon.pt() < muonPtCut)
209  continue;
210  if (fabs(muon.eta()) > muonEtaCut)
211  continue;
212 
213  double ecalisopf = -1.0;
214  if (EcalPFClusterIsoMap.isValid()) {
215  ecalisopf = (*EcalPFClusterIsoMap)[muonRef];
216  }
217 
218  double hcalisopf = -1.0;
219  if (HcalPFClusterIsoMap.isValid()) {
220  hcalisopf = (*HcalPFClusterIsoMap)[muonRef];
221  }
222 
223  double trackiso = -1.0;
224  if (TrackIsoMap.isValid()) {
225  trackiso = (*TrackIsoMap)[muonRef];
226  }
227 
228  vtxInd.reserve(vtxMuPair.size());
229  for (unsigned int i = 0; i < vtxMuPair.size(); i++) {
230  float dr2_1 = reco::deltaR2(((vtxMuPair[i]).first), muon);
231  float dr2_2 = reco::deltaR2(((vtxMuPair[i]).second), muon);
232  if ((dr2_1 < minDR2) || (dr2_2 < minDR2))
233  vtxInd.push_back(i);
234  }
235 
236  outMuons->emplace_back(muon.pt(),
237  muon.eta(),
238  muon.phi(),
239  muon.mass(),
240  recoMuonType,
241  track->charge(),
242  normalizedChi2,
243  ecalisopf,
244  hcalisopf,
245  trackiso,
246  validStandAloneMuonHits,
247  matchedStandAloneMuonStations,
248  nRecoMuonValidMuonHits,
249  nRecoMuonChambers,
250  nRecoMuonChambersCSCorDT,
251  nRecoMuonMatches,
252  nRecoMuonMatchedStations,
253  nRecoMuonExpectedMatchedStations,
254  recoMuonStationMask,
255  nRecoMuonMatchedRPCLayers,
256  recoMuonRPClayerMask,
257  track->hitPattern().numberOfValidPixelHits(),
258  track->hitPattern().numberOfValidStripHits(),
259  track->hitPattern().pixelLayersWithMeasurement(),
260  track->hitPattern().trackerLayersWithMeasurement(),
261  track->chi2(),
262  track->ndof(),
263  track->dxy(),
264  track->dz(),
265  track->qoverp(),
266  track->lambda(),
267  track->pt(),
268  track->phi(),
269  track->eta(),
270  track->dxyError(),
271  track->dzError(),
272  track->qoverpError(),
273  track->lambdaError(),
274  track->phiError(),
275  track->dsz(),
276  track->dszError(),
277  track->covariance(0, 1),
278  track->covariance(0, 2),
279  track->covariance(0, 3),
280  track->covariance(0, 4),
281  track->covariance(1, 2),
282  track->covariance(1, 3),
283  track->covariance(1, 4),
284  track->covariance(2, 3),
285  track->covariance(2, 4),
286  track->covariance(3, 4),
287  track->vx(),
288  track->vy(),
289  track->vz(),
290  vtxInd,
291  track->hitPattern().run3ScoutingHitPatternPOD());
292  vtxInd.clear();
293  }
294 
295  // Put output
296  iEvent.put(std::move(outMuons));
297  iEvent.put(std::move(dispVertices), "displacedVtx");
298 }
edm::getRef
helper::MatcherGetRef< C >::ref_type getRef(const Handle< C > &c, size_t k)
Definition: getRef.h:28
run3scouting_cff.nRecoMuonChambersCSCorDT
nRecoMuonChambersCSCorDT
Definition: run3scouting_cff.py:81
mps_fire.i
i
Definition: mps_fire.py:429
Run3ScoutingVertexCollection
std::vector< Run3ScoutingVertex > Run3ScoutingVertexCollection
Definition: Run3ScoutingVertex.h:84
HLTScoutingMuonProducer::MuonCollection_
const edm::EDGetTokenT< reco::MuonCollection > MuonCollection_
Definition: HLTScoutingMuonProducer.h:65
HLT_2023v12_cff.TrackIsoMap
TrackIsoMap
Definition: HLT_2023v12_cff.py:17580
HLTScoutingMuonProducer::EcalPFClusterIsoMap_
const edm::EDGetTokenT< RecoChargedCandMap > EcalPFClusterIsoMap_
Definition: HLTScoutingMuonProducer.h:68
run3scouting_cff.nRecoMuonExpectedMatchedStations
nRecoMuonExpectedMatchedStations
Definition: run3scouting_cff.py:84
HLTScoutingMuonProducer::HcalPFClusterIsoMap_
const edm::EDGetTokenT< RecoChargedCandMap > HcalPFClusterIsoMap_
Definition: HLTScoutingMuonProducer.h:69
run3scouting_cff.nRecoMuonMatchedStations
nRecoMuonMatchedStations
Definition: run3scouting_cff.py:83
HLT_2023v12_cff.HcalPFClusterIsoMap
HcalPFClusterIsoMap
Definition: HLT_2023v12_cff.py:17579
TrackCollections2monitor_cff.normalizedChi2
normalizedChi2
Definition: TrackCollections2monitor_cff.py:266
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
run3scouting_cff.nRecoMuonMatchedRPCLayers
nRecoMuonMatchedRPCLayers
Definition: run3scouting_cff.py:86
run3scouting_cff.nRecoMuonMatches
nRecoMuonMatches
Definition: run3scouting_cff.py:82
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
ALPAKA_ACCELERATOR_NAMESPACE::dR2
ALPAKA_FN_ACC static ALPAKA_FN_INLINE float dR2(Position4 pos1, Position4 pos2)
Definition: PFClusterSoAProducerKernel.dev.cc:31
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:541
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
run3scouting_cff.recoMuonStationMask
recoMuonStationMask
Definition: run3scouting_cff.py:85
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
run3scouting_cff.recoMuonRPClayerMask
recoMuonRPClayerMask
Definition: run3scouting_cff.py:87
edm
HLT enums.
Definition: AlignableModifier.h:19
run3scouting_cff.vtxInd
vtxInd
Definition: run3scouting_cff.py:162
HLT_2023v12_cff.EcalPFClusterIsoMap
EcalPFClusterIsoMap
Definition: HLT_2023v12_cff.py:17578
run3scouting_cff.nRecoMuonChambers
nRecoMuonChambers
Definition: run3scouting_cff.py:80
eostools.move
def move(src, dest)
Definition: eostools.py:511
HLT_2023v12_cff.displacedvertexCollection
displacedvertexCollection
Definition: HLT_2023v12_cff.py:17574
HLT_2023v12_cff.track
track
Definition: HLT_2023v12_cff.py:9627

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