CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
DiMuonHistograms Class Reference

#include <DiMuonHistograms.h>

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

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 DiMuonHistograms (const edm::ParameterSet &pset)
 
virtual ~DiMuonHistograms ()
 
- 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 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< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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)
 

Private Attributes

int etaBBin
 
int etaBin
 
double etaBMax
 
double etaBMin
 
double EtaCutMax
 
double EtaCutMin
 
int etaEBin
 
double etaECMax
 
double etaECMin
 
std::string EtaName [3]
 
int etaOvlpBin
 
std::vector< MonitorElement * > GlbGlbMuon_HM
 
std::vector< MonitorElement * > GlbGlbMuon_LM
 
double HighMassMax
 
double HighMassMin
 
std::vector< MonitorElement * > LooseLooseMuon
 
double LowMassMax
 
double LowMassMin
 
std::vector< MonitorElement * > MediumMediumMuon
 
std::string metname
 
int nGlbGlb
 
int nLooseLoose
 
int nMediumMedium
 
int nTightTight
 
edm::ParameterSet parameters
 
std::vector< MonitorElement * > SoftSoftMuon
 
std::vector< MonitorElement * > StaTrkMuon_HM
 
std::vector< MonitorElement * > StaTrkMuon_LM
 
MonitorElementtest
 
edm::EDGetTokenT< reco::BeamSpottheBeamSpotLabel_
 
std::string theFolder
 
edm::EDGetTokenT< edm::View< reco::Muon > > theMuonCollectionLabel_
 
MuonServiceProxytheService
 
edm::EDGetTokenT< reco::VertexCollectiontheVertexLabel_
 
std::vector< MonitorElement * > TightTightMuon
 
std::vector< MonitorElement * > TrkTrkMuon_HM
 
std::vector< MonitorElement * > TrkTrkMuon_LM
 

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, GlobalCacheLuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCacheRunContext
 
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::NoCacheglobalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr< dqmDetails::NoCacheglobalBeginRunSummary (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)
 
- 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)
 

Detailed Description

Class DiMuonHistograms

DQM monitoring for dimuon mass

Author: S.Folgueras, U. Oviedo

Definition at line 33 of file DiMuonHistograms.h.

Constructor & Destructor Documentation

DiMuonHistograms::DiMuonHistograms ( const edm::ParameterSet pset)

Definition at line 32 of file DiMuonHistograms.cc.

References conversionPostprocessing_cfi::etaBin, and metProducer_cfi::parameters.

32  {
33 
34  // initialise parameters:
35  parameters = pSet;
36 
37  // counter
38  nTightTight = 0;
39  nMediumMedium = 0;
40  nLooseLoose = 0;
41  nGlbGlb = 0;
42 
43  // declare consumes:
44  theMuonCollectionLabel_ = consumes<edm::View<reco::Muon> > (parameters.getParameter<edm::InputTag>("MuonCollection"));
45  theVertexLabel_ = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
46 
47  theBeamSpotLabel_ = mayConsume<reco::BeamSpot> (parameters.getParameter<edm::InputTag>("BeamSpotLabel"));
48 
49  etaBin = parameters.getParameter<int>("etaBin");
50  etaBBin = parameters.getParameter<int>("etaBBin");
51  etaEBin = parameters.getParameter<int>("etaEBin");
52 
53  etaBMin = parameters.getParameter<double>("etaBMin");
54  etaBMax = parameters.getParameter<double>("etaBMax");
55  etaECMin = parameters.getParameter<double>("etaECMin");
56  etaECMax = parameters.getParameter<double>("etaECMax");
57 
58  LowMassMin = parameters.getParameter<double>("LowMassMin");
59  LowMassMax = parameters.getParameter<double>("LowMassMax");
60  HighMassMin = parameters.getParameter<double>("HighMassMin");
61  HighMassMax = parameters.getParameter<double>("HighMassMax");
62 
63  theFolder = parameters.getParameter<string>("folder");
64 }
T getParameter(std::string const &) const
edm::ParameterSet parameters
edm::EDGetTokenT< edm::View< reco::Muon > > theMuonCollectionLabel_
edm::EDGetTokenT< reco::BeamSpot > theBeamSpotLabel_
edm::EDGetTokenT< reco::VertexCollection > theVertexLabel_
std::string theFolder
DiMuonHistograms::~DiMuonHistograms ( )
virtual

Definition at line 66 of file DiMuonHistograms.cc.

66 { }

Member Function Documentation

void DiMuonHistograms::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
override

Definition at line 98 of file DiMuonHistograms.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), ALCARECOTkAlJpsiMuMu_cff::charge, relativeConstraints::error, edm::Event::getByToken(), muon::isLooseMuon(), muon::isMediumMuon(), muon::isSoftMuon(), muon::isTightMuon(), edm::HandleBase::isValid(), LogTrace, metname, electronCleaner_cfi::muons, reco::BeamSpot::position(), position, reco::BeamSpot::sigmaZ(), and badGlobalMuonTaggersAOD_cff::vtx.

98  {
99 
100  LogTrace(metname)<<"[DiMuonHistograms] Analyze the mu in different eta regions";
102  iEvent.getByToken(theMuonCollectionLabel_, muons);
103 
104  // =================================================================================
105  // Look for the Primary Vertex (and use the BeamSpot instead, if you can't find it):
106  reco::Vertex::Point posVtx;
107  reco::Vertex::Error errVtx;
108  unsigned int theIndexOfThePrimaryVertex = 999.;
109 
111  iEvent.getByToken(theVertexLabel_, vertex);
112  if (vertex.isValid()){
113  for (unsigned int ind=0; ind<vertex->size(); ++ind) {
114  if ( (*vertex)[ind].isValid() && !((*vertex)[ind].isFake()) ) {
115  theIndexOfThePrimaryVertex = ind;
116  break;
117  }
118  }
119  }
120 
121  if (theIndexOfThePrimaryVertex<100) {
122  posVtx = ((*vertex)[theIndexOfThePrimaryVertex]).position();
123  errVtx = ((*vertex)[theIndexOfThePrimaryVertex]).error();
124  }
125  else {
126  LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";
127 
128  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
129  iEvent.getByToken(theBeamSpotLabel_,recoBeamSpotHandle);
130  reco::BeamSpot bs = *recoBeamSpotHandle;
131 
132  posVtx = bs.position();
133  errVtx(0,0) = bs.BeamWidthX();
134  errVtx(1,1) = bs.BeamWidthY();
135  errVtx(2,2) = bs.sigmaZ();
136  }
137 
138  const reco::Vertex vtx(posVtx,errVtx);
139 
140  if(!muons.isValid()) return;
141 
142  // Loop on muon collection
143  TLorentzVector Mu1, Mu2;
144  float charge = 99.;
145  float InvMass = -99.;
146 
147  //Eta regions
148  double EtaCutMin[] = {0, etaBMin, etaECMin};
149  double EtaCutMax[] = {2.4, etaBMax, etaECMax};
150 
151  for (edm::View<reco::Muon>::const_iterator muon1 = muons->begin(); muon1 != muons->end(); ++muon1){
152  LogTrace(metname)<<"[DiMuonHistograms] loop over 1st muon"<<endl;
153 
154  // Loop on second muons to fill invariant mass plots
155  for (edm::View<reco::Muon>::const_iterator muon2 = muon1; muon2 != muons->end(); ++muon2){
156  LogTrace(metname)<<"[DiMuonHistograms] loop over 2nd muon"<<endl;
157  if (muon1==muon2) continue;
158 
159  // Global-Global Muon
160  if (muon1->isGlobalMuon() && muon2->isGlobalMuon()) {
161  LogTrace(metname)<<"[DiMuonHistograms] Glb-Glb pair"<<endl;
162  reco::TrackRef recoCombinedGlbTrack1 = muon1->combinedMuon();
163  reco::TrackRef recoCombinedGlbTrack2 = muon2->combinedMuon();
164  Mu1.SetPxPyPzE(recoCombinedGlbTrack1->px(), recoCombinedGlbTrack1->py(),recoCombinedGlbTrack1->pz(), recoCombinedGlbTrack1->p());
165  Mu2.SetPxPyPzE(recoCombinedGlbTrack2->px(), recoCombinedGlbTrack2->py(),recoCombinedGlbTrack2->pz(), recoCombinedGlbTrack2->p());
166 
167  charge = recoCombinedGlbTrack1->charge()*recoCombinedGlbTrack2->charge();
168  if (charge < 0) {
169  InvMass = (Mu1+Mu2).M();
170  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
171  if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] &&
172  fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
173  if (InvMass < LowMassMax) GlbGlbMuon_LM[iEtaRegion]->Fill(InvMass);
174  if (InvMass > HighMassMin) GlbGlbMuon_HM[iEtaRegion]->Fill(InvMass);
175  }
176  }
177  }
178  // Also Tight-Tight Muon Selection
179  if ( muon::isTightMuon(*muon1, vtx) &&
180  muon::isTightMuon(*muon2, vtx) ) {
181  test->Fill(InvMass);
182  LogTrace(metname)<<"[DiMuonHistograms] Tight-Tight pair"<<endl;
183  if (charge < 0){
184  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
185  if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] &&
186  fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
187  if (InvMass > 55. && InvMass < 125.) TightTightMuon[iEtaRegion]->Fill(InvMass);
188  }
189  }
190  }
191  }
192  // Also Medium-Medium Muon Selection
193  if ( muon::isMediumMuon(*muon1) &&
194  muon::isMediumMuon(*muon2) ) {
195  test->Fill(InvMass);
196  LogTrace(metname)<<"[DiMuonHistograms] Medium-Medium pair"<<endl;
197  if (charge < 0){
198  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
199  if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] &&
200  fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
201  if (InvMass > 55. && InvMass < 125.) MediumMediumMuon[iEtaRegion]->Fill(InvMass);
202  }
203  }
204  }
205  }
206  // Also Loose-Loose Muon Selection
207  if ( muon::isLooseMuon(*muon1) &&
208  muon::isLooseMuon(*muon2) ) {
209  test->Fill(InvMass);
210  LogTrace(metname)<<"[DiMuonHistograms] Loose-Loose pair"<<endl;
211  if (charge < 0){
212  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
213  if(fabs(recoCombinedGlbTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack1->eta())<EtaCutMax[iEtaRegion] &&
214  fabs(recoCombinedGlbTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoCombinedGlbTrack2->eta())<EtaCutMax[iEtaRegion]){
215  if (InvMass > 55. && InvMass < 125.) LooseLooseMuon[iEtaRegion]->Fill(InvMass);
216  }
217  }
218  }
219  }
220  }
221 
222  // Now check for STA-TRK
223  if (muon2->isStandAloneMuon() && muon1->isTrackerMuon()) {
224  LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
225  reco::TrackRef recoStaTrack = muon2->standAloneMuon();
226  reco::TrackRef recoTrack = muon1->track();
227  Mu2.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(),recoStaTrack->pz(), recoStaTrack->p());
228  Mu1.SetPxPyPzE(recoTrack->px(), recoTrack->py(),recoTrack->pz(), recoTrack->p());
229 
230  charge = recoStaTrack->charge()*recoTrack->charge();
231  if (charge < 0) {
232  InvMass = (Mu1+Mu2).M();
233  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
234  if(fabs(recoStaTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoStaTrack->eta())<EtaCutMax[iEtaRegion] &&
235  fabs(recoTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack->eta())<EtaCutMax[iEtaRegion]){
236  if (InvMass < LowMassMax) StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
237  if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
238  }
239  }
240  }
241  }
242  if (muon1->isStandAloneMuon() && muon2->isTrackerMuon()) {
243  LogTrace(metname)<<"[DiMuonHistograms] STA-Trk pair"<<endl;
244  reco::TrackRef recoStaTrack = muon1->standAloneMuon();
245  reco::TrackRef recoTrack = muon2->track();
246  Mu1.SetPxPyPzE(recoStaTrack->px(), recoStaTrack->py(),recoStaTrack->pz(), recoStaTrack->p());
247  Mu2.SetPxPyPzE(recoTrack->px(), recoTrack->py(),recoTrack->pz(), recoTrack->p());
248 
249  charge = recoStaTrack->charge()*recoTrack->charge();
250  if (charge < 0) {
251  InvMass = (Mu1+Mu2).M();
252  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
253  if(fabs(recoStaTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoStaTrack->eta())<EtaCutMax[iEtaRegion] &&
254  fabs(recoTrack->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack->eta())<EtaCutMax[iEtaRegion]){
255  if (InvMass < LowMassMax) StaTrkMuon_LM[iEtaRegion]->Fill(InvMass);
256  if (InvMass > HighMassMin) StaTrkMuon_HM[iEtaRegion]->Fill(InvMass);
257  }
258  }
259  }
260  }
261 
262  // TRK-TRK dimuon
263  if (muon1->isTrackerMuon() && muon2->isTrackerMuon()) {
264  LogTrace(metname)<<"[DiMuonHistograms] Trk-Trk dimuon pair"<<endl;
265  reco::TrackRef recoTrack2 = muon2->track();
266  reco::TrackRef recoTrack1 = muon1->track();
267  Mu2.SetPxPyPzE(recoTrack2->px(), recoTrack2->py(),recoTrack2->pz(), recoTrack2->p());
268  Mu1.SetPxPyPzE(recoTrack1->px(), recoTrack1->py(),recoTrack1->pz(), recoTrack1->p());
269 
270  charge = recoTrack1->charge()*recoTrack2->charge();
271  if (charge < 0) {
272  InvMass = (Mu1+Mu2).M();
273  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
274  if(fabs(recoTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack1->eta())<EtaCutMax[iEtaRegion] &&
275  fabs(recoTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack2->eta())<EtaCutMax[iEtaRegion]){
276  if (InvMass < LowMassMax) TrkTrkMuon_LM[iEtaRegion]->Fill(InvMass);
277  if (InvMass > HighMassMin) TrkTrkMuon_HM[iEtaRegion]->Fill(InvMass);
278  }
279  }
280  }
281 
282 
283  LogTrace(metname)<<"[DiMuonHistograms] Soft-Soft pair"<<endl;
284 
285  if (muon::isSoftMuon(*muon1, vtx) &&
286  muon::isSoftMuon(*muon2, vtx) ) {
287 
288  if (charge < 0) {
289  InvMass = (Mu1+Mu2).M();
290  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
291 
292  if(fabs(recoTrack1->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack1->eta())<EtaCutMax[iEtaRegion] &&
293  fabs(recoTrack2->eta())>EtaCutMin[iEtaRegion] && fabs(recoTrack2->eta())<EtaCutMax[iEtaRegion]){
294  SoftSoftMuon[iEtaRegion]->Fill(InvMass);
295  }
296  }
297  }
298  }
299  }
300  } //muon2
301  } //Muon1
302 
303 
304 }
std::string metname
std::vector< MonitorElement * > MediumMediumMuon
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
std::vector< MonitorElement * > GlbGlbMuon_LM
math::Error< dimension >::type Error
covariance error matrix (3x3)
Definition: Vertex.h:43
bool isMediumMuon(const reco::Muon &)
bool isLooseMuon(const reco::Muon &)
bool isSoftMuon(const reco::Muon &, const reco::Vertex &)
edm::EDGetTokenT< edm::View< reco::Muon > > theMuonCollectionLabel_
edm::EDGetTokenT< reco::BeamSpot > theBeamSpotLabel_
std::vector< MonitorElement * > SoftSoftMuon
math::XYZPoint Point
point in the space
Definition: Vertex.h:39
double BeamWidthX() const
beam width X
Definition: BeamSpot.h:86
std::vector< MonitorElement * > TrkTrkMuon_LM
std::vector< MonitorElement * > StaTrkMuon_LM
edm::EDGetTokenT< reco::VertexCollection > theVertexLabel_
bool isValid() const
Definition: HandleBase.h:74
#define LogTrace(id)
std::vector< MonitorElement * > TrkTrkMuon_HM
double sigmaZ() const
sigma z
Definition: BeamSpot.h:80
double BeamWidthY() const
beam width Y
Definition: BeamSpot.h:88
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
static int position[264][3]
Definition: ReadPGInfo.cc:509
std::vector< MonitorElement * > StaTrkMuon_HM
std::vector< MonitorElement * > GlbGlbMuon_HM
const Point & position() const
position
Definition: BeamSpot.h:62
std::vector< MonitorElement * > TightTightMuon
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
std::vector< MonitorElement * > LooseLooseMuon
void DiMuonHistograms::bookHistograms ( DQMStore::IBooker ibooker,
edm::Run const &  ,
edm::EventSetup const &   
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 68 of file DiMuonHistograms.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), conversionPostprocessing_cfi::etaBin, and DQMStore::IBooker::setCurrentFolder().

70  {
71 
72  ibooker.cd();
73  ibooker.setCurrentFolder(theFolder);
74 
75  int nBin[3] = {etaBin,etaBBin,etaEBin};
76  EtaName[0] = ""; EtaName[1] = "_Barrel"; EtaName[2] = "_EndCap";
77  test = ibooker.book1D("test","InvMass_{Tight,Tight}",100, 0., 200.);
78  for (unsigned int iEtaRegion=0; iEtaRegion<3; iEtaRegion++){
79 
80  GlbGlbMuon_LM.push_back(ibooker.book1D("GlbGlbMuon_LM"+EtaName[iEtaRegion],"InvMass_{GLB,GLB}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
81  TrkTrkMuon_LM.push_back(ibooker.book1D("TrkTrkMuon_LM"+EtaName[iEtaRegion],"InvMass_{TRK,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
82  StaTrkMuon_LM.push_back(ibooker.book1D("StaTrkMuon_LM"+EtaName[iEtaRegion],"InvMass_{STA,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
83 
84  GlbGlbMuon_HM.push_back(ibooker.book1D("GlbGlbMuon_HM"+EtaName[iEtaRegion],"InvMass_{GLB,GLB}"+EtaName[iEtaRegion],nBin[iEtaRegion], HighMassMin, HighMassMax));
85  TrkTrkMuon_HM.push_back(ibooker.book1D("TrkTrkMuon_HM"+EtaName[iEtaRegion],"InvMass_{TRK,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], HighMassMin, HighMassMax));
86  StaTrkMuon_HM.push_back(ibooker.book1D("StaTrkMuon_HM"+EtaName[iEtaRegion],"InvMass_{STA,TRK}"+EtaName[iEtaRegion],nBin[iEtaRegion], HighMassMin, HighMassMax));
87 
88  // arround the Z peak
89  TightTightMuon.push_back(ibooker.book1D("TightTightMuon"+EtaName[iEtaRegion],"InvMass_{Tight,Tight}"+EtaName[iEtaRegion], nBin[iEtaRegion], LowMassMin, LowMassMax));
90  MediumMediumMuon.push_back(ibooker.book1D("MediumMediumMuon"+EtaName[iEtaRegion],"InvMass_{Medium,Medium}"+EtaName[iEtaRegion],nBin[iEtaRegion], LowMassMin, LowMassMax));
91  LooseLooseMuon.push_back(ibooker.book1D("LooseLooseMuon"+EtaName[iEtaRegion],"InvMass_{Loose,Loose}"+EtaName[iEtaRegion], nBin[iEtaRegion], LowMassMin, LowMassMax));
92 
93  // low-mass resonances
94  SoftSoftMuon.push_back(ibooker.book1D("SoftSoftMuon"+EtaName[iEtaRegion],"InvMass_{Soft,Soft}"+EtaName[iEtaRegion],nBin[iEtaRegion], 0.0, 55.0));
95  }
96 }
std::vector< MonitorElement * > MediumMediumMuon
void cd(void)
Definition: DQMStore.cc:269
std::vector< MonitorElement * > GlbGlbMuon_LM
std::string EtaName[3]
std::vector< MonitorElement * > SoftSoftMuon
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
std::vector< MonitorElement * > TrkTrkMuon_LM
std::vector< MonitorElement * > StaTrkMuon_LM
std::vector< MonitorElement * > TrkTrkMuon_HM
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
std::string theFolder
std::vector< MonitorElement * > StaTrkMuon_HM
std::vector< MonitorElement * > GlbGlbMuon_HM
std::vector< MonitorElement * > TightTightMuon
std::vector< MonitorElement * > LooseLooseMuon

Member Data Documentation

int DiMuonHistograms::etaBBin
private

Definition at line 54 of file DiMuonHistograms.h.

int DiMuonHistograms::etaBin
private

Definition at line 53 of file DiMuonHistograms.h.

double DiMuonHistograms::etaBMax
private

Definition at line 64 of file DiMuonHistograms.h.

double DiMuonHistograms::etaBMin
private

Definition at line 63 of file DiMuonHistograms.h.

double DiMuonHistograms::EtaCutMax
private

Definition at line 62 of file DiMuonHistograms.h.

double DiMuonHistograms::EtaCutMin
private

Definition at line 61 of file DiMuonHistograms.h.

int DiMuonHistograms::etaEBin
private

Definition at line 55 of file DiMuonHistograms.h.

double DiMuonHistograms::etaECMax
private

Definition at line 66 of file DiMuonHistograms.h.

double DiMuonHistograms::etaECMin
private

Definition at line 65 of file DiMuonHistograms.h.

std::string DiMuonHistograms::EtaName[3]
private

Definition at line 59 of file DiMuonHistograms.h.

int DiMuonHistograms::etaOvlpBin
private

Definition at line 56 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_HM
private

Definition at line 76 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::GlbGlbMuon_LM
private

Definition at line 75 of file DiMuonHistograms.h.

double DiMuonHistograms::HighMassMax
private

Definition at line 73 of file DiMuonHistograms.h.

double DiMuonHistograms::HighMassMin
private

Definition at line 72 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::LooseLooseMuon
private

Definition at line 82 of file DiMuonHistograms.h.

double DiMuonHistograms::LowMassMax
private

Definition at line 71 of file DiMuonHistograms.h.

double DiMuonHistograms::LowMassMin
private

Definition at line 70 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::MediumMediumMuon
private

Definition at line 83 of file DiMuonHistograms.h.

std::string DiMuonHistograms::metname
private

Definition at line 50 of file DiMuonHistograms.h.

int DiMuonHistograms::nGlbGlb
private

Definition at line 99 of file DiMuonHistograms.h.

int DiMuonHistograms::nLooseLoose
private

Definition at line 98 of file DiMuonHistograms.h.

int DiMuonHistograms::nMediumMedium
private

Definition at line 97 of file DiMuonHistograms.h.

int DiMuonHistograms::nTightTight
private

Definition at line 96 of file DiMuonHistograms.h.

edm::ParameterSet DiMuonHistograms::parameters
private
std::vector<MonitorElement*> DiMuonHistograms::SoftSoftMuon
private

Definition at line 85 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_HM
private

Definition at line 78 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::StaTrkMuon_LM
private

Definition at line 77 of file DiMuonHistograms.h.

MonitorElement* DiMuonHistograms::test
private
edm::EDGetTokenT<reco::BeamSpot> DiMuonHistograms::theBeamSpotLabel_
private

Definition at line 92 of file DiMuonHistograms.h.

std::string DiMuonHistograms::theFolder
private

Definition at line 94 of file DiMuonHistograms.h.

edm::EDGetTokenT<edm::View<reco::Muon> > DiMuonHistograms::theMuonCollectionLabel_
private

Definition at line 90 of file DiMuonHistograms.h.

MuonServiceProxy* DiMuonHistograms::theService
private

Definition at line 46 of file DiMuonHistograms.h.

edm::EDGetTokenT<reco::VertexCollection> DiMuonHistograms::theVertexLabel_
private

Definition at line 91 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::TightTightMuon
private

Definition at line 84 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_HM
private

Definition at line 80 of file DiMuonHistograms.h.

std::vector<MonitorElement*> DiMuonHistograms::TrkTrkMuon_LM
private

Definition at line 79 of file DiMuonHistograms.h.