CMS 3D CMS Logo
List of all members | Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes
L1TMuonDQMOffline Class Reference

#include <L1TMuonDQMOffline.h>

Inheritance diagram for L1TMuonDQMOffline:
DQMEDAnalyzer edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > > edm::stream::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

 L1TMuonDQMOffline (const edm::ParameterSet &ps)
 
virtual ~L1TMuonDQMOffline ()
 
- Public Member Functions inherited from DQMEDAnalyzer
virtual void beginRun (edm::Run const &, edm::EventSetup const &) final
 
virtual void beginStream (edm::StreamID id) final
 
virtual void bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &)=0
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer (void)
 
virtual void endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
virtual void endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
uint32_t streamId () const
 
- Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
 EDAnalyzer ()=default
 
- Public Member Functions inherited from edm::stream::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDAnalyzerBase ()
 
- 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
 
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::vector< ModuleDescription const * > &modules, 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 updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Protected Member Functions

virtual void beginLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c) override
 
virtual void bookControlHistos (DQMStore::IBooker &)
 
virtual void bookEfficiencyHistos (DQMStore::IBooker &ibooker, int ptCut)
 
virtual void bookHistograms (DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &iSetup) override
 
virtual void dqmBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override
 
virtual void dqmEndLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
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 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)
 

Private Member Functions

void analyze (const edm::Event &e, const edm::EventSetup &c) override
 
void getMuonGmtPairs (edm::Handle< l1t::MuonBxCollection > &gmtCands)
 
const reco::Vertex getPrimaryVertex (edm::Handle< reco::VertexCollection > &vertex, edm::Handle< reco::BeamSpot > &beamSpot)
 
void getProbeMuons (edm::Handle< edm::TriggerResults > &trigResults, edm::Handle< trigger::TriggerEvent > &trigEvent)
 
void getTightMuons (edm::Handle< reco::MuonCollection > &muons, const reco::Vertex &vertex)
 
bool matchHlt (edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon *mu)
 

Private Attributes

edm::ESHandle< MagneticField > m_BField
 
edm::EDGetTokenT< reco::BeamSpot > m_BsInputTag
 
std::map< std::string, MonitorElement * > m_ControlHistos
 
std::map< int, std::map< std::string, MonitorElement * > > m_EfficiencyHistos
 
std::vector< double > m_effVsEtaBins
 
std::vector< double > m_effVsPhiBins
 
std::vector< double > m_effVsPtBins
 
edm::EDGetTokenT< l1t::MuonBxCollection > m_GmtInputTag
 
std::vector< int > m_GmtPtCuts
 
std::string m_HistFolder
 
HLTConfigProvider m_hltConfig
 
BXVector< l1t::Muon > m_L1tL1tMuons
 
std::vector< l1t::MuonBxCollection > m_L1tMuons
 
float m_MaxGmtMuonDR
 
float m_MaxHltMuonDR
 
float m_MaxMuonEta
 
std::vector< MuonGmtPair > m_MuonGmtPairs
 
edm::EDGetTokenT< reco::MuonCollection > m_MuonInputTag
 
std::vector< const reco::Muon * > m_ProbeMuons
 
edm::ESHandle< Propagator > m_propagatorAlong
 
edm::ESHandle< Propagator > m_propagatorOpposite
 
std::vector< reco::MuonCollection > m_RecoMuons
 
std::vector< reco::Muon > m_RecoRecoMuons
 
std::vector< const reco::Muon * > m_TightMuons
 
std::vector< int > m_trigIndices
 
edm::EDGetTokenT< trigger::TriggerEvent > m_trigInputTag
 
std::vector< std::string > m_trigNames
 
std::string m_trigProcess
 
edm::EDGetTokenT< edm::TriggerResults > m_trigProcess_token
 
bool m_verbose
 
edm::EDGetTokenT< reco::VertexCollection > m_VtxInputTag
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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
 
- Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache > globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr< dqmDetails::NoCache > globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
 
static void globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
 
static void globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
 
- Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Detailed Description

Definition at line 113 of file L1TMuonDQMOffline.h.

Constructor & Destructor Documentation

L1TMuonDQMOffline::L1TMuonDQMOffline ( const edm::ParameterSet &  ps)

Definition at line 106 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, m_MaxGmtMuonDR, m_MaxHltMuonDR, m_MaxMuonEta, and m_verbose.

106  :
107  m_verbose(ps.getUntrackedParameter<bool>("verbose")),
108  m_HistFolder(ps.getUntrackedParameter<string>("histFolder")),
109  m_GmtPtCuts(ps.getUntrackedParameter< vector<int> >("gmtPtCuts")),
110  m_MuonInputTag(consumes<reco::MuonCollection>(ps.getUntrackedParameter<InputTag>("muonInputTag"))),
111  m_GmtInputTag(consumes<l1t::MuonBxCollection>(ps.getUntrackedParameter<InputTag>("gmtInputTag"))),
112  m_VtxInputTag(consumes<VertexCollection>(ps.getUntrackedParameter<InputTag>("vtxInputTag"))),
113  m_BsInputTag(consumes<BeamSpot>(ps.getUntrackedParameter<InputTag>("bsInputTag"))),
114  m_trigInputTag(consumes<trigger::TriggerEvent>(ps.getUntrackedParameter<InputTag>("trigInputTag"))),
115  m_trigProcess(ps.getUntrackedParameter<string>("trigProcess")),
116  m_trigProcess_token(consumes<edm::TriggerResults>(ps.getUntrackedParameter<InputTag>("trigProcess_token"))),
117  m_trigNames(ps.getUntrackedParameter<vector<string> >("triggerNames")),
118  m_effVsPtBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPtBins")),
119  m_effVsPhiBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsPhiBins")),
120  m_effVsEtaBins(ps.getUntrackedParameter<std::vector<double>>("efficiencyVsEtaBins"))
121 {
122  if (m_verbose) cout << "[L1TMuonDQMOffline:] ____________ Storage initialization ____________ " << endl;
123  // CB do we need them from cfi?
124  m_MaxMuonEta = 2.4;
125  m_MaxGmtMuonDR = 0.7;
126  m_MaxHltMuonDR = 0.1;
127  // CB ignored at present
128  //m_MinMuonDR = 1.2;
129 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
L1TMuonDQMOffline::m_BsInputTag
edm::EDGetTokenT< reco::BeamSpot > m_BsInputTag
Definition: L1TMuonDQMOffline.h:168
L1TMuonDQMOffline::m_MuonInputTag
edm::EDGetTokenT< reco::MuonCollection > m_MuonInputTag
Definition: L1TMuonDQMOffline.h:165
L1TMuonDQMOffline::m_trigProcess_token
edm::EDGetTokenT< edm::TriggerResults > m_trigProcess_token
Definition: L1TMuonDQMOffline.h:171
L1TMuonDQMOffline::m_effVsPhiBins
std::vector< double > m_effVsPhiBins
Definition: L1TMuonDQMOffline.h:174
L1TMuonDQMOffline::m_HistFolder
std::string m_HistFolder
Definition: L1TMuonDQMOffline.h:163
L1TMuonDQMOffline::m_trigNames
std::vector< std::string > m_trigNames
Definition: L1TMuonDQMOffline.h:172
L1TMuonDQMOffline::m_trigProcess
std::string m_trigProcess
Definition: L1TMuonDQMOffline.h:170
L1TMuonDQMOffline::m_GmtInputTag
edm::EDGetTokenT< l1t::MuonBxCollection > m_GmtInputTag
Definition: L1TMuonDQMOffline.h:166
L1TMuonDQMOffline::m_effVsPtBins
std::vector< double > m_effVsPtBins
Definition: L1TMuonDQMOffline.h:173
L1TMuonDQMOffline::m_GmtPtCuts
std::vector< int > m_GmtPtCuts
Definition: L1TMuonDQMOffline.h:164
L1TMuonDQMOffline::m_effVsEtaBins
std::vector< double > m_effVsEtaBins
Definition: L1TMuonDQMOffline.h:175
L1TMuonDQMOffline::m_VtxInputTag
edm::EDGetTokenT< reco::VertexCollection > m_VtxInputTag
Definition: L1TMuonDQMOffline.h:167
L1TMuonDQMOffline::m_trigInputTag
edm::EDGetTokenT< trigger::TriggerEvent > m_trigInputTag
Definition: L1TMuonDQMOffline.h:169
gather_cfg.cout
cout
Definition: gather_cfg.py:145
L1TMuonDQMOffline::~L1TMuonDQMOffline ( )
virtual

Definition at line 132 of file L1TMuonDQMOffline.cc.

132 { }

Member Function Documentation

void L1TMuonDQMOffline::analyze ( const edm::Event &  e,
const edm::EventSetup &  c 
)
overrideprivate

Definition at line 181 of file L1TMuonDQMOffline.cc.

References ecalDrivenElectronSeedsParameters_cff::beamSpot, BXVector< T >::begin(), gather_cfg::cout, BXVector< T >::end(), PVValHelper::eta, edm::EventSetup::get(), edm::Event::getByToken(), getMuonGmtPairs(), getPrimaryVertex(), getProbeMuons(), getTightMuons(), m_BField, m_BsInputTag, m_EfficiencyHistos, m_GmtInputTag, m_GmtPtCuts, m_MaxMuonEta, m_MuonGmtPairs, m_MuonInputTag, m_propagatorAlong, m_propagatorOpposite, m_trigInputTag, m_trigProcess_token, m_verbose, m_VtxInputTag, RPCpg::mu, electronCleaner_cfi::muons, phi, impactParameterTagInfos_cfi::primaryVertex, and EnergyCorrector::pt.

181  {
182 
183  Handle<reco::MuonCollection> muons;
184  iEvent.getByToken(m_MuonInputTag, muons);
185 
186  Handle<BeamSpot> beamSpot;
187  iEvent.getByToken(m_BsInputTag, beamSpot);
188 
189  Handle<VertexCollection> vertex;
190  iEvent.getByToken(m_VtxInputTag, vertex);
191 
192  Handle<l1t::MuonBxCollection> gmtCands;
193  iEvent.getByToken(m_GmtInputTag,gmtCands);
194 
195  Handle<edm::TriggerResults> trigResults;
196  iEvent.getByToken(m_trigProcess_token,trigResults);
197 
198  edm::Handle<trigger::TriggerEvent> trigEvent;
199  iEvent.getByToken(m_trigInputTag,trigEvent);
200 
201  eventSetup.get<IdealMagneticFieldRecord>().get(m_BField);
202  eventSetup.get<TrackingComponentsRecord>().get("PropagatorWithMaterial",m_propagatorAlong);
203  eventSetup.get<TrackingComponentsRecord>().get("PropagatorWithMaterialOpposite",m_propagatorOpposite);
204 
205  const Vertex primaryVertex = getPrimaryVertex(vertex,beamSpot);
206 
207  getTightMuons(muons,primaryVertex);
208  getProbeMuons(trigResults,trigEvent); // CB add flag to run on orthogonal datasets (no T&P)
209  getMuonGmtPairs(gmtCands);
210 
211 // MuonCollection::const_iterator muonIt = muons->begin();
212 // MuonCollection::const_iterator muonEnd = muons->end();
213 
214  vector<l1t::Muon> gmtContainer;// = gmtCands->getRecord(0).getGMTCands();
215 
216  for (auto mu = gmtCands->begin(0); mu != gmtCands->end(0); ++mu) {
217  gmtContainer.push_back(*mu);
218  }
219 
220 // vector<l1t::Muon>::const_iterator gmtIt = gmtContainer.begin();
221 // vector<l1t::Muon>::const_iterator gmtEnd = gmtContainer.end();
222 
223  if (m_verbose) cout << "[L1TMuonDQMOffline:] Computing efficiencies" << endl;
224 
225  vector<MuonGmtPair>::const_iterator muonGmtPairsIt = m_MuonGmtPairs.begin();
226  vector<MuonGmtPair>::const_iterator muonGmtPairsEnd = m_MuonGmtPairs.end();
227 
228  for(; muonGmtPairsIt!=muonGmtPairsEnd; ++muonGmtPairsIt) {
229  float eta = muonGmtPairsIt->eta();
230  float phi = muonGmtPairsIt->phi();
231  float pt = muonGmtPairsIt->pt();
232 
233  float gmtPt = muonGmtPairsIt->gmtPt();
234  int qual = muonGmtPairsIt->gmtQual();
235 
236  vector<int>::const_iterator gmtPtCutsIt = m_GmtPtCuts.begin();
237  vector<int>::const_iterator gmtPtCutsEnd = m_GmtPtCuts.end();
238 
239  for (; gmtPtCutsIt!=gmtPtCutsEnd; ++ gmtPtCutsIt) {
240  int gmtPtCut = (*gmtPtCutsIt);
241  bool gmtAboveCut = (gmtPt > gmtPtCut);
242 
243  stringstream ptCutToTag; ptCutToTag << gmtPtCut;
244  string ptTag = ptCutToTag.str();
245 
246  if (fabs(eta) < m_MaxMuonEta) {
247  m_EfficiencyHistos[gmtPtCut]["EffvsPt_" + ptTag + "_Den"]->Fill(pt);
248  m_EfficiencyHistos[gmtPtCut]["EffvsPt_OPEN_" + ptTag + "_Den"]->Fill(pt);
249  m_EfficiencyHistos[gmtPtCut]["EffvsPt_DOUBLE_" + ptTag + "_Den"]->Fill(pt);
250  m_EfficiencyHistos[gmtPtCut]["EffvsPt_SINGLE_" + ptTag + "_Den"]->Fill(pt);
251 
252  if (gmtAboveCut) {
253  m_EfficiencyHistos[gmtPtCut]["EffvsPt_" + ptTag + "_Num"]->Fill(pt);
254 
255  if (qual >= 4) m_EfficiencyHistos[gmtPtCut]["EffvsPt_OPEN_" + ptTag + "_Num"]->Fill(pt);
256  if (qual >= 8) m_EfficiencyHistos[gmtPtCut]["EffvsPt_DOUBLE_" + ptTag + "_Num"]->Fill(pt);
257  if (qual >= 12) m_EfficiencyHistos[gmtPtCut]["EffvsPt_SINGLE_" + ptTag + "_Num"]->Fill(pt);
258  }
259 
260  // efficiency in eta/phi at plateau
261  if (pt > 1.25*gmtPtCut) { // efficiency in eta/phi at plateau
262 
263  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_" + ptTag + "_Den"]->Fill(phi);
264  m_EfficiencyHistos[gmtPtCut]["EffvsEta_" + ptTag + "_Den"]->Fill(eta);
265 
266  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_OPEN_" + ptTag + "_Den"]->Fill(phi);
267  m_EfficiencyHistos[gmtPtCut]["EffvsEta_OPEN_" + ptTag + "_Den"]->Fill(eta);
268 
269  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_DOUBLE_" + ptTag + "_Den"]->Fill(phi);
270  m_EfficiencyHistos[gmtPtCut]["EffvsEta_DOUBLE_" + ptTag + "_Den"]->Fill(eta);
271 
272  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_SINGLE_" + ptTag + "_Den"]->Fill(phi);
273  m_EfficiencyHistos[gmtPtCut]["EffvsEta_SINGLE_" + ptTag + "_Den"]->Fill(eta);
274 
275  if (gmtAboveCut) {
276  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_" + ptTag + "_Num"]->Fill(phi);
277  m_EfficiencyHistos[gmtPtCut]["EffvsEta_" + ptTag + "_Num"]->Fill(eta);
278 
279  if (qual >= 4) {
280  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_OPEN_" + ptTag + "_Num"]->Fill(phi);
281  m_EfficiencyHistos[gmtPtCut]["EffvsEta_OPEN_" + ptTag + "_Num"]->Fill(eta);
282  }
283  if (qual >= 8) {
284  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_DOUBLE_" + ptTag + "_Num"]->Fill(phi);
285  m_EfficiencyHistos[gmtPtCut]["EffvsEta_DOUBLE_" + ptTag + "_Num"]->Fill(eta);
286  }
287  if (qual >= 12) {
288  m_EfficiencyHistos[gmtPtCut]["EffvsPhi_SINGLE_" + ptTag + "_Num"]->Fill(phi);
289  m_EfficiencyHistos[gmtPtCut]["EffvsEta_SINGLE_" + ptTag + "_Num"]->Fill(eta);
290  }
291  }
292  }
293  }
294  }
295  }
296  if (m_verbose) cout << "[L1TMuonDQMOffline:] Computation finished" << endl;
297 }
BXVector::end
const_iterator end(int bx) const
L1TMuonDQMOffline::m_propagatorOpposite
edm::ESHandle< Propagator > m_propagatorOpposite
Definition: L1TMuonDQMOffline.h:145
L1TMuonDQMOffline::m_BsInputTag
edm::EDGetTokenT< reco::BeamSpot > m_BsInputTag
Definition: L1TMuonDQMOffline.h:168
L1TMuonDQMOffline::m_MuonInputTag
edm::EDGetTokenT< reco::MuonCollection > m_MuonInputTag
Definition: L1TMuonDQMOffline.h:165
L1TMuonDQMOffline::m_trigProcess_token
edm::EDGetTokenT< edm::TriggerResults > m_trigProcess_token
Definition: L1TMuonDQMOffline.h:171
ecalDrivenElectronSeedsParameters_cff.beamSpot
beamSpot
Definition: ecalDrivenElectronSeedsParameters_cff.py:11
L1TMuonDQMOffline::getProbeMuons
void getProbeMuons(edm::Handle< edm::TriggerResults > &trigResults, edm::Handle< trigger::TriggerEvent > &trigEvent)
Definition: L1TMuonDQMOffline.cc:469
iEvent
int iEvent
Definition: GenABIO.cc:230
L1TMuonDQMOffline::m_GmtInputTag
edm::EDGetTokenT< l1t::MuonBxCollection > m_GmtInputTag
Definition: L1TMuonDQMOffline.h:166
L1TMuonDQMOffline::m_EfficiencyHistos
std::map< int, std::map< std::string, MonitorElement * > > m_EfficiencyHistos
Definition: L1TMuonDQMOffline.h:148
RPCpg::mu
const int mu
Definition: Constants.h:22
L1TMuonDQMOffline::m_MuonGmtPairs
std::vector< MuonGmtPair > m_MuonGmtPairs
Definition: L1TMuonDQMOffline.h:154
L1TMuonDQMOffline::m_BField
edm::ESHandle< MagneticField > m_BField
Definition: L1TMuonDQMOffline.h:143
L1TMuonDQMOffline::m_propagatorAlong
edm::ESHandle< Propagator > m_propagatorAlong
Definition: L1TMuonDQMOffline.h:144
L1TMuonDQMOffline::getMuonGmtPairs
void getMuonGmtPairs(edm::Handle< l1t::MuonBxCollection > &gmtCands)
Definition: L1TMuonDQMOffline.cc:507
impactParameterTagInfos_cfi.primaryVertex
primaryVertex
Definition: impactParameterTagInfos_cfi.py:5
L1TMuonDQMOffline::m_GmtPtCuts
std::vector< int > m_GmtPtCuts
Definition: L1TMuonDQMOffline.h:164
electronCleaner_cfi.muons
muons
Definition: electronCleaner_cfi.py:12
L1TMuonDQMOffline::getPrimaryVertex
const reco::Vertex getPrimaryVertex(edm::Handle< reco::VertexCollection > &vertex, edm::Handle< reco::BeamSpot > &beamSpot)
Definition: L1TMuonDQMOffline.cc:421
L1TMuonDQMOffline::m_VtxInputTag
edm::EDGetTokenT< reco::VertexCollection > m_VtxInputTag
Definition: L1TMuonDQMOffline.h:167
L1TMuonDQMOffline::m_trigInputTag
edm::EDGetTokenT< trigger::TriggerEvent > m_trigInputTag
Definition: L1TMuonDQMOffline.h:169
gather_cfg.cout
cout
Definition: gather_cfg.py:145
BXVector::begin
const_iterator begin(int bx) const
L1TMuonDQMOffline::getTightMuons
void getTightMuons(edm::Handle< reco::MuonCollection > &muons, const reco::Vertex &vertex)
Definition: L1TMuonDQMOffline.cc:453
void L1TMuonDQMOffline::beginLuminosityBlock ( edm::LuminosityBlock const &  lumiBlock,
edm::EventSetup const &  c 
)
overrideprotectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 171 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), and m_verbose.

171  {
172  if(m_verbose) cout << "[L1TMuonDQMOffline:] Called beginLuminosityBlock at LS=" << lumiBlock.id().luminosityBlock() << endl;
173 }
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void L1TMuonDQMOffline::bookControlHistos ( DQMStore::IBooker &  ibooker)
protectedvirtual

Definition at line 300 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), gather_cfg::cout, m_ControlHistos, m_EfficiencyHistos, m_HistFolder, m_verbose, dataset::name, Pi, and DQMStore::IBooker::setCurrentFolder().

Referenced by bookHistograms().

300  {
301  if(m_verbose) cout << "[L1TMuonDQMOffline:] Booking Control Plot Histos" << endl;
302 
303  ibooker.setCurrentFolder(m_HistFolder+"/control_variables");
304 
305  string name = "MuonGmtDeltaR";
306  m_ControlHistos[name] = ibooker.book1D(name.c_str(),name.c_str(),25.,0.,2.5);
307 
308  name = "NTightVsAll";
309  m_ControlHistos[name] = ibooker.book2D(name.c_str(),name.c_str(),5,-0.5,4.5,5,-0.5,4.5);
310 
311  name = "NProbesVsTight";
312  m_ControlHistos[name] = ibooker.book2D(name.c_str(),name.c_str(),5,-0.5,4.5,5,-0.5,4.5);
313 
315  string name1 = "TagMuonEta_Histo";
316  m_EfficiencyHistos[0][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
317  string name2 = "TagMuonPhi_Histo";
318  m_EfficiencyHistos[0][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
319  string name3 = "TagMuonPt_Histo";
320  m_EfficiencyHistos[0][name3] = ibooker.book1D(name3.c_str(),name3.c_str(),50,0.,100.);
321 //*****
322  name1 = "ProbeMuonEta_Histo";
323  m_EfficiencyHistos[0][name1] = ibooker.book1D(name1.c_str(),name1.c_str(),50,-2.5,2.5);
324  name2 = "ProbeMuonPhi_Histo";
325  m_EfficiencyHistos[0][name2] = ibooker.book1D(name2.c_str(),name2.c_str(),24,-TMath::Pi(),TMath::Pi());
326  name3 = "ProbeMuonPt_Histo";
327  m_EfficiencyHistos[0][name3] = ibooker.book1D(name3.c_str(),name3.c_str(),50,0.,100.);
328 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
L1TMuonDQMOffline::m_HistFolder
std::string m_HistFolder
Definition: L1TMuonDQMOffline.h:163
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
L1TMuonDQMOffline::m_EfficiencyHistos
std::map< int, std::map< std::string, MonitorElement * > > m_EfficiencyHistos
Definition: L1TMuonDQMOffline.h:148
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
L1TMuonDQMOffline::m_ControlHistos
std::map< std::string, MonitorElement * > m_ControlHistos
Definition: L1TMuonDQMOffline.h:149
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void L1TMuonDQMOffline::bookEfficiencyHistos ( DQMStore::IBooker &  ibooker,
int  ptCut 
)
protectedvirtual

Definition at line 331 of file L1TMuonDQMOffline.cc.

References DQMStore::IBooker::book1D(), gather_cfg::cout, m_EfficiencyHistos, m_effVsEtaBins, m_effVsPhiBins, m_effVsPtBins, m_HistFolder, m_verbose, HiGenCleaner_cff::ptCut, and DQMStore::IBooker::setCurrentFolder().

Referenced by bookHistograms().

331  {
332  if(m_verbose) cout << "[L1TMuonDQMOffline:] Booking Efficiency Plot Histos for pt cut = " << ptCut << endl;
333 
334  stringstream ptCutToTag; ptCutToTag << ptCut;
335  string ptTag = ptCutToTag.str();
336 
337  ibooker.setCurrentFolder(m_HistFolder+"/numerators_and_denominators");
338 
339  std::vector<float> effVsPtBins(m_effVsPtBins.begin(), m_effVsPtBins.end());
340  int nEffVsPtBins = effVsPtBins.size() - 1;
341  float* ptBinsArray = &(effVsPtBins[0]);
342 
343  string name1 = "EffvsPt_" + ptTag + "_Den";
344  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
345  string name2 = "EffvsPt_" + ptTag + "_Num";
346  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
347 
348  name1 = "EffvsPt_OPEN_" + ptTag + "_Den";
349  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
350  name2 = "EffvsPt_OPEN_" + ptTag + "_Num";
351  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
352 
353  name1 = "EffvsPt_DOUBLE_" + ptTag + "_Den";
354  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
355  name2 = "EffvsPt_DOUBLE_" + ptTag + "_Num";
356  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
357 
358  name1 = "EffvsPt_SINGLE_" + ptTag + "_Den";
359  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPtBins, ptBinsArray);
360  name2 = "EffvsPt_SINGLE_" + ptTag + "_Num";
361  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPtBins, ptBinsArray);
362 
364 
365  std::vector<float> effVsPhiBins(m_effVsPhiBins.begin(), m_effVsPhiBins.end());
366  int nEffVsPhiBins = effVsPhiBins.size() - 1;
367  float* phiBinsArray = &(effVsPhiBins[0]);
368 
369  std::vector<float> effVsEtaBins(m_effVsEtaBins.begin(), m_effVsEtaBins.end());
370  int nEffVsEtaBins = effVsEtaBins.size() - 1;
371  float* etaBinsArray = &(effVsEtaBins[0]);
372 
373  name1 = "EffvsPhi_" + ptTag + "_Den";
374  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
375  name2 = "EffvsPhi_" + ptTag + "_Num";
376  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
377 
378  name1 = "EffvsEta_" + ptTag + "_Den";
379  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
380  name2 = "EffvsEta_" + ptTag + "_Num";
381  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
382 
384 
385  name1 = "EffvsPhi_OPEN_" + ptTag + "_Den";
386  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
387  name2 = "EffvsPhi_OPEN_" + ptTag + "_Num";
388  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
389 
390  name1 = "EffvsEta_OPEN_" + ptTag + "_Den";
391  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
392  name2 = "EffvsEta_OPEN_" + ptTag + "_Num";
393  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
394 
396 
397  name1 = "EffvsPhi_DOUBLE_" + ptTag + "_Den";
398  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
399  name2 = "EffvsPhi_DOUBLE_" + ptTag + "_Num";
400  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
401 
402  name1 = "EffvsEta_DOUBLE_" + ptTag + "_Den";
403  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
404  name2 = "EffvsEta_DOUBLE_" + ptTag + "_Num";
405  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
406 
408 
409  name1 = "EffvsPhi_SINGLE_" + ptTag + "_Den";
410  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsPhiBins, phiBinsArray);
411  name2 = "EffvsPhi_SINGLE_" + ptTag + "_Num";
412  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsPhiBins, phiBinsArray);
413 
414  name1 = "EffvsEta_SINGLE_" + ptTag + "_Den";
415  m_EfficiencyHistos[ptCut][name1] = ibooker.book1D(name1.c_str(), name1.c_str(), nEffVsEtaBins, etaBinsArray);
416  name2 = "EffvsEta_SINGLE_" + ptTag + "_Num";
417  m_EfficiencyHistos[ptCut][name2] = ibooker.book1D(name2.c_str(), name2.c_str(), nEffVsEtaBins, etaBinsArray);
418 }
HiGenCleaner_cff.ptCut
ptCut
Definition: HiGenCleaner_cff.py:12
L1TMuonDQMOffline::m_effVsPhiBins
std::vector< double > m_effVsPhiBins
Definition: L1TMuonDQMOffline.h:174
L1TMuonDQMOffline::m_HistFolder
std::string m_HistFolder
Definition: L1TMuonDQMOffline.h:163
L1TMuonDQMOffline::m_effVsPtBins
std::vector< double > m_effVsPtBins
Definition: L1TMuonDQMOffline.h:173
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
L1TMuonDQMOffline::m_EfficiencyHistos
std::map< int, std::map< std::string, MonitorElement * > > m_EfficiencyHistos
Definition: L1TMuonDQMOffline.h:148
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
L1TMuonDQMOffline::m_effVsEtaBins
std::vector< double > m_effVsEtaBins
Definition: L1TMuonDQMOffline.h:175
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void L1TMuonDQMOffline::bookHistograms ( DQMStore::IBooker &  ibooker,
const edm::Run &  run,
const edm::EventSetup &  iSetup 
)
overrideprotectedvirtual

Definition at line 141 of file L1TMuonDQMOffline.cc.

References bookControlHistos(), bookEfficiencyHistos(), gather_cfg::cout, createfilelist::int, m_GmtPtCuts, m_hltConfig, m_trigIndices, m_trigNames, m_verbose, HLTConfigProvider::size(), and HLTConfigProvider::triggerName().

141  {
142  //book histos
143  bookControlHistos(ibooker);
144  vector<int>::const_iterator gmtPtCutsIt = m_GmtPtCuts.begin();
145  vector<int>::const_iterator gmtPtCutsEnd = m_GmtPtCuts.end();
146 
147  for (; gmtPtCutsIt!=gmtPtCutsEnd; ++ gmtPtCutsIt) {
148  bookEfficiencyHistos(ibooker, (*gmtPtCutsIt));
149  }
150 
151  vector<string>::const_iterator trigNamesIt = m_trigNames.begin();
152  vector<string>::const_iterator trigNamesEnd = m_trigNames.end();
153 
154  for (; trigNamesIt!=trigNamesEnd; ++trigNamesIt) {
155  TString tNameTmp = TString(*trigNamesIt); // use TString as it handles regex
156  TRegexp tNamePattern = TRegexp(tNameTmp,true);
157  int tIndex = -1;
158 
159  for (unsigned ipath = 0; ipath < m_hltConfig.size(); ++ipath) {
160  TString tmpName = TString(m_hltConfig.triggerName(ipath));
161  if (tmpName.Contains(tNamePattern)) {
162  tIndex = int(ipath);
163  m_trigIndices.push_back(tIndex);
164  }
165  }
166  if (tIndex < 0 && m_verbose) cout << "[L1TMuonDQMOffline:] Warning: Could not find trigger " << (*trigNamesIt) << endl;
167  }
168 }
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:76
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:91
L1TMuonDQMOffline::bookEfficiencyHistos
virtual void bookEfficiencyHistos(DQMStore::IBooker &ibooker, int ptCut)
Definition: L1TMuonDQMOffline.cc:331
L1TMuonDQMOffline::m_trigNames
std::vector< std::string > m_trigNames
Definition: L1TMuonDQMOffline.h:172
L1TMuonDQMOffline::bookControlHistos
virtual void bookControlHistos(DQMStore::IBooker &)
Definition: L1TMuonDQMOffline.cc:300
L1TMuonDQMOffline::m_GmtPtCuts
std::vector< int > m_GmtPtCuts
Definition: L1TMuonDQMOffline.h:164
L1TMuonDQMOffline::m_hltConfig
HLTConfigProvider m_hltConfig
Definition: L1TMuonDQMOffline.h:141
gather_cfg.cout
cout
Definition: gather_cfg.py:145
L1TMuonDQMOffline::m_trigIndices
std::vector< int > m_trigIndices
Definition: L1TMuonDQMOffline.h:177
void L1TMuonDQMOffline::dqmBeginRun ( const edm::Run &  run,
const edm::EventSetup &  iSetup 
)
overrideprotectedvirtual

Definition at line 134 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, HLTConfigProvider::init(), m_hltConfig, m_trigProcess, and m_verbose.

134  {
135  if (m_verbose) cout << "[L1TMuonDQMOffline:] Called beginRun." << endl;
136  bool changed = true;
137  m_hltConfig.init(run,iSetup,m_trigProcess,changed);
138 }
L1TMuonDQMOffline::m_trigProcess
std::string m_trigProcess
Definition: L1TMuonDQMOffline.h:170
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
Definition: HLTConfigProvider.cc:42
L1TMuonDQMOffline::m_hltConfig
HLTConfigProvider m_hltConfig
Definition: L1TMuonDQMOffline.h:141
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void L1TMuonDQMOffline::dqmEndLuminosityBlock ( edm::LuminosityBlock const &  lumiBlock,
edm::EventSetup const &  c 
)
protectedvirtual

Definition at line 176 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), and m_verbose.

176  {
177  if(m_verbose) cout << "[L1TMuonDQMOffline:] Called endLuminosityBlock at LS=" << lumiBlock.id().luminosityBlock() << endl;
178 }
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void L1TMuonDQMOffline::getMuonGmtPairs ( edm::Handle< l1t::MuonBxCollection > &  gmtCands)
private

Definition at line 507 of file L1TMuonDQMOffline.cc.

References BXVector< T >::begin(), gather_cfg::cout, MuonGmtPair::dR(), BXVector< T >::end(), PVValHelper::eta, m_ControlHistos, m_EfficiencyHistos, m_MaxGmtMuonDR, m_MuonGmtPairs, m_ProbeMuons, m_verbose, RPCpg::mu, phi, and EnergyCorrector::pt.

Referenced by analyze().

507  {
508 
509  m_MuonGmtPairs.clear();
510  if (m_verbose) cout << "[L1TMuonDQMOffline:] Getting muon GMT pairs" << endl;
511 
512  vector<const reco::Muon*>::const_iterator probeMuIt = m_ProbeMuons.begin();
513  vector<const reco::Muon*>::const_iterator probeMuEnd = m_ProbeMuons.end();
514  vector<l1t::Muon> gmtContainer;
515 
516  for (auto mu = gmtCands->begin(0); mu != gmtCands->end(0); ++mu) {
517  gmtContainer.push_back(*mu);
518  }
519 
520  vector<l1t::Muon>::const_iterator gmtIt;
521  vector<l1t::Muon>::const_iterator gmtEnd = gmtContainer.end();
522 
523  for (; probeMuIt!=probeMuEnd; ++probeMuIt) {
524  float eta = (*probeMuIt)->eta();
525  float phi = (*probeMuIt)->phi();
526  float pt = (*probeMuIt)->pt();
527 
528  m_EfficiencyHistos[0]["ProbeMuonEta_Histo"]->Fill(eta);
529  m_EfficiencyHistos[0]["ProbeMuonPhi_Histo"]->Fill(phi);
530  m_EfficiencyHistos[0]["ProbeMuonPt_Histo"]->Fill(pt);
531 
532  MuonGmtPair pairBestCand((*probeMuIt),0);
533 // pairBestCand.propagate(m_BField,m_propagatorAlong,m_propagatorOpposite);
534  gmtIt = gmtContainer.begin();
535 
536  for(; gmtIt!=gmtEnd; ++gmtIt) {
537  MuonGmtPair pairTmpCand((*probeMuIt),&(*gmtIt));
538 // pairTmpCand.propagate(m_BField,m_propagatorAlong,m_propagatorOpposite);
539 
540  if ( (pairTmpCand.dR() < m_MaxGmtMuonDR) && (pairTmpCand.dR() < pairBestCand.dR() ) ) {
541  pairBestCand = pairTmpCand;
542  }
543 
544  }
545  m_MuonGmtPairs.push_back(pairBestCand);
546  m_ControlHistos["MuonGmtDeltaR"]->Fill(pairBestCand.dR());
547  }
548 }
BXVector::end
const_iterator end(int bx) const
L1TMuonDQMOffline::m_ProbeMuons
std::vector< const reco::Muon * > m_ProbeMuons
Definition: L1TMuonDQMOffline.h:153
L1TMuonDQMOffline::m_EfficiencyHistos
std::map< int, std::map< std::string, MonitorElement * > > m_EfficiencyHistos
Definition: L1TMuonDQMOffline.h:148
RPCpg::mu
const int mu
Definition: Constants.h:22
L1TMuonDQMOffline::m_MuonGmtPairs
std::vector< MuonGmtPair > m_MuonGmtPairs
Definition: L1TMuonDQMOffline.h:154
L1TMuonDQMOffline::m_ControlHistos
std::map< std::string, MonitorElement * > m_ControlHistos
Definition: L1TMuonDQMOffline.h:149
gather_cfg.cout
cout
Definition: gather_cfg.py:145
BXVector::begin
const_iterator begin(int bx) const
const reco::Vertex L1TMuonDQMOffline::getPrimaryVertex ( edm::Handle< reco::VertexCollection > &  vertex,
edm::Handle< reco::BeamSpot > &  beamSpot 
)
private

Definition at line 421 of file L1TMuonDQMOffline.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), edm::HandleBase::isValid(), reco::BeamSpot::position(), impactParameterTagInfos_cfi::primaryVertex, and reco::BeamSpot::sigmaZ().

Referenced by analyze().

422  {
423  Vertex::Point posVtx;
424  Vertex::Error errVtx;
425 
426  bool hasPrimaryVertex = false;
427 
428  if (vertex.isValid()) {
429  vector<Vertex>::const_iterator vertexIt = vertex->begin();
430  vector<Vertex>::const_iterator vertexEnd = vertex->end();
431 
432  for (;vertexIt!=vertexEnd;++vertexIt) {
433  if (vertexIt->isValid() && !vertexIt->isFake()) {
434  posVtx = vertexIt->position();
435  errVtx = vertexIt->error();
436  hasPrimaryVertex = true;
437  break;
438  }
439  }
440  }
441 
442  if ( !hasPrimaryVertex ) {
443  posVtx = beamSpot->position();
444  errVtx(0,0) = beamSpot->BeamWidthX();
445  errVtx(1,1) = beamSpot->BeamWidthY();
446  errVtx(2,2) = beamSpot->sigmaZ();
447  }
448  const Vertex primaryVertex(posVtx,errVtx);
449  return primaryVertex;
450 }
reco::Vertex::Error
math::Error< dimension >::type Error
covariance error matrix (3x3)
Definition: Vertex.h:43
reco::Vertex::Point
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
reco::BeamSpot::BeamWidthX
double BeamWidthX() const
beam width X
Definition: BeamSpot.h:86
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
impactParameterTagInfos_cfi.primaryVertex
primaryVertex
Definition: impactParameterTagInfos_cfi.py:5
reco::BeamSpot::sigmaZ
double sigmaZ() const
sigma z
Definition: BeamSpot.h:80
reco::BeamSpot::BeamWidthY
double BeamWidthY() const
beam width Y
Definition: BeamSpot.h:88
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
void L1TMuonDQMOffline::getProbeMuons ( edm::Handle< edm::TriggerResults > &  trigResults,
edm::Handle< trigger::TriggerEvent > &  trigEvent 
)
private

Definition at line 469 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, PVValHelper::eta, m_ControlHistos, m_EfficiencyHistos, m_ProbeMuons, m_TightMuons, m_verbose, matchHlt(), phi, EnergyCorrector::pt, and mathSSE::sqrt().

Referenced by analyze().

470  {
471 
472  if (m_verbose) cout << "[L1TMuonDQMOffline:] getting probe muons" << endl;
473  m_ProbeMuons.clear();
474 
475  vector<const reco::Muon*>::const_iterator probeCandIt = m_TightMuons.begin();
476 // vector<const reco::Muon*>::const_iterator probeMuIt = m_ProbeMuons.begin();
477  vector<const reco::Muon*>::const_iterator tightMuonsEnd = m_TightMuons.end();
478 
479  for (; probeCandIt!=tightMuonsEnd; ++probeCandIt) {
480  bool tagHasTrig = false;
481  vector<const reco::Muon*>::const_iterator tagCandIt = m_TightMuons.begin();
482  float deltar = 0.;
483 
484  for (; tagCandIt!=tightMuonsEnd; ++tagCandIt) {
485  float eta = (*tagCandIt)->eta();
486  float phi = (*tagCandIt)->phi();
487  float pt = (*tagCandIt)->pt();
488  float dEta = eta - (*probeCandIt)->eta();
489  float dPhi = phi - (*probeCandIt)->phi();
490 
491  deltar = sqrt(dEta*dEta + dPhi*dPhi);
492 
493  if ( (*tagCandIt) == (*probeCandIt) || (deltar<0.7) ) continue; // CB has a little bias for closed-by muons
494  tagHasTrig |= matchHlt(trigEvent,(*tagCandIt));
495  if (tagHasTrig) {
496  m_EfficiencyHistos[0]["TagMuonEta_Histo"]->Fill(eta);
497  m_EfficiencyHistos[0]["TagMuonPhi_Histo"]->Fill(phi);
498  m_EfficiencyHistos[0]["TagMuonPt_Histo"]->Fill(pt);
499  }
500  }
501  if (tagHasTrig) m_ProbeMuons.push_back((*probeCandIt));
502  }
503  m_ControlHistos["NProbesVsTight"]->Fill(m_TightMuons.size(),m_ProbeMuons.size());
504 }
L1TMuonDQMOffline::matchHlt
bool matchHlt(edm::Handle< trigger::TriggerEvent > &triggerEvent, const reco::Muon *mu)
Definition: L1TMuonDQMOffline.cc:551
L1TMuonDQMOffline::m_ProbeMuons
std::vector< const reco::Muon * > m_ProbeMuons
Definition: L1TMuonDQMOffline.h:153
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
L1TMuonDQMOffline::m_EfficiencyHistos
std::map< int, std::map< std::string, MonitorElement * > > m_EfficiencyHistos
Definition: L1TMuonDQMOffline.h:148
deltar
Definition: deltar.py:1
L1TMuonDQMOffline::m_TightMuons
std::vector< const reco::Muon * > m_TightMuons
Definition: L1TMuonDQMOffline.h:152
L1TMuonDQMOffline::m_ControlHistos
std::map< std::string, MonitorElement * > m_ControlHistos
Definition: L1TMuonDQMOffline.h:149
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void L1TMuonDQMOffline::getTightMuons ( edm::Handle< reco::MuonCollection > &  muons,
const reco::Vertex &  vertex 
)
private

Definition at line 453 of file L1TMuonDQMOffline.cc.

References gather_cfg::cout, muon::isTightMuon(), m_ControlHistos, m_TightMuons, and m_verbose.

Referenced by analyze().

453  {
454 
455  if (m_verbose) cout << "[L1TMuonDQMOffline:] Getting tight muons" << endl;
456  m_TightMuons.clear();
457  MuonCollection::const_iterator muonIt = muons->begin();
458  MuonCollection::const_iterator muonEnd = muons->end();
459 
460  for(; muonIt!=muonEnd; ++muonIt) {
461  if (muon::isTightMuon((*muonIt), vertex)) {
462  m_TightMuons.push_back(&(*muonIt));
463  }
464  }
465  m_ControlHistos["NTightVsAll"]->Fill(muons->size(),m_TightMuons.size());
466 }
L1TMuonDQMOffline::m_TightMuons
std::vector< const reco::Muon * > m_TightMuons
Definition: L1TMuonDQMOffline.h:152
L1TMuonDQMOffline::m_ControlHistos
std::map< std::string, MonitorElement * > m_ControlHistos
Definition: L1TMuonDQMOffline.h:149
gather_cfg.cout
cout
Definition: gather_cfg.py:145
muon::isTightMuon
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:837
bool L1TMuonDQMOffline::matchHlt ( edm::Handle< trigger::TriggerEvent > &  triggerEvent,
const reco::Muon *  mu 
)
private

Definition at line 551 of file L1TMuonDQMOffline.cc.

References DEFINE_FWK_MODULE, deltaR(), trigger::TriggerEvent::filterIds(), trigger::TriggerEvent::filterIndex(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::getObjects(), m_hltConfig, m_MaxHltMuonDR, m_trigIndices, m_trigProcess, HLTConfigProvider::moduleLabels(), HLTConfigProvider::size(), and TriggerAnalyzer::trigObjs.

Referenced by getProbeMuons().

551  {
552 
553  double matchDeltaR = 9999;
554 
555  TriggerObjectCollection trigObjs = triggerEvent->getObjects();
556 
557  vector<int>::const_iterator trigIndexIt = m_trigIndices.begin();
558  vector<int>::const_iterator trigIndexEnd = m_trigIndices.end();
559 
560  for(; trigIndexIt!=trigIndexEnd; ++trigIndexIt) {
561  const vector<string> moduleLabels(m_hltConfig.moduleLabels(*trigIndexIt));
562  const unsigned moduleIndex = m_hltConfig.size((*trigIndexIt))-2;
563  const unsigned hltFilterIndex = triggerEvent->filterIndex(InputTag(moduleLabels[moduleIndex],"",m_trigProcess));
564 
565  if (hltFilterIndex < triggerEvent->sizeFilters()) {
566  const Keys triggerKeys(triggerEvent->filterKeys(hltFilterIndex));
567  const Vids triggerVids(triggerEvent->filterIds(hltFilterIndex));
568  const unsigned nTriggers = triggerVids.size();
569  for (size_t iTrig = 0; iTrig < nTriggers; ++iTrig) {
570  const TriggerObject trigObject = trigObjs[triggerKeys[iTrig]];
571  double dRtmp = deltaR((*mu),trigObject);
572  if (dRtmp < matchDeltaR) matchDeltaR = dRtmp;
573  }
574  }
575  }
576  return (matchDeltaR < m_MaxHltMuonDR);
577 }
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:76
trigger::TriggerEvent::filterKeys
const Keys & filterKeys(trigger::size_type index) const
Definition: TriggerEvent.h:111
trigger::TriggerEvent::filterIndex
trigger::size_type filterIndex(const edm::InputTag &filterTag) const
find index of filter in data-member vector from filter tag
Definition: TriggerEvent.h:123
L1TMuonDQMOffline::m_trigProcess
std::string m_trigProcess
Definition: L1TMuonDQMOffline.h:170
trigger::TriggerEvent::filterIds
const Vids & filterIds(trigger::size_type index) const
Definition: TriggerEvent.h:110
trigger::TriggerEvent::getObjects
const TriggerObjectCollection & getObjects() const
Definition: TriggerEvent.h:98
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:101
deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:81
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:20
L1TMuonDQMOffline::m_hltConfig
HLTConfigProvider m_hltConfig
Definition: L1TMuonDQMOffline.h:141
trigger::Vids
std::vector< int > Vids
Definition: TriggerTypeDefs.h:22
L1TMuonDQMOffline::m_trigIndices
std::vector< int > m_trigIndices
Definition: L1TMuonDQMOffline.h:177

Member Data Documentation

edm::ESHandle<MagneticField> L1TMuonDQMOffline::m_BField
private

Definition at line 143 of file L1TMuonDQMOffline.h.

Referenced by analyze().

edm::EDGetTokenT<reco::BeamSpot> L1TMuonDQMOffline::m_BsInputTag
private

Definition at line 168 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::map<std::string, MonitorElement*> L1TMuonDQMOffline::m_ControlHistos
private

Definition at line 149 of file L1TMuonDQMOffline.h.

Referenced by bookControlHistos(), getMuonGmtPairs(), getProbeMuons(), and getTightMuons().

std::map<int, std::map<std::string, MonitorElement*> > L1TMuonDQMOffline::m_EfficiencyHistos
private

Definition at line 148 of file L1TMuonDQMOffline.h.

Referenced by analyze(), bookControlHistos(), bookEfficiencyHistos(), getMuonGmtPairs(), and getProbeMuons().

std::vector<double> L1TMuonDQMOffline::m_effVsEtaBins
private

Definition at line 175 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos().

std::vector<double> L1TMuonDQMOffline::m_effVsPhiBins
private

Definition at line 174 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos().

std::vector<double> L1TMuonDQMOffline::m_effVsPtBins
private

Definition at line 173 of file L1TMuonDQMOffline.h.

Referenced by bookEfficiencyHistos().

edm::EDGetTokenT<l1t::MuonBxCollection> L1TMuonDQMOffline::m_GmtInputTag
private

Definition at line 166 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::vector<int> L1TMuonDQMOffline::m_GmtPtCuts
private

Definition at line 164 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and bookHistograms().

std::string L1TMuonDQMOffline::m_HistFolder
private

Definition at line 163 of file L1TMuonDQMOffline.h.

Referenced by bookControlHistos(), and bookEfficiencyHistos().

HLTConfigProvider L1TMuonDQMOffline::m_hltConfig
private

Definition at line 141 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms(), dqmBeginRun(), and matchHlt().

BXVector<l1t::Muon> L1TMuonDQMOffline::m_L1tL1tMuons
private

Definition at line 159 of file L1TMuonDQMOffline.h.

std::vector<l1t::MuonBxCollection> L1TMuonDQMOffline::m_L1tMuons
private

Definition at line 157 of file L1TMuonDQMOffline.h.

float L1TMuonDQMOffline::m_MaxGmtMuonDR
private

Definition at line 180 of file L1TMuonDQMOffline.h.

Referenced by getMuonGmtPairs(), and L1TMuonDQMOffline().

float L1TMuonDQMOffline::m_MaxHltMuonDR
private

Definition at line 181 of file L1TMuonDQMOffline.h.

Referenced by L1TMuonDQMOffline(), and matchHlt().

float L1TMuonDQMOffline::m_MaxMuonEta
private

Definition at line 179 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and L1TMuonDQMOffline().

std::vector<MuonGmtPair> L1TMuonDQMOffline::m_MuonGmtPairs
private

Definition at line 154 of file L1TMuonDQMOffline.h.

Referenced by analyze(), and getMuonGmtPairs().

edm::EDGetTokenT<reco::MuonCollection> L1TMuonDQMOffline::m_MuonInputTag
private

Definition at line 165 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::vector<const reco::Muon*> L1TMuonDQMOffline::m_ProbeMuons
private

Definition at line 153 of file L1TMuonDQMOffline.h.

Referenced by getMuonGmtPairs(), and getProbeMuons().

edm::ESHandle<Propagator> L1TMuonDQMOffline::m_propagatorAlong
private

Definition at line 144 of file L1TMuonDQMOffline.h.

Referenced by analyze().

edm::ESHandle<Propagator> L1TMuonDQMOffline::m_propagatorOpposite
private

Definition at line 145 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::vector<reco::MuonCollection> L1TMuonDQMOffline::m_RecoMuons
private

Definition at line 156 of file L1TMuonDQMOffline.h.

std::vector<reco::Muon> L1TMuonDQMOffline::m_RecoRecoMuons
private

Definition at line 158 of file L1TMuonDQMOffline.h.

std::vector<const reco::Muon*> L1TMuonDQMOffline::m_TightMuons
private

Definition at line 152 of file L1TMuonDQMOffline.h.

Referenced by getProbeMuons(), and getTightMuons().

std::vector<int> L1TMuonDQMOffline::m_trigIndices
private

Definition at line 177 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms(), and matchHlt().

edm::EDGetTokenT<trigger::TriggerEvent> L1TMuonDQMOffline::m_trigInputTag
private

Definition at line 169 of file L1TMuonDQMOffline.h.

Referenced by analyze().

std::vector<std::string> L1TMuonDQMOffline::m_trigNames
private

Definition at line 172 of file L1TMuonDQMOffline.h.

Referenced by bookHistograms().

std::string L1TMuonDQMOffline::m_trigProcess
private

Definition at line 170 of file L1TMuonDQMOffline.h.

Referenced by dqmBeginRun(), and matchHlt().

edm::EDGetTokenT<edm::TriggerResults> L1TMuonDQMOffline::m_trigProcess_token
private

Definition at line 171 of file L1TMuonDQMOffline.h.

Referenced by analyze().

bool L1TMuonDQMOffline::m_verbose
private

Definition at line 162 of file L1TMuonDQMOffline.h.

Referenced by analyze(), beginLuminosityBlock(), bookControlHistos(), bookEfficiencyHistos(), bookHistograms(), dqmBeginRun(), dqmEndLuminosityBlock(), getMuonGmtPairs(), getProbeMuons(), getTightMuons(), and L1TMuonDQMOffline().

edm::EDGetTokenT<reco::VertexCollection> L1TMuonDQMOffline::m_VtxInputTag
private

Definition at line 167 of file L1TMuonDQMOffline.h.

Referenced by analyze().

Generated for CMSSW Reference Manual by doxygen 1.8.11